整层水汽通量散度量级问题

左手疆环

诚心求助各位大师~~：
      我最近在论坛里看了很多关于求解水汽通量的问题，自己借鉴了前辈们的程序，但是计算出来的量级还是有些问题，麻烦高手们指点一下。
     我用的是fnl1*1°资料。
     首先我计算了比湿，利用的公式是es=(6.112*exp(17.67*(tmpprs.1-273.15)/(tmpprs.1-29.65)))，qs=(0.62197*es/(prs-0.378*es))， q=rhprs.2*qs/100.0。
     然后计算了单层水汽qu=q*ugrdprs.2/9.8*1000,qv=q*vgrdprs.2/9.8*1000，单位是g/(s*cm*hpa)。
    最后我计算了整层的水汽通量和散度。问题就是vint这个函数。我定义了quall=vint(pressfc.2(z=1)/100,qu.1,100)*9.8/1000，原因是我想得到整层水汽通量的量级是kg/(m*s)，*9.8是因为计算单层的qu已经除过9.8了，vint本身又除了一遍所以乘回来。/1000是因为本来qu单位是g/(s*cm*hpa)，积分后转化为kg/(m*s)需要除以10，但是如果vint单除以10的话我计算得到的量级特别大，超过1e2(我查阅文献整层水汽一般是1e2 kg/(m*s)的量级)。后来我再翻阅了关于vint函数，发现它的计算涉及f/g，f=100，所以我就认为vint函数又多乘了个100，因此我就把原本的/10改为了/1000。这样量级才吻合。
     可是计算整层水汽散度的时候问题就来了，我直接用的 hdivg(quall,qvall) ，但画图发现量级是 1e-4 kg/(m**2*s)，比一般文献里的 1e-5 kg/(m**2*s) 大了1个量级。我想了很久也没想明白，有可能就是vint函数那里搞错了，可是如果不那样的话又不知道怎样做，希望各位大师耐心看完我的问题，能帮我解答一下哈，谢谢！！！

river

我觉得吧vint积分是一个内部过程，和你之前除没除9.8没什么关系吧，还有就是从g转换到kg就是除以1000，这是肯定的，积分还能让国际标准的g和kg的换算关系发生变化呢。所以你就老老实实积分然后除以1000看看怎么样

左手疆环

river 发表于 2013-11-2 13:05

                               
登录/注册后可看大图

我觉得吧vint积分是一个内部过程，和你之前除没除9.8没什么关系吧，还有就是从g转换到kg就是除以1000，这是 ...

恩  号 我试试哈  谢了哈~~~~

yuejun0026

vint的函数计算里面它用的公式就是100/g*sum(expr * Delta(level))，所以在这个vint里面就已经除以了g了

yuejun0026

你可以看看这个上说的http://www.iges.org/grads/gadoc/gradfuncvint.html

道狼52203

{:eb315:}

leeol

yuejun0026 发表于 2013-12-3 20:31
vint的函数计算里面它用的公式就是100/g*sum(expr * Delta(level))，所以在这个vint里面就已经除以了g了

那请问如LZ情况，使用vint计算整层水汽通量时，vint(pressfc.2(z=1)/100,qu.1,100)*9.8，这个后面的*9.8就不要了？

yuejun0026

leeol 发表于 2014-3-25 20:31
那请问如LZ情况，使用vint计算整层水汽通量时，vint(pressfc.2(z=1)/100,qu.1,100)*9.8，这个后面的*9.8 ...

是的   公式是这样给的

leeol

yuejun0026 发表于 2014-3-27 17:10
是的   公式是这样给的

那还是要请教，LZ在之前计算单层水汽qu=q*ugrdprs.2/9.8*1000,qv=q*vgrdprs.2/9.8*1000，除以了一个9.8，vint后不需要乘回来？

yuejun0026

leeol 发表于 2014-3-28 10:13
那还是要请教，LZ在之前计算单层水汽qu=q*ugrdprs.2/9.8*1000,qv=q*vgrdprs.2/9.8*1000，除以了一个9.8， ...

额  这个是楼主弄的程序  我也不清楚  反正vint很难用   我一般都不用这个