grads编程求助

凉茶

现有1948-2011年ncep月平均湿度的nc文件
现要求1958年至2011年冬季年平均湿度 得出grd文件
取值范围 70E-140E  0-60N  8层（300 400 500 600 700 850 925 1000）
有点急  谢谢各位大神啦

凉茶

求一个gs文件

topmad

明确的告诉你，求不到

自己动手哇，有不会的再来有针对性的提问，伸手党不可取

凉茶

topmad 发表于 2012-5-8 20:03

                               
登录/注册后可看大图

明确的告诉你，求不到

自己动手哇，有不会的再来有针对性的提问，伸手党不可取

我也不想做个伸手党啊 我自己做了不过好像出了点问题 自己调了好久都调不过来 所以就来了
你能给我留个qq么  帮我看看怎么样啊

topmad

自己调不通，把代码发上来，让懂的人帮你看看哪里出问题了

凉茶

topmad 发表于 2012-5-8 20:22

                               
登录/注册后可看大图

自己调不通，把代码发上来，让懂的人帮你看看哪里出问题了

'reinit'
'sdfopen d:\lunwen\yw\shum.mon.mean.nc'
'set gxout fwrite'
'set fwrite d:\lunwen\yw\win1.grd'
'set x 70 140'
'set y 0 60'
i=1
while(i<=54)
m=(i-1)*12
'set t 132'
'define winave=ave(shum,t='m',t='m'+2)'
'set lev 1000'
'd winave'
'set lev 925'
'd winave'
'set lev 850'
'd winave'
'set lev 700'
'd winave'
'set lev 600'
'd winave'
'set lev 500'
'd winave'
'set lev 400'
'd winave'
'set lev 300'
'd winave'
i=i+1
endwhile
'disable fwrite'
;
帮忙看看哈 谢啦

善人/jw

你提出来后，哪里有问题呢

兰溪之水

层次设置放到define前面去，还有t='m+2'

冰枫

层次也可以做个循环，用z表示

冰点

这里的经纬度x,y应该改为lon,lat……‘set t ’132这句应该是说设置时间为1958年12月吧，不过后面那句一写，这句好像无效了，不如直接将'set t '132去掉，而把求平均那句改为ave( ,t='m+132',t='m+134')……