请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 7215|回复: 11

[分享资料] 求助:怎么用GrADS设置高度变量

[复制链接]

新浪微博达人勋

发表于 2013-11-6 20:33:43 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
我现在用的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

新浪微博达人勋

0
早起挑战累计收入
发表于 2013-11-6 20:37:02 | 显示全部楼层
楼主恳求你好好先看基础知识
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-7 07:32:34 | 显示全部楼层
说清楚问题,你想要什么,高度循环有好多方法,你想的这个是最麻烦的。而且不是所有能在gs里运行的命令就能直接输入黑屏窗口运行的
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-7 09:42:32 | 显示全部楼层
z=1
while(z<=17)
'set z 'z''
'd something'
z=z+1
endwhile
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-11-7 17:57:45 | 显示全部楼层

我想用一个纬向风的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,什么都不明白净问一些特傻的问题,见笑了~
还是感谢你的回答!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-11-7 17:59:26 | 显示全部楼层
边缘人糊涂乐 发表于 2013-11-7 09:42
z=1
while(z

谢谢你的回答!菜鸟倍受鼓励~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-7 18:36:39 | 显示全部楼层
yrovl 发表于 2013-11-7 17:57
我想用一个纬向风的NC数据计算全球的质量流函数,这个数据包括17层,一个前辈给了一个计算质量流函数的参 ...

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

新浪微博达人勋

 楼主| 发表于 2013-11-7 19:41:21 | 显示全部楼层
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
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-11-7 19:44:37 | 显示全部楼层
额。。。图片好像没显示
错误提示是说:
syntax error:invalid operand
'psi1'not a variable or function name
error ocurred at column 1
DISPLAY error:Invalid expression
   expression=psi1
.
.
.

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-7 22:15:03 | 显示全部楼层
yrovl 发表于 2013-11-7 19:41
这是我的gs文件
'reinit'
'sdfopen c:\V\vwnd.mon.ltm.nc'

看的稀里糊涂的,水平有限,等别人来看吧
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表