求大神帮忙看下在计算位涡过程中出现的问题

7795216223 · 发表于 2014-4-15 20:21:12

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

您需要登录才可以下载或查看，没有帐号？立即注册

x

'set lon 50 160'
'set lat 10 70'
'set z 1'
'define prs=100000'
'define qv=SPFHprs'
'define e=prs*qv/(0.62197+qv)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(TMPprs)-log(e)-4.805)'
'undefine e'
'define theta=TMPprs*pow((100000/prs),(0.2854*(1.0-0.28*qv)))'
'define thse=theta*exp(((3376./tlcl)-2.54)*qv*(1.0+0.81*qv))'
'd thse'
这是程序中的一部分，为保证正确分层进行计算，可还会出现问题？

lqouc · 发表于 2014-4-15 20:29:10

报错就是说你取对数的时候出现了负数，仔细检查你的公式是不是温度项应该为热力学温度。

7795216223 · 发表于 2014-4-15 20:34:05

lqouc 发表于 2014-4-15 20:29
报错就是说你取对数的时候出现了负数，仔细检查你的公式是不是温度项应该为热力学温度。

TMPprs是Temperature [K]

lqouc · 发表于 2014-4-15 20:50:51

7795216223 发表于 2014-4-15 20:34
TMPprs是Temperature [K]

那热力学温度怎么会是负数呢？
仔细检查一下原始数据吧。还有就是ctl的缺省设置。

7795216223 · 发表于 2014-4-15 21:21:08

lqouc 发表于 2014-4-15 20:50
那热力学温度怎么会是负数呢？
仔细检查一下原始数据吧。还有就是ctl的缺省设置。

嗯，谢谢了!

7795216223 · 发表于 2014-4-16 10:30:40

7795216223 发表于 2014-4-15 21:21
嗯，谢谢了!

看了一下数据，发现比湿存在<0的情况，不知道咋么处理。咋办啊？

		自动登录	找回密码
密码			立即注册