wrf模式输出结果画假相当位温的gs脚本（含公式说明和出处）

东风急流

转模式的同学常会分析假相当位温这个变量，但模式不能直接输出这个变量，大气物理学课本中的计算公式又较为复杂（需要计算抬升凝结高度处温度）。现在结合大气物理学课本中公式和中国气象局发布的《中华人民共和国气象行业标准》中的公式来计算假相当位温（家园中资料也很多，但较为杂乱），并对每一步做了详细说明，现在将主要部分贴出，如果有不严谨的地方还请大家指正。
1、在运行ARWpost的时候需要输出tc和td（露点温度）这两个变量,最好再输出theta(位温)
2、运行此脚本：
*计算饱和水汽压，TC是摄氏温标的温度(大气物理学P21,公式2.2.9a)
'es=6.1078*exp(17.27*tc/(tc+273.16-35.86))'
*计算饱和比湿(大气物理学P20,公式2.2.5)
'qs=0.622*es/(lev-0.378*es)'
*计算比湿            
'qa=RH*qs/100.0'

*计算假相当位温（不知道位温时）(中华人民共和国气象行业标准---南海夏季风监测指标P10)
'thse=(tc+273.15)*exp(0.28586*log(1000.0/lev)+2500.0*qa/(338.52-0.24*(tc+273)+1.24*td))'
*计算假相当位温（已知位温时）
*'thse=theta*exp(2500.0*qa/(338.52-0.24*(tc+273)+1.24*td))'
；
3、最终输出结果如图:台风中心高温高湿

极地冷高压

位温脚本搞的可以

cch513

愿意分享的人都是好人！

九】

楼主是个大好人啊

ww_wen

谢谢楼主分享！

问题少女

楼主您好，十分感谢您的分享，每个公式都有出处，辛苦楼主了。
看完楼主的帖子，我有如下困惑：
“*计算假相当位温（不知道位温时）(中华人民共和国气象行业标准---南海夏季风监测指标P10)
'thse=(tc+273.15)*exp(0.28586*log(1000.0/lev)+2500.0*qa/(338.52-0.24*(tc+273)+1.24*td))'
*计算假相当位温（已知位温时）
*'thse=theta*exp(2500.0*qa/(338.52-0.24*(tc+273)+1.24*td))'
”
这一部分，我看了那个文档，其中是这样说明的：
C:\Users\Administrator\Desktop\QQ图片20161229122228.png
T和Td的单位是℃，楼主的公式里边，分母下面(tc+273)单位是K，td单位是℃，应该不能相加吧？

问题少女

刚刚的图片没有上传，文档中的说明如下：
其中T为850hPa温度，单位：℃；q为850hPa比湿,单位: 克/克；Td为露点温度，

phot839

这个脚本很实用

初翠翠

感谢楼主的分享

angelaPP

谢谢楼主！