也谈ncl计算 假相当位温

活在当下de

参考了好几个帖子的计算方法，总算算出来了，代码如下：

pre_levels = 850
tmp_ave    = 10.
rhum_ave   = 10^((7.5*6)/(237.3+6)) / 10^((7.5*10)/(237.3+10)) *100

es = where( tmp_ave.lt.0, 6.112*exp(17.67* tmp_ave /( tmp_ave +243.5)),\
                           6.11*10^(7.5* tmp_ave /( tmp_ave +237.3)) )
q  = rhum_ave * (0.62197*es/( pre_levels - 0.378*es))/100.
e  = pre_levels * q/(0.62197+q)+10^(-10)
tlcl  = 55.0+2840.0/(3.5*log(tmp_ave +273.15)-log(e)-4.805)
theta = ( tmp_ave +273.15)*(1000. / pre_levels )^(0.2854*(1.0-0.28*q))
eqt   = theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q))   ;假相当位温


主要考虑了 es 在0℃的分界，结果与《气象常用参数和物理量查算表》中的结果一致：

343625321

楼主，你画的850hPa的，prs是固定的。但是如果想画剖面图怎么办呢?

活在当下de

哇哦，感谢清风的贡献！{:5_235:}

shiyan101

感谢楼主的贡献非常有用

指甲钳

请问楼主，tmp_ave    = 10.是固定的吗？还是按照不同的资料计算？

活在当下de

指甲钳 发表于 2015-2-6 07:41
请问楼主，tmp_ave    = 10.是固定的吗？还是按照不同的资料计算？

不是固定的，程序中算的事 850hPa、10℃ 时的假相当位温

超台十八

这些公式，找 wmo 文件比较有谱。

slc2005011020

感谢楼主
我是学习的

993099565

赞！！！！！！！！！！

sunn_huia

这些公式，找 wmo 文件比较有谱。