求助水汽通量散度的量级问题

lacie翼

话说 不是要除以g么

尽头的尽头

269984288 发表于 2011-8-9 09:23

                               
登录/注册后可看大图

好吧，我再想想，谢谢了哈

嗯，哈哈，要是问题解决了记得告诉大家，我们也学习下

269984288

lacie翼 发表于 2011-8-9 09:37

                               
登录/注册后可看大图

话说 不是要除以g么

话说不是已经除了g么？{:3_44:}

269984288

尽头的尽头 发表于 2011-8-9 11:56

                               
登录/注册后可看大图

嗯，哈哈，要是问题解决了记得告诉大家，我们也学习下

好的~~{:3_51:}

bluexlm

100要除，因为fnl资料里rh是几十的量级；1000要乘，因为比湿要求g/kg，最后量级是e-5没错，有些文献是这个量级，当然也有-6到-8的，可能他们错了。
另外比湿（饱和）的计算公式五花八门：
1.meteoinfo的
es = 6.112*DataMath.Exp(17.2693822*(air-273.16)/(air-29.65))
qs = 0.62197*es/(prs-0.378*es)
2.下面这个也常见，最大区别就是29.65改为35.86
'define es=6.11*exp(17.26*(tt-273.16)/(tt-35.86))'
'define qs=0.622*es/(lev-0.378*es)'
3.mofangbao的，写的mixr，其实就是比湿
#'tc=(tmpprs-273.16)'
'td=tc-((14.55+0.114*tc)*(1-0.01*rh) + pow((2.5+0.007*tc)*(1-0.01*rh),3) + (15.9+0.37*tc)*pow((1-0.01*rh),14))'
'vapr=6.112*exp((17.67*td)/(td+243.5))'
'e=vapr*1.001+(lev-100)/900*0.0034'
'define mixr=0.62137*(e/(lev-e))*1000'
好像还有，林林总总，还有更让人头晕的么。我也用fnl的2m比湿验证了一下，发现气象学书上这个公式可行，这个~~~实在有点晕
******验证比湿计算，气象学书上公式*******
*'define tc=tmp2m-273.16'
*饱和水汽压hpa
*'define es=6.11*pow(10,(7.45*tc/(235+tc)))'
*水汽压hpa
*'define e=es*rh2m/100'
*比湿g/g;kg/kg
*'define q=0.622*e*100/pressfc'
*'d q'
*'d spfh2m'

269984288

bluexlm 发表于 2011-8-24 13:01

                               
登录/注册后可看大图

100要除，因为fnl资料里rh是几十的量级；1000要乘，因为比湿要求g/kg，最后量级是e-5没错，有些文献是这个量 ...

哇哇~~~谢谢啦~~~热心人士啊~~~{:3_48:}{:3_62:}，还总结了那么多先，膜拜……

bluexlm

刚看了何立富的一篇文章，里面水汽通量散度是-6量级，还有‘数值预报产品的种类和含义

’也说-7~-6，这是为啥
另:
*水汽通量单位g/(cm*hpa*s)，矢量
*水汽通量散度单位g/(cm^2*hpa*s)
*整层水汽通量单位g/(cm*s)
*整层水汽通量散度单位g/(cm^2*s)

bluexlm

验证了一下，meteoinfo、mofangbao和气象学书上公式同fnl2m比湿差异不大，用35.86则比较大了

269984288

bluexlm 发表于 2011-8-25 19:20

                               
登录/注册后可看大图

刚看了何立富的一篇文章，里面水汽通量散度是-6量级，还有‘数值预报产品的种类和含义

’也说-7~-6，这 ...

{:3_61:}我也不知道啊~~~反正是按天原定义来算出来的……

bluexlm

这个是我的图