[Thread Prev][Thread Next][Index]
Re: a problem with FFTA
Hi everyone,
Just a PS to my note yesterday. Using monthly data requires
some caution, as the units of "month" on different monthly
axes do not necessarily mean the same thing. So one might
get a "drift" when regridding, especially noticeable with
long time series, and which might have an impact when doing
time series analysis like an FFT.
An approach that sidesteps this risk is
DEFINE AXIS/units=days/T=<lo>:<hi>/npoints=`T[gt=sst,t=<lo>:<hi>@ngd]` regt
In the example from the post, this would be:
LET ntim = T[gt=sst,T=01-JAN-1960:01-DEC-1997@ngd]
DEFINE AXIS/units=days/T=01-JAN-1960:01-DEC-1997/NPOINTS=`ntim` regt
Note the units are days, but the length of the time step
is still a "month". We nowuse npoints equal time steps.
Be careful that the endpoints match points on the original axis.
Or you could use use index limits and do it this way:
SHOW GRID sst ! see that the time axis has 456 points
LET tt = T[gt=sst]
DEFINE SYMBOL trng "`tt[L=1]`:`tt[L=456]`"
DEFINE AXIS/UNITS=`tt,RETURN=tunit`/T=($trng)/NPOINTS=`tt[t=($trng)@ngd]` regt
Ansley Manke
[Thread Prev][Thread Next][Index]
Dept of Commerce /
NOAA /
OAR /
PMEL /
TMAP
Contact Us | Privacy Policy | Disclaimer | Accessibility Statement