<div dir="ltr"><div><div><div>For the units section of a vars entry for netCDF files, you need to input the dimension ordering in the netcdf file. This is easily obtained by the order of the dimensions in the header entry for each array in the output from ncdump. As you have shown, for q, the order is time, a vertical dimension, then lat and lon. There are special one-character codes that you need to use in a Grads control file. It's pretty straightforward:<br></div>t for time<br></div>z for vertical dimension<br></div><div>y for south-north dimension<br></div><div>x for west-east dimension<br></div><div>With that, it should make sense that what you want to put in place of '99' is t,z,y,x. I don't know if spaces are a concern, but to be save, don't include any.<br><br></div><div>Jeff Duda<br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Nov 13, 2015 at 1:53 PM, Luis Blacutt <span dir="ltr"><<a href="mailto:luis.blacutt@gmail.com" target="_blank">luis.blacutt@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 Grads users<div><br></div><div>I have a problem generating a ctl file to read 11 years of data downloaded from ECWMF. I can read individually with sdfopen, however, when I tried to create a ctl file I think I have problems declaring the variable. My ctl file looks like this</div><div><div><br></div><div>DSET ^q-%<a href="http://y4.nc" target="_blank">y4.nc</a></div><div>UNDEF -999.99</div><div>OPTIONS little_endian yrev zrev template</div><div>DTYPE netcdf</div><div>XDEF 111 linear 275.000000 0.500000</div><div>YDEF 135 linear -55.000000 0.500000</div><div>TDEF 16072 linear 00:00Z01jan1999 6hr<br></div><div>zdef 11 levels 1000.000000 925.000000 850.000000 700.000000 500.000000 300.000000 250.000000 200.000000 100.000000 20.000000 10.000000</div><div>vars 1<br></div><div>q 11 99 Specific humidity [kg kg**-1]</div><div>endvars</div></div><div><br></div><div>I read that I can not use 99 for my variable definition, however I'm clueless about what to put in that line</div><div><br></div><div>The ncdump -h from an individual netcdf file looks like this</div><div><div>netcdf q-1999 {</div><div>dimensions:</div><div><span style="white-space:pre-wrap"> </span>longitude = 111 ;</div><div><span style="white-space:pre-wrap"> </span>latitude = 135 ;</div><div><span style="white-space:pre-wrap"> </span>level = 11 ;</div><div><span style="white-space:pre-wrap"> </span>time = UNLIMITED ; // (1460 currently)</div><div>variables:</div><div><span style="white-space:pre-wrap"> </span>float longitude(longitude) ;</div><div><span style="white-space:pre-wrap"> </span>longitude:units = "degrees_east" ;</div><div><span style="white-space:pre-wrap"> </span>longitude:long_name = "longitude" ;</div><div><span style="white-space:pre-wrap"> </span>float latitude(latitude) ;</div><div><span style="white-space:pre-wrap"> </span>latitude:units = "degrees_north" ;</div><div><span style="white-space:pre-wrap"> </span>latitude:long_name = "latitude" ;</div><div><span style="white-space:pre-wrap"> </span>int level(level) ;</div><div><span style="white-space:pre-wrap"> </span>level:units = "millibars" ;</div><div><span style="white-space:pre-wrap"> </span>level:long_name = "pressure_level" ;</div><div><span style="white-space:pre-wrap"> </span>int time(time) ;</div><div><span style="white-space:pre-wrap"> </span>time:units = "hours since 1900-01-01 00:00:0.0" ;</div><div><span style="white-space:pre-wrap"> </span>time:long_name = "time" ;</div><div><span style="white-space:pre-wrap"> </span>time:calendar = "gregorian" ;</div><div><span style="white-space:pre-wrap"> </span>short q(time, level, latitude, longitude) ;</div><div><span style="white-space:pre-wrap"> </span>q:scale_factor = 3.60802682733904e-07 ;</div><div><span style="white-space:pre-wrap"> </span>q:add_offset = 0.0118165291556829 ;</div><div><span style="white-space:pre-wrap"> </span>q:_FillValue = -32767s ;</div><div><span style="white-space:pre-wrap"> </span>q:missing_value = -32767s ;</div><div><span style="white-space:pre-wrap"> </span>q:units = "kg kg**-1" ;</div><div><span style="white-space:pre-wrap"> </span>q:long_name = "Specific humidity" ;</div><div><span style="white-space:pre-wrap"> </span>q:standard_name = "specific_humidity" ;</div><div><br></div><div>// global attributes:</div><div><span style="white-space:pre-wrap"> </span>:Conventions = "CF-1.0" ;</div><div><span style="white-space:pre-wrap"> </span>:history = "2015-11-03 16:49:08 GMT by grib_to_netcdf-1.13.1: grib_to_netcdf /data/data04/scratch/netcdf-atls20-a562cefde8a29a7288fa0b8b7f9413f7-AHRBQo.target -o /data/data04/scratch/netcdf-atls20-a562cefde8a29a7288fa0b8b7f9413f7-ZNlmj2.nc -utime" ;</div><div>}</div></div><div><br></div><div>Best regards</div><span class="HOEnZb"><font color="#888888"><div>Luis</div></font></span></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">Jeff Duda<br>Graduate research assistant<br>University of Oklahoma School of Meteorology<br>Center for Analysis and Prediction of Storms<br></div>
</div>