[Thread Prev][Thread Next][Index]
Presently available EOF external functions
Hello Ferreters -
This seems to be an appropriate time to say a few things about the
EOF external functions available in V5.n of Ferret. These are
EOF_SPACE, EOF_STAT, and EOF_TFUNC. See the V5 documentation for
descriptions of the use of these.
These EOF routines were originally designed for a special purpose:
to compute EOFs on gappy time series. Gaps in the time series
introduce special complications which require a much more difficult
(and CPU-consuming) code to deal with. Briefly, it is possible to
make estimates of the EOFs and time functions (PCs) that minimize
the expected square error in a least squares sense. (Ask me for
details; the algebra is not easy to write down in an e-mail). I
wrote these routines in fortran about 15 years ago, based on
unpublished notes by Dudley Chelton.
These EOF routines are now in ferret because someone asked the PMEL
ferret team to make them into external functions. A more systematic
approach would have made the more usual (and MUCH faster) type of
EOF routines external functions first, but the squeaky wheel gets
the grease, so this is the way it turned out. We are now working
on making better EOF external functions available.
So while the present routines definitely work, they are not ideal for
probably the great majority of EOF work that most of us encounter.
For complete time series, they are far slower than a usual EOF
routine (roughly nx*nt times slower). Lots of needless computation
is being done if there are no gaps. On the other hand, if you DO
have a gappy time series, then these routines are very useful and
hard to come by otherwise.
Billy K
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
William S. Kessler
NOAA / Pacific Marine Environmental Laboratory
7600 Sand Point Way NE
Seattle Wa 98115 USA
Tel: 206-526-6221 o__ ____
Fax: 206-526-6744 _,>/'_ -----
E-mail: kessler@pmel.noaa.gov (_) \(_) ------
Home page: http://www.pmel.noaa.gov/~kessler
[Thread Prev][Thread Next][Index]
Dept of Commerce /
NOAA /
OAR /
PMEL /
TMAP
Contact Us | Privacy Policy | Disclaimer | Accessibility Statement