求助：怎么用GrADS设置高度变量

yrovl

我现在用的NC数据数据共有17层，我想在开始的设置中设置每一层的高度，在gs文件中写到
lev.1=1000;lev.2=925;lev.3=850......直至lev.17=10
可是因为gs文件运行总出错，我就把gs文件中的命令一句一句的输到黑屏窗口中，结果输到lev.1=1000时，就会提示我unknown command:lev.1=1000
但或直接输lev1=1000就没有错误。可是我后面要做高度上的循环，只能用lev.1=1000的形式，不知道该怎么办，
grads菜鸟，请高人指点迷津！感激不尽~

mofangbao

楼主恳求你好好先看基础知识

river

说清楚问题，你想要什么，高度循环有好多方法，你想的这个是最麻烦的。而且不是所有能在gs里运行的命令就能直接输入黑屏窗口运行的

边缘人糊涂乐

z=1
while(z<=17)
'set z 'z''
'd something'
z=z+1
endwhile

yrovl

river 发表于 2013-11-7 07:32

                               
登录/注册后可看大图

说清楚问题，你想要什么，高度循环有好多方法，你想的这个是最麻烦的。而且不是所有能在gs里运行的命令就能 ...

我想用一个纬向风的NC数据计算全球的质量流函数，这个数据包括17层，一个前辈给了一个计算质量流函数的参考答案，在里面他把这17层分别设置了出来：
'reinit'
'sdfopen c:\gddzy\vwnd.mon.ltm.nc'

lev.1=1000 ; lev.2=925 ; lev.3=850 ; lev.4=700  ; lev.5=600  ; lev.6=500
lev.7=400  ; lev.8=300 ; lev.9=250 ; lev.10=200 ; lev.11=150 ; lev.12=100
lev.13=70  ; lev.14=50 ; lev.15=30 ; lev.16=20  ; lev.17=10

'set gxout fwrite'
'set fwrite c:\gddzy\psai.mon.ltm.dat'

a=6.37e6
pi=3.14159
g=9.8
我不明白为什么经纬度不用一个点一个点的设置，而高度就要这样一个高度层一个高度层地赋值。
后来，我进行了一下修改后，我的gs文件就不能运行出来了，所以就想在黑屏上一句一句输入我的命令来查错了，然后就出现了我昨天发的问题。
第一次用grads，什么都不明白净问一些特傻的问题，见笑了~
还是感谢你的回答！

yrovl

边缘人糊涂乐 发表于 2013-11-7 09:42

                               
登录/注册后可看大图

z=1
while(z

谢谢你的回答！菜鸟倍受鼓励~

river

yrovl 发表于 2013-11-7 17:57

                               
登录/注册后可看大图

我想用一个纬向风的NC数据计算全球的质量流函数，这个数据包括17层，一个前辈给了一个计算质量流函数的参 ...

不是问题傻不傻，是没问到点子上去。还有你只贴了一部分gs，别人看了可能该是一头雾水，来回答问题的人找不到解决问题的切入点，给人的感觉就是你留了一道作业题，要求看贴的人来直接给你一个答案。所以你把你修改过的gs全部贴上来好了，你修改的是哪几句，这样别人才可以根据你的叙述和你的gs来帮你找问题

yrovl

river 发表于 2013-11-7 18:36

                               
登录/注册后可看大图

不是问题傻不傻，是没问到点子上去。还有你只贴了一部分gs，别人看了可能该是一头雾水，来回答问题的人找 ...

这是我的gs文件
'reinit'
'sdfopen c:\V\vwnd.mon.ltm.nc'
lev.1=1000; lev.2=925;lev.3=850;lev.4=700;lev.5=600;lev.6=500
lev.7=400; lev.8=300; lev.9=250; lev.10=200; lev.11=150; lev.12=100;
lev.13=70; lev.14=50; lev.15=30; lev.16=20; lev.17=10

'set gxout fwrite'
'set fwrite c:\V\psi.dat'
a=6.37e6
pi=3.14159
g=9.78


zz=1
while(zz<=16)
  temp=zz+1
  dp.zz=(lev.temp-lev.zz)*100
  zz=zz+1
endwhile

'set x 1 144'
'set y 1 73'

tt=1
while(t<=12)
   'set z 1 17'
   'set tt'tt
   'define Vave=ave(vwnd,lon=0,lon=360,-b)'

    zz=17
    'set z'zz
    'define psi'zz'=0'

     zz=16
     while(zz>=1)
      temp1=zz+1
      'define Xishu=2*'pi'*'a'*cos(lat*'pi'/180)/'g'
      'psi'zz'=psi'temp1'-Xishu*Vave*'dp.zz''
       zz=zz-1
      endwhile
tt=tt+1
endwhile

   zz=1
    while(zz<=17)
     'd psi'zz
    zz=zz+1
    endwhile
'disable fwrite'
这个是质量流函数的计算公式
[v]=g/(2πa cos&#8289;&#8709; )  &#8706;φ/&#8706;p
其中[v]是全球纬向平均风，φ就是质量流函数
运行gs文件以后，黑屏上就总提示这个错误
c:\gddzy\55.png

yrovl

额。。。图片好像没显示
错误提示是说：
syntax error:invalid operand
'psi1'not a variable or function name
error ocurred at column 1
DISPLAY error:Invalid expression
   expression=psi1
.
.
.

river

yrovl 发表于 2013-11-7 19:41

                               
登录/注册后可看大图

这是我的gs文件
'reinit'
'sdfopen c:\V\vwnd.mon.ltm.nc'

看的稀里糊涂的，水平有限，等别人来看吧