Hi
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.
Cheers
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:tim.johns@xxxxxxxxxxxxxxxx" ;
: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 ;
}
