GRADS求湿位涡出错，求大神指导！

故乡月(环)

脚本如下：
'reinit'
'sdfopen D:\air.2014.nc'
'sdfopen D:\rhum.2014.nc'
'sdfopen D:\uwnd.2014.nc'
'sdfopen D:\vwnd.2014.nc'


'set mpdset cnworld'
'draw map'
'set lat 0 50'
'set lon 100 120'


#'set ylevs 995 900 800 700 600 500 400 300'
'set lev 1000 100'
'set zlog on'


'set t 880'

*****计算假相当位温*****
'set grads off'
'define prs=lev'
'define es=6.112*exp(17.67*（air-273.15)/(air-29.65))'
'define qs=0.62197*es/(prs-0.378*es)'
'define qv=rhum.2*qs/100'
'define e=prs*qv/(0.62197+qv)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(air)-log(e)-4.775)'

'define theta=air*pow((1000./prs),(0.2854*(1.0-0.28*qv)))'
'define eqt=theta*exp(((3376./tlcl)-2.54)*qv*(1.0+0.81*qv))'

*****计算湿位涡3*****
'define vo=hcurl(uwnd.3,vwnd.4)'
'define f=2*7.292*sin(lat*3.14159/180.0)*0.00001'
'define g=9.8'
'define dp=100*(lev(z-1)-lev(z+1))'
'define deqt=eqt(z-1)-eqt(z+1)'
'define du=uwnd.3(z-1)-uwnd.3(z+1)'
'define dv=vwnd.4(z-1)-vwnd.4(z+1)'
'define dx=2.0*6370949.0*cos(lat*3.14159/180.0)*3.14159/180.0'
'define dy=2.0*6370949.0*3.14159/180.0'
'define dtx=cdiff(eqt,x)'
'define dty=cdiff(eqt,y)'
'define pv1=-g*(vo+f)*deqt/dp'
'define pv2=g*((dv/dp)*(dtx/dx)-(du/dp)*(dty/dy))'
'define pv=pv1+pv2'

'set lon 107'
'set lat 10 40'
'set gxout contour'
'set csmooth on'
'set cthick 7'
'set clopts -1 -1 0.10'
'd pv1*10e5'
*'d pv2*10e6'
*'d pv*10e5'
*'print'
*'disable print'
'q time'
res=subwrd(result,3)
'draw title 'res'    'pv

'printim D:\0606.jpg white'
运行的时候出现问题，grads卡着没有反应，grads提示的信息如下：

传说中的谁

grads不是告诉你括号"("有问题么

user_qxjy

应该是这行有语法错误：'define es=6.112*exp(17.67*（air-273.15)/(air-29.65))' ；其中“*（”可能是中文状态下输入的左括号，应改为西文状态下输入

river

括号的问题，注意英文状态下的才行

故乡月(环)

谢谢大家！一个小失误，搞了我整整一天，终于弄出来了，除了那个符号的问题外，rhum数据只有8层，玛德没有注意到，害的我积分的时候整死了！

当年明月

楼主在吗？想问问您最后是怎么解决的？改成括号了貌似也不行、、{:eb303:}

213de88

嗯 受教了

fyy123

想问一下为啥在网上其他地方看到的是
dp=100*(lev(z+1)-lev(z-1))'
'define deqt=eqt(z+1)-eqt(z-1)'
'define du=uwnd.3(z+1)-uwnd.3(z-1)'
'define dv=vwnd.4(z+1)-vwnd.4(z-1)'
不知道哪个是对的？
菜鸟求答