[Thread Prev][Thread Next][Index]
Re: [ferret_users] seasonal climatologies
- To: Izidine Pinto <izidinep@xxxxxxxxx>
- Subject: Re: [ferret_users] seasonal climatologies
- From: Peter Szabo <szabpet83@xxxxxxxxx>
- Date: Wed, 16 Jun 2010 10:35:34 +0200
- Cc: ferret user ferret user <ferret_users@xxxxxxxx>
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:mime-version:received:received:in-reply-to :references:date:message-id:subject:from:to:cc:content-type; bh=/fXHlxuJBPfo2ZnxeqwwPbfOUzCDknTpBNVSUmGNI6o=; b=so4Ku2x3DDoZkmjTuomQjSmo9OIMgFzvvRoqB5Xx1vL8mW/4fJbyUgP5opGJizCCGQ Ur1GhckYrnju+luQ60KsXJ6Ywd2ZfxB/4zaCms3eCFP+6ZPZ3rS/1qDcH+bAKeEFxPmC SZG4j4vEIXi70Hapv8r/wiybMBZhGlPSin5b0=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=mime-version:in-reply-to:references:date:message-id:subject:from:to :cc:content-type; b=AfuqAnF8dCQxsL7F6JSyzWqcmc9EPC9UfiEo9cJnzbL4LCRXRt2oi0/02/tLrmbTu5 dsd/YR+ZE0TZkQTOyCqQ4+AcW0ZdZkugXlr2bohmRJrrErrP8awDZsF7jJqa+I54dQ8M ZY0QDR2avjrKUam4/SUyZ/SmAH4I7ccIPEVpc=
- In-reply-to: <AANLkTinWJipAgbCe0qYYCpMdnqP6HzKy1-n_pG0-EQ9U@xxxxxxxxxxxxxx>
- References: <AANLkTinWJipAgbCe0qYYCpMdnqP6HzKy1-n_pG0-EQ9U@xxxxxxxxxxxxxx>
- Sender: owner-ferret_users@xxxxxxxx
Dear Izidine,
after you have created a 12-month climatology, you wil have J, F, ... D (l=1:12)
let pre_monthlycl = precip[d=1,gt=month_reg@mod] ! 12-month climatology
with the running smoothing (averaging) you will get an average for JFM, MAM, ... ,OND (l=1:12, but only l=2:10 will exist)
let pre_3monthav = pre_monthlycl[l=@sbx:3] ! running 3-month average
and this will select MAM, JJA and SON (l=1:4, but only l=:2:4 will exist)
save/file=
pre_seasonal_alt.nc pre_3monthav[l=1:10:3] ! save only DJF,MAM,JJA,SON
Consequently, they will differ.
Here is my demonstration:
yes? let t2mmod=t2m[l=1:12] !this equals to pre_monthlycl
yes? let t2msbx=t2mmod[l=@sbx:3] !~pre_3monthav
yes? let t2m3month=t2msbx[l=1:10:3] !~pre_3monthav[l=1:10:3]
listings:
T2MMOD T2MSBX
1: -3.000 ....
2: -5.500 -4.867
3: -6.100 -4.833
4: -2.900 -3.600
5: -1.800 -2.967
6: -4.200 -1.700
7: 0.900 -1.033
8: 0.200 -0.500
9: -2.600 -1.133
10: -1.000 -1.667
11: -1.400 -1.433
12: -1.900 ....
T2M3MONTH
1: ....
2: -3.600
3: -1.033
4: -1.667
Peter
On Tue, Jun 15, 2010 at 11:58 AM, Izidine Pinto
<izidinep@xxxxxxxxx> wrote:
Dear ferret users
I am comparing two seasonal means of two files
file 1 contains daily data from 1jan1981 to 31dec1999
and file 2 contains daily data only for 'all' DJF for the same period
1> ./daily.ppt.1981.1999.era40.nc
name title I J K L
PRECIP TOTAL PRECIPITATION RATE 1:127 1:110 1:1 1:6939
2> ./ppt.DJF.daily.1981_1999.era40.nc (default)
name title I J K L
PRECIP TOTAL PRECIPITATION RATE 1:127 1:110 1:1 1:1683
let pre_monthlycl = precip[d=1,gt=month_reg@mod] ! 12-month climatology
let pre_3monthav = pre_monthlycl[l=@sbx:3] ! running 3-month average
save/file=pre_seasonal_alt.nc pre_3monthav[l=1:10:3] ! save only DJF,MAM,JJA,SON
shade pre_3monthav[l=1] ! DJF average
shade PRECIP[d=2,l=1:1683@ave] ! DJF average
I was expecting to have the same result
Hope someone can help me to understand why the results differs
Thanks in advance
--
I.P.
"In these matters the only certainty is that nothing is certain". Pliny the Elder (23 AD - 79 AD)
[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce /
NOAA /
OAR /
PMEL /
TMAP
Privacy Policy | Disclaimer | Accessibility Statement