Re: [ferret_users] saving seasonal values on to a continuos time axis

[Steeven Paul Yerraguntla <steevenpaul@xxxxxxxxx>]

Hi nuncio,
      In my previous mail, there is a confusion in variable naming. So please ignore that, and plz. check the following
        use sst.nc
        let month=l[gt=temp]
        let DJF = if (mod(month,12) le 2) OR (mod(month,12) le 2) then 1

        let MAM = if (mod(month,12) ge 3) AND (mod(month,12) le 5) then 1
        let JJA = if (mod(month,12) ge 6) AND (mod(month,12) le 8) then 1
        let SON = if (mod(month,12) ge 9) AND (mod(month,12) le 11) then 1

        let winter_sst=temp*DJF
        let spring_sst=temp*MAM
        let summer_sst=temp*JJA
        let fall_sst  =temp*SON

    let sst_compr_djf = COMPRESSL(winter_sst)
        list/x=60/y=15/l=1:10  sst_compr_djf

        PAUSE

       ! 1870-2010 --> 140 years --> 140*3=420 data points (time axis)
     

    !defining to abstract axes one to put the data on the time axis

    define axis/t=1:420:1 t420

    let sst_djf = winter_compr_djf[gt=t420@ASN]
    list/x=60/y=15/t=1:10 sst_djf   !(or we can save the data or plot)
    
    ! Similarly we can extract MAM, JJA, SON.    Please take care in case of DJF values extraction and check before you use the extracted data.


regards,
- Paul.

On Thu, Oct 13, 2011 at 4:18 PM, Steeven Paul Yerraguntla <steevenpaul@xxxxxxxxx> wrote:

Hai  Nuncio,
          How are you?  I usually extract months, the way you want, as follows.This jnl file i saw somewhere in ferret threads provided by Jaison (i think). 

        use sst.nc
        let month=l[gt=temp]
        let DJF = if (mod(month,12) le 2) OR (mod(month,12) le 2) then 1
        let MAM = if (mod(month,12) ge 3) AND (mod(month,12) le 5) then 1
        let JJA = if (mod(month,12) ge 6) AND (mod(month,12) le 8) then 1
        let SON = if (mod(month,12) ge 9) AND (mod(month,12) le 11) then 1
        let winter_sst=temp*DJF
        let spring_sst=temp*MAM
        let summer_sst=temp*JJA
        let fall_sst  =temp*SON

    let sst_compr = COMPRESSL(wintersst)
        list/x=60/y=15/l=1:10  sst_compr
        PAUSE

       ! 1870-2010 --> 140 years --> 140*3=420 data points (time axis)
    

    !defining to abstract axes one to put the data on the time axis

    define axis/t=1:420:1 t420
    let sst_djf = winter_sst[gt=t420@ASN]
    list/x=60/y=15/t=1:10 sst_djf   !(or we can save the data or plot)
   
    ! we can extract the MAM, JJA, SON same as above. Please take care in case of DJF values extraction and check before you use the extracted data.

regards,
- Paul.

On Thu, Oct 13, 2011 at 2:40 PM, <nuncio@xxxxxxxxx> wrote:

Dear Ferreters
    I like to save all the JJA,SON,DJF values of a 30 year long time series of global SST to seperate files on to  a continuous time axis.  My data begins at 1980 january and ends at 2010 dec
I need to save 30*3 values continuously, ie after I make each file, t axis must run from 1 to 90, with out bad values in between.

Best regards
nuncio


Nuncio Murukesh, Ph.D
Research Scientist.
National Centre for Antarctic and Ocean Research
Headland Sada, Vasco Da Gama
Goa, 403804
India.
Ph: Cell-91 9890357423
   off- 91 832 2525428

--

Steeven Paul Y                              
Research Fellow
National Institute of Oceanography(NIO)
Regional Center, Visakhapatnam (RCV)                 
176 Lawsons Bay Colony
Visakhapatnam - 530017

    

--

Steeven Paul Y                              
Research Fellow
National Institute of Oceanography(NIO)
Regional Center, Visakhapatnam (RCV)                 
176 Lawsons Bay Colony
Visakhapatnam - 530017