ncl计算假相当位温和水汽通量

气象妹子 · 发表于 2015-11-2 10:38:03

学习学习，，最近正在算

cqlqlq · 发表于 2016-4-13 12:02:12

学习！

{:5_231:}

cqlqlq · 发表于 2016-4-13 12:02:20

{:5_201:}

simon7 · 发表于 2016-6-4 19:19:28

楼主以上为计算水汽通量的，水汽通量散度如何计算啊？ncl中散度公式如何编写，请教了

tureqiuyao · 发表于 2016-12-6 13:46:12

b1= where(tk .gt. 263.0, P-0.278*exp(17.26*(tk-273.16)/(tk-35.86)),\
P-0.278*exp(21.87*(tk-273.16)/(tk-7.66)))
请问此处为何取0.278，要计算饱和比湿，根据公式q=εe/p-0.378e,分母不是应该取0.378吗？求大神解答啊！

tcbyle · 发表于 2017-7-21 08:35:06

学习啦！

御风者 · 发表于 2017-10-25 21:05:15

楼主的这段程序存在问题
  a1 = where(tk .gt. 263.0, 0.622*6.11*exp(17.26*(tk-273.16)/(tk-35.86)), \
                  0.622*6.11*exp(21.87*(tk-273.16)/(tk-7.66))) ;
   b1= where(tk .gt. 263.0, P-0.278*exp(17.26*(tk-273.16)/(tk-35.86)),\
               P-0.278*exp(21.87*(tk-273.16)/(tk-7.66))) ;
   qs1=a1/b1
饱和比温公式为 qs = 0.622*es/p-0.378es，所以上面程序中b1中的0.278 应该换成 0.378*6.11=2.31。

ql455020268 · 发表于 2018-10-5 03:56:33

tureqiuyao 发表于 2016-12-6 13:46
b1= where(tk .gt. 263.0, P-0.278*exp(17.26*(tk-273.16)/(tk-35.86)),\
P-0.278* ...

我看了一下，楼主的a/b不应该是饱和混合比吗，应该是es/(p-es)吧，感觉楼主程序计算的是饱和比湿了

ql455020268 · 发表于 2018-10-5 03:56:36

tureqiuyao 发表于 2016-12-6 13:46
b1= where(tk .gt. 263.0, P-0.278*exp(17.26*(tk-273.16)/(tk-35.86)),\
P-0.278* ...

我看了一下，楼主的a/b不应该是饱和混合比吗，应该是es/(p-es)吧，感觉楼主程序计算的是饱和比湿了

ql455020268 · 发表于 2018-10-5 03:58:32

御风者发表于 2017-10-25 21:05
楼主的这段程序存在问题
a1 = where(tk .gt. 263.0, 0.622*6.11*exp(17.26*(tk-273.16)/(tk-35.86)), \
...

应该是饱和混合比吧，a1/b1,正确的是es/(p-es),p为气压，es为饱和水汽压

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