[gradsusr] netcdf file has two times defined for the same timestep

Charles Seman Charles.Seman at noaa.gov
Tue Jun 28 15:08:58 EDT 2011


Roberto,

Try extracting the variable of interest into a separate file (if you 
have access to NCO utilities, use "ncks")...

I hope this helps,
Chuck

Roberto Mera wrote:
> 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
> 
> 
> 
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr

-- 

Please note that Charles.Seman at noaa.gov should be considered my NOAA
email address, not cjs at gfdl.noaa.gov.

********************************************************************
  Charles Seman                                Charles.Seman at noaa.gov
  U.S. Department of Commerce / NOAA / OAR
  Geophysical Fluid Dynamics Laboratory         voice: (609) 452-6547
  201 Forrestal Road                              fax: (609) 987-5063
  Princeton, NJ  08540-6649            http://www.gfdl.noaa.gov/~cjs/
********************************************************************

"The contents of this message are mine personally and do not reflect any
official or unofficial position of the United States Federal Government,
the United States Department of Commerce, or NOAA."



More information about the gradsusr mailing list