[Thread Prev][Thread Next][Index]

Re: [ferret_users] selective climatology



Hi Timothee,
There isn't an & syntax in Ferret, but you can do this kind of thing with a "mask"

! Define a variable with the time coordinates, hours since t0 or whatever units it is
let tt = t[gt=hanch]
! Locate the edges of your time segments
let t1 = tt[t=1-dec-1972]
let t2 = tt[t=1-feb-1973]
let t3 = tt[t=1-aug-1976]
let t4 = tt[t=1-jul-1977]

! Define the mask. There is an implied "else missing"
let tmask = if ( (tt gt t1) and (tt lt t2) ) or ( (tt GT t3) and (tt LT t4) ) then 1

! This variable has the value of hatch during your El Nino times, missing otherwise
let hanch_NINO = hatch* tmask

!(try a plot a time series of hatch_NINO to see the mask work:
! plot/x=/y= hatch_NINO)
and proceed as you have done in your script


Timothée Brochier wrote:

Hello,

I want to do a climatology over the warm events of El nino,
how can I select these times interval?

I tried this way, but it doesnt work

yes? let hanch_NINO= hanch[((T=1-dec-1972:1-feb-1973) & (T=1-aug-1976:1-jul-1977))]
yes? let hanch_NINO_clim = hanch_NINO[gt=seasonal_reg@mod,d=1]
yes? fill/l=4 hanch_NINO_clim
**ERROR: command syntax: ((T=1-DEC-1972:1-FEB-1973) & (T=1-AUG-1976:1-JUL-1977)

(later I'll agregate all the warm episodes, until today, about 10...)

many thanks

Timothee B.


Este mensaje ha sido analizado por MailScanner
en busca de virus y otros contenidos peligrosos,
y se considera que esta limpio.
MailScanner agradece por el apoyo.


[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP

Contact Us | Privacy Policy | Disclaimer | Accessibility Statement