<div dir="ltr"><div><div><div>   MFC = Moisture Flux Convergence<br><br>   MFC_advect = -( u*(dq/dx)+v*(dq/dy) )    ; advection term<br>   MFC_conv    = -q*( (du/dx)+  (dv/dy) )      ; con(div)-vergence<br><br>   MFC = MFC_advect + MFC_conv<br><br></div><div><br></div><div><br><br>-----------------------------<br><br></div>In shea_util.ncl there are several functions. These are *NOT* supported.<br></div>I put them there for people to look at and modify.  They require global grids<br></div><div>ordered South-&gt;North because I use spherical harmonic functions. <br></div><div>Documentation is in the functions only. <br></div><div><div><br><pre>%&gt; less $NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl <br><br>--------------------<br>undef(&quot;vertint&quot;)<br>function vertint (X:numeric, dp:numeric, con[1]:numeric \<br>                 ,long_name:string, units:string, iopt:integer)<br>; <br>; Perform vertical integration:  dim_sum_n(X*dp,lind)<br></pre>------------------------------<br>undef(&quot;advect_gradient&quot;)<br>function advect_gradient (u:numeric, v:numeric, x:numeric \<br>                         ,gridType[1]:integer, con[1]:numeric \<br>                         ,long_name:string, units:string, iopt[1]:integer)<br>;<br>; linear advection of a scalar quantity: UV.GRADIENT(X)<br>;                                      : U*(dX/dlon) + V*(dX/dlat)<br>;<br>; Requires:<br>;    (1) grid be global cuz spherical harmonics are used<br>;    (2) required dimension order ([time,[lev,]]lat,lon)<br>;    (3) the input grids *must* be ordered South==&gt;North<br><br></div><div>Divergence can be calculated via spherical harmonics<br><br></div><div>Good luck<br></div><div><br><div><br></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Dec 18, 2014 at 6:48 PM, tms_l <span dir="ltr">&lt;<a href="mailto:littithomas@gmail.com" target="_blank">littithomas@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div>hi,<br><br></div>does any one having any idea of calculating vertically integrated moisture flux convergence at different pressure levels , using ncl?<br><br></div>thanks in advance<br><br></div>thomas<br></div>
</blockquote></div><br></div>