<div dir="ltr"><div><div><div><div><div><div><div><div><div><div><div><div><div>Lyndon,<br></div>Absolutely. I calclulate MFC frequently. Use cdiff. But before I give the full explanation, I should note that your formula already gives the full MFC field. hdivg(q*u,q*v) should include the advection term. You can verify this by comparing that formula to those below that I will give you.<br><br></div>The trick with this is to define your grid spacing before doing anything else. While your model integration may use a fixed grid spacing (whether in degrees or linear distance or otherwise), if you use a control file to display the data, and if it includes XDEF and YDEF entries, then the grid on which your data are calculated may differ from that on which the model was integrated, so you can't use your scalar grid spacing value in that case. You can define your data grid spacing using trigonometry. The equations to calculate that are<br><br></div>dx = Re * cos(lat) * cdiff(lon,x)<br></div>dy = Re * cdiff(lat,y)<br><br></div>In the above equations, Re is the radius of the earth, and you need to make sure to convert between radians and degrees. Grads functions assume inputs are in radians and also output in radians. However, the lat and lon fields are in degrees, so convert. Both formulas are also instances of the arc length formula for a circle, which assumes the angle is in radians. So basically, convert everything to radians in those equations.<br><br></div>Once you have those terms, you have all the fields you need (assuming you have mixing ratio and wind components defined on the same grid already):<br><br></div>advection term = -u * cdiff(q,x) / dx - v * cdiff(q,y) / dy<br></div>convergence term = q*hdivg(u,v)<br><br></div>To verify that the advection term is already contained in hdivg(q*u,q*v), it should suffice to display <br></div>hdivg(q*u,q*v) - q * hdivg(u,v)<br></div>This field should match the advection term field above. You can also display <br>hdivg(q*u,q*v) - q * hdivg(u,v) - [ -u * cdiff(q,x) / dx - v * cdiff(q,y) / dy] and it should be 0.<br><br></div>I might be wrong about this since the divergence operator in spherical coordinates includes an additional term to represent convergence for purely meridional wind, so you may want to crack open a dynamics textbook and find what that term is and include it in here somewhere. Either way, I hope this helps.<br><br></div>Jeff Duda<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Jan 27, 2017 at 11:35 PM, Lyndon Mark Olaguera <span dir="ltr"><<a href="mailto:olagueralyndonmark429@gmail.com" target="_blank">olagueralyndonmark429@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">Dear All,<div><br></div><div>I'm trying to calculate the moisture flux convergence using the following commands:</div><div><br></div><div><div>*At 850 mb </div><div><br></div><div>'set dfile 1'</div><div>'define q=shum(z=3)'</div><div><br></div><div>'set dfile 2'</div><div>'define u=uwnd(z=3)'</div><div><br></div><div>'set dfile 3'</div><div>'define v=vwnd(z=3)'</div><div><br></div><div>'color -0.2 0.2 0.02 -kind blue->white->red'<br></div><div>'define mconv=(1)*hdivg(u*q,v*q)*1e6'</div><div><br></div><div>According to this link, I have to add the advection term.</div><div><br></div><div><a href="http://www.spc.noaa.gov/publications/banacos/mfc-sls.pdf" target="_blank">http://www.spc.noaa.gov/<wbr>publications/banacos/mfc-sls.<wbr>pdf</a></div><div><br></div><div><img src="cid:ii_159e390675f5b588" alt="Inline image 1" style="margin-right:0px" height="134" width="406"><br></div><div><br></div><div>Question:</div><div><br></div><div>1. Is there a way to do this in grads?</div><div>2. Or is there a function in GRADS that can calculate partial derivatives?</div><div><br></div><div>I'll appreciate any help.</div><div><br></div><div>Best,</div><div><br></div><div><div class="m_274993068877275854gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><font face="comic sans ms, sans-serif"><b><u>Lyndon Mark P. Olaguera</u></b></font><div><br></div></div></div></div></div></div></div></div></div></div>
</div></div>
<br>______________________________<wbr>_________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/<wbr>listinfo/gradsusr</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Jeff Duda<br>Post-doctoral research fellow<br>University of Oklahoma School of Meteorology<br></div></div></div></div></div>
</div>