去气候态的相关

末端

下面的是我写的gs，不知道为什么老是出错。jul_ty.ctl是我提取的30的七月份的台风个数去气候态的结果，如果用7月份SST去气候态的结果就不会错，jun_anomaly.ctl是30年的六月份SST去气候态的结果，其他月份的都不行。
'reinit'
'set gxout fwrite'
'set fwrite d:\test\test052301.grd'
'open e:\wmf01\month\jul_ty.ctl'
'set t 1 30'
'index=anomaly'
'close 1'
'open e:\wmf01\month\jun_anomaly.ctl'
'set gxout shaded'
'set x 1 316'
'set y 1 98'
'set z 1'
'set t 1'
'd tcorr(index,tanomaly,t=1,t=30)'
'disable fwrite'
'close 1'
'reinit'
jul_ty.ctl如下
dset e:\wmf01\month\jul_ty.grd
undef 9.9e+08
xdef 1 linear 1 1
ydef 1 linear 1 1
zdef 1 levels 1000
tdef 30 linear jul1979 1yr
vars 1
anomaly 1 99 unit:
endvars
jun_anomaly.ctl如下

dset e:\wmf01\month\jun_anomaly.grd
undef -9.99E+08
xdef 316 linear 102.25 0.5
ydef 98 linear -18.25 0.5
zdef 1 levels 5
tdef 30 linear jun1979 1yr
vars 1
tanomaly 1 99 unit:deg
endvars
jul_anomaly.ctl如下
dset e:\wmf01\month\jul_anomaly.grd
undef -9.99E+08
xdef 316 linear 102.25 0.5
ydef 98 linear -18.25 0.5
zdef 1 levels 5
tdef 30 linear jul1979 1yr
vars 1
tanomaly 1 99 unit:deg
endvars
不知道我的什么地方出错了

做个霸气的木头

你的脚本是为了提取数据，但是中途加了句‘set gxout shaded’，会画出图形来的，文件就不会输出，你把这行去了，还有，不同的月份天数不一样，你注意下！