求助脚本运行出错

tszhaoweiwen

'reinit'
'open d:\1004\fnl.ctl'
'set t 1'
'set lon 108 120'
'set lat 20 25'
'set lev 850'
'define T850=tmpprs-273.16'
'define RH=rhprs'
'define Td850=T850-((14.55+0.114*T850)*(1-0.01*rh) + pow((2.5+0.007*T850)*(1-0.01*rh),3) + (15.9+0.117*T850)*pow((1-0.01*rh),14))'
'define F850=2*mag(UGRDprs,VGRDprs)'

'set lev 500'
'define T500=tmpprs-273.16'
'define F500=2*mag(UGRDprs,VGRDprs)'

'define I=12*Td850+20(T850+Td850-2*T500-49)+2*F850+F500'
'set cint 2'
'set cthick 6'
'd I'
'printim d:/I2.png  white'

运行脚本一直出错，求修改

river

提示什么错误了？提问要说清楚，尽可能的提供线索。

'define I=12*Td850+20(T850+Td850-2*T500-49)+2*F850+F500'
没有运算符号？

tszhaoweiwen

river 发表于 2015-10-26 17:26
提示什么错误了？提问要说清楚，尽可能的提供线索。

'define I=12*Td850+20(T850+Td850-2*T500-49)+2*F8 ...

真的是漏了个*号  谢谢你

tszhaoweiwen

river 发表于 2015-10-26 17:26
提示什么错误了？提问要说清楚，尽可能的提供线索。

'define I=12*Td850+20(T850+Td850-2*T500-49)+2*F8 ...

想用grasd画一个 强天气威胁指数SWEAT图，表达式为：
I=12D+20（TT-49）+2F850+F500+125（S+0.2）
其中I代表SWEAT；D＝850hPa露点温度(℃)，若D是负数，此项为0；
TT＝（T+Td）850-2T500，若TT小于49，则20(TT－49)项等于0；
grads可以实现这个IF语句吗？脚本也要‘’括起来吗？

麻瓜520

tszhaoweiwen 发表于 2015-10-26 19:45
想用grasd画一个 强天气威胁指数SWEAT图，表达式为：
I=12D+20（TT-49）+2F850+F500+125（S+0.2）
其 ...

循环语句和判断语句不需要加''

river

grads实现你说的这个当然没问题了，但是不是所有语句都需要用单引号括起来啊，主要是看是不是变量，一般的define的变量都是需要括起来的