ncl计算水汽平衡方程

柚柚子

最近在算水汽平衡方程moisture budget遇到了一些问题请教一下大家！
根据水汽平衡方程，降水和蒸发之差等于水汽通量散度的垂直积分：
其中，P为降水，E蒸发，g重力加速度，Pw水的密度，u水平风场，q比湿，p气压。
等式左边的单位均能换算成mm/s，因此等式右边的单位也应该是mm/s,但由气象家园里ncl计算整层水汽通量散度的经验贴[url=http://bbs.06climate.com/NCL计算单层和整层水汽通量、水汽通量散度的单位和量级%20http://bbs.06climate.com/forum.p ... mp;fromuid=49064%20(出处:%20气象家园)]NCL计算单层和整层水汽通量、水汽通量散度的单位和量级 http://bbs.06climate.com/forum.p ... 5&fromuid=49064 (出处: 气象家园)[/url]：vibeta(uv2dv_cfd(q*u/g))*10000000→10e-5kg/(m2*s) 在没有除以水的密度1000kg/m3的情况下水汽通量散度垂直积分的单位就已经是mm/s （kg/(m2*s)实际上是等于mm/s）。这是一个疑问。
根据方程自己推算了一下单位：

发现等式左右两边单位好像是不相等的？是不是我没算对呀？

用wrfout数据计算水汽平衡方程，蒸发E用的是潜热通量（单位W/m2）LH*0.4*(10^-6)转换为蒸发量（单位mm/s）,水汽通量散度根据上面气象家园的帖子计算1000hPa-150hPa垂直积分，等式两边差异很大。

有没有算过的大佬讲解一下？谢谢！