<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">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">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">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">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">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>