[Thread Prev][Thread Next][Index]
[ferret_users] regriding from discontinuous time axis
Hi Ferreters
I would like to regrid a variable on a new time axis. The problem is
that the source grid is discontinuous (with gaps of several days) ; the
tranformation fills the gap with interpolated or averaged values
(depending of the transformation used). Is it possible to avoid this
behaviour (I mean fill the gap with missing values) ?
Here is an example :
! Discontinuous time axis from January 1 to January 12 (no values from 3
to 10)
Let time = {0,0.5,1.0,1.5,2.0,10.0,10.5,11.0,11.5,12.0}
Define Axis /T /From /T0=1-JAN-2000 /Units=days my_axis = time
Let var = TSEQUENCE( RANDU( time ) )
Let fake_var = T[GT=my_axis]
Let var_regrid = RESHAPE( var, fake_var )
! New time axis
Define Axis /T=01-JAN-2000:13-JAN-2000:1 /Units=days /T0=1-JAN-2000
my_axis_1dAV
Plot /Line /Sym=18 var_regrid
! Now each day has a good value
Plot /Over /Line /Sym=20 /Co=2 var_regrid[ GT=my_axis_1dAV@AVE ]
Thanks !
Bruno
[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce /
NOAA /
OAR /
PMEL /
Ferret
Privacy Policy | Disclaimer | Accessibility Statement