- 积分
- 7778
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-6-10
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
在执行脚本时报错一下内容:
Warning from LOG: Data has 39823 values <= zero
These were set to the undefined value
我的脚本内容是这样的:
'reinit'
'open c:/test2/echhxa50.ctl '
'open c:/test2/echrxa85.ctl '
'open c:/test2/echtxa85.ctl '
'define tk=tmp850mb.3' ;*某高度层的开氏温度K
'define tc=tk-273.16' ;*某高度层的摄氏温度C
'define rh=rh850mb.2' ;*某高度层的相对湿度Relative humidity%
'define prs=hgt500mb.1' ;*获得某层高度的气压
******************************************************************************************
*求饱和水汽压Tetens经验公式
*水面es,tk开氏温度,tc摄氏温度
if(tk>273.16)
'define es=6.1078*exp(17.2693882*tc/(tk-35.86))'
endif
*冰面es,tk开氏温度,tc摄氏温度
if(tk<=273.16)
'define es=6.1078*exp(21.8745584*tc/(tk-7.66))'
endif
*饱和比湿
'define qs=0.622*es/(prs-0.378*es)'
*用相对湿度等求比湿
'q=rh*qs/100'
*水汽压
'e=prs*q/(0.622+q)'
*凝结高度的绝对温度,tk起始面上绝对温度K,
'define tlcl=55.0+2840.0/(3.5*log(tk)-log(e)-4.805)'
*求假相当位温Bolton公式,se为开氏温度K
'define theta=tk*pow((1000/prs),(0.2854*(1.0-0.28*q)))'
'define se=theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q))'
******************************************************************************************
'd se'
求助这个错误怎么改啊?
|
|