grads计算位势高度差，脚本出错！

看云淡风清

我想做的是计算两个区域的位势高度差，资料为hgt.mon.mean.nc，首先我想提取1961-2014年冬季的数据，然后再计算两个区域的位势高度差，求大神帮我看看，谢谢！
提取冬季的gs：
'reinit'
'sdfopen  E:\hgt.mon.mean.nc'
'set fwrite E:\hgt\hgt_winter_1961_2014.grd'
'set gxout fwrite'
'set x 1 144'
'set y 1 73'
'set lev 500'
i=168
while(i<=806)
'd ave(hgt,t='i',t='i+2')'
i=i+12
endwhile
'disable fwrite'
;
配置ctl：
dset E:\hgt\hgt_winter_1961_2014.grd
undef  -9.96921e+36
title monthly mean hgt from the NCEP Reanalysis
xdef  144 linear  0  2.5
ydef  73 linear  -90 2.5
zdef  1 linear 1 1
tdef  54 linear 1961 1yr
vars 1
hgt 1 99  monthly mean geopotential height
endvars
计算位势高度差:
'reinit'
'open  E:\hgt\hgt_winter_1961_2014.ctl'
'set x 1'
'set y 1'
'set z 1'
'set t 1  54'
'define ghta=tloop(aave(hgt,lon=55,lon=70,lat=40,lat=55))'
'define ghtb=tloop(aave(hgt,lon=130,lon=140,lat=35,lat=45))'
'define ghta1=ghta-ave(ghta,t=1,t=54)'
'define ghtb1=ghtb-ave(ghtb,t=1,t=54)'
'd (ghta1-ghtb1)'
;

huyuhengwan

哪一个开始出错的？fwrite就出错了？

river

错误提示？

妖怪别跑

？？？？？？？？

看云淡风清

huyuhengwan 发表于 2016-7-26 12:34
哪一个开始出错的？fwrite就出错了？

显示的是ctl的那个时间出错，该怎么设啊

看云淡风清

river 发表于 2016-7-26 12:45
错误提示？

错误提示是ctl的时间设置有问题，我第1步提取了每年冬季平均的位势高度，那个.grd文件的ctl该怎么写啊

造字仓颉

我也希望学习这方面的东西

风之精灵

ave(ghta,t=1,t=54)请重新定义一个新的变量，并使用modify命令使其具备与ghta相同的时间序列，再计算ghta1。同理，ghtb1的计算也如此修改。

river

看云淡风清 发表于 2016-7-26 15:03
错误提示是ctl的时间设置有问题，我第1步提取了每年冬季平均的位势高度，那个.grd文件的ctl该怎么写啊

先检查你提取出来的资料是不是正常，看看出图情况怎么样

赤道雨

来学习的？？？！！！！！！！！