<div dir="ltr"><div>Change:  <br></div><div>   s1   = static_stability (pc, temp, 0, 0)</div>To:<br><div><div>   s1  = static_stability (pc, temp, 1, 0)</div><div><br></div><div>Or, replace the 'conform' and the above 's1' use:<br></div><div><br></div><div>   s1  = static_stability (p , temp, 1, 0) <br></div><div><br></div><div>=====</div><div>%> ncl heatsource.ncl</div><div><br></div><div><br></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Sep 6, 2019 at 6:25 AM Lyndz via ncl-talk <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr" class="gmail-m_4438535427870773800m_-3431812809093330190gmail_signature"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Dear NCL experts,</div><div><br></div><div><b>[Details]</b></div><div>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 (<a href="https://renqlsysu.github.io/2019/02/01/apparent_heat_source/" target="_blank">https://renqlsysu.github.io/2019/02/01/apparent_heat_source/</a>)</div><div><br></div><div>I have uploaded my files in my dropbox for verification:</div><div><br></div><div><a href="https://www.dropbox.com/sh/jx6gr9huan5v7lp/AADDr0pULF3cmMgP7LMVQVCma?dl=0" target="_blank">https://www.dropbox.com/sh/jx6gr9huan5v7lp/AADDr0pULF3cmMgP7LMVQVCma?dl=0</a><br></div><div><br></div><div><br></div><div><b>[Problem]</b></div><div>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 <b><i>-Inf</i></b> for all grid points.</div><div><br></div><div>I would like to ask for suggestions on how to compute this correctly in NCL.</div><div><br></div><div>Sincerely,</div><div><br></div><div>Lyndz</div><div><br></div></div></div></div></div></div></div></div></div></div></div>
_______________________________________________<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></blockquote></div>