grads运行gs按回车就消失！

气象局

我在学着水汽输送场的grads程序，我借用了群论坛里的程序，自己做了些修改，由于自己是初学者可能会犯些低级错误，请大家见谅，有错的地方希望朋友能帮我指出来，感激不尽！下面是自己做了些修改的程序，请大家帮忙！
* This is a script for displaying moisture convergence
* Written by Michael Maxwell
*   
* rh    = Relative Humidity in %
* t     = Temp at *set level in degrees Kelvin
* tc    = Temp in degrees C
* td    = Dewpoint at *set level in degrees C
* e     = Vapor pressure
* mixr  = Mixing ratio
* u     = U-wind in m/s
* v     = V-wind in m/s
* mconv = moisture convergence/divergence. convergence is positive and divergence is
negative.
'reinit'
'sdfopen c:\zl\t.nc'
'sdfopen c:\zl\u.nc'
'sdfopen c:\zl\v.nc'
'sdfopen c:\zl\sr.nc'
'set lev 700'
'set lat 20 35'
'set lon 100 110'
'set t 64'
'set mpdset cnworld'
'set xlopts 1 4 0.15'
'set ylopts 1 4 0.15'
'set grads off'
'set timelab off'
'set grid off'
*‘tc=(tmpprs-273.16)'
'td=t.1-((14.55+0.114*t.1)*(1-0.01*q.4) + pow((2.5+0.007*t.1)*(1-0.01*q.4),3) +
(15.9+0.37*t.1)*pow((1-0.01*q.4),14))'
'vapr=6.112*exp((17.67*td)/(td+243.5))'
'e=vapr*1.001+(lev-100)/900*0.0034'
'define mixr=0.62137*(e/(lev-e))*1000/9.8'
通过此gs文件算出的水汽通量的单位是g/(cm.s),一般为 十几到几十
'define qx=u.2*mixr'
'define qy=v.3*mixr'
'define mconv=hdivg(qx,qy)*1e6'
'enable print c:\tx\sq.gmf'
'set arrscl 1 500'
'set lon 124'
'set cmin 30'
'set gxout shaded'
'd mag(qx,qy)'
'set cmin 30'
'set gxout contour'
'd mag(qx,qy)'
'd theta'
'd skip(u,8);skip(v,8)'
'print'
'disable print'
'printim c:/tx/qfs.gif white'
在grads一运行，按回车，grads就消失了，前提是grads能运行其他的gs文件，我觉得就是这个程序出问题了，请大家帮忙看看！

才仔

*tc=(tmpprs-273.16)'   *是打错了？？

river

如果是gs问题的话，一般不会导致grads消失的。所以感觉还是你的grads有问题。换台电脑运行试一下

liutaoz21

*tc=(tmpprs-273.16)'这句的星号“*“应该是撇“ ' “吧！

气象局

才仔 发表于 2013-6-19 22:03

                               
登录/注册后可看大图

*tc=(tmpprs-273.16)'   *是打错了？？

没有打错，这是原文件里的。我也不明白，但是它是注释我觉得也就没影响了！

气象局

river 发表于 2013-6-19 22:28

                               
登录/注册后可看大图

如果是gs问题的话，一般不会导致grads消失的。所以感觉还是你的grads有问题。换台电脑运行试一下

但是我的grads能运行其他的gs文件哦！我交换着操作，运行其他的gs文件没问题，就是运行这个gs文件就出问题。

气象局

liutaoz21 发表于 2013-6-20 03:37

                               
登录/注册后可看大图

*tc=(tmpprs-273.16)'这句的星号“*“应该是撇“ ' “吧！

原文件是#我查了很多，找不到在gs中#号表示什么，所以我就改为了*，" ' "是我弄丢了，我已经补上了，但是我认为*是注释哦，不影响运行哦！所以我认为问题可能不是这个造成的！

气象局

谢谢大家的参与，今早花了两个小时终于琢磨出来了，问题已解决了！主要问题是在那个露点温度的公式太长了，应该加上续行符，续行符为&。大家可以试一试哈！再次感谢大家！
！

才仔

气象局 发表于 2013-6-20 09:44

                               
登录/注册后可看大图

谢谢大家的参与，今早花了两个小时终于琢磨出来了，问题已解决了！主要问题是在那个露点温度的公式太长了， ...

在哪加，td那行最后面？？？

气象局

才仔 发表于 2013-6-20 12:24

                               
登录/注册后可看大图

在哪加，td那行最后面？？？

是的哦！