<div dir="ltr"><div><div dir="ltr" class="m_-5172389275171024825gmail_signature" data-smartmail="gmail_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 Sir Dennis and NCL experts,</div><div><br></div><div>Thank you for sending the script about this.</div><div><br></div><div>I still have one question about the units of the following line:</div><div><br></div><div><br></div><div> Q1 = Cp*dTdt - Cp*(omega*ss-advT)<br> Q1@long_name = "Total Diabatic Heating as the Apparent Heat Source"<br> Q1@units = "K/s" <br></div><div><br></div><div>From dimensional analysis:</div><div><br></div><div>Cp = J/K*kg</div><div>dTdt = K/s</div><div>omega = Pa/s</div><div>ss = K/Pa</div><div>advT = K/s</div><div><br></div><div>Q1=(J/Kkg)*K/s - (J/Kkg)*(Pa/s*(K/Pa) - K/s)</div><div><br></div><div>Therefore the unit of Q1 should be J/kg/s (Joules per kg per second)?</div><div><br></div><div>I would like the unit of Q1 to be in W/m2. I will still need to integrate it vertically.</div><div><br></div><div>I'll appreciate any help about this.</div><div><br></div><div>Sincerely,</div><div>Lyndz</div><div><br></div><div><br></div><div><br></div></div></div></div></div></div></div></div></div></div></div><br><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">---------- Forwarded message ---------<br>From: <strong class="gmail_sendername" dir="auto">Dennis Shea</strong> <span dir="auto"><<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>></span><br>Date: Mon, Sep 9, 2019 at 12:17 PM<br>Subject: Re: [ncl-talk] Computing apparent heat source from reanalysis data<br>To: Lyndz <<a href="mailto:olagueralyndonmark429@gmail.com" target="_blank">olagueralyndonmark429@gmail.com</a>><br></div><br><br><div dir="ltr"><div>Yes. All arrays MUST have the same orientation.</div><div>Please look at what I sent. <br></div><div>=====<br></div><div>You must be careful with units.</div><div><br></div><div>Please ... questions only to ncl-talk. <br></div><div>===============</div><div><div>The attached is a sketch of calculating Q1 and Q2. It is not tested!</div><div><br></div><div>You can use this as a template of how to proceed.</div><div>Again, you must be careful with units. <br></div><div><br></div><div>That is your responsibility.<br></div><div><br></div></div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Sep 8, 2019 at 9:04 PM Lyndz <<a href="mailto:olagueralyndonmark429@gmail.com" target="_blank">olagueralyndonmark429@gmail.com</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"><font face="arial, sans-serif">Dear Sir Dennis, </font><div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif">Thank you for your help.<br></font><div><br></div><div><span style="font-family:arial,sans-serif">I have another question about the equation</span><span style="font-family:arial,sans-serif"><font color="#ffffff"> </font><span style="font-size:12px;white-space:pre-wrap"><font color="#000000">Q</font><font color="#000000"> = Cp*dTdt - Cp*(omega*s1+advTcfd).</font></span></span><br></div><div><span style="background-color:rgb(255,255,255)"><span style="font-size:12px;white-space:pre-wrap"><font face="arial, sans-serif" color="#000000"><br></font></span></span></div><div><span style="background-color:rgb(255,255,255)"><span style="font-size:12px;white-space:pre-wrap"><font face="arial, sans-serif" color="#000000">I inverted the orientation for the latitude (N-S) when I computed <i>advTcfd </i>using the ::-1. </font></span></span></div><div><span style="background-color:rgb(255,255,255)"><span style="font-size:12px;white-space:pre-wrap"><font color="#000000"><font face="arial, sans-serif">Do I also need to do this when computing dTdt? </font></font></span></span></div><div><span style="background-color:rgb(255,255,255)"><span style="font-size:12px;white-space:pre-wrap"><font color="#000000"><font face="arial, sans-serif"><br></font></font></span></span></div><div><span style="background-color:rgb(255,255,255)"><span style="font-size:12px;white-space:pre-wrap"><font color="#000000"><font face="arial, sans-serif"><br></font></font></span></span></div><div><span style="background-color:rgb(255,255,255)"><span style="font-size:12px;white-space:pre-wrap"><font color="#000000"><font face="arial, sans-serif">Sincerely,</font></font></span></span></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Sep 9, 2019 at 11:58 AM Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@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>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" target="_blank">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="m_-5172389275171024825m_-7723108804868556284gmail-m_-3595094406077490079gmail-m_-8571885415037845368gmail-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>
</blockquote></div>
</blockquote></div>
</div></div>