Re: [ferret_users] regridding of time axis

["Ansley C. Manke" <ansley.b.manke@xxxxxxxx>]

Hi Steffie,

This would be similar to what is described here for creating a yearly time series; using either @AVE or another regridding transformation
../../faq/averaging-to-regrid-a-time-series

Except that you will want to define a seasonal axis instead of a yearly one. You can base one axis on the definition of an existing axis, taking every nth point, so starting with a monthly axis is an easy way to do this. If your incoming data already has a well-defined time axis, you can base it on that. Here's how that looks, using the dataset monthly_navy_winds that comes with the Ferret distribution:

yes? use monthly_navy_winds
yes? show grid uwnd
GRID GDN1
name axis # pts start end subset
FNOCX LONGITUDE 144mr 20E 17.5E(377.5) full
FNOCY LATITUDE 73 r 90S 90N full
normal Z
TIME TIME 132 r 16-JAN-1982 20:00 17-DEC-1992 03:30 full

! Use the coordinate values of the time axis from the file to define a new time axis
yes? let tt = t[gt=uwnd]
yes? define axis/t/t0="`uwnd,return=t0`"/units="`uwnd,return=tunits`" time4 = tt[L=2:144:3]
!-> define axis/t/t0="14-JAN-1980 14:00:00"/units="hour" time4 = tt[L=2:144:3]

! The original data:
yes? list/t=1-jan-1982:1-jan-1983/x=180/y=0 uwnd
VARIABLE : ZONAL WIND (M/S)
FILENAME : monthly_navy_winds.cdf
FILEPATH : /home/users/tmap/ferret/linux/fer_dsets/data/
SUBSET : 12 points (TIME)
LONGITUDE: 180E
LATITUDE : 0
180E
65
16-JAN-1982 / 1: -2.738
16-FEB-1982 / 2: -2.331
18-MAR-1982 / 3: -3.251
18-APR-1982 / 4: -3.124
18-MAY-1982 / 5: -1.818
18-JUN-1982 / 6: -2.007
18-JUL-1982 / 7: -0.367
17-AUG-1982 / 8: -0.900
17-SEP-1982 / 9: 3.223
17-OCT-1982 / 10: 2.167
17-NOV-1982 / 11: 2.875
17-DEC-1982 / 12: -0.007

! The data on the seasonal axis:
yes? list/t=1-jan-1982:1-jan-1983/x=180/y=0 uwnd[gt=time4@ave]
VARIABLE : ZONAL WIND (M/S)
regrid: 2192 hour on T@AVE
FILENAME : monthly_navy_winds.cdf
FILEPATH : /home/users/tmap/ferret/linux/fer_dsets/data/
SUBSET : 4 points (TIME)
LONGITUDE: 180E
LATITUDE : 0
180E
65
16-FEB-1982 / 1: -2.773
18-MAY-1982 / 2: -2.316
17-AUG-1982 / 3: 0.652
17-NOV-1982 / 4: 1.678


If you are reading data from a text file, then you'd first define the monthly time axis (See DEFINE AXIS) and a grid for the data, and set up to read the data as described here, and then define the seasonal axis as above.

-Ansley

On 11/19/2017 10:25 PM, Steffie Tom wrote:

hii

We have data from 1979-2015, which consists of various parameters like temperature etc. we need to create seasonal data i.e. say for summer jan-feb-march for all years i.e. separately for 1979, then jan-feb-march for 1980 etc. i.e in one file we want 111 months.

First of all how to achieve this ?

After that how to achieve the seasonal average of the same data above ?

thanks