- 积分
- 50
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-5-4
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 monkey小可 于 2015-4-13 22:05 编辑
各位版友好!
小弟最近用grads用6小时资料计算月平均的结果,用的gs文件如下:
'reinit'
'sdfopen D:\Paper\0413/LU1990_2040-2060.nc'
'set fwrite D:\Paper\0413/LU1990-2040-2060.grd'
'set gxout fwrite'
data0=0
year=2040
while(year<=2060)
month=1
while(month<=12)
if(month=1|month=3|month=5|month=7|month=8|month=10|month=12)
data=124
else
if(month=4|month=6|month=9|month=11)
data=120
else
if(month=2&(year=2020|year=2024|year=2028|year=2032|year=2036|year=2040))
data=116
else
if(month=2)
data=112
endif
endif
endif
endif
data0=data0+data
t1=data0-data+1
t3=data0
say t1 t3
'set x 1 116'
'set y 1 102'
'set z 1'
'set t 1'
'define cc=ave(t2,t='t1',t='t3')'
'd cc'
*'d ave(rainc+rainnc,t='t1',t='t2')'
month=month+1
endwhile
year=year+1
endwhile
'disable fwrite'
'reinit'
但是得到的结果画出来的图却是这个样子的
1
contouring: -2.5e+38 to 2.5e+38 interval 5e+37
图的ctl文件是这样的:
dset D:\Paper\0413\LU1990-2040-2060.grd
options byteswapped
undef -5.66684e+22
title OUTPUT FROM WRF V3.6.1 MODEL
pdef 116 102 lcc 36.500 102.000 58.500 51.500 60.00000 30.00000 102.00000 50000.000 50000.000
xdef 426 linear 53.98895 0.22522523
ydef 236 linear 8.76476 0.22522523
zdef 1 linear 1 1
tdef 252 linear 00Z01JAN2020 1mo
VARS 1
T2 1 0 TEMP at 2 M (K)
*RAINC 1 0 ACCUMULATED TOTAL CUMULUS PRECIPITATION (mm)
ENDVARS
我已经被这个问题困扰了好久,以前也遇到过这种问题,感觉像是缺测值参与了计算,不知道是不是修改一下缺测值就ok呢?还是gs文件里面出了问题呢?请各位版友不吝赐教!谢谢啦
|
|