<div dir="ltr"><div><div><div>Corey,<br></div>I would tend to think the uncommented ZDEF entry should have worked. But the commented one should have also, unless that line exceeds 255 characters, as specified in the rules for the ZDEF entry of a control file (see here: <a href="http://cola.gmu.edu/grads/gadoc/descriptorfile.html#ZDEF">http://cola.gmu.edu/grads/gadoc/descriptorfile.html#ZDEF</a>).<br><br></div>Did you mess around with adding &quot;options pascals&quot;? Also, what is your dtype? It might help if you attached your control file text for troubleshooting.<br><br></div>Jeff<br><div><br><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, May 29, 2016 at 11:25 AM, Corey Gabriel <span dir="ltr">&lt;<a href="mailto:cjgabriel7@gmail.com" target="_blank">cjgabriel7@gmail.com</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 Jeff.<div><br></div><div>Thanks.  Fortunately I was able to use a couple of other tools to perform the calculations </div><div><br></div><div>I did try plotting the data with an arbitrary set of vertical levels.  That leads to an error message stating that ZDEF is now an invalid description file record.  Even though I was able to handle this elsewhere, I would still like to understand what I am doing wrong with GrADS.  I have copied the dimension records from my data descriptor file below.  To the ZDEF entries reflect what you mean by <span style="font-size:12.8px">plotting the data with an arbitrary set of </span><span style="font-size:12.8px">vertical levels?</span></div><div><br></div><div><div><br></div><div>XDEF lon 288 linear 0 1.25</div><div>YDEF lat 192 linear -90 0.94240837 </div><div>ZDEF 26 levels 1000 970 930 850 800 700 600 500 420 370 300 270 230 200 160 140 120 100 85 70 50 40 25 15 8 4 (arbitrary z axis values)</div><div>*ZDEF 26 levels 992.556 970.554 929.648 867.160 787.702 696.796 600.524 510.455 433.895 368.817 313.501 266.481 226.513 192.539 163.662 139.115 118.250 100.514 85.439 70.059 53.114 37.230 23.944625 13.967214 7.388 3.544 (<b>values matching the sigma-pressure levels)</b></div><div>TDEF time 840 linear 01jan2020 1mo</div><div><br></div><div>Thanks again!</div><div>Corey</div><div><br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, May 29, 2016 at 12:00 PM,  <span dir="ltr">&lt;<a href="mailto:gradsusr-request@gradsusr.org" target="_blank">gradsusr-request@gradsusr.org</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Send gradsusr mailing list submissions to<br>
        <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
        <a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
or, via email, send a message with subject or body &#39;help&#39; to<br>
        <a href="mailto:gradsusr-request@gradsusr.org" target="_blank">gradsusr-request@gradsusr.org</a><br>
<br>
You can reach the person managing the list at<br>
        <a href="mailto:gradsusr-owner@gradsusr.org" target="_blank">gradsusr-owner@gradsusr.org</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than &quot;Re: Contents of gradsusr digest...&quot;<br>
<br>
<br>
Today&#39;s Topics:<br>
<br>
   1. Re: Atmosphere_hybrid_sigma_pressure_coordinate (Jeff Duda)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Sun, 29 May 2016 00:16:50 -0500<br>
From: Jeff Duda &lt;<a href="mailto:jeffduda319@gmail.com" target="_blank">jeffduda319@gmail.com</a>&gt;<br>
Subject: Re: [gradsusr] Atmosphere_hybrid_sigma_pressure_coordinate<br>
To: GrADS Users Forum &lt;<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>&gt;<br>
Message-ID:<br>
        &lt;<a href="mailto:CAAig09BEnghZGYfm3kKgX-JWs9xyPSHBoFk%2B9P%2BFcsXs92he4A@mail.gmail.com" target="_blank">CAAig09BEnghZGYfm3kKgX-JWs9xyPSHBoFk+9P+FcsXs92he4A@mail.gmail.com</a>&gt;<br>
Content-Type: text/plain; charset=&quot;utf-8&quot;<br>
<br>
Have you actually tried just plotting the data with an arbitrary set of<br>
vertical levels? The values you specify for ZDEF in your control file do<br>
not have to exactly match the vertical levels of your data. What really<br>
matters is how many vertical levels there are.<br>
<br>
Now, if you want to start performing calculations such as the spatial<br>
average of Q at a specific level like 850 mb, that&#39;s a bit more tricky, but<br>
as long as you know the structure of the vertical levels, you should be<br>
able to calculate it.<br>
<br>
Jeff Duda<br>
<br>
On Thu, May 26, 2016 at 12:12 PM, Corey Gabriel &lt;<a href="mailto:cjgabriel7@gmail.com" target="_blank">cjgabriel7@gmail.com</a>&gt;<br>
wrote:<br>
<br>
&gt; Dear gradsusr community,<br>
&gt;<br>
&gt; I am attempting to create a zonal mean plot of a variable (Q), which has<br>
&gt; 26 levels.<br>
&gt;<br>
&gt; The vertical coordinates are &quot;hybrid sigma pressure.&quot;  I have created a<br>
&gt; data descriptor file and entered the vertical coordinates in the following<br>
&gt; manner: ZDEF 26 levels L1 L2 ... L26.<br>
&gt;<br>
&gt; I have reviewed the GrADS documentation, and it seems that several<br>
&gt; examples suggest that the vertical coordinates should be in pure pressure<br>
&gt; coordinates and that GrADS may not read these hybrid coordinates.<br>
&gt;<br>
&gt; The relevant ncdump infomration is provided below.  Just below the<br>
&gt; dimensions, you can see that the formula from which the sigma coordinates<br>
&gt; are calculated from actual pressure and surface pressure are provided.<br>
&gt;<br>
&gt; Has anyone attempted to plot hybrid sigma pressure coordinates in GrADS<br>
&gt; before?  If so, would someone be able to provide an explanation of how to<br>
&gt; convert the vertical coordinate into something GrADS will read?<br>
&gt;<br>
&gt; Thanks so much!<br>
&gt; Corey<br>
&gt;<br>
&gt; ncdump -c <a href="http://q.corey30.001.cam.2020-01-2089-12.nc" rel="noreferrer" target="_blank">q.corey30.001.cam.2020-01-2089-12.nc</a><br>
&gt;<br>
&gt; netcdf q.corey30.001.cam.2020-01-2089-12 {<br>
&gt;<br>
&gt; dimensions:<br>
&gt;<br>
&gt; lev = 26 ;<br>
&gt;<br>
&gt; time = UNLIMITED ; // (840 currently)<br>
&gt;<br>
&gt; lat = 192 ;<br>
&gt;<br>
&gt; lon = 288 ;<br>
&gt;<br>
&gt; nbnd = 2 ;<br>
&gt;<br>
&gt; variables:<br>
&gt;<br>
&gt; double lev(lev) ;<br>
&gt;<br>
&gt; lev:long_name = &quot;hybrid level at midpoints (1000*(A+B))&quot; ;<br>
&gt;<br>
&gt; lev:units = &quot;level&quot; ;<br>
&gt;<br>
&gt; lev:positive = &quot;down&quot; ;<br>
&gt;<br>
&gt; lev:standard_name = &quot;atmosphere_hybrid_sigma_pressure_coordinate&quot; ;<br>
&gt;<br>
&gt; lev:formula_terms = &quot;a: hyam b: hybm p0: P0 ps: PS&quot; ;<br>
&gt;<br>
&gt; double time(time) ;<br>
&gt;<br>
&gt; time:long_name = &quot;time&quot; ;<br>
&gt;<br>
&gt; time:units = &quot;days since 2004-01-01 00:00:00&quot; ;<br>
&gt;<br>
&gt; time:calendar = &quot;noleap&quot; ;<br>
&gt;<br>
&gt; time:bounds = &quot;time_bnds&quot; ;<br>
&gt;<br>
&gt; double lat(lat) ;<br>
&gt;<br>
&gt; lat:long_name = &quot;latitude&quot; ;<br>
&gt;<br>
&gt; lat:units = &quot;degrees_north&quot; ;<br>
&gt;<br>
&gt; double lon(lon) ;<br>
&gt;<br>
&gt; lon:long_name = &quot;longitude&quot; ;<br>
&gt;<br>
&gt; lon:units = &quot;degrees_east&quot; ;<br>
&gt;<br>
&gt; double time_bnds(time, nbnd) ;<br>
&gt;<br>
&gt; time_bnds:long_name = &quot;time interval endpoints&quot; ;<br>
&gt;<br>
&gt; float Q(time, lev, lat, lon) ;<br>
&gt;<br>
&gt; Q:mdims = 1 ;<br>
&gt;<br>
&gt; Q:units = &quot;kg/kg&quot; ;<br>
&gt;<br>
&gt; Q:long_name = &quot;Specific humidity&quot; ;<br>
&gt;<br>
&gt; Q:cell_methods = &quot;time: mean&quot; ;<br>
&gt;<br>
&gt;<br>
&gt; // global attributes:<br>
&gt;<br>
&gt; :Conventions = &quot;CF-1.0&quot; ;<br>
&gt;<br>
&gt; :source = &quot;CAM&quot; ;<br>
&gt;<br>
&gt; :case = &quot;corey30&quot; ;<br>
&gt;<br>
&gt; :title = &quot;UNSET&quot; ;<br>
&gt;<br>
&gt; :logname = &quot;cgabriel&quot; ;<br>
&gt;<br>
&gt; :host = &quot;ys0861&quot; ;<br>
&gt;<br>
&gt; :Version = &quot;$Name$&quot; ;<br>
&gt;<br>
&gt; :revision_Id = &quot;$Id$&quot; ;<br>
&gt;<br>
&gt; :initial_file =<br>
&gt; &quot;/glade/p/acd/tilmes/inputdata/init/ccmi/b.e11.TSREFC2.f19.f19.ccmi23.1deg_26L.001&quot;<br>
&gt; ;<br>
&gt;<br>
&gt; :topography_file = &quot;/glade/p/cesmdata/cseg/inputdata/atm/cam/topo/<br>
&gt; <a href="http://USGS-gtopo30_0.9x1.25_remap_c051027.nc" rel="noreferrer" target="_blank">USGS-gtopo30_0.9x1.25_remap_c051027.nc</a>&quot;<br>
&gt;<br>
&gt; _______________________________________________<br>
&gt; gradsusr mailing list<br>
&gt; <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
&gt;<br>
&gt;<br>
<br>
<br>
--<br>
Jeff Duda<br>
Graduate research assistant<br>
University of Oklahoma School of Meteorology<br>
Center for Analysis and Prediction of Storms<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20160529/1d67c2cb/attachment-0001.html" rel="noreferrer" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20160529/1d67c2cb/attachment-0001.html</a><br>
<br>
------------------------------<br>
<br>
_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br>
<br>
End of gradsusr Digest, Vol 75, Issue 29<br>
****************************************<br>
</blockquote></div><br></div>
<br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">Jeff Duda<br>Graduate research assistant<br>University of Oklahoma School of Meteorology<br>Center for Analysis and Prediction of Storms<br></div>
</div>