[ncl-talk] CAPE and CIN calculation

Bill Ladwig ladwig at ucar.edu
Wed Feb 17 12:10:52 MST 2016


The bug is due to incorrect units for PSFC, which should be in hPa.  So, in
the WRFUserARW.ncl script, for the line you mentioned above:

cape = wrf_cape_3d( P, tk, QV, z, HGT, PSFC, True )

should be:

cape = wrf_cape_3d( P, tk, QV, z, HGT, PSFC/100., True )

See Dave's reply for how to view the wrf_cape_3d C/Fortran source code.

Bill


On Wed, Feb 17, 2016 at 11:14 AM, Dan Chen <dchen at ucar.edu> wrote:

> Hi,
>
> I'm wondering if below email is successfully sent to ncl-talk list? If so,
> can anyone help?
>
> Or can you let me know how the CIN and CAPE are calculated in cape_3d
> function? I checked the  code in
> $NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl, and found the line
> "cape = wrf_cape_3d( P, tk, QV, z, HGT, PSFC, True )", but I can't find the
> code where *wrf_cape_3d* is located.  can you help to check?
>
> Many thanks,
> Dan
>
>
>
>
> On 2/12/16 5:27 PM, Dan Chen wrote:
>
> Hi,
>
> I found there is a bug fix of CAPE_3d calculation in 6.4.0. But I can't
> find the v6.4 by module load ncl on yellowstone. It seemed like the v6.3 is
> the latest.
>
> can you provide me the path of v6.4.0 on yellowstone and how to use it?
>
> Many thanks,
> Dan
>
>
>
>
> _______________________________________________
> 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/20160217/499701fc/attachment.html 


More information about the ncl-talk mailing list