[ncl-talk] Computing apparent heat source from reanalysis data

[Dennis Shea]

Change:
   s1   = static_stability (pc, temp, 0, 0)
To:
   s1  = static_stability (pc, temp, 1, 0)

Or, replace the 'conform' and the above 's1' use:

   s1  = static_stability (p , temp, 1, 0)

=====
%> ncl heatsource.ncl



On Fri, Sep 6, 2019 at 6:25 AM Lyndz via ncl-talk <ncl-talk at ucar.edu> wrote:

> Dear NCL experts,
>
> *[Details]*
> I am calculating the apparent heat source from NCEP R2 reanalysis data
> following the equation of Yanai et al. (1973) as in the attached image (
> https://renqlsysu.github.io/2019/02/01/apparent_heat_source/)
>
> I have uploaded my files in my dropbox for verification:
>
> https://www.dropbox.com/sh/jx6gr9huan5v7lp/AADDr0pULF3cmMgP7LMVQVCma?dl=0
>
>
> *[Problem]*
> Attached in this email is the script that I am currently working on. It
> runs fine but the values of the apparent heat source is *-Inf* for all
> grid points.
>
> I would like to ask for suggestions on how to compute this correctly in
> NCL.
>
> Sincerely,
>
> Lyndz
>
> _______________________________________________
> 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/20190908/868b1ada/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: heatsource.ncl
Type: application/octet-stream
Size: 3655 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190908/868b1ada/attachment.obj>