grads循环问题 坐标轴为时间维度轴 有没有朋友愿意耐心的帮我看看 有点复杂

superman828

我的资料是ncep 40年的逐月平均高度场 H‘*是H’=H(月平均)-H（40年平均）；H'*=ave(H',lon=0,lon=360)
我只能大概把程序大概写成这样子，大神帮我看一样我哪些地方有错误
还是时间坐标轴需要设置吗 怎么设置？
万分感谢
介于程序下载扣金钱 那我就复制到下面咯
'reinit'
'sdfopen E:\hjw\hgt.mon.mean.nc'
'set lon 0 360'
'set lat 30 60'
'set lev 500'
'define aa=ave(hgt,t=1,t=480)'
mo=1
while(mo<=12)
'define bb=0'
i=mo
while(i<=480)
'set t 'i''
'define bb=hgt-aa'
'define bb=ave(bb,lon=0,lon=360)'
'define cc=bb*bb'
i=i+12
endwhile
'define cc=cc/40'
'define cc=ave(cc,lon=0,lon=360)'
'define cc=cc**0.5'
'set gxout contour'
'd cc'
mo=mo+1
endwhile
'printim e:\Ia.png white'

superman828

图中的变量就是公式需要计算的Ia

superman828

e:\a.png

治愈茶包

前来学习学习