<div dir="ltr"><div>Sometime ago, I  created a test diabatic heating function. I  can not remember why I did not make it public. Possibly because doing so required extensive testing, writing documentation and support. My guess is that I just did not have the time.  Sorry about that!</div><div>======</div><div>See:  <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>  <a href="https://www.ncl.ucar.edu/Applications/Scripts/Q1Q2_yanai_1.ncl">Q1Q2_yanai.ncl</a></div><div><br></div><div>Read the text for this example.  The 'Q1' term represents the total diabatic heating.</div><div>*Maybe* this will be helpful.    Sorry, I can not offer more support.</div><div>----<br></div><div>Also: see Walter Hannah's page: "<b>A Comparison of Methods for Estimating Diabatic Heating</b>"</div><div><br></div><div><a href="http://hannahlab.org/a-comparison-of-methods-for-estimating-diabatic-heating/"><b>http://hannahlab.org/a-comparison-of-methods-for-estimating-diabatic-heating/</b></a></div><div><br></div><div>He used NCL scripts. I am not sure if these are still available.</div><div><br></div><div>=======================</div><div><br></div><div>Good Luck</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Apr 28, 2022 at 10:23 AM Lyndz via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu">ncl-talk@mailman.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"><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>I saw this function/tutorial for calculating the Potential vorticity from the NCL webpage.</div><div><br></div><div><a href="https://www.ncl.ucar.edu/Applications/pot_vort.shtml#:~:text=There%20are%20two%20potential%20vorticity,static%20stability%20and%20potential%20temperature" target="_blank">https://www.ncl.ucar.edu/Applications/pot_vort.shtml#:~:text=There%20are%20two%20potential%20vorticity,static%20stability%20and%20potential%20temperature</a>.<br></div><div><br></div><div>I would like to compute for the Potential Vorticity Tendency (PVT) on an isobaric level and global rectilinear grid reanalysis.</div><div><br></div><div>According to the attached file, PVT can be decomposed into </div><div>PVT= B* [ horizontal Advection + vertical advection + diabatic heating + friction ]</div><div><br></div><div>I can compute for the horizontal advection using this: <a href="https://www.ncl.ucar.edu/Document/Functions/Contributed/advect_variable.shtml" target="_blank">https://www.ncl.ucar.edu/Document/Functions/Contributed/advect_variable.shtml</a></div><div><br></div><div>My problem is how to calculate the diabatic heating here. Any suggestions on how I can do this in NCL?</div><div><br></div><div>-Lyndz</div><div><br></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@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div>