<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style>
</head>
<body class='hmmessage'><div dir='ltr'>
Arlindo&nbsp;thanks for your suggestions. I'll try them.<div>Hyacinth</div><div><br><div>&gt; <br>&gt; Message: 5<br>&gt; Date: Sat, 17 Sep 2011 10:12:57 -0400<br>&gt; From: Arlindo da Silva &lt;dasilva@alum.mit.edu&gt;<br>&gt; Subject: Re: [gradsusr] Moisture flux convergence/divergence<br>&gt; To: GrADS Users Forum &lt;gradsusr@gradsusr.org&gt;<br>&gt; Message-ID:<br>&gt;         &lt;CALz7xZf5MaHEvP1bX1Z93ysYNg5xVCGOvBdBaM5aHTkEs6+HSw@mail.gmail.com&gt;<br>&gt; Content-Type: text/plain; charset="iso-8859-1"<br>&gt; <br>&gt; 2011/9/16 Hyacinth Nnamchi &lt;hyacinth.1@hotmail.com&gt;<br>&gt; <br>&gt; &gt;  Jeff,<br>&gt; &gt;<br>&gt; &gt; I seems my question wasn't well framed, and so you dealt with the typo. I<br>&gt; &gt; have attached two Figs showing what the calculations below will give and<br>&gt; &gt; what exactly I'll like to do.<br>&gt; &gt;<br>&gt; &gt;<br>&gt; &gt;From the attached plot it appears that you want to plot the divergent<br>&gt; component of the mass flux, as given by the Helmholtz decomposition:<br>&gt; <br>&gt;            http://en.wikipedia.org/wiki/Helmholtz_decomposition<br>&gt; <br>&gt; The first step is to compute the velocity potential<br>&gt; <br>&gt;    ga-&gt; chi = fish_chi(zon,mer)<br>&gt; <br>&gt; (As Roger mentions, make sure to average in the vertical first if this is<br>&gt; what you want.) Next, the divergent component of the mass flux is related to<br>&gt; the gradient of chi; see this for more information:<br>&gt; <br>&gt;  Calculating the Divergent<br>&gt; Wind&lt;http://opengrads.org/doc/udxt/fish/fish.html#calculating_the_divergent_wind&gt;<br>&gt; <br>&gt; Notice that you will need an opengrads build in order to use the fish<br>&gt; extension; more info here:<br>&gt; <br>&gt;     http://opengrads.org/wiki/index.php?title=Installing_the_OpenGrADS_Bundle<br>&gt; <br>&gt;  Good luck,<br>&gt; <br>&gt;    Arlindo<br>&gt; <br>&gt; <br>&gt; <br>&gt; <br>&gt; <br>&gt; <br>&gt; &gt; Thanks,<br>&gt; &gt;<br>&gt; &gt; Hyacinth<br>&gt; &gt;<br>&gt; &gt;<br>&gt; &gt; &gt; Date: Thu, 15 S
ep 2011 17:34:50 -0500&gt;<br>&gt; &gt; From: Jeffrey Duda &lt;jdduda@iastate.edu&gt;&gt;<br>&gt; &gt;  Subject: Re: [gradsusr] Moisture flux convergence/divergence<br>&gt; &gt; &gt; To: GrADS Users Forum &lt;gradsusr@gradsusr.org&gt;<br>&gt; &gt; &gt; Message-ID:<br>&gt; &gt; &gt; &lt;CAP9LQ-VNkxcu-+ey_QkZU7f4MHf+VwJDvav=Dci3ns+x+pjG8Q@mail.gmail.com&gt;<br>&gt; &gt; &gt; Content-Type: text/plain; charset="iso-8859-1"<br>&gt; &gt;<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; Set the graphics output to vector ('set gxout vector'), then 'd zon;mer'<br>&gt; &gt; &gt; (note the semicolon instead of a comma).<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; Jeff Duda<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; 2011/9/15 Hyacinth Nnamchi &lt;hyacinth.1@hotmail.com&gt;<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; &gt; Hello,<br>&gt; &gt; &gt; &gt;<br>&gt; &gt; &gt; &gt; This is not directly about GrADS but since most of address similar!<br>&gt; &gt; &gt; &gt; problems, I guess someone on this forum will be able to answer my<br>&gt; &gt; question.<br>&gt; &gt; &gt; &gt;<br>&gt; &gt; &gt; &gt; I understand that to calculate moisture transport, given zonal,<br>&gt; &gt; meridional<br>&gt; &gt; &gt; &gt; winds and specific humidity, can simply do:<br>&gt; &gt; &gt; &gt;<br>&gt; &gt; &gt; &gt; 'define zon=uwnd.1*shum.3'<br>&gt; &gt; &gt; &gt; 'define mer=vwnd.2*shum.3'<br>&gt; &gt; &gt; &gt; 'd zon,mer'<br>&gt; &gt; &gt; &gt;<br>&gt; &gt; &gt; &gt; But, I want to plot vectors indicating field of net moisture flux (all<br>&gt; &gt; arrows converging/diverging) instead of just transport.<br>&gt; &gt; &gt; &gt;<br>&gt; &gt; &gt; &gt; Thanks for any ideas.<br>&gt; &gt; &gt; &gt;<br>&gt; &gt; &gt; &gt; Hyacinth<br>&gt; &gt; &gt; &gt;<br>&gt; &gt; &gt; &gt;<br>&gt; &gt; &gt; &gt; __<br>&gt; &gt;<br>&gt; &gt; _______________________________________________<br>&gt; &gt; gradsusr mailing list<br>&gt; &gt; gradsusr@gradsusr.org<br>&gt; &gt; http://gradsusr.org/mailman/listinfo/gradsusr<br>&gt; &gt;<br>&gt; &gt;<br>&gt; <br>&gt;&gt; ***************************************
*<br></div></div>                                               </div></body>
</html>