爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 3526|回复: 1

[分享资料] 报错 Data has 39823 values <= zero These were set to the undefined value

[复制链接]

新浪微博达人勋

发表于 2019-1-29 13:30:26 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

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'

求助这个错误怎么改啊?



密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2019-1-31 09:09:02 | 显示全部楼层
没人帮我...只能自己来了,通过尝试,发现最后加上
'printim c:\test2\thetase51.png white'
'c'
;
就可以了,也就是把图输出到文件中,同时检查一下公式的逻辑关系、参数设置的错误就可以了
一开始想的是,把set和输出的语句都忽略掉,先让grads能出图以后再调图,现在看来直接就把调图语句加上不容易出错,当然还是要明白每句话说的是啥。
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表