GRADS处理NCEP资料的问题

东篱

计算NCEP位势高度场逐月平均hgt.mon.mean.nc中,夏季（6-8月）600高度的逐年平均，并写入文件f:\data\600.hgt.ave.dat中。编写了一个.gs文件，但是运行出错，请各位帮忙检查下问题出在哪里，谢谢
.gs:
'reinit'
'sdfopen f:\data\hgt.mon.mean.nc'
'set gxout fwrite'
'set fwrite f:\data\600.hgt.ave.dat'

'set lon 60 110'
'set lat 20 50'
'set lev 600'

tt=1;ist=0;ied=0

while(tt<=780)

        ist=tt+5;ied=tt+7
        'd ave(hgt, t='ist',t='ied')'
tt=tt+12

endwhile

'disable fwrite'
;


错误提示：
No hardcopy metafile open
All files closed; all defined objects released
All GrADS attributes have been reinitialized
Syntax Error
  Error occurred on line 14
  In file f:\data\600.hgt.ave.1.gs

x1n3r

是不是'set gxout fwrite'和'set fwrite f:\data\600.hgt.ave.dat'两句写反了，应该颠倒下位置。。

river

错误提示已经定位了，变量错误发生在第14行。楼主你那两句 ist=tt+5;ied=tt+7和 'd ave(hgt, t='ist',t='ied')'前面的空白的是怎么加进去的，是不是用Tab键了？

传说中的谁

每一句顶格写，不要像其他语言代码那样写的有层次感

冰点

tt=1
while(tt<=780)   
   'set t 'tt
   'd ave(hgt, t='tt'+5,t='tt'+7)'
tt=tt+12
endwhile

这样写应该可以，好像大家都这样写……
编辑的时候我把   ‘set t ’tt这一句放到while下面，保存之后它总是跑是上面那一行，不知道什么情况。

东篱

x1n3r 发表于 2013-5-23 06:25

                               
登录/注册后可看大图

是不是'set gxout fwrite'和'set fwrite f:\data\600.hgt.ave.dat'两句写反了，应该颠倒下位置。。

谢谢回复，不过这个位置没有问题。我已经解决了。

东篱

river 发表于 2013-5-23 07:19

                               
登录/注册后可看大图

错误提示已经定位了，变量错误发生在第14行。楼主你那两句 ist=tt+5;ied=tt+7和 'd ave(hgt, t='ist',t='ie ...

谢谢，真的是这样的。问题已经解决了。

东篱

传说中的谁 发表于 2013-5-23 08:53

                               
登录/注册后可看大图

每一句顶格写，不要像其他语言代码那样写的有层次感

谢谢，是这样的，以后会注意了。这已经是你第二次帮我解决问题了，太感谢了

东篱

冰点 发表于 2013-5-23 11:30

                               
登录/注册后可看大图

tt=1
while(tt

谢谢你的回复，主要问题是顶格问题，已经解决了。非常感谢。

东篱

东篱在这里感谢各位的帮忙了，气象论坛真的是个好地方。已经帮我解决了好多问题了，还可以学到不少东西