Re: [ferret] climatological day

[Ansley Manke <Ansley.B.Manke@xxxxxxxx>]

Hi Patrick,
Perhaps somebody will have a more elegant idea, but it seems that it is
required that a calendar axis be defined with a month, day, and year;
and you need the axis to be a calendar axis in order to regrid your data
to it and get the modulo days.   You can at least customize the labels on
the plot, so that the date associated with the climatological day won't be visible.

yes? def axis/T="1-JAN-0000:00:00:00":"1-JAN-0000:23:30:00":0.5/units=hours/modulo/calendar=NOLEAP
my_day
yes? let all_feb = my_variable[t="1-feb-1987:00:00":28-feb-1987:23:30"]
yes? let clim_feb = all_feb[gt=my_day@MOD]

yes? plot/set clim_feb
yes? ppl txlint, 2, 0     ! Label every second hour, and do not label the day
yes? ppl plot

Patrick Brockmann wrote:

> Hi all,
>
> Is there a way to define a climatological day ?
>
> def axis/T="1-JAN-0000:00:00:00":"1-JAN-0000:23:30:00":0.5/units=hours/modulo/calendar=NOLEAP my_day
> could do the thing but is refrenced to a specific day, here for exemple the 1st of january.
>
> Some info from my variable
> yes? show axis MY_TAXIS
>  name       axis              # pts   start                end
>  MY_TAXIS  TIME             17520 r   01-JAN-1987 00:00    31-DEC-1987 23:30
> T0 = 15-JAN-1901
> CALENDAR = NOLEAP
>    Axis span (to cell edges) = 8760
>
> My purpose is to calculate the climatological day from a selected month
>
> Thanks for any help
> Patrick
>
> --
>
> mailto:Patrick.Brockmann@ipsl.jussieu.fr
> 01.44.27.21.10 --> IPSL (Jussieu, Tour 26, 4eme)
> 01.69.08.40.09 --> LSCE (Orme des merissiers 706, Bureau 3)
> http://www.ipsl.jussieu.fr/~brocksce/