<div dir="ltr"><div><div>Have you actually tried just plotting the data with an arbitrary set of vertical levels? The values you specify for ZDEF in your control file do not have to exactly match the vertical levels of your data. What really matters is how many vertical levels there are.<br><br></div>Now, if you want to start performing calculations such as the spatial average of Q at a specific level like 850 mb, that's a bit more tricky, but as long as you know the structure of the vertical levels, you should be able to calculate it.<br><br></div>Jeff Duda<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, May 26, 2016 at 12:12 PM, Corey Gabriel <span dir="ltr"><<a href="mailto:cjgabriel7@gmail.com" target="_blank">cjgabriel7@gmail.com</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">Dear gradsusr community,<div><br></div><div>I am attempting to create a zonal mean plot of a variable (Q), which has 26 levels. </div><div><br></div><div>The vertical coordinates are "hybrid sigma pressure." I have created a data descriptor file and entered the vertical coordinates in the following manner: ZDEF 26 levels L1 L2 ... L26.</div><div><br></div><div>I have reviewed the GrADS documentation, and it seems that several examples suggest that the vertical coordinates should be in pure pressure coordinates and that GrADS may not read these hybrid coordinates.</div><div><br></div><div>The relevant ncdump infomration is provided below. Just below the dimensions, you can see that the formula from which the sigma coordinates are calculated from actual pressure and surface pressure are provided.</div><div><br></div><div>Has anyone attempted to plot hybrid sigma pressure coordinates in GrADS before? If so, would someone be able to provide an explanation of how to convert the vertical coordinate into something GrADS will read?</div><div><br></div><div>Thanks so much!</div><div>Corey</div><div><br></div><div>ncdump -c <a href="http://q.corey30.001.cam.2020-01-2089-12.nc" target="_blank">q.corey30.001.cam.2020-01-2089-12.nc</a><br></div><div>
<p><span>netcdf q.corey30.001.cam.2020-01-2089-12 {</span></p>
<p><span>dimensions:</span></p>
<p><span><span> </span>lev = 26 ;</span></p>
<p><span><span> </span>time = UNLIMITED ; // (840 currently)</span></p>
<p><span><span> </span>lat = 192 ;</span></p>
<p><span><span> </span>lon = 288 ;</span></p>
<p><span><span> </span>nbnd = 2 ;</span></p>
<p><span>variables:</span></p>
<p><span><span> </span>double lev(lev) ;</span></p>
<p><span><span> </span><span> </span>lev:long_name = "hybrid level at midpoints (1000*(A+B))" ;</span></p>
<p><span><span> </span><span> </span>lev:units = "level" ;</span></p>
<p><span><span> </span><span> </span>lev:positive = "down" ;</span></p>
<p><span><span> </span><span> </span>lev:standard_name = "atmosphere_hybrid_sigma_pressure_coordinate" ;</span></p>
<p><span><span> </span><span> </span>lev:formula_terms = "a: hyam b: hybm p0: P0 ps: PS" ;</span></p>
<p><span><span> </span>double time(time) ;</span></p>
<p><span><span> </span><span> </span>time:long_name = "time" ;</span></p>
<p><span><span> </span><span> </span>time:units = "days since 2004-01-01 00:00:00" ;</span></p>
<p><span><span> </span><span> </span>time:calendar = "noleap" ;</span></p>
<p><span><span> </span><span> </span>time:bounds = "time_bnds" ;</span></p>
<p><span><span> </span>double lat(lat) ;</span></p>
<p><span><span> </span><span> </span>lat:long_name = "latitude" ;</span></p>
<p><span><span> </span><span> </span>lat:units = "degrees_north" ;</span></p>
<p><span><span> </span>double lon(lon) ;</span></p>
<p><span><span> </span><span> </span>lon:long_name = "longitude" ;</span></p>
<p><span><span> </span><span> </span>lon:units = "degrees_east" ;</span></p>
<p><span><span> </span>double time_bnds(time, nbnd) ;</span></p>
<p><span><span> </span><span> </span>time_bnds:long_name = "time interval endpoints" ;</span></p>
<p><span><span> </span>float Q(time, lev, lat, lon) ;</span></p>
<p><span><span> </span><span> </span>Q:mdims = 1 ;</span></p>
<p><span><span> </span><span> </span>Q:units = "kg/kg" ;</span></p>
<p><span><span> </span><span> </span>Q:long_name = "Specific humidity" ;</span></p>
<p><span><span> </span><span> </span>Q:cell_methods = "time: mean" ;</span></p>
<p><span></span><br></p>
<p><span>// global attributes:</span></p>
<p><span><span> </span><span> </span>:Conventions = "CF-1.0" ;</span></p>
<p><span><span> </span><span> </span>:source = "CAM" ;</span></p>
<p><span><span> </span><span> </span>:case = "corey30" ;</span></p>
<p><span><span> </span><span> </span>:title = "UNSET" ;</span></p>
<p><span><span> </span><span> </span>:logname = "cgabriel" ;</span></p>
<p><span><span> </span><span> </span>:host = "ys0861" ;</span></p>
<p><span><span> </span><span> </span>:Version = "$Name$" ;</span></p>
<p><span><span> </span><span> </span>:revision_Id = "$Id$" ;</span></p>
<p><span><span> </span><span> </span>:initial_file = "/glade/p/acd/tilmes/inputdata/init/ccmi/b.e11.TSREFC2.f19.f19.ccmi23.1deg_26L.001" ;</span></p>
<p><span><span> </span><span> </span>:topography_file = "/glade/p/cesmdata/cseg/inputdata/atm/cam/topo/<a href="http://USGS-gtopo30_0.9x1.25_remap_c051027.nc" target="_blank">USGS-gtopo30_0.9x1.25_remap_c051027.nc</a>"</span></p></div></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>