<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Well, since a bug was just discovered when using xdfopen to read CAM output, I suggest you make small mods to your descriptor and use 'open' instead. All you should have to do is remove the dimension names from your XDEF, YDEF, ZDEF, and TDEF statements, and add UNDEF. <div>--Jennifer</div><div><br><div><br><div><br><div><div>On Feb 10, 2011, at 11:24 AM, Steven Chan wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div>Hello All:<br><br>I have been struggling to open a certain netCDF file with xdfopen.<br><br>GRADS version used: 2.0.a9 & 2.0.a9.oga.1<br>Error message:<br><blockquote type="cite">ga-> xdfopen fqd000_nch.cam2.h1.2005-07-01-00000.xdf<br></blockquote><blockquote type="cite">Scanning Descriptor File: fqd000_nch.cam2.h1.2005-07-01-00000.xdf<br></blockquote><blockquote type="cite">gadsdf: SDF file does not have any non-coordinate variables.<br></blockquote><blockquote type="cite">SDF Descriptor file fqd000_nch.cam2.h1.2005-07-01-00000.xdf was not <br></blockquote><blockquote type="cite">successfully opened & parsed.<br></blockquote><br>I am using the following "xdf control file":<br><blockquote type="cite">DSET ^fqd000_nch.cam2.h1.2005-07-01-00000.nc<br></blockquote><blockquote type="cite">TITLE fqd000_nch.cam2.h1.2005-07-01-00000.nc<br></blockquote><blockquote type="cite">*<br></blockquote><blockquote type="cite">* Vertical = Hybrid<br></blockquote><blockquote type="cite">*<br></blockquote><blockquote type="cite">OPTIONS 365_day_calendar<br></blockquote><blockquote type="cite">DTYPE netcdf<br></blockquote><blockquote type="cite">XDEF lon 1152 linear 0.00 0.3125<br></blockquote><blockquote type="cite">YDEF lat 768 linear -90.00 0.2346805737<br></blockquote><blockquote type="cite">ZDEF lev 30 levels<br></blockquote><blockquote type="cite">992.556095123291 976.325407391414 957.485479535535 936.198398470879<br></blockquote><blockquote type="cite">912.64454694464 887.020248919726 859.53476652503 820.858368650079<br></blockquote><blockquote type="cite">763.404481112957 691.389430314303 609.778694808483 524.687174707651<br></blockquote><blockquote type="cite">445.992574095726 379.100903868675 322.241902351379 273.910816758871<br></blockquote><blockquote type="cite">232.828618958592 197.908086702227 168.225079774857 142.994038760662<br></blockquote><blockquote type="cite">121.547240763903 103.317126631737 87.8212302923203 72.0124505460262<br></blockquote><blockquote type="cite">54.5954797416925 38.2682997733355 24.6122200042009 14.3566322512925<br></blockquote><blockquote type="cite">7.59481964632869 3.64346569404006<br></blockquote><blockquote type="cite">TDEF time 1 linear 00Z01JUL2005 3hr<br></blockquote><blockquote type="cite">*<br></blockquote><blockquote type="cite">VARS 11<br></blockquote><blockquote type="cite">OMEGA=>omega 30 t,z,y,x Vertical velocity (pressure) (Pa/s)<br></blockquote><blockquote type="cite">PBLH=>pblh 0 t,y,x PBL Height (m)<br></blockquote><blockquote type="cite">PRECC=>precc 0 t,y,x Convective precipitation rate (liq + ice) <br></blockquote><blockquote type="cite">(m/s)<br></blockquote><blockquote type="cite">PRECL=>precl 0 t,y,x Large-scale (stable) precipitation rate <br></blockquote><blockquote type="cite">(liq + ice) (m/s)<br></blockquote><blockquote type="cite">PRECSH=>precsh 0 t,y,x Shallow Convection precipitation rate (m/s)<br></blockquote><blockquote type="cite">PS=>ps 0 t,y,x Surface pressure (Pa)<br></blockquote><blockquote type="cite">Q=>q 30 t,z,y,x Specific humidity (kg/kg)<br></blockquote><blockquote type="cite">QFLX=>qflx 0 t,y,x Surface water flux (kg/m2/s)<br></blockquote><blockquote type="cite">T=>tmpprs 30 t,z,y,x Temperature (K)<br></blockquote><blockquote type="cite">U=>ugrdprs 30 t,z,y,x Zonal Wind (m/s)<br></blockquote><blockquote type="cite">V=>vgrdprs 30 t,z,y,x Meridonal Wind (m/s)<br></blockquote><blockquote type="cite">ENDVARS<br></blockquote><br>A netCDF file with the following ncdump information:<br><blockquote type="cite">netcdf fqd000_nch.cam2.h1.2005-07-01-00000 {<br></blockquote><blockquote type="cite">dimensions:<br></blockquote><blockquote type="cite"> lat = 768 ;<br></blockquote><blockquote type="cite"> lon = 1152 ;<br></blockquote><blockquote type="cite"> slat = 767 ;<br></blockquote><blockquote type="cite"> slon = 1152 ;<br></blockquote><blockquote type="cite"> lev = 30 ;<br></blockquote><blockquote type="cite"> ilev = 31 ;<br></blockquote><blockquote type="cite"> isccp_prs = 7 ;<br></blockquote><blockquote type="cite"> isccp_tau = 7 ;<br></blockquote><blockquote type="cite"> isccp_prstau = 49 ;<br></blockquote><blockquote type="cite"> time = UNLIMITED ; // (1 currently)<br></blockquote><blockquote type="cite"> tbnd = 2 ;<br></blockquote><blockquote type="cite"> chars = 8 ;<br></blockquote><blockquote type="cite">variables:<br></blockquote><blockquote type="cite"> double P0 ;<br></blockquote><blockquote type="cite"> P0:long_name = "reference pressure" ;<br></blockquote><blockquote type="cite"> P0:units = "Pa" ;<br></blockquote><blockquote type="cite"> double lat(lat) ;<br></blockquote><blockquote type="cite"> lat:long_name = "latitude" ;<br></blockquote><blockquote type="cite"> lat:units = "degrees_north" ;<br></blockquote><blockquote type="cite"> double lon(lon) ;<br></blockquote><blockquote type="cite"> lon:long_name = "longitude" ;<br></blockquote><blockquote type="cite"> lon:units = "degrees_east" ;<br></blockquote><blockquote type="cite"> double slat(slat) ;<br></blockquote><blockquote type="cite"> slat:long_name = "staggered latitude" ;<br></blockquote><blockquote type="cite"> slat:units = "degrees_north" ;<br></blockquote><blockquote type="cite"> double slon(slon) ;<br></blockquote><blockquote type="cite"> slon:long_name = "staggered longitude" ;<br></blockquote><blockquote type="cite"> slon:units = "degrees_east" ;<br></blockquote><blockquote type="cite"> double w_stag(slat) ;<br></blockquote><blockquote type="cite"> w_stag:long_name = "staggered latitude weights" ;<br></blockquote><blockquote type="cite"> double lev(lev) ;<br></blockquote><blockquote type="cite"> lev:long_name = "hybrid level at midpoints (1000*(A+B))" ;<br></blockquote><blockquote type="cite"> lev:units = "level" ;<br></blockquote><blockquote type="cite"> lev:positive = "down" ;<br></blockquote><blockquote type="cite"> lev:standard_name = <br></blockquote><blockquote type="cite">"atmosphere_hybrid_sigma_pressure_coordinate" ;<br></blockquote><blockquote type="cite"> lev:formula_terms = "a: hyam b: hybm p0: P0 ps: PS" ;<br></blockquote><blockquote type="cite"> double ilev(ilev) ;<br></blockquote><blockquote type="cite"> ilev:long_name = "hybrid level at interfaces <br></blockquote><blockquote type="cite">(1000*(A+B))" ;<br></blockquote><blockquote type="cite"> ilev:units = "level" ;<br></blockquote><blockquote type="cite"> ilev:positive = "down" ;<br></blockquote><blockquote type="cite"> ilev:standard_name = <br></blockquote><blockquote type="cite">"atmosphere_hybrid_sigma_pressure_coordinate" ;<br></blockquote><blockquote type="cite"> ilev:formula_terms = "a: hyai b: hybi p0: P0 ps: PS" ;<br></blockquote><blockquote type="cite"> double isccp_prs(isccp_prs) ;<br></blockquote><blockquote type="cite"> isccp_prs:long_name = "Mean ISCCP pressure" ;<br></blockquote><blockquote type="cite"> isccp_prs:units = "mb" ;<br></blockquote><blockquote type="cite"> isccp_prs:isccp_prs_bnds = 0., 180., 310., 440., 560., <br></blockquote><blockquote type="cite">680., 800., 1000. ;<br></blockquote><blockquote type="cite"> double isccp_tau(isccp_tau) ;<br></blockquote><blockquote type="cite"> isccp_tau:long_name = "Mean ISCCP optical depth" ;<br></blockquote><blockquote type="cite"> isccp_tau:units = "unitless" ;<br></blockquote><blockquote type="cite"> isccp_tau:isccp_tau_bnds = 0., 0.3, 1.3, 3.6, 9.4, <br></blockquote><blockquote type="cite">23., 60., 379. ;<br></blockquote><blockquote type="cite"> double isccp_prstau(isccp_prstau) ;<br></blockquote><blockquote type="cite"> isccp_prstau:long_name = "Mean pressure (mb).mean <br></blockquote><blockquote type="cite">optical depth (unitless)/1000" ;<br></blockquote><blockquote type="cite"> isccp_prstau:units = "mixed" ;<br></blockquote><blockquote type="cite"> double time(time) ;<br></blockquote><blockquote type="cite"> time:long_name = "time" ;<br></blockquote><blockquote type="cite"> time:units = "days since 2005-01-01 00:00:00" ;<br></blockquote><blockquote type="cite"> time:calendar = "noleap" ;<br></blockquote><blockquote type="cite"> time:bounds = "time_bnds" ;<br></blockquote><blockquote type="cite"> double time_bnds(time, tbnd) ;<br></blockquote><blockquote type="cite"> time_bnds:long_name = "time interval endpoints" ;<br></blockquote><blockquote type="cite"> char date_written(time, chars) ;<br></blockquote><blockquote type="cite"> char time_written(time, chars) ;<br></blockquote><blockquote type="cite"> int ntrm ;<br></blockquote><blockquote type="cite"> ntrm:long_name = "spectral truncation parameter M" ;<br></blockquote><blockquote type="cite"> int ntrn ;<br></blockquote><blockquote type="cite"> ntrn:long_name = "spectral truncation parameter N" ;<br></blockquote><blockquote type="cite"> int ntrk ;<br></blockquote><blockquote type="cite"> ntrk:long_name = "spectral truncation parameter K" ;<br></blockquote><blockquote type="cite"> int ndbase ;<br></blockquote><blockquote type="cite"> ndbase:long_name = "base day" ;<br></blockquote><blockquote type="cite"> int nsbase ;<br></blockquote><blockquote type="cite"> nsbase:long_name = "seconds of base day" ;<br></blockquote><blockquote type="cite"> int nbdate ;<br></blockquote><blockquote type="cite"> nbdate:long_name = "base date (YYYYMMDD)" ;<br></blockquote><blockquote type="cite"> int nbsec ;<br></blockquote><blockquote type="cite"> nbsec:long_name = "seconds of base date" ;<br></blockquote><blockquote type="cite"> int mdt ;<br></blockquote><blockquote type="cite"> mdt:long_name = "timestep" ;<br></blockquote><blockquote type="cite"> mdt:units = "s" ;<br></blockquote><blockquote type="cite"> int nlon(lat) ;<br></blockquote><blockquote type="cite"> nlon:long_name = "number of longitudes" ;<br></blockquote><blockquote type="cite"> int wnummax(lat) ;<br></blockquote><blockquote type="cite"> wnummax:long_name = "cutoff Fourier wavenumber" ;<br></blockquote><blockquote type="cite"> double hyai(ilev) ;<br></blockquote><blockquote type="cite"> hyai:long_name = "hybrid A coefficient at layer <br></blockquote><blockquote type="cite">interfaces" ;<br></blockquote><blockquote type="cite"> double hybi(ilev) ;<br></blockquote><blockquote type="cite"> hybi:long_name = "hybrid B coefficient at layer <br></blockquote><blockquote type="cite">interfaces" ;<br></blockquote><blockquote type="cite"> double hyam(lev) ;<br></blockquote><blockquote type="cite"> hyam:long_name = "hybrid A coefficient at layer <br></blockquote><blockquote type="cite">midpoints" ;<br></blockquote><blockquote type="cite"> double hybm(lev) ;<br></blockquote><blockquote type="cite"> hybm:long_name = "hybrid B coefficient at layer <br></blockquote><blockquote type="cite">midpoints" ;<br></blockquote><blockquote type="cite"> double gw(lat) ;<br></blockquote><blockquote type="cite"> gw:long_name = "gauss weights" ;<br></blockquote><blockquote type="cite"> int ndcur(time) ;<br></blockquote><blockquote type="cite"> ndcur:long_name = "current day (from base day)" ;<br></blockquote><blockquote type="cite"> int nscur(time) ;<br></blockquote><blockquote type="cite"> nscur:long_name = "current seconds of current day" ;<br></blockquote><blockquote type="cite"> int date(time) ;<br></blockquote><blockquote type="cite"> date:long_name = "current date (YYYYMMDD)" ;<br></blockquote><blockquote type="cite"> double co2vmr(time) ;<br></blockquote><blockquote type="cite"> co2vmr:long_name = "co2 volume mixing ratio" ;<br></blockquote><blockquote type="cite"> double ch4vmr(time) ;<br></blockquote><blockquote type="cite"> ch4vmr:long_name = "ch4 volume mixing ratio" ;<br></blockquote><blockquote type="cite"> double n2ovmr(time) ;<br></blockquote><blockquote type="cite"> n2ovmr:long_name = "n2o volume mixing ratio" ;<br></blockquote><blockquote type="cite"> double f11vmr(time) ;<br></blockquote><blockquote type="cite"> f11vmr:long_name = "f11 volume mixing ratio" ;<br></blockquote><blockquote type="cite"> double f12vmr(time) ;<br></blockquote><blockquote type="cite"> f12vmr:long_name = "f12 volume mixing ratio" ;<br></blockquote><blockquote type="cite"> double sol_tsi(time) ;<br></blockquote><blockquote type="cite"> sol_tsi:long_name = "total solar irradiance" ;<br></blockquote><blockquote type="cite"> sol_tsi:units = "W/m2" ;<br></blockquote><blockquote type="cite"> int datesec(time) ;<br></blockquote><blockquote type="cite"> datesec:long_name = "current seconds of current date" ;<br></blockquote><blockquote type="cite"> int nsteph(time) ;<br></blockquote><blockquote type="cite"> nsteph:long_name = "current timestep" ;<br></blockquote><blockquote type="cite"> float OMEGA(time, lev, lat, lon) ;<br></blockquote><blockquote type="cite"> OMEGA:units = "Pa/s" ;<br></blockquote><blockquote type="cite"> OMEGA:long_name = "Vertical velocity (pressure)" ;<br></blockquote><blockquote type="cite"> float PBLH(time, lat, lon) ;<br></blockquote><blockquote type="cite"> PBLH:units = "m" ;<br></blockquote><blockquote type="cite"> PBLH:long_name = "PBL height" ;<br></blockquote><blockquote type="cite"> float PRECC(time, lat, lon) ;<br></blockquote><blockquote type="cite"> PRECC:units = "m/s" ;<br></blockquote><blockquote type="cite"> PRECC:long_name = "Convective precipitation rate (liq <br></blockquote><blockquote type="cite">+ ice)" ;<br></blockquote><blockquote type="cite"> float PRECL(time, lat, lon) ;<br></blockquote><blockquote type="cite"> PRECL:units = "m/s" ;<br></blockquote><blockquote type="cite"> PRECL:long_name = "Large-scale (stable) precipitation <br></blockquote><blockquote type="cite">rate (liq + ice)" ;<br></blockquote><blockquote type="cite"> float PRECSH(time, lat, lon) ;<br></blockquote><blockquote type="cite"> PRECSH:units = "m/s" ;<br></blockquote><blockquote type="cite"> PRECSH:long_name = "Shallow Convection precipitation <br></blockquote><blockquote type="cite">rate" ;<br></blockquote><blockquote type="cite"> float PS(time, lat, lon) ;<br></blockquote><blockquote type="cite"> PS:units = "Pa" ;<br></blockquote><blockquote type="cite"> PS:long_name = "Surface pressure" ;<br></blockquote><blockquote type="cite"> float Q(time, lev, lat, lon) ;<br></blockquote><blockquote type="cite"> Q:units = "kg/kg" ;<br></blockquote><blockquote type="cite"> Q:long_name = "Specific humidity" ;<br></blockquote><blockquote type="cite"> float QFLX(time, lat, lon) ;<br></blockquote><blockquote type="cite"> QFLX:units = "kg/m2/s" ;<br></blockquote><blockquote type="cite"> QFLX:long_name = "Surface water flux" ;<br></blockquote><blockquote type="cite"> float T(time, lev, lat, lon) ;<br></blockquote><blockquote type="cite"> T:units = "K" ;<br></blockquote><blockquote type="cite"> T:long_name = "Temperature" ;<br></blockquote><blockquote type="cite"> float U(time, lev, lat, lon) ;<br></blockquote><blockquote type="cite"> U:units = "m/s" ;<br></blockquote><blockquote type="cite"> U:long_name = "Zonal wind" ;<br></blockquote><blockquote type="cite"> float V(time, lev, lat, lon) ;<br></blockquote><blockquote type="cite"> V:units = "m/s" ;<br></blockquote><blockquote type="cite"> V:long_name = "Meridional wind" ;<br></blockquote><blockquote type="cite">// global attributes:<br></blockquote><blockquote type="cite"> :Conventions = "CF-1.0" ;<br></blockquote><blockquote type="cite"> :source = "CAM" ;<br></blockquote><blockquote type="cite"> :case = "fqd000_nch" ;<br></blockquote><blockquote type="cite"> :title = "UNSET" ;<br></blockquote><blockquote type="cite"> :logname = "juliob" ;<br></blockquote><blockquote type="cite"> :host = "jaguarpf-batch8" ;<br></blockquote><blockquote type="cite"> :Version = "$Name$" ;<br></blockquote><blockquote type="cite"> :revision_Id = "$Id$" ;<br></blockquote><blockquote type="cite"> :initial_file = <br></blockquote><blockquote type="cite">"/tmp/work/juliob/init/qdLL05.cam2.i.1998-01-01-00000.nc" ;<br></blockquote><blockquote type="cite"> :topography_file = <br></blockquote><blockquote type="cite">"/tmp/proj/ccsm/inputdata/atm/cam/topo/USGS_gtopo30_0.23x0.31_remap_c061107.nc" <br></blockquote><blockquote type="cite">;<br></blockquote><blockquote type="cite">}<br></blockquote><br>If I remove VARS/ENDVARS from the xdf, the file can open, but the <br>metadata information is erroneous and data cannot be displayed:<br><blockquote type="cite">Scanning Descriptor File: fqd000_nch.cam2.h1.2005-07-01-00000.xdf<br></blockquote><blockquote type="cite">SDF file fqd000_nch.cam2.h1.2005-07-01-00000.nc is open as file 1<br></blockquote><blockquote type="cite">LON set to 0 360<br></blockquote><blockquote type="cite">LAT set to -90 90<br></blockquote><blockquote type="cite">LEV set to 992.556 992.556<br></blockquote><blockquote type="cite">Time values set: 2005:7:1:0 2005:7:1:0<br></blockquote><blockquote type="cite">E set to 1 1<br></blockquote><blockquote type="cite">ga-> q file<br></blockquote><blockquote type="cite">File 1 : fqd000_nch.cam2.h1.2005-07-01-00000.nc<br></blockquote><blockquote type="cite"> Descriptor: fqd000_nch.cam2.h1.2005-07-01-00000.xdf<br></blockquote><blockquote type="cite"> Binary: fqd000_nch.cam2.h1.2005-07-01-00000.nc<br></blockquote><blockquote type="cite"> Type = Gridded<br></blockquote><blockquote type="cite"> Xsize = 1152 Ysize = 768 Zsize = 30 Tsize = 1 Esize = 1<br></blockquote><blockquote type="cite"> Number of Variables = 11<br></blockquote><blockquote type="cite"> omega 30 t,z,y,x Vertical velocity (pressure)<br></blockquote><blockquote type="cite"> omega 0 t,y,x PBL height<br></blockquote><blockquote type="cite"> omega 0 t,y,x Convective precipitation rate (liq + ice)<br></blockquote><blockquote type="cite"> omega 0 t,y,x Large-scale (stable) precipitation rate (liq + ice)<br></blockquote><blockquote type="cite"> omega 0 t,y,x Shallow Convection precipitation rate<br></blockquote><blockquote type="cite"> omega 0 t,y,x Surface pressure<br></blockquote><blockquote type="cite"> omega 30 t,z,y,x Specific humidity<br></blockquote><blockquote type="cite"> omega 0 t,y,x Surface water flux<br></blockquote><blockquote type="cite"> omega 30 t,z,y,x Temperature<br></blockquote><blockquote type="cite"> omega 30 t,z,y,x Zonal wind<br></blockquote><blockquote type="cite"> omega 30 t,z,y,x Meridional wind<br></blockquote><br>I am fairly sure I am doing xdfopen correctly, but I may be there <br>something silly that I am just not realizing.<br>Thank you ahead for the responses.<br><br>Steven Chan<br><br>_______________________________________________<br>gradsusr mailing list<br><a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>http://gradsusr.org/mailman/listinfo/gradsusr<br></div></blockquote></div><br><div apple-content-edited="true"> <span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: auto; -khtml-text-decorations-in-effect: none; text-indent: 0px; -apple-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><div style="word-wrap: break-word; -khtml-nbsp-mode: space; -khtml-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: auto; -khtml-text-decorations-in-effect: none; text-indent: 0px; -apple-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: auto; -khtml-text-decorations-in-effect: none; text-indent: 0px; -apple-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><div>--</div><div>Jennifer M. Adams</div><div>IGES/COLA</div><div>4041 Powder Mill Road, Suite 302</div><div>Calverton, MD 20705</div><div><a href="mailto:jma@cola.iges.org">jma@cola.iges.org</a></div><div><br class="khtml-block-placeholder"></div><br class="Apple-interchange-newline"></span></span></div></span> </div><br></div></div></div></body></html>