<div dir="ltr"><div><div><div><div><div><div><div>You had print statements in the code. You should have looked at the output.<br><br>p  = f-&gt;lv_ISBL3                                 ; hPa<br>print(p)<br>t<br><br></div>You had<br><br></div>[snip]<br>     p = int2flt(f-&gt;lv_ISBL3)             ; &lt;=== hPa<br>    psfc = f-&gt;PRES_252_SFC        ; &quot;ground&quot; pressure =&gt; pressure  (Pa)<br>        ; change units and calculate needed variables<br>    psfc = psfc/100<br>    psfc@units = &quot;hPa&quot;<br><br>    p = p/100                                   ; NO!  why<br>    p@units = &quot;hPa&quot;                       ; p was in hPa<br></div>[snip]<br><br></div>Note: there is no need for:  <br>    p = int2flt(f-&gt;lv_ISBL3)  <br></div>the following is fine<br></div>    p = f-&gt;lv_ISBL3<br><br></div>See attachment<br><div><div><div><div><div><br></div></div></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Jan 29, 2016 at 3:12 PM, Michelle Serino <span dir="ltr">&lt;<a href="mailto:mmserino@tamu.edu" target="_blank">mmserino@tamu.edu</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Dennis,<div><br></div><div>Thank you for getting back to me; your advice helped!  I am now trying to adapt your advice to a separate script I have for the 25-km RUC2-Analysis, and I receive the error:</div><div><br></div><div><p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo">fatal:Le: Number of dimensions do not match, can&#39;t continue</p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo">fatal:[&quot;Execute.c&quot;:8575]:Execute: Error occurred at or near line 97 in file sounding25.ncl</p><p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><br></p><p style="margin:0px;line-height:normal"><font face="arial, helvetica, sans-serif">I did my best to make sure that I converted the variable names to those of the 25-km format and that I didn&#39;t make any other simple mistakes.  I could easily have overlooked something, though.  Again, my script is attached, and the direct link to the 25-km data is below. Thank you for your continued help!</font></p><p style="margin:0px;line-height:normal"><font face="arial, helvetica, sans-serif"><br></font></p><p style="margin:0px;line-height:normal"><font face="arial, helvetica, sans-serif"><a href="https://www.dropbox.com/s/fhq2s4o8vp0yiol/ruc2anl_252_20090605_1200_000.grb?dl=0" target="_blank">https://www.dropbox.com/s/fhq2s4o8vp0yiol/ruc2anl_252_20090605_1200_000.grb?dl=0</a><br></font></p><p style="margin:0px;line-height:normal"><font face="arial, helvetica, sans-serif"><br></font></p><p style="margin:0px;line-height:normal"><font face="arial, helvetica, sans-serif">Best regards,</font></p><p style="margin:0px;line-height:normal"><font face="arial, helvetica, sans-serif"><br></font></p><p style="margin:0px;line-height:normal"><font face="arial, helvetica, sans-serif">Michelle Serino</font></p></div></div><div class="gmail_extra"><span class=""><br clear="all"><div><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr">Michelle Serino<div>M.S. Candidate</div><div><span style="font-size:small">Texas A&amp;M University</span><br></div><div>Atmospheric Sciences Dept.</div><div>Eller O&amp;M, Suite 1017</div><div>College Station, TX 77843<br></div><div><a href="tel:%28484%29%20769-0949" value="+14847690949" target="_blank">(484) 769-0949</a></div></div></div></div></div></div></div></div></div></div>
<br></span><div><div class="h5"><div class="gmail_quote">On Thu, Jan 28, 2016 at 8:45 PM, Dennis Shea <span dir="ltr">&lt;<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>&gt;</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>A key is  the &#39;level_type&#39; attribute associated with the pressure variables (PRES_). There are numerous PRES_ variables. In particular,<br><br>      float PRES_P0_L1_GLC0 ( ygrid_0, xgrid_0 )<br>         center :       US National Weather Service - NCEP (WMC)<br>         production_status :    Operational products<br>         long_name :    Pressure<br>         units :        Pa<br>         _FillValue :   1e+20<br>         coordinates :  gridlat_0 gridlon_0<br>         grid_type :    Lambert Conformal can be secant or tangent, conical or bipolar<br>         parameter_discipline_and_category :    Meteorological products, Mass<br>         parameter_template_discipline_category_number :        ( 0, 0, 3, 0 )<br>         level_type :   Ground or water surface   &lt;*******************************************<br>         level :         0<br>         forecast_time :        0<br>         forecast_time_units :  hours<br>         initial_time : 06/05/2009 (12:00)<br><br></div>It is a guess that this is the surface pressure (psfc). All pressure levels &gt; psfc must be extrapolated. [ *Extrapolation is always dangerous!] <br><br></div>The &#39;ind&#39; function can be used to get all &#39;valid&#39; levels [ p&lt;psfc ]. Try attached.<br><br></div>Cheers<br><div><div><div><div><br></div></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div>On Thu, Jan 28, 2016 at 11:47 AM, Michelle Serino <span dir="ltr">&lt;<a href="mailto:mmserino@tamu.edu" target="_blank">mmserino@tamu.edu</a>&gt;</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><div dir="ltr"><div>Good afternoon,</div><div><br></div><div>I am plotting RUC2-Analysis data on a skew-T.  Because this is analysis data, it seems that the data are always interpolated to the 1000-mb surface, regardless of the chosen location.  For example, in a high elevation location like Denver, the data extend down to 1000 mb, which is unrealistic.  This creates incorrect surface values, which results in a very incorrect value of surface CAPE.</div><div><br></div><div>I am wondering if, while still plotting on a skew-T, there is a way to specify the data that is plotted.  Or, if there is another, more realistic, way around this.  With a real sounding, I would not have this problem, so I am sure there is a way to fix this for the analysis data.  I have attached the script I am using, along with its plot and an observed sounding for the same location (Denver).  The direct link to the analysis data I am using is provided below.  I appreciate any help!</div><div><br></div><div><a href="https://www.dropbox.com/s/qpe96hmhh9hsi0h/ruc2anl_130_20090605_1200_000.grb2?dl=0" target="_blank">https://www.dropbox.com/s/qpe96hmhh9hsi0h/ruc2anl_130_20090605_1200_000.grb2?dl=0</a><br></div><div><br></div><div>Best regards,</div><div><br></div><div>Michelle Serino</div><span><font color="#888888"><div><br></div><div><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr">Michelle Serino<div>M.S. Candidate</div><div><span style="font-size:small">Texas A&amp;M University</span><br></div><div>Atmospheric Sciences Dept.</div><div>Eller O&amp;M, Suite 1017</div><div>College Station, TX 77843<br></div><div><a href="tel:%28484%29%20769-0949" value="+14847690949" target="_blank">(484) 769-0949</a></div></div></div></div></div></div></div></div></div></div>
</font></span></div>
<br></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><br>
<br></blockquote></div><br></div>
</blockquote></div><br></div></div></div>
</blockquote></div><br></div>