<div dir="ltr"><div><a href="https://www.ncl.ucar.edu/Applications/wind.shtml"><b>https://www.ncl.ucar.edu/Applications/wind.shtml</b></a></div><div><br></div><div>See bottom example.</div><div><p>
This example uses daily mean data from NOAA/OAR/ESRL PSD, Boulder, Colorado, USA.  Specifically, 
<a href="https://www.esrl.noaa.gov/psd/data/gridded/data.ncep.reanalysis2.html"><b>NCEP_Reanalysis 2</b></a>
daily mean data spanning 3-9 April 1993. 
</p><p>
This is a <b>test script</b>. It is not fully tested. The returned quantities are:
</p><pre>         q1 =   (dT/dt) - [omega*static_stability - V.grad(T)]  ; K/day
         q2 = -[(dH/dt) +  V.grad(H) +omega*(dH/dp)]            ; g/(kg-Pa)</pre></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Sep 13, 2019 at 4:28 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>Dear NCL-experts,</div><div><br></div><div>*This is related to my previous emails here.</div><div><br></div><div><b>[Details]</b></div><div><b><br></b></div><div>I am calculating the heat sources and sinks. My input files can be downloaded from these links:</div><div><br></div><div><b>uwnd:</b> <a href="https://www.dropbox.com/s/9esjzhe9q9bodow/lag-lead_uwnd.nc?dl=0" target="_blank">https://www.dropbox.com/s/9esjzhe9q9bodow/lag-lead_uwnd.nc?dl=0</a></div><div><b>vwnd:</b> <a href="https://www.dropbox.com/s/o4nf0vhv2u9c80v/lag-lead_vwnd.nc?dl=0" target="_blank">https://www.dropbox.com/s/o4nf0vhv2u9c80v/lag-lead_vwnd.nc?dl=0</a></div><div><b>omega:</b> <a href="https://www.dropbox.com/s/kgoe8sf7z7peq1a/lag-lead_omega.nc?dl=0" target="_blank">https://www.dropbox.com/s/kgoe8sf7z7peq1a/lag-lead_omega.nc?dl=0</a></div><div><b>mixing ratio: </b><a href="https://www.dropbox.com/s/u1w0hg8fvoecqhr/lag-lead_mixratio.nc?dl=0" target="_blank">https://www.dropbox.com/s/u1w0hg8fvoecqhr/lag-lead_mixratio.nc?dl=0</a></div><div><b>temp:</b> <a href="https://www.dropbox.com/s/9wxf5lbzjsrfukq/lag-lead_temp.nc?dl=0" target="_blank">https://www.dropbox.com/s/9wxf5lbzjsrfukq/lag-lead_temp.nc?dl=0</a></div><div><br></div><div>Attached in this email is the script that I am using.</div><div><br></div><div><br></div><div><b>[Problems]</b></div><div><b><br></b></div><div><b>[1] </b>I am getting -Inf for Q2. The error occurs when calculating <b>dHdp</b>. I can't figure out why this gives -Inf. The calculation is straightforward as:</div><div>   p2 = conform(H,p,1)</div>   dHdp = center_finite_diff_n (H,p2,False,1,0)<div><br></div><div><b>[2] </b>The values of Q1 is very large (min=-62600.10594520718   max=81387.3026968562)</div><div>I was expecting this to be less than 1000. </div><div><br></div><div>Any suggestions on how I can do this correctly in NCL?</div><div><br></div><div>I'll appreciate any help....</div><div><br></div><div>Sincerely,</div><div><div>Lyndz</div><div><b><br></b></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>