<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->North because I use spherical harmonic functions. <br></div><div>Documentation is in the functions only. <br></div><div><div><br><pre>%> less $NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl <br><br>--------------------<br>undef("vertint")<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("advect_gradient")<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==>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"><<a href="mailto:littithomas@gmail.com" target="_blank">littithomas@gmail.com</a>></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>