请教各位问题 谢谢大家回答

river

jeffggg 发表于 2013-5-6 11:52

                               
登录/注册后可看大图

这个是不等间隔的  循环的话 原来有232个数 要和这8个数做差 怎么实现呢

两层循环啊，一层循环年，一层循环33--40旬。你只要图的话，直接画图就行了吧，没必要生成资料再配ctl画图。
'reinit'
'sdfopen I:\CMAP\cmap\3\cmap\precip.pentad.mean.nc'
'set x 1 144'
'set y 1 72'
'set lev 0 '
'set t 33 40'
'A=ave(precip,t+0,t=2084,73)'
i=1
while(i<=29)
j=33
while(j<=40)
tt=(i-1)*73+j
'set t 'tt''
'define B=precip-A'
'd B'
'printim I:\CMAP\cmap\3\cmap\'tt'jupingNC.bmp white'
'c'
j=j+1
endwhile
i=i+1
endwhile
;
grads的基础还是不行啊，也不太会变通。还得继续努力啊

river

jeffggg 发表于 2013-5-6 11:52

                               
登录/注册后可看大图

这个是不等间隔的  循环的话 原来有232个数 要和这8个数做差 怎么实现呢

上面的能够出图，但是我不确定求出来正好是几个距平，因为有涉及到前面隐式循环求平均的，不知道能不能和下面距平对上，你自己验证吧