大神帮我看看我计算假相当位温的代码是什么问题，感激不尽。

杨玲玲1014

'reinit'
'sdfopen e:\thse\thse.nc'
'enable print e:\thse\1.gmf'
'set grads off'
'set grid off'
'set lon 110 122'
'set lat 28 34'
'set mpdset cnworld'
'set map 1 1 1'
'set vpage 0 11 0 8.5'
'set parea 0.8 10.8 1 8.3'
'set xlopts 1 4 0.20'
'set ylopts 1 4 0.20'
'set xlint 5'
'set ylint -2'
'set csmooth on'
'set z 1'
'set t 1'
'set grads off'
'define es=6.112*exp(17.67*(t-273.15)/(t-35))'
'define e=850*q/(0.62197+q)+1e(-10)'
'define tlcl=55+2840.0/(3.5*log(t)-log(e)-4.805)'
'define theta=t*pow((1000/850),(0.2854*(1.0-0.28*q)))'
'define eqt=theta*exp(((3376.0/tlcl)-2.54)*q*(1.0+0.81q))'
'set ccolor 1'
'set gxout contour'
'set grads off'
'set cint 4'
'd eqt-273.15'
'set mpdset cnriver'
'set map 1 1 10'
'draw map'
'print'
'disable print'
'reinit

lqouc

无图，无ctl，无报错，无问题描述。
只能说850hPa上eqt的计算公式没错，其他的猜不出。

ps：grads的问题发到fortran版了。

杨玲玲1014

直接用的是nc格式的数据，所以没有写ctl文件，[img=29,30]file:///C:/Users/Administrator/AppData/Roaming/Tencent/Users/2286280986/QQ/WinTemp/RichOle/XD2WBUWD[%7D%6035_SFTX]FAN9.png[/img]

river

*------------------------------------------------------------------
*  GrADS script:  the.gs
*
*  Script purpose:  Computes equivalent potential temperature (theta-e) in K.
*                   Uses the definition of theta-e from Bolton (1980).
*
*  Author:  George H. Bryan, Penn State:  bryan@essc.psu.edu
*
*  Last modified:  10 July 2000
*
*  Assumes:  t = temperature (K)
*            prs = pressure (Pa)
*            q = water vapor mixing ratio (unitless)
*-----------------------------------------------------------------
'define prs=lev'
'define es=100*(6.112*exp(17.67*(t-273.15)/(t-29.65)))'
"define q=rh*(0.62197*es/(prs-es))/100."
"define e=(prs/100.)*q/(0.62197+q)+1e-10"
"define tlcl=55.0+2840.0/(3.5*log(t)-log(e)-4.805)"
"undefine e"
"define theta=t*pow((100000./prs),(0.2854*(1.0-0.28*q)))"
"define the=theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q))"
"undefine tlcl"
"undefine theta"

杨玲玲1014

错误代码是：Expected operator or '>'
Error ocurred at column 20
DEFINE error: Invadid expression.
0 [sig] opengrads 3408 open _stackdumpfile: dumping stack trace to opengrads.exe.stackdump

lqouc

杨玲玲1014 发表于 2015-9-28 22:03
错误代码是：Expected operator or '>'
Error ocurred at column 20
DEFINE error: Invadid expression.
...

第20行看起来并没有什么错。。。。。感觉没少东西
重启下电脑或者换个电脑也这样？
不行就试试四楼给的代码，自己去掉有rh的那句，应该是差不多的。

杨玲玲1014

谢谢，终于改好咯

catgrace_z

杨玲玲1014 发表于 2015-9-29 14:54
谢谢，终于改好咯

你好，请问你是怎么解决的？我也出现了类似的问题