odd ut_calendar behavior

From: Rich Neale <rneale_at_nyahnyahspammersnyahnyah>
Date: Mon Dec 13 2010 - 14:06:57 MST

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