- 积分
- 2139
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-4-15
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
我用的是CMORPH降水融合资料,共24个时次,时间间隔为一小时。在run *.gs时总会报错non-numeric args to numeric operation。以下是我的ctl和gs文件,请大神帮忙看一下问题在哪里,谢谢~
ctl文件:
DSET F:\main\CMORPH\SEVP_CLI_CHN_MERGE_FY2_PRE_HOUR_GRID_0.10-%y4%m2%d2%h2.grd
UNDEF -999.0
OPTIONS little_endian template
TITLE China Hourly Merged Precipitation Analysis
xdef 700 linear 70.05 0.10
ydef 440 linear 15.05 0.10
ZDEF 1 LEVELS 1
TDEF 24 LINEAR 00Z31Mar2013 1hr
VARS 2
crain 1 00 CH01 combined analysis (mm/Hour)
gsamp 1 00 CH02 gauge numbers
ENDVARS
gs文件:
'reinit''open F:\main\CHN_PRCP_HOUR_MERG_DISPLAY_0.1deg.lnx.ctl''set grads off''set grid off''set mpdset cnworld''set font 5''set lat 17.9 20.5''set lon 108.2 111.5''set z 1''set gxout shaded''set xlopts 1 5 0.2''set ylopts 1 5 0.2''set xlint 0.5''set cmin 0.01''set t 1''r1=crain''set t 2''r2=crain''set t 3''r3=crain'(...省略中间几个时次的设置,都差不多)'set t 23''r23=crain''set t 24''r24=crain' 'set clevs 4 8 10 12 14 16 18 20 22 24''province-basemap hain smth9(r1+r2+r3+r4+r5+r6+r7+r8+r9+r10+r11+r12+r13+r14+r15+r16+r17)''q w2xy 110.37 20.00' x=subwrd(result,3) y=subwrd(result,6) 'writehz 'x+0.1' 'y-0.1' 琼山 1 2.3 1 1.0 0.5 0 20' 这是提示错误的第87行'draw mark 2 'x' 'y' 0.06''q w2xy 110.25 20.00' x=subwrd(result,3) y=subwrd(result,6) 'writehz 'x-0.2' 'y-0.3' 海口 1 2.3 1 1 0.5 0 20''draw mark 2 'x' 'y' 0.06''q w2xy 108.67 19.1' x=subwrd(result,3) y=subwrd(result,6) 'writehz 'x+0.1' 'y-0.1' 东方 1 2.3 1 1 0.5 0 20''draw mark 2 'x' 'y' 0.06''q w2xy 109.68 19.91' x=subwrd(result,3) y=subwrd(result,6) 'writehz 'x-0.1' 'y-0.1' 林高 1 2.3 1 1 0.5 0 20''draw mark 2 'x' 'y' 0.06'(还有其他很多站点,节省空间我就不贴了,都差不多)'cbar_interp 1 1 1 9.5 4.2''set strsiz 0.18''draw string 9.65 6.7 mm''printim F:\main\1a.png white'; run这个gs文件,会出现如下的错误提示:
|
-
|