changing x-axis into time-axis

[Patrick Joeckel <joeckel@xxxxxxxxxxxxxxxxx>]

Dear Ferreters,

with the following script I try to generate a "climatology" from scattered
data. I have ascii-data in the form  "day lat value", where "day" is the day of
the year (1...365) and "lat" is the latitude. Using SCATTER2GRIDLAPLACE_XY(...),
it is possible to achieve a gridded dataset on a x-y grid (see script below). 
However since I would
like to compare with other data, that is on a time-y grid, I would like to
transfer the x-axis (which is the day number) to a ferret time-axis. Is this
possible (e.g using SAMPLEI(..) , RESHAPE(...) etc. ). Any suggestions are
appreciated !

Yours, sincerely

                 Patrick Joeckel

!---------------------------------------------------
! open the file
file/var="day,lat,value" data.txt

! define grid for objective analysis
DEFINE AXIS/x=1:365:1 xax
DEFINE AXIS/y=-90:90:5/unit=lat yax

! gridding of scattered data
LET data_pre = scatter2gridlaplace_xy (day, lat, value, x[gx=xax],
y[gy=yax], 5., 5) 

! this result looks OK, however it is on a latitude-x - grid
shade data_pre   

! define the time axis and grid 
! finally the data should be transferred to this grid
DEFINE AXIS/t="01-Jan-1993":"31-Dec-1993"/npoints=365 tax
DEFINE GRID/y=yax/t=tax gyt

! ??? how to proceed ???

! ------------------------------------------------


-- 
---------------------------------------------------------------
Patrick Joeckel
Max Planck Institute for Chemistry
Joh.-J. Becher Weg 27         
55128 Mainz                   phone:++49-6131-305452      
Germany                       fax  :++49-6131-305436     
                             e-mail:joeckel@mpch-mainz.mpg.de
                         http://www.mpch-mainz.mpg.de/~joeckel
---------------------------------------------------------------
Never trust an operating system you don't have sources for.