<div dir="ltr">Regarding point 1: I see that now.  I am not sure what I was thinking before, but thanks for the clarification.<div><br></div><div>Regarding point 2: I am not surprised that CAPE is calculated like that, making the reported value that of surface-based CAPE.  It seems to be incorrect only because the surface values for T and Td are incorrect.  Any ideas why that may be?  I realize that this is an analysis, but I expected the RUC to be a little more accurate/consistent.</div><div><br></div><div>Thank you,</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&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>(484) 769-0949</div></div></div></div></div></div></div></div></div></div>
<br><div class="gmail_quote">On Fri, Nov 6, 2015 at 3:49 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">[1]<br>
ncl_filedump is a command line utility. It operates like  &#39;ncdump -h ...&#39;<br>
<br>
It is not meant for use within a script<br>
<br>
---<br>
The same information is available within the script via<br>
<span class=""><br>
    dir  = &quot;/atmomounts/home/grad/mserino/Desktop/&quot;<br>
    fil  = &quot;ruc2anl_252_20090605_2000_000.grb&quot;<br>
<br>
; force a &#39;time&#39; dimension<br>
<br>
    setfileoption(&quot;grb&quot;,&quot;SingleElementDimensions&quot;,&quot;Initial_time&quot;)<br>
<br>
    f     = addfile(dir+fil,&quot;r&quot;)<br>
</span>    print(f)                              ; same info as:  ncl_fildump<br>
-itime  rucanl_252 ....grb<br>
<br>
[2]<br>
Well, it seems the software used has chosen a low level to begin the<br>
ascent of the parcel.<br>
The area enclosed by the sounding (black line and dashed red line) is<br>
quite large.<br>
 Why this occurred, I do not know.<br>
