Re: [ferret_users] complete seasonal cycle

["william s. kessler" <william.s.kessler@xxxxxxxx>]

Hi Pierre -

If the axis of "rain" is modulo (wraps around the year), then plot  
more than 1 year, while restricting the plot limits:

plot/set/l=0:13 rain     ! plot 14 months
ppl time,w000001010000,w000012312359   ! restrict the plot axes to  
one year
ppl txlint,1,0     ! label only months, not the year
ppl plot

The PPL TIME statement acts like /HLIMITS/VLIMITS for non-time axes;  
it restricts the plot but not the data plotted.
It's syntax is PPL TIME,wyyyymmddhhmm, wyyyymmddhhmm

To do January to January, I think this would work:

plot/set/l=1:13 rain     ! plot 13 months
ppl time,w000001150000,w000101150000   ! restrict the plot axes to 15- 
jan-0000 to 15-jan-0001
ppl txlint,1,0     ! label only months, not the year
ppl plot

This kind of structure works just as well for FILL plots, filling in  
the shading from 1 Jan through 31 Dec.

Billy K



On 25Apr 2008, at 11:13 AM, Pierre Sepulchre wrote:

> Dear Ferreters,
>
> I am having a little problem plotting seasonal cycles.
> Let's say I have a classical 12-timestep variable and i want to  
> plot the variable monthly values for one gridpoint :
>
> plot/x=-50/y=20 rain[d=1] is working fine.
>
>
> But, what I would like to get is a plot with the last value on the  
> right being equal to the first on the left, i.e a plot starting  
> with january and finishing with january (13 values instead of 12,  
> just so that it's easier to see the transition between december and  
> january).
>
> I tried this trick :
>
> define axis/t=1:12:1/modulo mytimeaxis
> plot/x=-50/y=20/hlim=1:13 rain[d=1,gt=mytimeaxis@asn]
>
> but it didn't work.
> If anyone can help, i would really appreciate it ^_^
>
> Thanks,
> Pierre
>
> <pierrese.vcf>