[gradsusr] netcdf file has two times defined for the same timestep
Roberto Mera
rjmera at coas.oregonstate.edu
Tue Jun 28 14:42:27 EDT 2011
Hi. I'm trying to map variables from a netcdf file that has two separate
defined times for the exact same timestep (1 month). Using sdfopen
displays only two of the variables where there are 22. Using a ctl file
it maps out a constant value 0. I've isolated the problem as follows:
there is a time0 and time1 defined in the netcdf and two of the
variables have the former as its time but the rest correspond to the
latter. Below is the ncdump of the file. I'm using v2.0.a9.oga.1 on a
unix platform. Please advice.
netcdf zzzrlargescalenov98.nc
{
dimensions:
time0 = UNLIMITED ; // (1 currently)
time1 = 1 ;
z0 = 1 ;
z1 = 1 ;
z2 = 1 ;
z3 = 1 ;
z4 = 3 ;
z5 = 1 ;
z6 = 1 ;
latitude0 = 145 ;
longitude0 = 192 ;
latitude1 = 144 ;
longitude1 = 192 ;
variables:
float time0(time0) ;
time0:standard_name = "time" ;
time0:axis = "T" ;
time0:calendar = "360_day" ;
time0:units = "days since 1959-12-01 00:00:00" ;
time0:meaning_period = "720 hours" ;
time0:meaning_sampling_frequency = "3 hours" ;
float time1(time1) ;
time1:standard_name = "time" ;
time1:axis = "T" ;
time1:calendar = "360_day" ;
time1:units = "days since 1959-12-01 00:00:00" ;
time1:meaning_period = "720 hours" ;
time1:meaning_sampling_frequency = "1 hours" ;
float z0(z0) ;
z0:standard_name = "height" ;
z0:units = "m" ;
z0:direction = "up" ;
float z1(z1) ;
float z2(z2) ;
z2:standard_name = "atmosphere_hybrid_sigma_pressure_coordinate" ;
z2:direction = "down" ;
float z3(z3) ;
z3:standard_name = "height" ;
z3:units = "m" ;
z3:direction = "up" ;
float z4(z4) ;
z4:standard_name = "air_pressure" ;
z4:units = "hPa" ;
z4:direction = "down" ;
float z5(z5) ;
z5:standard_name = "air_pressure" ;
z5:units = "hPa" ;
z5:direction = "down" ;
float z6(z6) ;
z6:standard_name = "height" ;
z6:units = "m" ;
z6:direction = "up" ;
float latitude0(latitude0) ;
latitude0:standard_name = "latitude" ;
latitude0:units = "degrees_north" ;
latitude0:axis = "Y" ;
float longitude0(longitude0) ;
longitude0:standard_name = "longitude" ;
longitude0:units = "degrees_east" ;
longitude0:axis = "X" ;
float latitude1(latitude1) ;
latitude1:standard_name = "latitude" ;
latitude1:units = "degrees_north" ;
latitude1:axis = "Y" ;
float longitude1(longitude1) ;
longitude1:standard_name = "longitude" ;
longitude1:units = "degrees_east" ;
longitude1:axis = "X" ;
float field186(time0, z0, latitude0, longitude0) ;
field186:_FillValue = -1.073742e+09f ;
field186:stash_item = "201" ;
field186:stash_section = "1" ;
field186:field_code = "186" ;
field186:standard_name = "surface_net_downward_shortwave_flux" ;
field186:long_name = "NET DOWN SURFACE SW FLUX: SW TS ONLY" ;
field186:units = "W m-2" ;
field186:cell_method = "time: mean " ;
float field201(time0, z1, latitude0, longitude0) ;
field201:_FillValue = -1.073742e+09f ;
field201:stash_item = "208" ;
field201:stash_section = "1" ;
field201:field_code = "201" ;
field201:standard_name = "toa_outgoing_shortwave_flux" ;
field201:long_name = "OUTGOING SW RAD FLUX (TOA)" ;
field201:units = "W m-2" ;
field201:cell_method = "time: mean " ;
float field207(time0, z1, latitude0, longitude0) ;
field207:_FillValue = -1.073742e+09f ;
field207:stash_item = "209" ;
field207:stash_section = "1" ;
field207:field_code = "207" ;
field207:standard_name =
"toa_outgoing_shortwave_flux_assuming_clear_sky" ;
field207:long_name = "CLEAR-SKY (II) UPWARD SW FLUX (TOA)" ;
field207:units = "W m-2" ;
field207:cell_method = "time: mean " ;
float field187(time0, z0, latitude0, longitude0) ;
field187:_FillValue = -1.073742e+09f ;
field187:stash_item = "201" ;
field187:stash_section = "2" ;
field187:field_code = "187" ;
field187:standard_name = "surface_net_downward_longwave_flux" ;
field187:long_name = "NET DOWN SURFACE LW RAD FLUX" ;
field187:units = "W m-2" ;
field187:cell_method = "time: mean " ;
float field206(time0, z1, latitude0, longitude0) ;
field206:_FillValue = -1.073742e+09f ;
field206:stash_item = "205" ;
field206:stash_section = "2" ;
field206:field_code = "206" ;
field206:standard_name = "toa_net_upward_longwave_flux" ;
field206:long_name = "OUTGOING LW RAD FLUX (TOA)" ;
field206:units = "W m-2" ;
field206:cell_method = "time: mean " ;
float field210(time0, z1, latitude0, longitude0) ;
field210:_FillValue = -1.073742e+09f ;
field210:stash_item = "206" ;
field210:stash_section = "2" ;
field210:field_code = "210" ;
field210:standard_name =
"toa_net_upward_longwave_flux_assuming_clear_sky" ;
field210:long_name = "CLEAR-SKY (II) UPWARD LW FLUX (TOA)" ;
field210:units = "W m-2" ;
field210:cell_method = "time: mean " ;
float field178(time1, z2, latitude0, longitude0) ;
field178:_FillValue = -1.073742e+09f ;
field178:stash_item = "217" ;
field178:stash_section = "3" ;
field178:field_code = "178" ;
field178:standard_name = "surface_upward_sensible_heat_flux" ;
field178:long_name = "SURFACE B.LAYER HEAT FLUXES W/M2" ;
field178:units = "W m-2" ;
field178:cell_method = "time: mean " ;
float field56(time1, z3, latitude1, longitude1) ;
field56:_FillValue = -1.073742e+09f ;
field56:stash_item = "225" ;
field56:stash_section = "3" ;
field56:field_code = "56" ;
field56:standard_name = "eastward_wind" ;
field56:long_name = "10 METRE WIND U-COMP" ;
field56:units = "m s-1" ;
field56:cell_method = "time: mean " ;
float field57(time1, z3, latitude1, longitude1) ;
field57:_FillValue = -1.073742e+09f ;
field57:stash_item = "226" ;
field57:stash_section = "3" ;
field57:field_code = "57" ;
field57:standard_name = "northward_wind" ;
field57:long_name = "10 METRE WIND V-COMP" ;
field57:units = "m s-1" ;
field57:cell_method = "time: mean " ;
float field180(time1, z0, latitude0, longitude0) ;
field180:_FillValue = -1.073742e+09f ;
field180:stash_item = "234" ;
field180:stash_section = "3" ;
field180:field_code = "180" ;
field180:standard_name = "surface_upward_latent_heat_flux" ;
field180:long_name = "SURFACE LATENT HEAT FLUX W/M2" ;
field180:units = "W m-2" ;
field180:cell_method = "time: mean " ;
float field16(time1, z3, latitude0, longitude0) ;
field16:_FillValue = -1.073742e+09f ;
field16:stash_item = "236" ;
field16:stash_section = "3" ;
field16:field_code = "16" ;
field16:standard_name = "air_temperature" ;
field16:long_name = "TEMPERATURE AT 1.5M" ;
field16:units = "K" ;
field16:cell_method = "time: mean " ;
float field88(time1, z3, latitude0, longitude0) ;
field88:_FillValue = -1.073742e+09f ;
field88:stash_item = "245" ;
field88:stash_section = "3" ;
field88:field_code = "88" ;
field88:standard_name = "relative_humidity" ;
field88:long_name = "RELATIVE HUMIDITY AT 1.5M" ;
field88:units = "1" ;
field88:cell_method = "time: mean " ;
float field17(time1, z3, latitude0, longitude0) ;
field17:_FillValue = -1.073742e+09f ;
field17:stash_item = "250" ;
field17:stash_section = "3" ;
field17:field_code = "17" ;
field17:standard_name = "dew_point_temperature" ;
field17:long_name = "DEWPOINT AT 1.5M" ;
field17:units = "K" ;
field17:cell_method = "time: mean " ;
float field90(time1, z0, latitude0, longitude0) ;
field90:_FillValue = -1.073742e+09f ;
field90:stash_item = "216" ;
field90:stash_section = "5" ;
field90:field_code = "90" ;
field90:standard_name = "precipitation_flux" ;
field90:long_name = "TOTAL PRECIPITATION RATE KG/M2/S" ;
field90:units = "kg m-2 s-1" ;
field90:cell_method = "time: mean " ;
float field106(time1, z0, latitude0, longitude0) ;
field106:_FillValue = -1.073742e+09f ;
field106:stash_item = "208" ;
field106:stash_section = "8" ;
field106:field_code = "106" ;
field106:standard_name = "liquid_water_content_of_soil_layer" ;
field106:long_name = "SOIL MOISTURE CONTENT" ;
field106:units = "kg m-2" ;
field106:cell_method = "time: mean " ;
float field56_1(time1, z4, latitude1, longitude1) ;
field56_1:_FillValue = -1.073742e+09f ;
field56_1:stash_item = "201" ;
field56_1:stash_section = "15" ;
field56_1:field_code = "56" ;
field56_1:standard_name = "eastward_wind" ;
field56_1:long_name = "U COMPNT OF WIND ON PRESSURE LEVELS" ;
field56_1:units = "m s-1" ;
field56_1:cell_method = "time: mean " ;
float field57_1(time1, z4, latitude1, longitude1) ;
field57_1:_FillValue = -1.073742e+09f ;
field57_1:stash_item = "202" ;
field57_1:stash_section = "15" ;
field57_1:field_code = "57" ;
field57_1:standard_name = "northward_wind" ;
field57_1:long_name = "V COMPNT OF WIND ON PRESSURE LEVELS" ;
field57_1:units = "m s-1" ;
field57_1:cell_method = "time: mean " ;
float field1(time1, z5, latitude0, longitude0) ;
field1:_FillValue = -1.073742e+09f ;
field1:stash_item = "202" ;
field1:stash_section = "16" ;
field1:field_code = "1" ;
field1:standard_name = "geopotential_height" ;
field1:long_name = "GEOPOTENTIAL HEIGHT: PRESSURE LEVELS" ;
field1:units = "m" ;
field1:cell_method = "time: mean " ;
float field16_1(time1, z5, latitude0, longitude0) ;
field16_1:_FillValue = -1.073742e+09f ;
field16_1:stash_item = "203" ;
field16_1:stash_section = "16" ;
field16_1:field_code = "16" ;
field16_1:standard_name = "air_temperature_at_pressure_level" ;
field16_1:long_name = "TEMPERATURE ON PRESSURE LEVELS" ;
field16_1:units = "K" ;
field16_1:cell_method = "time: mean " ;
float field88_1(time1, z5, latitude0, longitude0) ;
field88_1:_FillValue = -1.073742e+09f ;
field88_1:stash_item = "204" ;
field88_1:stash_section = "16" ;
field88_1:field_code = "88" ;
field88_1:standard_name = "relative_humidity_at_pressure_level" ;
field88_1:long_name = "RELATIVE HUMIDITY ON PRESSURE LEVELS" ;
field88_1:units = "1" ;
field88_1:cell_method = "time: mean " ;
float field8(time1, z6, latitude0, longitude0) ;
field8:_FillValue = -1.073742e+09f ;
field8:stash_item = "222" ;
field8:stash_section = "16" ;
field8:field_code = "8" ;
field8:standard_name = "air_pressure_at_sea_level" ;
field8:long_name = "PRESSURE AT MEAN SEA LEVEL" ;
field8:units = "Pa" ;
field8:cell_method = "time: mean " ;
float field16_2(time1, z0, latitude0, longitude0) ;
field16_2:_FillValue = -1.073742e+09f ;
field16_2:stash_item = "24" ;
field16_2:stash_section = "0" ;
field16_2:field_code = "16" ;
field16_2:standard_name = "surface_temperature" ;
field16_2:long_name = "SURFACE TEMPERATURE AFTER TIMESTEP" ;
field16_2:units = "K" ;
field16_2:cell_method = "time: mean " ;
// global attributes:
:Conventions = "CF-1.3" ;
:alpham = 0.5f ;
}
Appending that, here is the ctl file:
DSET ^zzzrlargescalenov98.nc
dtype netcdf
TITLE Pacific northwest
UNDEF -99999
xdef 192 linear 0 1.875
ydef 145 linear -90 1.25
ZDEF 1 LINEAR 1 1
TDEF 1 LINEAR 0Z01OCT2006 1MO
VARS 3
field186=>field186 1 T,Z,Y,X surface_downwelling_shortwave_flux (W m-2)
field201=>field201 1 T,Z,Y,X precipitation_flux (kg/m2/s)
field16_2=>field16_2 1 T,Z,Y,X surface_temperature (K)
ENDVARS
-Robert
More information about the gradsusr
mailing list