<p dir="ltr">Hi,<br>
Thanks a lot :-)<br>
This is exactly what I am looking for.</p>
<p dir="ltr">Best,<br>
/M</p>
<div class="gmail_quote">On 11 Jul 2016 04:12, "Dennis Shea" <<a href="mailto:shea@ucar.edu">shea@ucar.edu</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div>Well, NCL has several standard integration functions:<br><br><a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/simpne.shtml" target="_blank">http://www.ncl.ucar.edu/Document/Functions/Built-in/simpne.shtml</a><br><a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/simpeq.shtml" target="_blank">http://www.ncl.ucar.edu/Document/Functions/Built-in/simpeq.shtml</a><br><a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/ftcurvi.shtml" target="_blank">http://www.ncl.ucar.edu/Document/Functions/Built-in/ftcurvi.shtml</a><br><br>===<br></div>However, I suggest<b> </b>weighting the values by the appropriate layer thicknesses.<br><b><br> </b>vertical_average = SUM[q*dp]/SUM[dp]<br> vertical_integral = SUM[q*dp]<br><br></div> LWP = SUM[q*dp]/gravity ; g/m^2<br><br>===<br><br></div><div>'dp' can be calculated via:<br></div><div> <a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/dpres_hybrid_ccm.shtml" target="_blank">http://www.ncl.ucar.edu/Document/Functions/Built-in/dpres_hybrid_ccm.shtml</a><br> <a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/dpres_plevel.shtml" target="_blank">http://www.ncl.ucar.edu/Document/Functions/Built-in/dpres_plevel.shtml</a><br><br>===<br><br></div><div>The 6.4.0 NCL has<br><br> <a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/wgt_vertical_n.shtml" target="_blank">http://www.ncl.ucar.edu/Document/Functions/Contributed/wgt_vertical_n.shtml</a><br><br></div><div>The *beta* 6.4.0 can be downloaded from:<br><br> <a href="http://www.cgd.ucar.edu/~shea/contributed.ncl_beta_640" target="_blank">http://www.cgd.ucar.edu/~shea/contributed.ncl_beta_640</a><br><br></div><div>To use it:<br><br></div><div> load "./contributed.ncl_beta_640<br><br></div><div>or<br><br> load "/where_ever/you_put_it/contributed.ncl_beta_640</div><div>===<br></div>Another, definition of LWP is here<br> <a href="https://climatedataguide.ucar.edu/climate-data/liquid-water-path-overview" target="_blank">https://climatedataguide.ucar.edu/climate-data/liquid-water-path-overview</a><br><div> <div>
<div><br></div><div>Good Luck<br></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Jul 9, 2016 at 4:50 PM, Marston Johnston <span dir="ltr"><<a href="mailto:shejo284@gmail.com" target="_blank">shejo284@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">Hi,<div><br></div><div>I'm wondering if NCL has a function for calculating the vertically integrated ice/liquid water path? I did this a long time ago in python using a trapezoid rule but I can't seem to find something similar in NCL. If anyone has a code to do this, I would appreciate some tips in how to do it in NCL.</div><span><font color="#888888"><div><br></div><div>/M</div><div><br></div><div><br>
</div></font></span></div>
<br>_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</blockquote></div>