[ncl-talk] Ice/liquid water path

[Dennis Shea]

Well, NCL has several standard integration functions:

http://www.ncl.ucar.edu/Document/Functions/Built-in/simpne.shtml
http://www.ncl.ucar.edu/Document/Functions/Built-in/simpeq.shtml
http://www.ncl.ucar.edu/Document/Functions/Built-in/ftcurvi.shtml

===
However, I suggest weighting the values by the appropriate layer
thicknesses.

    vertical_average  = SUM[q*dp]/SUM[dp]
    vertical_integral   = SUM[q*dp]

     LWP = SUM[q*dp]/gravity    ; g/m^2

===

'dp' can be calculated via:

http://www.ncl.ucar.edu/Document/Functions/Built-in/dpres_hybrid_ccm.shtml

http://www.ncl.ucar.edu/Document/Functions/Built-in/dpres_plevel.shtml

===

The 6.4.0 NCL has


http://www.ncl.ucar.edu/Document/Functions/Contributed/wgt_vertical_n.shtml

The *beta* 6.4.0 can be downloaded from:

       http://www.cgd.ucar.edu/~shea/contributed.ncl_beta_640

To use it:

   load "./contributed.ncl_beta_640

or

   load "/where_ever/you_put_it/contributed.ncl_beta_640
===
Another, definition of LWP is here
  https://climatedataguide.ucar.edu/climate-data/liquid-water-path-overview

Good Luck

On Sat, Jul 9, 2016 at 4:50 PM, Marston Johnston <shejo284 at gmail.com> wrote:

> Hi,
>
> 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.
>
> /M
>
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160710/d7dc73c3/attachment.html