[Thread Prev][Thread Next][Index]

[ferret_users] calendar names lowercase, monthly noleap output



Hi,

Just a note on a problem of calendar names  that are written with uppercase.
It causes troubles with CDO (climate data operators)
http://www.mpimet.mpg.de/cdo

And indeed the CF convention mentions that one should use "noleap" and not "NOLEAP".
Same for other calendar names

Tested with last ferret v6.61
yes? def axis/t="01-JAN-2010":"31-DEC-2011":1/units="days"/cal="noleap" mytaxis
yes? save/clobber/file=file_noleap.nc t[gt=mytaxis]

netcdf file_noleap {
dimensions:
       MYTAXIS = UNLIMITED ; // (730 currently)
variables:
       double MYTAXIS(MYTAXIS) ;
               MYTAXIS:units = "DAYS since 1901-01-15 00:00:00" ;
               MYTAXIS:axis = "T" ;
               MYTAXIS:calendar = "NOLEAP" ;

By the way, it should be nice to find a way to define easily a
time axis with days units in a noleap calendar for monthy output
For now, I use the following lines

let year_ref=($01)
let year_start=($02)
let year_end=($03)

let nbyears=year_end-year_start+1
let nbdays={31,28,31,30,31,30,31,31,30,31,30,31}

DEFINE AXIS/UNITS="days"/T0="01-JAN-`year_ref`"/EDGES/CAL="noleap" monthaxis = `(year_start-year_ref)*365`+xcat(0,xsequence(365*(j[j=1:`nbyears`]-1)+nbdays[i=@rsum]))
show axis monthaxis

Patrick



[Thread Prev][Thread Next][Index]

Contact Us
Dept of Commerce / NOAA / OAR / PMEL / TMAP

Privacy Policy | Disclaimer | Accessibility Statement