对各个等压级的比湿求和问题？

hustlh

比湿shum数据等压级有5级（set z 1 5），我想把各个等级的比湿相加求和，然后输出；
写脚本过程中出了问题，请高手帮我看看

'set z 1 5'
'set t 1 365'
'set gxout fwrite'
'set fwrite g:\sqhl.grd'
'define a=shum(z=1)'
'define b=shum(z=2)'
'define c=shum(z=3)'
'define d=shum(z=4)'
'define e=shum(z=5)'
'set z 1'
'define sqhl=sum(0.5*a+b+c+d+0.5*e)'
i=1
while(i<=365)
'set t 'i''
'd sqhl'
i=i+1
endwhile
'disable fwrite'
'reinit'


我在想是不是define有问题，但是我又不知道怎么改正，希望大家帮帮指正！

ps：很急~~~~    拜托各位了

river

  你这个整个不理解grads的自带函数用法啊，sum和ave 、vint、hdivg之类的一样，都是有确定的用法的。不适随随便就能用的。自己去看grads官方文档，用法和注意事项都很详细。
  在gs求和和fortran里面是一样啊，直接加就行了，干吗还要弄个sum函数啊。
  楼主grads基础不好。grads和gs就相当与fortran和for f90的关系是差不多的，grads是用来运行代码的，gs就是你编的代码。如果你在做fortran程序的时候出错了，你是怎么找错的呢？肯定首先是fortran程序的错误提示是吧，那grads也一样，出了错有错误提示，那就是改错的线索。可是楼主提问怎么不贴错误提示呢？

river

还有一点就是，没必要提前定义比湿的高度，在求和的过程中固定就可以。
'define sqhl=shum(z=1)+shum(z=2)+shum(z=3)+shum(z=4)+shum(z=5)'
还有就是这一句要放到循环里

hustlh

river 发表于 2013-12-5 17:41

                               
登录/注册后可看大图

还有一点就是，没必要提前定义比湿的高度，在求和的过程中固定就可以。
'define sqhl=shum(z=1)+shum(z=2) ...

好的，清楚了。感谢river了！