假相当位温脚本

virl · 发表于 2015-3-15 22:16:50

谢谢

{:eb502:}

xiaocaoqiqiao · 发表于 2015-6-9 21:52:19

Jaxon 发表于 2014-12-15 18:05
问下楼主这是什么问题呢？资料分辨率应该不影响吧？

我也出现这样的问题应该是t0或者e出现负值，导致log(t0)或log(e)不能求对数，请问你最后解决了吗？

虾米go · 发表于 2015-7-20 21:06:25

也正在找公式的出处

虾米go · 发表于 2015-7-20 21:26:31

雪お轩漠发表于 2015-3-15 10:50
根据北大版《大气物理学》P20
我觉得程序22行的
'define e=prs*qq/(0.62197+qq)+1e-10'

是啊，我也感觉不应该加吧

四叶草 · 发表于 2015-8-18 16:26:13

楼主，谢谢你的分享！但其中有个公式本人感觉有点疑惑。如下：

'definee=prs*qq/(0.62197+qq)+1e-10' （1）

见《大气物理学》P20，e和r的公式自己推导一下就出来了

本人看了《大气物理课本》p20，发现楼主定义变量（1）应该是用下图红框内的公式，公式内的r表示的是混合比，而楼主定义变量中的qq应该是比湿，本人认为套用的公式应该是下图中公式（2.2.5），水汽压与比湿的关系应该是：e=q*p/(0.622+0.378*q)

《东亚副热带季风雨带建立特征及其降水性质分析》文中里的假相当位温计算公式中，比湿恰恰用r表示，两者间会不会搞混了。如下图：

四叶草 · 发表于 2015-8-18 17:00:05

四叶草发表于 2015-8-18 16:26
楼主，谢谢你的分享！但其中有个公式本人感觉有点疑惑。如下：
'definee=prs*qq/(0.62197+qq)+1e-10' （1 ...

疑惑解决了！由于通常大气中的混合比和比湿都小于0.04，因此可以认为q≈r，故楼主的计算方法可取。下面是本人利用'define e=prs*qq/(0.62197+qq)+1e-10'（1）和 'define e=q*p/(0.622+0.378*q)' （2）出的图，图像显示几乎一样。

（1）计算的结果

（2）计算的结果

dzn · 发表于 2015-9-22 14:23:49

不错的帖子，

银酱赛高 · 发表于 2015-10-16 15:59:13

楼主，请教一下，如果我要画假相当位温的垂直剖面图，需要用到不同的es计算公式吗？（因为es经验公式大多有条件限制，不知道在不同高度层次是否可用同一条公式？）

chengongyan · 发表于 2015-11-7 12:10:06

我用楼主的帖子已经画出来了，跟MICAPS的值差不多，谢谢楼主大公无私的分享

sunflower · 发表于 2016-1-22 09:10:32

楼主给力，感谢谢分享

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

[分享资料] 假相当位温脚本