[Thread Prev][Thread Next][Index]

[ferret_users] defining common time axis for 2 datasets with different origin/end



Hi Ferret users,

I would like to compare (bias and correlation maps) 2 monthly precipitation datasets.

The first one, D1, starts on 1rst january 1901 and ends on 1rst december 2006, so 106 years and 1272 months: L varies from 1 to 1272
 01-JAN-1901 /    1:    60.
 31-JAN-1901 /    2:    10.
 03-MAR-1901 /    3:    31.
....
 01-OCT-2006 / 1270:    88.
 01-NOV-2006 / 1271:     9.
 01-DEC-2006 / 1272:    16.


The second one, D2, starts at 16 jan 1951 and ends at 16 december 2007, so 57 years and 684 values: L varies from 1 to 684
 16-JAN-1951 /   1:   21.4
 15-FEB-1951 /   2:    3.3
 18-MAR-1951 /   3:    7.0
...
 16-OCT-2007 / 682:   42.4
 16-NOV-2007 / 683:    2.3
 16-DEC-2007 / 684:    0.0

In order to compare D1 and D2 over the same period, I defined a common time axis/grid over 1951-2006 period such as : define axis/t="1-jan-1951":"31-dec-2006":30.436875/units=day axd ; define grid/t=axd gtd
and regrid the two datasets on this axis.

let p1 = pm[d=D1, gt=gtd]
let p2 = pm[d=D2, gt=gtd]

But when I list these new variables on one point new datasets, I have the following output:
                     P1     P2
01-JAN-1951 /   1:    4.0   ....
31-JAN-1951 /   2:   36.5   12.3
02-MAR-1951 /   3:   33.2    5.1
...
31-OCT-2006 / 671:   11.2   11.3
01-DEC-2006 / 672:   16.8    9.1
31-DEC-2006 / 673:   ....    8.7

Consequently, it's impossible to compare these 2 regridded datasets.
Then, how can I define the correct time axis (i.e. identical L values) wich will give consistent regrided datasets ?

Thanks in advance.

--
Francois DELCLAUX
------------------------------------------------------------
UMR HydroSciences Montpellier
Universite Montpellier II - Place Eugene Bataillon
Case courrier MSE
34095  Montpellier Cedex 5      FRANCE
http://www.hydrosciences.fr/
mailto: delclaux@xxxxxxxxxxxxxxxxxxx
Tel : (33) (0)4 67 14 90 11      Fax : (33) (0)4 67 14 47 74
------------------------------------------------------------



[Thread Prev][Thread Next][Index]

Contact Us
Dept of Commerce / NOAA / OAR / PMEL / TMAP

Privacy Policy | Disclaimer | Accessibility Statement