<div dir="ltr">My apologies. After switching back and forth between my 13-km script and the 25-km script, I was mixed up. Thanks for your help.<div><br></div><div>Best regards,</div><div><br></div><div>Michelle</div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature"><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&M University</span><br></div><div>Atmospheric Sciences Dept.</div><div>Eller O&M, Suite 1017</div><div>College Station, TX 77843<br></div><div>(484) 769-0949</div></div></div></div></div></div></div></div></div></div>
<br><div class="gmail_quote">On Sun, Jan 31, 2016 at 8:30 AM, 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><div><div><div><div>You had print statements in the code. You should have looked at the output.<br><br>p = f->lv_ISBL3 ; hPa<br>print(p)<br>t<br><br></div>You had<br><br></div>[snip]<br> p = int2flt(f->lv_ISBL3) ; <=== hPa<br> psfc = f->PRES_252_SFC ; "ground" pressure => pressure (Pa)<br> ; change units and calculate needed variables<br> psfc = psfc/100<br> psfc@units = "hPa"<br><br> p = p/100 ; NO! why<br> p@units = "hPa" ; p was in hPa<br></div>[snip]<br><br></div>Note: there is no need for: <br> p = int2flt(f->lv_ISBL3) <br></div>the following is fine<br></div> p = f->lv_ISBL3<br><br></div>See attachment<br><div><div><div><div><div><br></div></div></div></div></div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Jan 29, 2016 at 3:12 PM, Michelle Serino <span dir="ltr"><<a href="mailto:mmserino@tamu.edu" target="_blank">mmserino@tamu.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">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't continue</p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo">fatal:["Execute.c":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'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><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&M University</span><br></div><div>Atmospheric Sciences Dept.</div><div>Eller O&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><div class="gmail_quote">On Thu, Jan 28, 2016 at 8:45 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>A key is the 'level_type' 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 <*******************************************<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 > psfc must be extrapolated. [ *Extrapolation is always dangerous!] <br><br></div>The 'ind' function can be used to get all 'valid' levels [ p<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"><<a href="mailto:mmserino@tamu.edu" target="_blank">mmserino@tamu.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><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&M University</span><br></div><div>Atmospheric Sciences Dept.</div><div>Eller O&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>
</div></div></blockquote></div><br></div>