[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