Hi Emilie,
What trouble did you have when you tried using @MOD? That is really
what you
want to use, to average all the January values, then all the February
values, and so on, returning a 12-month result.
yes? use CMAP.1979-2
yes? use climatological_axes
yes? cancel data climatological_axes ! the axes will still be defined
yes? let clim_rain = rain1[d=1,x=10,y=0,gt=month_reg@mod]
yes? list clim_rain
Emilie Vanvyve wrote:
Hello dear
Ferret users,
I have the CMAP monthly data from 1979 to 2002 and would like to get
the mean values for every January, February..., and December. I looked
at the online manual and the FAQ and discovered the existence of the
"climatological_axes". I tried it, could not use the @MOD (as advised)
and so I used the @ASN.
When I check the results, Ferret just took my first year and copied it
repeatively across the next years, so that the year 1980 in the new
climatological variable = 1979 values of the original variable, the
year 1981 = 1979, and so on and so forth.
Example: I focus on one point (around 10°N 0°E).
yes? use CMAP.1979-2002.nc
yes? show data
currently SET data sets:
1> ./CMAP.1979-2002.nc (default)
name title I J K L
RAIN1 precipitation (gauge, satellite 1:144 1:72 ... 1:288
ERROR1 relative error for rain1 1:144 1:72 ... 1:288
RAIN2 precipitation (gauge, satellite 1:144 1:72 ... 1:288
ERROR2 relative error for rain2 1:144 1:72 ... 1:288
yes? list rain1[d=1,l=1:24,x=10,y=0]
VARIABLE : precipitation (gauge, satellite, models)
(mm/day)
DATA SET : CPC Merged Analysis of Precipitation (CMAP)
FILENAME : CMAP.1979-2002.nc
SUBSET : 24 points (TIME)
CALENDAR : 360_DAY
LONGITUDE: 8.8E
LATITUDE : 1.2S
8.8E
4
30-JAN-1979 / 1: 4.870
30-FEB-1979 / 2: 5.780
30-MAR-1979 / 3: 4.190
30-APR-1979 / 4: 6.630
30-MAY-1979 / 5: 4.130
30-JUN-1979 / 6: 0.530
30-JUL-1979 / 7: 0.620
30-AUG-1979 / 8: 0.180
30-SEP-1979 / 9: 0.400
30-OCT-1979 / 10: 2.020
30-NOV-1979 / 11: 7.650
30-DEC-1979 / 12: 4.780
30-JAN-1980 / 13: 4.820
30-FEB-1980 / 14: 2.680
30-MAR-1980 / 15: 8.200
30-APR-1980 / 16: 4.220
30-MAY-1980 / 17: 2.240
30-JUN-1980 / 18: 0.530
30-JUL-1980 / 19: 0.800
30-AUG-1980 / 20: 0.290
30-SEP-1980 / 21: 2.660
30-OCT-1980 / 22: 8.630
30-NOV-1980 / 23: 7.130
30-DEC-1980 / 24: 3.550
yes? use climatological_axes
*** NOTE: regarding /opt/ferret/go/climatological_axes.cdf ...
*** NOTE: Climatological axes SEASONAL_REG, MONTH_REG, and
MONTH_IRREG defined
yes? let clim_rain = rain1[d=1,x=10,y=0,gt=month_reg@asn]
yes? list clim_rain[l=1:24]
VARIABLE :
RAIN1[D=CMAP.1979-2002,X=10,Y=0,GT=MONTH_REG@ASN]
SUBSET : 24 points (TIME)
LONGITUDE: 8.8E
LATITUDE : 1.2S
8.8E
4
16-JAN / 1: 4.870
15-FEB / 2: 5.780
17-MAR / 3: 4.190
16-APR / 4: 6.630
16-MAY / 5: 4.130
16-JUN / 6: 0.530
16-JUL / 7: 0.620
16-AUG / 8: 0.180
15-SEP / 9: 0.400
16-OCT / 10: 2.020
15-NOV / 11: 7.650
16-DEC / 12: 4.780
15-JAN / 13: 4.870
14-FEB / 14: 5.780
17-MAR / 15: 4.190
16-APR / 16: 6.630
17-MAY / 17: 4.130
16-JUN / 18: 0.530
17-JUL / 19: 0.620
16-AUG / 20: 0.180
15-SEP / 21: 0.400
16-OCT / 22: 2.020
15-NOV / 23: 7.650
16-DEC / 24: 4.780
So clim_rain(1979)
= rain(1979), clim_rain(1980) = rain(1979), ...!
Where am I wrong?
Thanks for any suggestion!
Emilie V.
____________________________________________________________
Emilie VANVYVE
Physicist, PhD student
Université catholique de
Louvain (UCL)
Institut d'astronomie et de géophysique G. Lemaître (ASTR)
Chemin du Cyclotron, 2
1348 Louvain-la-Neuve (Belgium)
Phone: +32-(0)10-473300 // Fax: +32-(0)10-474722
E-mail: vanvyve@astr.ucl.ac.be // Web: www.astr.ucl.ac.be
|