Re: [ferret_users] Plotting with a paleo time scale

[James Orr <James.Orr@xxxxxxxxxxxx>]

Olivier,

Here is a simple workaround if it is only a question of the y-axis label. 
Just replace your last command with the following 4 lines:

shade/nolab/set t2m[y=-2:2@ave,l=1:5000@sbx:250]
ppl ylab "Years before Present"
ppl xlab "Longitude"
ppl shade

James

On Wed, 10 Oct 2018, Olivier Marti wrote:

> Hello,
>
> I want to plot data using a paleo time scale : time is decreasing from a positive value to zero (present day).
>
> The following script gives more details : 
>
> USE"https://vesg.ipsl.upmc.fr/thredds/dodsC/store/p86mart/IPSLCM5A2/PROD/Holocene/TR5AS-Vlr01/ATM/Analyse/TS_YE/TR
> 5AS-Vlr01_20000101_79991231_1Y_t2m.nc"
>
> ! Model Years are from 2000 to 7999, which correspond to nothing :-(
> ! I can define a time axis starting at 1, which it's less confusing
> LET nt = `t2m_yearmean,RETURN=lsize`
> DEFINE AXIS/T=1:`nt`:1/UNITS="Years" YearFromStart
>
> ! But the labelling I want is Year before present (present being 1850)
> LET/TITLE="Years BP"   YearBP=6000-T[GT=YearFromStart]+1
>
> ! I can not define the following axis because it's reversed
> !DEFINE AXIS/T/UNITS="Years"/NAME="YearBP" YearBP
>
>
> LET t2m = t2m_yearmean[GT=YearFromStart@ASN]-273.15
>
>
> ! How can I use YearBP as an axis to have proper labelling 
> ! of the time axis on the following plot  ?
> SHADE  t2m[Y=-2:2@AVE,L=1:5000@SBX:250] 
>
>
> Any idea ?
>
> Thanks,
>
> Olivier
>
> -- 
> Dr. Olivier Marti  - mailto:olivier.marti@xxxxxxxxxxxx
> Laboratoire des Sciences du Climat et de l'Environnement
> Institut Pierre Simon Laplace
> Laboratoire CEA-CNRS-UVSQ  -  UMR 8212
> Tel : +33 1 69 08 77 27 - Mob : +33 6 45 36 43 74
> Fax : +33 1 69 08 30 73
> <http://www.lsce.ipsl.fr/Pisp/olivier.marti>

--
LSCE/IPSL, Laboratoire des Sciences du Climat et de l'Environnement
CEA-CNRS-UVSQ

LSCE/IPSL, CEA Saclay           http://www.ipsl.jussieu.fr/~jomce
Bat. 712 - Orme                 mailto:  James.Orr@xxxxxxxxxxxx
Point courrier 132
F-91191 Gif-sur-Yvette Cedex    Phone:   (33) (0)1 69 08 39 73
FRANCE                          Fax:     (33) (0)1 69 08 30 73