[gradsusr] A GrADS scripts that calculates moisture convergence

[Rashed Mahmood]

Hi Allen

Does "lev" in the script means a fixed level value (e.g., 850, 700, 500
etc)?

Thanks
Rashed

On Fri, Feb 11, 2011 at 9:53 PM, Teddy Allen <teddyallen at yahoo.com> wrote:

>   * This is a script for displaying moisture convergence
>
> * Written by Michael Maxwell
>
> *
>
> * rh = Relative Humidity in %
>
> * t = Temp at *set level in degrees Kelvin
>
> * tc = Temp in degrees C
>
> * td = Dewpoint at *set level in degrees C
>
> * e = Vapor pressure
>
> * mixr = Mixing ratio
>
> * u = U-wind in m/s
>
> * v = V-wind in m/s
>
> * mconv = moisture convergence/divergence. convergence is positive and
> divergence is negative.
>
> 'tc = (t-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.117*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'
>
> 'mixr = 0.62197*(e/(lev-e))*1000'
>
> 'mconv = (-1)*hdivg(u*mixr,v*mixr)*1e4'
>
> 'd mconv'
>
>
>
> http://www.teddyallen.com
>
> "To move quickly, go alone. To move far, go *together*."
>
> African Proverb
>
>
>  ------------------------------
> *From:* saiguran loisulie <saiguranloi at yahoo.com>
> *To:* GrADS Users <gradsusr at gradsusr.org>
> *Sent:* Fri, February 11, 2011 8:40:51 AM
> *Subject:* [gradsusr] A GrADS scripts that calculates moisture convergence
>
>   Hi,
>
> I wanted to download a script that calculates moisture convergence from
> ftp://grads.iges.org/grads/scripts/mconv.gs; but I failed due to access
> restrictions imposed here.
>
> Can someone email me this script or the similar or better one? I need it to
> analyse my data to study moisture dynamics.
>
> Thank you in advance for your cordial help.
>
> Saiguran Loisulie
> Climate Systems Analysis Group
> Environmental and Geographical Science Department
> University of Cape Town
> Private Bag X3
> Rondebosch
> 7701
> SOUTH AFRICA
>
> CELL 0748 068 230
>
> ------------------------------
> It's here! Your new message!
> Get new email alerts<http://us.rd.yahoo.com/evt=49938/*http://tools.search.yahoo.com/toolbar/features/mail/>with the free Yahoo!
> Toolbar.<http://us.rd.yahoo.com/evt=49938/*http://tools.search.yahoo.com/toolbar/features/mail/>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110226/08615b4d/attachment-0003.html