<div dir="ltr">Hi Dennis, <div>ok . it was the "tweener" capability that I found most useful in 6.3.0. I used it to calculate the 600-400 mb average temperature. This was useful for identifying "warm core" cyclones. Also to calculate environmental steering flow. The environmental steering flow would be calculated with different upper boundaries to accommodate shallow and deep tropical storms. The top pressure level of the averaging layer would almost always be .gt. min(plev) (i.e. below the upper-most pressure level). <div><br></div><div>I'm not sure what the solution would be for pathological cases, but the old version always worked for me. <br><div><br></div><div><br><div><br></div><div><br><div><br></div></div></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Apr 4, 2017 at 4:33 PM, Dennis Shea <span dir="ltr"><<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div>I'll look. dpres_plevel was changed to handle a 'pathological' case over the Himalayas.<br><br></div><div>That said ... I should have changed the documentation to:<br><br><em style="margin:0px;padding:0px;color:rgb(0,0,0);font-family:verdana,sans-serif;font-size:13.3333px;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:left;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255)">ptop</em><span style="color:rgb(0,0,0);font-family:verdana,sans-serif;font-size:13.3333px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:left;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);display:inline;float:none"></span><p style="margin:10px 0px 20px;padding:0px;color:rgb(51,51,51);font-family:verdana,sans-serif;font-weight:normal;font-size:13.3333px;line-height:12pt;font-variant-ligatures:normal;font-variant-caps:normal;font-style:normal;letter-spacing:normal;text-align:left;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255)" class="m_2271263875853642661gmail-indent">A scalar specifying the top of the column. Must be: (a) .le. min(plev) and (b) have the same units as<span class="m_2271263875853642661gmail-Apple-converted-space"> </span><em style="margin:0px;padding:0px">plev</em>.</p><p style="margin:10px 0px 20px;padding:0px;color:rgb(51,51,51);font-family:verdana,sans-serif;font-weight:normal;font-size:13.3333px;line-height:12pt;font-variant-ligatures:normal;font-variant-caps:normal;font-style:normal;letter-spacing:normal;text-align:left;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255)" class="m_2271263875853642661gmail-indent">---<br></p></div><div>'ptop' was always meant to be ptop .le. min(plevel)<br><br></div><div>It was never meant to be an arbitrary *tweener* level. Hence the name ... ptop<span class=""><br><br>plev = (/200., 800 /)<br>pbot = 800.<br></span>ptop = min(plev) ; should be <= min(plev) <br>dp = dpres_plevel(plev, pbot, ptop, 0)<span class=""><br>print("plev="+plev+" dp="+dp) <br><br></span>---<br>(0) plev=200 dp=300<br>(1) plev=800 dp=300<br>--<br>ptop = 50<br>dp = dpres_plevel(plev, pbot, ptop, 0)<span class=""><br>print("plev="+plev+" dp="+dp) <br><br></span>(0) plev=200 dp=450<br>(1) plev=800 dp=300<br><br></div><div><br></div><br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="h5">On Mon, Apr 3, 2017 at 5:23 PM, David Ahijevych <span dir="ltr"><<a href="mailto:ahijevyc@ucar.edu" target="_blank">ahijevyc@ucar.edu</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><div dir="ltr">Hi, <div><br></div><div>I enjoy using the dpres_plevel_Wrap function to calculate weighted-averages, but the function changed from 6.3.0 to 6.4.0. <div><br></div><div>Here is a simple script :</div><div><br></div><div><div>load "$NCARG_ROOT/lib/ncarg/nclscri<wbr>pts/csm/contributed.ncl"</div><div>plev = (/200., 800 /)</div><div>pbot = 800.</div><div>ptop = 250.</div><div>dp = dpres_plevel_Wrap(plev, pbot, ptop, 0)</div><div>print("plev="+plev+" dp="+dp)</div></div><div><br></div><div><br></div><div>In 6.3.0, the result is</div></div><div><p style="margin:0px;font-size:11px;line-height:normal;font-family:menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">(0)<span class="m_2271263875853642661m_7093132981062599207gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>plev=200 dp=250</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">(1)<span class="m_2271263875853642661m_7093132981062599207gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>plev=800 dp=300</span></p></div><div><span style="font-variant-ligatures:no-common-ligatures">In 6.4.0 the result is</span></div><div><span style="font-variant-ligatures:no-common-ligatures"><p style="margin:0px;font-size:11px;line-height:normal;font-family:menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">(0)<span class="m_2271263875853642661m_7093132981062599207gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>plev=200 dp=9.96921e+36</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">(1)<span class="m_2271263875853642661m_7093132981062599207gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>plev=800 dp=9.96921e+36</span></p><div><span style="font-variant-ligatures:no-common-ligatures"><br></span></div><div><span style="font-variant-ligatures:no-common-ligatures">I think the 6.3.0 way is correct. </span></div><div><span style="font-variant-ligatures:no-common-ligatures"><br></span></div><div><span style="font-variant-ligatures:no-common-ligatures">Dave</span></div><div><span style="font-variant-ligatures:no-common-ligatures"><br></span></div></span></div></div>
<br></div></div>______________________________<wbr>_________________<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/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</blockquote></div><br></div>