problem opening netcdf file

Jennifer Adams jma at COLA.IGES.ORG
Wed May 31 11:00:55 EDT 2006


Lee,
You also need to read the documentation on using NetCDF data files with
GrADS:
http://www.iges.org/grads/gadoc/SDFdescriptorfile.html

You might get xdfopen to work if you put an XDEF, YDEF, ZDEF, and TDEF
entries in your xdf-type descriptor. The YDEF   entry would have to
borrow a name from one of the other coordinate variables since none
exists for the Y dimension. Something like:
ydef lat 1 linear 0 1

Jennifer

On May 31, 2006, at 10:19 AM, Grenfell Lee wrote:

> Dear Grads enthusiasts,
>
> I have a problem opening a netcdf file in grads - I have tried both
> sdfopen and xdfopen.
>
> Details below.
>
> Many thanks,
>
> Lee
>
> Dr John Lee Grenfell
> German Space Agency (DLR)
> Institute for Planetary Research
> Rutherfordstr. 2
> 12489 Berlin
> Germany
> Tel: +49 30 67055-337
>
> My Version of Grads
> --------------------
>  o This is a 32-bit LITTLE ENDIAN machine version.
>  o Command line editing (readline) ENABLED.
>  o CIRES/CDC (http://www.cdc.noaa.gov) SDF/XDF interface ENABLED.
>      Use sdfopen/xdfopen to read NetCDF files.
>  o DTYPE netcdf is ENABLED; DTYPE hdfsds is DISABLED.
>  o OPeNDAP (a.k.a. DODS) gridded data interface DISABLED.
>  o OPeNDAP (a.k.a. DODS) station data interface DISABLED.
>  o PCMDI (http://www-pcmdi.llnl.gov) LATS interface ENABLED.
>      This version is configured to write GRIB and NetCDF files.
>  o DAO (http://dao.gsfc.nasa.gov) Athena Widget GUI DISABLED.
>  o NRL/DAO/PCMDI XA or ImageMagick Image Output DISABLED.
>  o printim command for direct png/gif output ENABLED.
>      (via the GD Library -- http://www.boutell.com/gd)
>
> My Platform
> -----------
> OTS44 Cluster Machine
>
>
> Description of Problem
> ----------------------
>
> When I try to open my netcdf file I get the error message:
>
> with sdfopen:
> ga-> sdfopen arch.10yr.base.nc
> Scanning self-describing file:  arch.10yr.base.nc
> SDF file has no discernable X coordinate.
>
> or, with xdfopen:
> ga-> xdfopen arch.10yr.base.nc
> Scanning Descriptor File:  arch.10yr.base.nc
> XDFOpen Error:  Unknown keyword in description file
>   --> The invalid description file record is:
>   --> cdf
>   The data file was not opened.
>
>
> Output from ncdump -c
> ---------------------
>
> netcdf arch.10yr.base {
> dimensions:
>         lat = 35 ;
>         lev = 121 ;
>         time = UNLIMITED ; // (147 currently)
> variables:
>         float latitudes(lat) ;
>                 latitudes:units = "angular degrees" ;
>         float levels(lev) ;
>                 levels:units = "kilometers" ;
>         float time(time) ;
>                 time:units = "days" ;
>         float temperature(time, lev, lat) ;
>                 temperature:units = "K" ;
>         float PotTemp(time, lev, lat) ;
>                 PotTemp:units = "K" ;
>         float Tcloud(time, lev, lat) ;
>                 Tcloud:units = "K" ;
>         float Tobs(time, lev, lat) ;
>                 Tobs:units = "K" ;
>         float totdens(time, lev, lat) ;
>                 totdens:units = "molec/cm3" ;
>         float u(time, lev, lat) ;
>                 u:units = "m/s" ;
>         float v(time, lev, lat) ;
>                 v:units = "m/s" ;
>         float w(time, lev, lat) ;
>                 w:units = "m/s" ;
>         float cf(time, lev, lat) ;
>                 cf:units = "s-3" ;
>         float chi(time, lev, lat) ;
>                 chi:units = "m2/s" ;
>         float chiM(time, lev, lat) ;
>                 chiM:units = "kg/m/s" ;
>         float q(time, lev, lat) ;
>                 q:units = "K/day" ;
>         float SolHeat(time, lev, lat) ;
>                 SolHeat:units = "K/day" ;
>         float ChemHeat(time, lev, lat) ;
>                 ChemHeat:units = "K/day" ;
>         float LatHeat(time, lev, lat) ;
>                 LatHeat:units = "K/day" ;
>         float qcool(time, lev, lat) ;
>                 qcool:units = "K/day" ;
>         float fx(time, lev, lat) ;
>                 fx:units = "m/s/day" ;
>         float kzz(time, lev, lat) ;
>                 kzz:units = "m2/s" ;
>         float fr(time, lev, lat) ;
>                 fr:units = "m/s/day" ;
>         float kyy(time, lev, lat) ;
>                 kyy:units = "m2/s" ;
>         float pw1_amp(time, lev, lat) ;
>                 pw1_amp:units = "m" ;
>         float pw2_amp(time, lev, lat) ;
>                 pw2_amp:units = "m" ;
>         float abar1(time, lev, lat) ;
>                 abar1:units = "cm2/cm3" ;
>         float abar2(time, lev, lat) ;
>                 abar2:units = "cm2/cm3" ;
>         float abar1pmc(time, lev, lat) ;
>                 abar1pmc:units = "cm2/cm3" ;
>         float abar2pmc(time, lev, lat) ;
>                 abar2pmc:units = "cm2/cm3" ;
>         float h2oTot(time, lev, lat) ;
>                 h2oTot:units = "vmr" ;
>         float h2oSat(time, lev, lat) ;
>                 h2oSat:units = "vmr" ;
>         float h2oCond(time, lev, lat) ;
>                 h2oCond:units = "vmr" ;
>         float h2oCond_sed(time, lev, lat) ;
>                 h2oCond_sed:units = "vmr" ;
>         float RH(time, lev, lat) ;
>                 RH:units = "%" ;
>         float SH(time, lev, lat) ;
>                 SH:units = "g h2o/kg air" ;
>         float n2o(time, lev, lat) ;
>                 n2o:units = "vmr" ;
>         float ch4(time, lev, lat) ;
>                 ch4:units = "vmr" ;
>         float h2o(time, lev, lat) ;
>                 h2o:units = "vmr" ;
>         float co(time, lev, lat) ;
>                 co:units = "vmr" ;
>         float hcl(time, lev, lat) ;
>                 hcl:units = "vmr" ;
>         float hbr(time, lev, lat) ;
>                 hbr:units = "vmr" ;
>         float co2(time, lev, lat) ;
>                 co2:units = "vmr" ;
>         float h2(time, lev, lat) ;
>                 h2:units = "vmr" ;
>         float ch2o(time, lev, lat) ;
>                 ch2o:units = "vmr" ;
>         float cfc10(time, lev, lat) ;
>                 cfc10:units = "vmr" ;
>         float cfc11(time, lev, lat) ;
>                 cfc11:units = "vmr" ;
>         float cfc12(time, lev, lat) ;
>                 cfc12:units = "vmr" ;
>         float cfc113(time, lev, lat) ;
>                 cfc113:units = "vmr" ;
>         float cfc114(time, lev, lat) ;
>                 cfc114:units = "vmr" ;
>         float cfc115(time, lev, lat) ;
>                 cfc115:units = "vmr" ;
>         float hcfc22(time, lev, lat) ;
>                 hcfc22:units = "vmr" ;
>         float ch3ccl3(time, lev, lat) ;
>                 ch3ccl3:units = "vmr" ;
>         float ch3cl(time, lev, lat) ;
>                 ch3cl:units = "vmr" ;
>         float ha1211(time, lev, lat) ;
>                 ha1211:units = "vmr" ;
>         float ha1301(time, lev, lat) ;
>                 ha1301:units = "vmr" ;
>         float hf(time, lev, lat) ;
>                 hf:units = "vmr" ;
>         float ch3br(time, lev, lat) ;
>                 ch3br:units = "vmr" ;
>         float chbr3(time, lev, lat) ;
>                 chbr3:units = "vmr" ;
>         float o3(time, lev, lat) ;
>                 o3:units = "vmr" ;
>         float o3p(time, lev, lat) ;
>                 o3p:units = "vmr" ;
>         float o1d(time, lev, lat) ;
>                 o1d:units = "vmr" ;
>         float h(time, lev, lat) ;
>                 h:units = "vmr" ;
>         float oh(time, lev, lat) ;
>                 oh:units = "vmr" ;
>         float ho2(time, lev, lat) ;
>                 ho2:units = "vmr" ;
>         float h2o2(time, lev, lat) ;
>                 h2o2:units = "vmr" ;
>         float cl(time, lev, lat) ;
>                 cl:units = "vmr" ;
>         float cl2(time, lev, lat) ;
>                 cl2:units = "vmr" ;
>         float clo(time, lev, lat) ;
>                 clo:units = "vmr" ;
>         float oclo(time, lev, lat) ;
>                 oclo:units = "vmr" ;
>         float hocl(time, lev, lat) ;
>                 hocl:units = "vmr" ;
>         float clono2(time, lev, lat) ;
>                 clono2:units = "vmr" ;
>         float clno2(time, lev, lat) ;
>                 clno2:units = "vmr" ;
>         float cl2o2(time, lev, lat) ;
>                 cl2o2:units = "vmr" ;
>         float n(time, lev, lat) ;
>                 n:units = "vmr" ;
>         float no(time, lev, lat) ;
>                 no:units = "vmr" ;
>         float no2(time, lev, lat) ;
>                 no2:units = "vmr" ;
>         float no3(time, lev, lat) ;
>                 no3:units = "vmr" ;
>         float hno3(time, lev, lat) ;
>                 hno3:units = "vmr" ;
>         float n2o5(time, lev, lat) ;
>                 n2o5:units = "vmr" ;
>         float ho2no2(time, lev, lat) ;
>                 ho2no2:units = "vmr" ;
>         float br(time, lev, lat) ;
>                 br:units = "vmr" ;
>         float bro(time, lev, lat) ;
>                 bro:units = "vmr" ;
>         float hobr(time, lev, lat) ;
>                 hobr:units = "vmr" ;
>         float brono2(time, lev, lat) ;
>                 brono2:units = "vmr" ;
>         float brcl(time, lev, lat) ;
>                 brcl:units = "vmr" ;
>         float ch3o2(time, lev, lat) ;
>                 ch3o2:units = "vmr" ;
>         float ch3ooh(time, lev, lat) ;
>                 ch3ooh:units = "vmr" ;
>         float o2(time, lev, lat) ;
>                 o2:units = "vmr" ;
>         float ClOx(time, lev, lat) ;
>                 ClOx:units = "vmr" ;
>         float Cly(time, lev, lat) ;
>                 Cly:units = "vmr" ;
>         float Bry(time, lev, lat) ;
>                 Bry:units = "vmr" ;
>         float NOx(time, lev, lat) ;
>                 NOx:units = "vmr" ;
>         float NOy(time, lev, lat) ;
>                 NOy:units = "vmr" ;
>         float totH(time, lev, lat) ;
>                 totH:units = "vmr" ;
>
> // global attributes:
>                 :lat = "latitudes" ;
>                 :lev = "levels" ;
>                 :model_name = "SOCRATES" ;
>                 :run_label = "10 year earth control on DLR cluster" ;
>                 :nlon = "diurnal_average" ;
> data:
>
>  time = 730366, 730391, 730416, 730441, 730466, 730491, 730516, 730541,
>     730566, 730591, 730616, 730641, 730666, 730691, 730716, 730741,
> 730766,
>     730791, 730816, 730841, 730866, 730891, 730916, 730941, 730966,
> 730991,
>     731016, 731041, 731066, 731091, 731116, 731141, 731166, 731191,
> 731216,
>     731241, 731266, 731291, 731316, 731341, 731366, 731391, 731416,
> 731441,
>     731466, 731491, 731516, 731541, 731566, 731591, 731616, 731641,
> 731666,
>     731691, 731716, 731741, 731766, 731791, 731816, 731841, 731866,
> 731891,
>     731916, 731941, 731966, 731991, 732016, 732041, 732066, 732091,
> 732116,
>     732141, 732166, 732191, 732216, 732241, 732266, 732291, 732316,
> 732341,
>     732366, 732391, 732416, 732441, 732466, 732491, 732516, 732541,
> 732566,
>     732591, 732616, 732641, 732666, 732691, 732716, 732741, 732766,
> 732791,
>     732816, 732841, 732866, 732891, 732916, 732941, 732966, 732991,
> 733016,
>     733041, 733066, 733091, 733116, 733141, 733166, 733191, 733216,
> 733241,
>     733266, 733291, 733316, 733341, 733366, 733391, 733416, 733441,
> 733466,
>     733491, 733516, 733541, 733566, 733591, 733616, 733641, 733666,
> 733691,
>     733716, 733741, 733766, 733791, 733816, 733841, 733866, 733891,
> 733916,
>     733941, 733966, 733991, 734016 ;
> }
>
>
--
Jennifer Miletta Adams
Center for Ocean-Land-Atmosphere Studies (COLA)
4041 Powder Mill Road, Suite 302
Calverton, MD 20705 USA
jma at cola.iges.org



More information about the gradsusr mailing list