Hi, I am getting strange behavior from ut_calendar when I run the
following on cgd systems for those interested). I am trying to test that
the date of a CAM file is January 2006.
inptr =
addfile("/project/convection/juliob/CAMruns/twdLL05/twdLL05.cam2.h1.2006-01-01-00000.nc","r")
print(ut_calendar(inptr->time, -3))
*Result
Variable: unnamed (return)
Type: integer
Total Size: 4 bytes
1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
Number Of Attributes: 1
calendar : noleap
(0) 2005123124
I would expect this to be 2006010100
Any idea why this is 24 at the end and not 00?
Below is an ncdump -h on the file.
The time units are a little odd
time:units = "days since 2005-01-01 12:00:00" ;
And I am wondering if ut_calendar is having difficulty with the '12'
Thanks!!
Rich
*********************************************
netcdf twdLL05.cam2.h1.2006-01-01-00000 {
dimensions:
lat = 96 ;
lon = 144 ;
slat = 95 ;
slon = 144 ;
lev = 30 ;
ilev = 31 ;
isccp_prs = 7 ;
isccp_tau = 7 ;
isccp_prstau = 49 ;
time = UNLIMITED ; // (1 currently)
tbnd = 2 ;
chars = 8 ;
variables:
double P0 ;
P0:long_name = "reference pressure" ;
P0:units = "Pa" ;
double lat(lat) ;
lat:long_name = "latitude" ;
lat:units = "degrees_north" ;
double lon(lon) ;
lon:long_name = "longitude" ;
lon:units = "degrees_east" ;
double slat(slat) ;
slat:long_name = "staggered latitude" ;
slat:units = "degrees_north" ;
double slon(slon) ;
slon:long_name = "staggered longitude" ;
slon:units = "degrees_east" ;
double w_stag(slat) ;
w_stag:long_name = "staggered latitude weights" ;
double lev(lev) ;
lev:long_name = "hybrid level at midpoints (1000*(A+B))" ;
lev:units = "level" ;
lev:positive = "down" ;
lev:standard_name =
"atmosphere_hybrid_sigma_pressure_coordinate" ;
lev:formula_terms = "a: hyam b: hybm p0: P0 ps: PS" ;
double ilev(ilev) ;
ilev:long_name = "hybrid level at interfaces
(1000*(A+B))" ;
ilev:units = "level" ;
ilev:positive = "down" ;
ilev:standard_name =
"atmosphere_hybrid_sigma_pressure_coordinate" ;
ilev:formula_terms = "a: hyai b: hybi p0: P0 ps: PS" ;
double isccp_prs(isccp_prs) ;
isccp_prs:long_name = "Mean ISCCP pressure" ;
isccp_prs:units = "mb" ;
isccp_prs:isccp_prs_bnds = 0., 180., 310., 440., 560.,
680., 800., 1000. ;
double isccp_tau(isccp_tau) ;
isccp_tau:long_name = "Mean ISCCP optical depth" ;
isccp_tau:units = "unitless" ;
isccp_tau:isccp_tau_bnds = 0., 0.3, 1.3, 3.6, 9.4, 23.,
60., 379. ;
double isccp_prstau(isccp_prstau) ;
isccp_prstau:long_name = "Mean pressure (mb).mean
optical depth (unitless)/1000" ;
isccp_prstau:units = "mixed" ;
double time(time) ;
time:long_name = "time" ;
time:units = "days since 2005-01-01 12:00:00" ;
time:calendar = "noleap" ;
time:bounds = "time_bnds" ;
double time_bnds(time, tbnd) ;
time_bnds:long_name = "time interval endpoints" ;
char date_written(time, chars) ;
char time_written(time, chars) ;
int ntrm ;
ntrm:long_name = "spectral truncation parameter M" ;
int ntrn ;
ntrn:long_name = "spectral truncation parameter N" ;
int ntrk ;
ntrk:long_name = "spectral truncation parameter K" ;
int ndbase ;
ndbase:long_name = "base day" ;
int nsbase ;
nsbase:long_name = "seconds of base day" ;
int nbdate ;
nbdate:long_name = "base date (YYYYMMDD)" ;
int nbsec ;
nbsec:long_name = "seconds of base date" ;
int mdt ;
mdt:long_name = "timestep" ;
mdt:units = "s" ;
int nlon(lat) ;
nlon:long_name = "number of longitudes" ;
int wnummax(lat) ;
wnummax:long_name = "cutoff Fourier wavenumber" ;
double hyai(ilev) ;
hyai:long_name = "hybrid A coefficient at layer
interfaces" ;
double hybi(ilev) ;
hybi:long_name = "hybrid B coefficient at layer
interfaces" ;
double hyam(lev) ;
hyam:long_name = "hybrid A coefficient at layer
midpoints" ;
double hybm(lev) ;
hybm:long_name = "hybrid B coefficient at layer
midpoints" ;
double gw(lat) ;
gw:long_name = "gauss weights" ;
int ndcur(time) ;
ndcur:long_name = "current day (from base day)" ;
int nscur(time) ;
nscur:long_name = "current seconds of current day" ;
int date(time) ;
date:long_name = "current date (YYYYMMDD)" ;
double co2vmr(time) ;
co2vmr:long_name = "co2 volume mixing ratio" ;
double ch4vmr(time) ;
ch4vmr:long_name = "ch4 volume mixing ratio" ;
double n2ovmr(time) ;
n2ovmr:long_name = "n2o volume mixing ratio" ;
double f11vmr(time) ;
f11vmr:long_name = "f11 volume mixing ratio" ;
double f12vmr(time) ;
f12vmr:long_name = "f12 volume mixing ratio" ;
double sol_tsi(time) ;
sol_tsi:long_name = "total solar irradiance" ;
sol_tsi:units = "W/m2" ;
int datesec(time) ;
datesec:long_name = "current seconds of current date" ;
int nsteph(time) ;
nsteph:long_name = "current timestep" ;
float OMEGA(time, lev, lat, lon) ;
OMEGA:units = "Pa/s" ;
OMEGA:long_name = "Vertical velocity (pressure)" ;
float PBLH(time, lat, lon) ;
PBLH:units = "m" ;
PBLH:long_name = "PBL height" ;
float PRECC(time, lat, lon) ;
PRECC:units = "m/s" ;
PRECC:long_name = "Convective precipitation rate (liq +
ice)" ;
float PRECL(time, lat, lon) ;
PRECL:units = "m/s" ;
PRECL:long_name = "Large-scale (stable) precipitation
rate (liq + ice)" ;
float PRECSH(time, lat, lon) ;
PRECSH:units = "m/s" ;
PRECSH:long_name = "Shallow Convection precipitation
rate" ;
float PS(time, lat, lon) ;
PS:units = "Pa" ;
PS:long_name = "Surface pressure" ;
float Q(time, lev, lat, lon) ;
Q:units = "kg/kg" ;
Q:long_name = "Specific humidity" ;
float T(time, lev, lat, lon) ;
T:units = "K" ;
T:long_name = "Temperature" ;
float U(time, lev, lat, lon) ;
U:units = "m/s" ;
U:long_name = "Zonal wind" ;
float V(time, lev, lat, lon) ;
V:units = "m/s" ;
V:long_name = "Meridional wind" ;
// global attributes:
:Conventions = "CF-1.0" ;
:source = "CAM" ;
:case = "twdLL05" ;
:title = "atm case twdLL" ;
:logname = "juliob" ;
:host = "be1005en.ucar.ed" ;
:Version = "$Name$" ;
:revision_Id = "$Id$" ;
:initial_file =
"/fs/cgd/csm/inputdata/atm/cam/inic/fv/cami-mam3_0000-01-01_1.9x2.5_L30_c090306.nc"
;
:topography_file =
"/fs/cgd/csm/inputdata/atm/cam/topo/USGS-gtopo30_1.9x2.5_remap_c050602.nc" ;
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Dec 13 14:07:23 2010
This archive was generated by hypermail 2.1.8 : Wed Dec 22 2010 - 16:10:23 MST