Re: [ferret_users] Just pick up a particular month from monthlydata...

[Rob Scott <rscott@xxxxxxxxxxxxxxxxxx>]

This is how I would do it :

Suppose your monthly variable is called "sst"


! mask out months you don't want:
let month = l[gt=sst]   !month since Jan 1900 so 1=Jan, 2=Feb,
13=Feb, etc.
                       
! make sure you sst grid starts on Jan.


let jjason = if (mod(month,12) gt 5) and (mod(month,12) lt 12) then 1
else 1/0
let aso = if (mod(month,12) gt 7) and (mod(month,12) lt 11) then 1 else
1/0

let sst_season = sst * aso

! regrid onto yearly axis
define
axis/t="01-JAN-1900":"31-DEC-2000":12/edges/units=month
taxis_an
   !show axis/t taxis_an
let sst_season_an = area_season[gt=taxis_an@ave]

! there you go!
! (though I'd like to know if there's a better way)

Hope that helps,
Rob





At 10:18 AM 2/10/2006, Seon Tae wrote:
> Hi all
>
> I have 100 year-monthly SST anomaly data averaged in NINO 3 region. I
> want to plot time series after picking up SST value of a particular month
> of each year.
> I mean, for example, I want to obtain time series plot for January for
> 100 years.
> I tried to find it in the Q and A archive on the website but i
> couldn't
> Could you please help me?
>
> Thank very much in advance.....
> -Seon Tae-

Rob Scott
Research Associate
Institute for Geophysics
University of Texas at Austin
4412 Spicewood Springs Rd.
Building 600
Austin, TX
78759

Office: +1-512-471-0375