Hi ferreters,
I have model outputs on a 360_day calendar (Roms output of scenarios with
IPSL forcings).
I have one file per month, for 20 years.
I would like to plot a variable (ex SST) from 01-Jan-1980 to
30-Dec-2000 for example (1 month after 1 month).
To do this, I define a time-axis at the same time of my first plot:
> plot/nolab/z=4/vlimits=11:21/color=black/line
temp_grid[d=3,x=-11:0@ave,y=41:48.5@ave,T="01-JAN-1980":"30-DEC-2000":24]
the first variable temp_grid is defined from 01-Jan-1980 to 30-Jan-1980
but Ferret answer me an error of time (see attached file):
indicating that the temp_grid variable is defined from October 2008 to February
2030.
If I check with the command
> show grid temp_grid
the time definition of temp_grid in January 1980 is confirmed.
And if I plot this variable without defining explicitly the time-axis in
the command:
> plot/nolab/z=4/vlimits=11:21/color=black/line
temp_grid[d=3,x=-11:0@ave,y=41:48.5@ave]
it is plotted in January 1980….
So it seems that an error is induced by the definition of the time-axis
in the “plot” command.
Is the 360_day calendar a problem?
The time axis and grid are previously defined as:
> define
axis/calendar=360_day/T="01-JAN-1980:12:00":"30-JAN-1980:12:00":24
taxis_01
> define grid/T=taxis_01/X=xaxis/Y=yaxis/Z=zdepth
gridBiscay01
And then the variable temp_grid:
> let/d=3/bad=-999 temp_z = IF
mask_rho[d=1] LT 1 THEN (-999) ELSE ZAXREPLACE(temp[d=3],depth,z[gz=zdepth])
> let/d=3 temp_grid =
temp_z[d=3,g=gridBiscay01@ASN]
Do you have any idea to correct this error, or to do it in another way?
Thank you very much!
_____________________________________________
Dr
AZTI - Tecnalia /
Marine research Division
Herrera Kaia Portualde z/g
20110 Pasaia (Gipuzkoa) -
Tel: 943 00 48 00 (ext. 890) - Fax: 946 572 555
e-mail: mchifflet@xxxxxxxxxxx
www.azti.es
; www.tecnalia.info
|
|
|
|
|
|
|
|
|
|
|
|
