This archive contains answers to questions sent to Unidata support through mid-2025. Note that the archive is no longer being updated. We provide the archive for reference; many of the answers presented here remain technically correct, even if somewhat outdated. For the most up-to-date information on the use of NSF Unidata software and data services, please consult the Software Documentation first.
Hi Harry, > I use the Fortran-77 library to access NetCDF files. don't ask me why > - I must have done something bad in a previous life. > > Anyway, I'm currently working with GCM data from CERA, in particular > the new ENSEMBLES E1 stabilization runs. > > All working fine, until I get to the HadGEM2 runs: > > 1. Program asks for the variable ID for 'lat' - gets '3'. > > 2. Program asks for the data from variable '3' - gets data for > time_bnds! > > If I do an ncdump, lat is the FOURTH variable. But the first is > 'height', which is scalar. Output from 'ncdump -h follows. > > Is the F77 library ignoring scalar variables? Program working fine > for all other models. I can't reproduce the problem you are seeing. When I create the netCDF file using ncgen and your CDL output, then run ncdump on the result, it shows everything is read correctly. A little program to print the variable ID for 'lat' returns 4, as it should: program kzd implicit none include 'netcdf.inc' character*(*) FILE_NAME parameter (FILE_NAME='HADGEM2_E1_1_MM_pr_1-1201.nc') integer ncid, varid, retval retval = nf_open(FILE_NAME, NF_NOWRITE, ncid) if (retval .ne. nf_noerr) call handle_err(retval) retval = nf_inq_varid(ncid, 'lat', varid) if (retval .ne. nf_noerr) call handle_err(retval) print *, 'lat id = ', varid retval = nf_close(ncid) if (retval .ne. nf_noerr) call handle_err(retval) end subroutine handle_err(errcode) implicit none include 'netcdf.inc' integer errcode print *, 'Error: ', nf_strerror(errcode) stop 2 end --Russ > Harry > Ian "Harry" Harris > Climatic Research Unit > School of Environmental Sciences > University of East Anglia > Norwich NR4 7TJ > United Kingdom > > ncdump -h HADGEM2_E1_1_MM_pr_1-1201.nc > netcdf HADGEM2_E1_1_MM_pr_1-1201 { > dimensions: > time = UNLIMITED ; // (1201 currently) > bound = 2 ; > lat = 145 ; > lon = 192 ; > variables: > float height ; > height:missing_value = 1.e+20f ; > height:name = "height" ; > height:positive = "up" ; > height:long_name = "height" ; > height:standard_name = "height" ; > height:units = "m" ; > height:axis = "Z" ; > float time(time) ; > time:bounds = "time_bnds" ; > time:long_name = "Time" ; > time:standard_name = "time" ; > time:units = "days since 1970-01-01 00:00:00" ; > time:calendar = "360_day" ; > time:axis = "T" ; > float time_bnds(time, bound) ; > float lat(lat) ; > lat:bounds = "lat_bnds" ; > lat:point_spacing = "even" ; > lat:long_name = "Latitude" ; > lat:standard_name = "latitude" ; > lat:units = "degrees_north" ; > lat:axis = "Y" ; > float lat_bnds(lat, bound) ; > float lon(lon) ; > lon:bounds = "lon_bnds" ; > lon:point_spacing = "even" ; > lon:long_name = "Longitude" ; > lon:standard_name = "longitude" ; > lon:axis = "X" ; > lon:units = "degrees_east" ; > lon:modulo = 360.f ; > lon:topology = "circular" ; > float lon_bnds(lon, bound) ; > float pr(time, lat, lon) ; > pr:comment = "Small negative data values generated > by the model are reset to zero" ; > pr:long_name = "Precipitation [time:mean; Mean over > one month sampling all time steps.]" ; > pr:standard_name = "precipitation_flux" ; > pr:missing_value = -1.073742e+09f ; > pr:coordinates = "height" ; > pr:units = "kg m-2 s-1" ; > pr:stash_section = 5 ; > pr:stash_item = 216 ; > pr:stash_model = 1 ; > pr:cell_methods = "time:mean" ; > pr:name = "Precipitation" ; > pr:lookup_source = "defaults (cdunifpp V0.11)" ; > pr:stash_mapping = "m01s05i216" ; > > // global attributes: > :Conventions = "CF-1.0" ; > :mean_type = "Mean over one month sampling all time > steps." ; > :title = "ENSEMBLES Stream Two:ense1_1:mm" ; > :institution = "Met Office Hadley Centre" ; > :source = "HadGEM2AO" ; > :contact = "Tim Johns:address@hidden" ; > :history = "Thu Feb 26 08:36:10 2009: /usr/local/ > badc/linux/suse10/nco/bin/ncrcat > METO_HadGEM2AO_ense1_1_mm_pr_199912010000-200912010000.nc > METO_HadGEM2AO_ense1_1_mm_pr_200912010000-201912010000.nc > METO_HadGEM2AO_ense1_1_mm_pr_201912010000-202912010000.nc > METO_HadGEM2AO_ense1_1_mm_pr_202912010000-203912010000.nc > METO_HadGEM2AO_ense1_1_mm_pr_203912010000-204912010000.nc > METO_HadGEM2AO_ense1_1_mm_pr_204912010000-205912010000.nc > METO_HadGEM2AO_ense1_1_mm_pr_205912010000-206912010000.nc > METO_HadGEM2AO_ense1_1_mm_pr_206912010000-207912010000.nc > METO_HadGEM2AO_ense1_1_mm_pr_207912010000-208912010000.nc > METO_HadGEM2AO_ense1_1_mm_pr_208912010000-209912010000.nc > METO_HadGEM2AO_ense1_1_mm_pr_209912010000-210912010000.nc > METO_HadGEM2AO_ense1_1_mm_pr.nc\n", > "Exported from MetO MaitreD/MASS data > archive to BADC in native PP format: Converted to netCDF at BADC > using BADC HadGEM Virtual Variable Suite version 0.2.1" ; > :comment = "Start from a point (Dec 1 1999) taken > from the end of the 20c3m_1 experiment and run to the end of the 21st > century. Anthropogenic forcings were changed through the run to > reflect the ENSEMBLES E1 scenario, i.e. changes in: CO2, CH4, N2O, > CFC11*, CFC12, Ozone, Vegetation Distribution, Sulphur-cycle > Emissions, Soot Emissions, Biomass Emissions. No variations in solar > radiation and volcanic aerosols. (In this experiment CFC11* or > equivalent-CFC11 is used to represent other halocarbon gases apart > from CFC12)." ; > :references = "Collins et al. 2008, Hadley Centre > Technical Note no. 74: Bellouin et al. 2007, Hadley Centre Technical > Note no. 73: Johns et al. 2006, J. Climate, 19(7), 1327-1353: Martin > et al. 2006, J. Climate, 19(7), 1274-1301: Ringer et al. 2006, J. > Climate, 19(7), 1302-1326: Stott et al. 2006, J. Climate, 19(12), > 2763-2782: Royer et al. 2008, Report No: D2A.0.4; Project: GOCE- > CT-2003-505539: Lowe et al. 2008, Submitted to EOS: McLaren et al. > 2006, J. Geophys. Res., 111, C12014" ; > :nco_openmp_thread_number = 1 ; > } > > _______________________________________________ > netcdfgroup mailing list > address@hidden > For list information or to unsubscribe, visit: > http://www.unidata.ucar.edu/mailing_lists/ > > Russ Rew UCAR Unidata Program address@hidden http://www.unidata.ucar.edu Ticket Details =================== Ticket ID: KZD-970215 Department: Support netCDF Priority: Normal Status: Closed