<span class="HOEnZb"><font color="#888888"><br>
D<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
<br>
On Fri, Nov 6, 2015 at 2:13 PM, Michelle Serino &lt;<a href="mailto:mmserino@tamu.edu">mmserino@tamu.edu</a>&gt; wrote:<br>
&gt; Thank you for getting back to me and for providing that script.  Also, thank<br>
&gt; you to Dave Brown for his input.  I have some follow-up questions/issues:<br>
&gt;<br>
&gt; 1) Why isn&#39;t &quot;ncl_filedump -itime ruc....grb&quot; in the script that you sent<br>
&gt; me?  Specifically, I would use a line like &quot;ncl_filedump -itime<br>
&gt; ruc2anl_252_20090605_2300_000.grb&quot; directly below the &quot;setfileoption&quot; line.<br>
&gt; Do I not need that line to fix the dimension problem you mentioned earlier?<br>
&gt; However, when I added that line to my code, I got a syntax error regarding<br>
&gt; the hyphenated portion of the line.  I typed it exactly as you provided it.<br>
&gt;<br>
&gt; 2) When I use my coordinates of interest (41.147, -104.801) in the script<br>
&gt; you sent, I still get inaccurate values of T, Td, and CAPE.  It is strange<br>
&gt; that the output seems ok for the coordinates you tried, but not for the one<br>
&gt; I need.<br>
&gt;<br>
&gt; Thanks again for your help.<br>
&gt;<br>
&gt; Michelle<br>
&gt;<br>
&gt; Michelle Serino<br>
&gt; M.S. Candidate<br>
&gt; Texas A&amp;M University<br>
&gt; Atmospheric Sciences Dept.<br>
&gt; Eller O&amp;M, Suite 1017<br>
&gt; College Station, TX 77843<br>
&gt; (484) 769-0949<br>
&gt;<br>
&gt; On Fri, Nov 6, 2015 at 2:52 PM, Dennis Shea &lt;<a href="mailto:shea@ucar.edu">shea@ucar.edu</a>&gt; wrote:<br>
&gt;&gt;<br>
&gt;&gt; The attached script produces the attached png files.<br>
&gt;&gt;<br>
&gt;&gt; Good luck<br>
&gt;&gt;<br>
&gt;&gt; On Tue, Nov 3, 2015 at 9:31 PM, Dennis Shea &lt;<a href="mailto:shea@ucar.edu">shea@ucar.edu</a>&gt; wrote:<br>
&gt;&gt; &gt; The different dimension names returned by &#39;ncl_filedump&#39; and your<br>
&gt;&gt; &gt; script is due to the fact that your script uses<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; ; force a &#39;time&#39; dimension<br>
&gt;&gt; &gt;     setfileoption(&quot;grb&quot;,&quot;SingleElementDimensions&quot;,&quot;Initial_time&quot;)<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; By default, ncl_filedump does not do this. However, the dimension<br>
&gt;&gt; &gt; names will match if you force the time dimension via<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; %&gt; ncl_filedump -itime ruc....grb<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; Credit to Dave Brown.<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; On Fri, Oct 23, 2015 at 11:07 AM, Michelle Serino &lt;<a href="mailto:mmserino@tamu.edu">mmserino@tamu.edu</a>&gt;<br>
&gt;&gt; &gt; wrote:<br>
&gt;&gt; &gt;&gt; Ok, I will install version 6.3.0.  Thank you so much for your prompt<br>
&gt;&gt; &gt;&gt; help!<br>
&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; Michelle Serino<br>
&gt;&gt; &gt;&gt; M.S. Candidate<br>
&gt;&gt; &gt;&gt; Texas A&amp;M University<br>
&gt;&gt; &gt;&gt; Atmospheric Sciences Dept.<br>
&gt;&gt; &gt;&gt; Eller O&amp;M, Suite 1017<br>
&gt;&gt; &gt;&gt; College Station, TX 77843<br>
&gt;&gt; &gt;&gt; (484) 769-0949<br>
&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; On Fri, Oct 23, 2015 at 10:40 AM, Dennis Shea &lt;<a href="mailto:shea@ucar.edu">shea@ucar.edu</a>&gt; wrote:<br>
&gt;&gt; &gt;&gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; offline<br>
&gt;&gt; &gt;&gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; [1] You should upgrade to 6.3.0<br>
&gt;&gt; &gt;&gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; [2]<br>
&gt;&gt; &gt;&gt;&gt; I have attached a file that contains the &#39;wind_direction&#39; and<br>
&gt;&gt; &gt;&gt;&gt; &#39;wind_speed&#39; functions.<br>
&gt;&gt; &gt;&gt;&gt; They may seem rather complicated. In fact, they are! However, to<br>
&gt;&gt; &gt;&gt;&gt; handles &#39;calm&#39; conditions and various dimension issues the NCL<br>
&gt;&gt; &gt;&gt;&gt; function must ve very robust.<br>
&gt;&gt; &gt;&gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; load &quot;./wind.ncl&quot;<br>
&gt;&gt; &gt;&gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; On Fri, Oct 23, 2015 at 9:29 AM, Michelle Serino &lt;<a href="mailto:mmserino@tamu.edu">mmserino@tamu.edu</a>&gt;<br>
&gt;&gt; &gt;&gt;&gt; wrote:<br>
&gt;&gt; &gt;&gt;&gt; &gt; I apologize for my vagueness.  The error I get is &quot;Undefined<br>
&gt;&gt; &gt;&gt;&gt; &gt; identifier:<br>
&gt;&gt; &gt;&gt;&gt; &gt; (wind_direction) is undefined, can&#39;t continue.&quot;  And I forgot to<br>
&gt;&gt; &gt;&gt;&gt; &gt; include<br>
&gt;&gt; &gt;&gt;&gt; &gt; the<br>
&gt;&gt; &gt;&gt;&gt; &gt; fact that I am using version 6.0.0, but I will use the equation in<br>
&gt;&gt; &gt;&gt;&gt; &gt; place<br>
&gt;&gt; &gt;&gt;&gt; &gt; of<br>
&gt;&gt; &gt;&gt;&gt; &gt; that function for now.<br>
&gt;&gt; &gt;&gt;&gt; &gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt; Michelle Serino<br>
&gt;&gt; &gt;&gt;&gt; &gt; M.S. Candidate<br>
&gt;&gt; &gt;&gt;&gt; &gt; Texas A&amp;M University<br>
&gt;&gt; &gt;&gt;&gt; &gt; Atmospheric Sciences Dept.<br>
&gt;&gt; &gt;&gt;&gt; &gt; Eller O&amp;M, Suite 1017<br>
&gt;&gt; &gt;&gt;&gt; &gt; College Station, TX 77843<br>
&gt;&gt; &gt;&gt;&gt; &gt; (484) 769-0949<br>
&gt;&gt; &gt;&gt;&gt; &gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt; On Fri, Oct 23, 2015 at 10:23 AM, Dennis Shea &lt;<a href="mailto:shea@ucar.edu">shea@ucar.edu</a>&gt; wrote:<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; I added the following print statement. The 1st 10 values for each<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; location are below.<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt;         wspd = sqrt(u^2 + v^2)<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt;         wdir   = wind_direction(u,v,0)<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt;         print(wspd+&quot;  &quot;+u+&quot;  &quot;+v+&quot;  &quot;+wdir)<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; (0)     location=(37.25,-107.75)  grid=(37.2053,-107.789)<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; (0)     17.4863  15.6  7.9  243.142<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; (1)     21.3026  19.4  8.8  245.601<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; (2)     29.0661  28  7.8  254.434<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; (3)     37.2979  35.7  10.8  253.168<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; (4)     39.3146  36  15.8  246.304<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; (5)     38.2027  33.2  18.9  240.348<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; (6)     36.0325  30.3  19.5  237.236<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; (7)     34.2118  27.9  19.8  234.638<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; (8)     33.1622  26.3  20.2  232.474<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; (9)     32.453  25.4  20.2  231.506<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; [snip]<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; (0)     location=(43.25,-101.75)  grid=(43.1875,-101.65)<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; (0)     17.4152  17.3  2  263.405<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; (1)     23.0235  22.8  -3.2  277.989<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; (2)     29.3097  29.1  -3.5  276.858<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; (3)     43.9915  43.8  -4.1  275.348<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; (4)     46.5897  46.4  -4.2  275.172<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; (5)     41.3241  41.2  -3.2  274.441<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; (6)     37.191  37.1  -2.6  274.009<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; (7)     34.8695  34.8  -2.2  273.617<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; (8)     32.9439  32.9  -1.7  272.958<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; (9)     31.1315  31.1  -1.4  272.578<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; =====<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; When you say error you never specified  the specific error message.<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; Was it something like &#39; wind_direction&#39; not found?<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; The &#39;wind_direction&#39; function  was introduced with NCL 6.1.0.<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; Perhaps<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; you have an old version of NCL. The current release is 6.3.0<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; wdir = atan2(u,v)*57.29579 + 180<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; On Fri, Oct 23, 2015 at 8:46 AM, Michelle Serino<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &lt;<a href="mailto:mmserino@tamu.edu">mmserino@tamu.edu</a>&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; wrote:<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; For reference, this error is regarding the line that says, &quot;wdir<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; =<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; f-&gt;wind_direction(u,v,0).&quot;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; Michelle Serino<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; M.S. Candidate<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; Texas A&amp;M University<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; Atmospheric Sciences Dept.<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; Eller O&amp;M, Suite 1017<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; College Station, TX 77843<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; (484) 769-0949<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; On Fri, Oct 23, 2015 at 9:39 AM, Michelle Serino<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; &lt;<a href="mailto:mmserino@tamu.edu">mmserino@tamu.edu</a>&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; wrote:<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt; Thank you for your quick reply!  To clarify, what did you<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt; accomplish<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt; by<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt; submitting the JIRA ticket?  Also, when I run your modified<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt; script,<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt; I<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt; am now<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt; getting an error that some vector subscripts need to be<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt; integers.<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt; Did<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt; you<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt; run into this as well?<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt; Thanks again,<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt; Michelle<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt; Michelle Serino<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt; M.S. Candidate<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt; Texas A&amp;M University<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt; Atmospheric Sciences Dept.<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt; Eller O&amp;M, Suite 1017<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt; College Station, TX 77843<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt; (484) 769-0949<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt; On Thu, Oct 22, 2015 at 6:55 PM, Dennis Shea &lt;<a href="mailto:shea@ucar.edu">shea@ucar.edu</a>&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt; wrote:<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; My speculation is that different grib tables are being<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; accessed.<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; JIRA-2297<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; On Thu, Oct 22, 2015 at 5:36 PM, Dennis Shea &lt;<a href="mailto:shea@ucar.edu">shea@ucar.edu</a>&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; wrote:<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt; There appears to be a naming error.<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt; When I do<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt; %&gt; ncl_filedump ruc2anl_252_20090605_2000_000.grb | less<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt; I get the variable names you use.<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt; ===<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt; When I execute your script I ran into errors. I inserted a<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt; f = addfile(...)<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt; print(f)                ; should match ncl_filedump BUT they<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt; don&#39;t<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt; ... NCL need to fix<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt; The variable names and dimension names are changed from the<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt; ncl_filedump.<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt; I will fill out a JIRA ticket.<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt; However, we can use the variable names from the &#39;print(f)&#39;,<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt; the<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt; script<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt; works correctly.<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt; ===<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt; See attached<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt; On Thu, Oct 22, 2015 at 4:55 PM, Michelle Serino<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt; &lt;<a href="mailto:mmserino@tamu.edu">mmserino@tamu.edu</a>&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt; wrote:<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; Good afternoon,<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; I am modifying one of NCAR&#39;s scripts (skewt_6.ncl) to plot a<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; skew-T<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; with<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; RUC2 Analysis data in grib1 format.  I have modified the<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; script<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; accordingly<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; to reflect the change in file type and to apply to the case<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; I am<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; studying,<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; but most everything else I have left as it was originally.<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; I am<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; receiving<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; an error that, after some extensive research and attempts, I<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; have<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; been<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; unable to fix, so I appreciate any help that anyone can<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; provide.<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; My code:<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl&quot;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl&quot;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/skewt_func.ncl&quot;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; begin<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; ; --- Read RUC GRIB file------------;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;     dir  = &quot;/atmomounts/home/grad/mserino/Desktop/&quot;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;     fil  = &quot;ruc2anl_252_20090605_2000_000.grb&quot;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;     ; force a &#39;time&#39; dimension<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; setfileoption(&quot;grb&quot;,&quot;SingleElementDimensions&quot;,&quot;Initial_time&quot;)<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;     f     = addfile(dir+fil,&quot;r&quot;)<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;     p     = int2flt(f-&gt;lv_ISBL2)               ;  ( lv_ISBL2<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; )<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;     time  = f-&gt;initial_time0_encoded  ; yyyymmddhh.hh_frac<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;     ; RUC grid point locations<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;     lat2d = f-&gt;gridlat_252              ; ( ygrid_0, xgrid_0<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; )<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;     lon2d = f-&gt;gridlon_252<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;     print(&quot;lat2d: min=&quot;+min(lat2d)+&quot;  ;  max=&quot;+max(lat2d))<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;     print(&quot;lon2d: min=&quot;+min(lon2d)+&quot;  ;  max=&quot;+max(lon2d))<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;     ;p     = p*0.01                    ; change units<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;     p@units = &quot;hPa&quot;                   ; skewT, mixhum_ptrh<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; use<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; mb<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; (hPa)<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;     ; --- Specify one or more locations<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;     lat   = (/ 37.25    , 43.25 /)<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;     lon   = (/-107.75 ,-101.75   /)<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;     npts  = dimsizes(lat)<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;     ; create plot(s)<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;     skewtOpts                 = True<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;     skewtOpts@DrawColAreaFill = True    ; default is False<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;     dataOpts   = True<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;     dataOpts@PrintZ = True<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;     do n=0,npts-1           ; loop over each grid pt<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;         ; find grid point nearest the user specified<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; location<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;         nm  = getind_latlon2d (lat2d,lon2d, lat(n), lon(n))<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;         nn  = nm(0,0)<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;         mm  = nm(0,1)<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;         print(&quot;location=(&quot;+lat(n)+&quot;,&quot;+lon(n)+&quot;)<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; grid=(&quot;+lat2d(nn,mm)+&quot;,&quot;+lon2d(nn,mm)+&quot;)&quot;)<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;         tk   = f-&gt;TMP_252_ISBL(0,:,nn,mm)<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;         z    = f-&gt;HGT_P0_252_ISBL(0,:,nn,mm)<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;         rh   = f-&gt;R_H_P0_252_ISBL(0,:,nn,mm)<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;         u    = f-&gt;U_GRD_252_ISBL(0,:,nn,mm)<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;         v    = f-&gt;V_GRD_252_ISBL(0,:,nn,mm)<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;         ; change units and calculate needed variables<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;         tc   = tk - 273.15<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;         tc@units= &quot;degC&quot;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;         q    = mixhum_ptrh (p, tk, rh, 2)<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;         q@units = &quot;kg/kg&quot;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;         tdc  = dewtemp_trh(tk,rh) - 273.15<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;         tdc@units = &quot;degC&quot;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;         wspd = sqrt(u^2 + v^2)<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;         wdir = wind_direction(u,v,0)<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;         itime= toint(time)<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;         skewtOpts@tiMainString    = &quot;RUC: &quot;+itime+&quot;:<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; (&quot;+lat(n)+&quot;,&quot;+lon(n)+&quot;)&quot;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;         ; each location will have a different file name<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;         wks  = gsn_open_wks (&quot;png&quot;,<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &quot;ruc2anl_skewt_&quot;+itime+&quot;_&quot;+sprinti(&quot;%0.3i&quot;,n))<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;         skewt_bkgd = skewT_BackGround (wks, skewtOpts)<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;         skewt_data = skewT_PlotData   (wks, skewt_bkgd,<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; p,tc,tdc,z,<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; wspd,wdir, dataOpts)<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;         draw (skewt_bkgd)<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;         draw (skewt_data)<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;         frame(wks)<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;     end do<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; end<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; The error I get regards the lines in red:<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; Either file (f) isn&#39;t defined or variable (lv_ISBL2) is not<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; a<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; variable<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; in<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; the file<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; The correct file, f, is defined because I am able to print<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; its<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; contents.<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; The variable &quot;lv_ISBL2&quot; is also a defined variable in the<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; file,<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; so<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; I<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; am not<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; sure why it is not recognized.  My code is attached, and the<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; RUC2<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; data<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; are<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; available at the link below.  Thank you again, and I look<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; forward<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; to<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; hearing<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; from you.<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; ruc2anl_252_20090605_2000_000.grb<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; Best regards,<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; Michelle Serino<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; M.S. Candidate<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; Texas A&amp;M University<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; Atmospheric Sciences Dept.<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; Eller O&amp;M, Suite 1017<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; College Station, TX 77843<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; (484) 769-0949<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; _______________________________________________<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; ncl-talk mailing list<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; <a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; List instructions, subscriber options, unsubscribe:<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;<br>
&gt;&gt; &gt;&gt;&gt; &gt;<br>
&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;<br>
&gt;<br>
&gt;<br>
</div></div></blockquote></div><br></div>