水汽通量中的比湿怎么求？

hxin0702

目前参考大家关于水汽通量gs文件。有下面几大疑惑，请高手解答啊.
准备文件NCEP的再分析日资料，air.2012.nc' ；rhum.2012.nc' ；uwnd.2012.nc' ；vwnd.2012.nc'；
（1）通过计算饱和水汽压es,饱和比湿qs，比湿q，风速V来计算水汽通量uvq
'define es=6.112*exp(17.2693822*(air-273.16)/(air-29.65))'              
'define qs=0.622*es/(lev-0.378*es)'
'define q=rhum.2*qs/100'
'define v=mag(uwnd.3,vwnd.4)'
'define uvq=v*q/9.8'
(2)第二种方法是通过计算饱和水汽压es,和下载的相对湿度rhum资料直接求比湿q，进而求水汽通量uvq
'define es=6.112*exp(17.2693822*(air-273.16)/(air-29.65))'  
'define q=rhum.2*es/100'
'define v=mag(uwnd.3,vwnd.4)'
'define uvq=v*q/(9.8*1000)'
以上两种方法主要分歧是比湿怎么求？我更偏向第一种，但是做出来的图有问题。还有，既然是用比湿资料，NCEP不是有比湿数据shum吗，怎么大家都不用呢？还请指教。

river

貌似一般都是先求qu qv 然后d mag(qu,qv)这样。你说的比湿的求法可以直接使用ncep的shum资料。你看一下这个帖子http://bbs.06climate.com/forum.php?mod=viewthread&tid=10552，就是这么算的，结果倒是没输出水汽通量，但是其他的出来的物理量好像没什么问题

youjianlaodao

有高人回答么？

blui

我也不知道怎么算水汽通量。。。我用的是再分析2资料，可是没有比湿。。

tiantian

a和b的值不能直接用的，温度在-15度以上按水面处理，在-40度以下按冰面处理，在二者之间的按冰水混合处理。

qixiangren34

楼主提出的这个问题请清风回答吧

wangxh0206

最近也正思考这个问题，求解！

beesfly

很纳闷，不太明白为何不直接用ncep1中的shum资料？

somnuszhouying

同样奇怪怎么不直接用nc中的比湿值呢？求解

赫小猫

nc资料中不是就有比湿的数据吗？