mconv.gs

[Rauscher Sara A]

Hi Jared,

I don't use the mconv.gs script, but this is what I use to obtain
vertically integrated moisture convergence:

set t 1
set lev 1000 300
set q = hus.1 *specific humidity in kg/kg
set u = ua.1  *u wind, m/s
set v = va.1  *v wind, m/s
define uq = u*q
define vq = v*q
define xconst=const(hus,1000,-a) *set bottom level to 1000 mb
set lev 1000
define vu1 = vint(xconst,uq,300)
define vv1 = vint(xconst,vq,300)
define div = hdivg(vu1,vv1)

I believe this gives units of kg/s; if you multiply by 86400 you should
get mm/day. I used specific humidity instead of mixing ratio, perhaps
this isn't correct, but you can easily convert from one to the other.

If I am wrong about the units or the method, I hope someone will correct
this.

good luck,
sara



Jared Bowden wrote:
> Could somebody confirm to me how to get vertically integrated horizontal
> moisture convergence using the mconv.gs script.  And what are the units of
> the final output.
>
> Thanks,
> Jared Bowden
> PhD. Candidate
> North Carolina State University
>

--

Dr. Sara A. Rauscher
Earth System Physics Section
The Abdus Salam International Centre for Theoretical Physics
Strada Costiera 11
34014 Trieste ITALY
ph: +39 040 2240 225 fax: +39 040 2240 449
email: srausche at ictp.it