Hi<div><br><div><div> The model estimate rainfall is come from WRF model ouput. But, I need to calculate the vertically integrated moisture convergence from model output u,v,and q.</div></div><div>Do you have a sample script for calculating the vertically integrated moisture convergence?</div>
<div>How do you add gravity and density term to the equation? How should I modified the script? Very appreciated for your input.</div><div><br></div><div>Chuan-Chi<br><br><div class="gmail_quote">2012/7/29 <span dir="ltr"><<a href="mailto:gradsusr-request@gradsusr.org" target="_blank">gradsusr-request@gradsusr.org</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Send gradsusr mailing list submissions to<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
or, via email, send a message with subject or body 'help' to<br>
<a href="mailto:gradsusr-request@gradsusr.org">gradsusr-request@gradsusr.org</a><br>
<br>
You can reach the person managing the list at<br>
<a href="mailto:gradsusr-owner@gradsusr.org">gradsusr-owner@gradsusr.org</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of gradsusr digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
1. Maxloc (Stefan Gofferje)<br>
2. Re: gradsusr Digest, Vol 29, Issue 55 (Chuan-Chi Tu)<br>
3. Re: gradsusr Digest, Vol 29, Issue 55 (Jeff Duda)<br>
4. Time interpolation (Stefan Gofferje)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Sat, 28 Jul 2012 19:37:20 +0300<br>
From: Stefan Gofferje <<a href="mailto:stefan@saakeskus.fi">stefan@saakeskus.fi</a>><br>
Subject: [gradsusr] Maxloc<br>
To: <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
Message-ID: <<a href="mailto:50141540.4080804@saakeskus.fi">50141540.4080804@saakeskus.fi</a>><br>
Content-Type: text/plain; charset=ISO-8859-1<br>
<br>
Hi,<br>
<br>
I was digging through the docs and found maxloc. From the description,<br>
it seems ideal from drawing "H" and "L" into a pressure chart but I<br>
somehow can't figure out how this works.<br>
Are there any examples anywhere?<br>
<br>
--Stefan<br>
<br>
--<br>
Suomen S??keskus - Hyv? s?? on makuasia<br>
<br>
Stefan Gofferje<br>
Haukantie 2 B 5<br>
FIN-37600 Valkeakoski<br>
GSM: <a href="tel:%2B358%20%2841%29%207290730" value="+358417290730">+358 (41) 7290730</a><br>
<br>
<a href="http://www.saakeskus.fi/" target="_blank">http://www.saakeskus.fi/</a><br>
<a href="mailto:stefan@saakeskus.fi">stefan@saakeskus.fi</a><br>
<br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Sat, 28 Jul 2012 21:30:17 -0600<br>
From: Chuan-Chi Tu <<a href="mailto:beth0509@gmail.com">beth0509@gmail.com</a>><br>
Subject: Re: [gradsusr] gradsusr Digest, Vol 29, Issue 55<br>
To: <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
Message-ID:<br>
<<a href="mailto:CAHkS%2BLaLwgF7w_2YsE6tavpbnC3iNzX8kmpOE%2BkQxY1icp-3cg@mail.gmail.com">CAHkS+LaLwgF7w_2YsE6tavpbnC3iNzX8kmpOE+kQxY1icp-3cg@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Hi,<br>
I use u, v in m/s.<br>
I use the vint function. I has the gravity term in it.<br>
<a href="http://www.iges.org/grads/gadoc/gradfuncvint.html" target="_blank">http://www.iges.org/grads/gadoc/gradfuncvint.html</a><br>
<br>
I am not sure if I add gravity term again. How do I add gravity and density<br>
term to the equation.<br>
<br>
The large value I get is about 120 mm/hr. The model simulated rain rate is<br>
50-90 mm/hr. These two terms<br>
should be close. Thanks for the help.<br>
<br>
Chuan-Chi<br>
<br>
2012/7/28 <<a href="mailto:gradsusr-request@gradsusr.org">gradsusr-request@gradsusr.org</a>><br>
<br>
> Send gradsusr mailing list submissions to<br>
> <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
><br>
> To subscribe or unsubscribe via the World Wide Web, visit<br>
> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
> or, via email, send a message with subject or body 'help' to<br>
> <a href="mailto:gradsusr-request@gradsusr.org">gradsusr-request@gradsusr.org</a><br>
><br>
> You can reach the person managing the list at<br>
> <a href="mailto:gradsusr-owner@gradsusr.org">gradsusr-owner@gradsusr.org</a><br>
><br>
> When replying, please edit your Subject line so it is more specific<br>
> than "Re: Contents of gradsusr digest..."<br>
><br>
><br>
> Today's Topics:<br>
><br>
> 1. Re: calculating vertically integrated moisture convergence<br>
> (Jeff Duda)<br>
><br>
><br>
> ----------------------------------------------------------------------<br>
><br>
> Message: 1<br>
> Date: Fri, 27 Jul 2012 11:59:44 -0500<br>
> From: Jeff Duda <<a href="mailto:jeffduda319@gmail.com">jeffduda319@gmail.com</a>><br>
> Subject: Re: [gradsusr] calculating vertically integrated moisture<br>
> convergence<br>
> To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>><br>
> Message-ID:<br>
> <<br>
> <a href="mailto:CAAig09DuYno-BOwJYY7Ss1zfOO-f76fweTQr3HFywFueiF%2B5jQ@mail.gmail.com">CAAig09DuYno-BOwJYY7Ss1zfOO-f76fweTQr3HFywFueiF+5jQ@mail.gmail.com</a>><br>
> Content-Type: text/plain; charset="iso-8859-1"<br>
><br>
> What units are your u and v in? If they're not in not m/s, that may be the<br>
> cause. Also, it looks more like you are attempting to compute precipitable<br>
> water. If so, I think you are missing the density of liquid water and the<br>
> gravitational constant term.<br>
><br>
> Jeff Duda<br>
><br>
> On Fri, Jul 27, 2012 at 10:08 AM, Chuan-Chi Tu <<a href="mailto:beth0509@gmail.com">beth0509@gmail.com</a>> wrote:<br>
><br>
> > Hi<br>
> ><br>
> > Can anyone help on calculating vertically integrated moisture<br>
> convergence?<br>
> ><br>
> > I use following script but the value looks too high<br>
> ><br>
> > q: specific humidity kg/kg<br>
> > PSFC:hPa<br>
> ><br>
> > set lev 1000 300<br>
> > define qu=q*u<br>
> > define qv=q*v<br>
> > define conv=-hdivg(qu,qv)<br>
> ><br>
> > set lev 1000<br>
> > d vint(PSFC,conv,300)*60*60<br>
> > *(mm/hr)<br>
> ><br>
> > Can anyone help on correcting it? Thank you for your help.<br>
> ><br>
> > _______________________________________________<br>
> > gradsusr mailing list<br>
> > <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
> > <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
> ><br>
> ><br>
><br>
><br>
> --<br>
> Jeff Duda<br>
> Graduate research assistant<br>
> University of Oklahoma School of Meteorology<br>
> Center for Analysis and Prediction of Storms<br>
> -------------- next part --------------<br>
> An HTML attachment was scrubbed...<br>
> URL:<br>
> <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20120727/c6d2110f/attachment-0001.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20120727/c6d2110f/attachment-0001.html</a><br>
><br>
> ------------------------------<br>
><br>
> _______________________________________________<br>
> gradsusr mailing list<br>
> <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
><br>
><br>
> End of gradsusr Digest, Vol 29, Issue 55<br>
> ****************************************<br>
><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20120728/b4329fd8/attachment-0001.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20120728/b4329fd8/attachment-0001.html</a><br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Sat, 28 Jul 2012 23:53:36 -0500<br>
From: Jeff Duda <<a href="mailto:jeffduda319@gmail.com">jeffduda319@gmail.com</a>><br>
Subject: Re: [gradsusr] gradsusr Digest, Vol 29, Issue 55<br>
To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>><br>
Message-ID:<br>
<CAAig09AHnEGK=<a href="mailto:bJbCK1Jndx2hD08L-nG2fkFDH-yaj-xiYyKRQ@mail.gmail.com">bJbCK1Jndx2hD08L-nG2fkFDH-yaj-xiYyKRQ@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
I'm not familiar with using that method to estimate rain rate. Sorry, I<br>
cannot help you further.<br>
<br>
On Sat, Jul 28, 2012 at 10:30 PM, Chuan-Chi Tu <<a href="mailto:beth0509@gmail.com">beth0509@gmail.com</a>> wrote:<br>
<br>
> Hi,<br>
> I use u, v in m/s.<br>
> I use the vint function. I has the gravity term in it.<br>
> <a href="http://www.iges.org/grads/gadoc/gradfuncvint.html" target="_blank">http://www.iges.org/grads/gadoc/gradfuncvint.html</a><br>
><br>
> I am not sure if I add gravity term again. How do I add gravity and<br>
> density term to the equation.<br>
><br>
> The large value I get is about 120 mm/hr. The model simulated rain rate is<br>
> 50-90 mm/hr. These two terms<br>
> should be close. Thanks for the help.<br>
><br>
> Chuan-Chi<br>
><br>
> 2012/7/28 <<a href="mailto:gradsusr-request@gradsusr.org">gradsusr-request@gradsusr.org</a>><br>
><br>
>> Send gradsusr mailing list submissions to<br>
>> <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
>><br>
>> To subscribe or unsubscribe via the World Wide Web, visit<br>
>> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
>> or, via email, send a message with subject or body 'help' to<br>
>> <a href="mailto:gradsusr-request@gradsusr.org">gradsusr-request@gradsusr.org</a><br>
>><br>
>> You can reach the person managing the list at<br>
>> <a href="mailto:gradsusr-owner@gradsusr.org">gradsusr-owner@gradsusr.org</a><br>
>><br>
>> When replying, please edit your Subject line so it is more specific<br>
>> than "Re: Contents of gradsusr digest..."<br>
>><br>
>><br>
>> Today's Topics:<br>
>><br>
>> 1. Re: calculating vertically integrated moisture convergence<br>
>> (Jeff Duda)<br>
>><br>
>><br>
>> ----------------------------------------------------------------------<br>
>><br>
>> Message: 1<br>
>> Date: Fri, 27 Jul 2012 11:59:44 -0500<br>
>> From: Jeff Duda <<a href="mailto:jeffduda319@gmail.com">jeffduda319@gmail.com</a>><br>
>> Subject: Re: [gradsusr] calculating vertically integrated moisture<br>
>> convergence<br>
>> To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>><br>
>> Message-ID:<br>
>> <<br>
>> <a href="mailto:CAAig09DuYno-BOwJYY7Ss1zfOO-f76fweTQr3HFywFueiF%2B5jQ@mail.gmail.com">CAAig09DuYno-BOwJYY7Ss1zfOO-f76fweTQr3HFywFueiF+5jQ@mail.gmail.com</a>><br>
>> Content-Type: text/plain; charset="iso-8859-1"<br>
>><br>
>> What units are your u and v in? If they're not in not m/s, that may be<br>
>> the<br>
>> cause. Also, it looks more like you are attempting to compute<br>
>> precipitable<br>
>> water. If so, I think you are missing the density of liquid water and the<br>
>> gravitational constant term.<br>
>><br>
>> Jeff Duda<br>
>><br>
>> On Fri, Jul 27, 2012 at 10:08 AM, Chuan-Chi Tu <<a href="mailto:beth0509@gmail.com">beth0509@gmail.com</a>><br>
>> wrote:<br>
>><br>
>> > Hi<br>
>> ><br>
>> > Can anyone help on calculating vertically integrated moisture<br>
>> convergence?<br>
>> ><br>
>> > I use following script but the value looks too high<br>
>> ><br>
>> > q: specific humidity kg/kg<br>
>> > PSFC:hPa<br>
>> ><br>
>> > set lev 1000 300<br>
>> > define qu=q*u<br>
>> > define qv=q*v<br>
>> > define conv=-hdivg(qu,qv)<br>
>> ><br>
>> > set lev 1000<br>
>> > d vint(PSFC,conv,300)*60*60<br>
>> > *(mm/hr)<br>
>> ><br>
>> > Can anyone help on correcting it? Thank you for your help.<br>
>> ><br>
>> > _______________________________________________<br>
>> > gradsusr mailing list<br>
>> > <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
>> > <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
>> ><br>
>> ><br>
>><br>
>><br>
>> --<br>
>> Jeff Duda<br>
>> Graduate research assistant<br>
>> University of Oklahoma School of Meteorology<br>
>> Center for Analysis and Prediction of Storms<br>
>> -------------- next part --------------<br>
>> An HTML attachment was scrubbed...<br>
>> URL:<br>
>> <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20120727/c6d2110f/attachment-0001.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20120727/c6d2110f/attachment-0001.html</a><br>
>><br>
>> ------------------------------<br>
>><br>
>> _______________________________________________<br>
>> gradsusr mailing list<br>
>> <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
>> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
>><br>
>><br>
>> End of gradsusr Digest, Vol 29, Issue 55<br>
>> ****************************************<br>
>><br>
><br>
><br>
> _______________________________________________<br>
> gradsusr mailing list<br>
> <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
><br>
><br>
<br>
<br>
--<br>
Jeff Duda<br>
Graduate research assistant<br>
University of Oklahoma School of Meteorology<br>
Center for Analysis and Prediction of Storms<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20120728/3fff224e/attachment-0001.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20120728/3fff224e/attachment-0001.html</a><br>
<br>
------------------------------<br>
<br>
Message: 4<br>
Date: Sun, 29 Jul 2012 15:37:26 +0300<br>
From: Stefan Gofferje <<a href="mailto:stefan@saakeskus.fi">stefan@saakeskus.fi</a>><br>
Subject: [gradsusr] Time interpolation<br>
To: <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
Message-ID: <<a href="mailto:50152E86.8040008@saakeskus.fi">50152E86.8040008@saakeskus.fi</a>><br>
Content-Type: text/plain; charset=ISO-8859-1<br>
<br>
Hi,<br>
<br>
is it possible to do time interpolation in map plots with grads?<br>
Obviously, when I use gxout line, some interpolation is done.<br>
I'd like to create movies/animations from gfs data but I'd like them to<br>
be a bit smoother than 3h-steps, let's say, 1h steps or even 30min-steps.<br>
<br>
How would I do this with grads?<br>
<br>
--Stefan<br>
<br>
--<br>
Suomen S??keskus - Hyv? s?? on makuasia<br>
<br>
Stefan Gofferje<br>
Haukantie 2 B 5<br>
FIN-37600 Valkeakoski<br>
GSM: <a href="tel:%2B358%20%2841%29%207290730" value="+358417290730">+358 (41) 7290730</a><br>
<br>
<a href="http://www.saakeskus.fi/" target="_blank">http://www.saakeskus.fi/</a><br>
<a href="mailto:stefan@saakeskus.fi">stefan@saakeskus.fi</a><br>
<br>
<br>
<br>
<br>
------------------------------<br>
<br>
_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br>
<br>
End of gradsusr Digest, Vol 29, Issue 56<br>
****************************************<br>
</blockquote></div><br></div></div>