T_TIDE
package (Pawlowicz et al.,
2002). It results from measurements made in 1975 at Tuktoyaktuk,
Northwest Territories, Canada.The data set is a test contains 1584 points, some of which have NA for sea-level height. The first point in the set is at time 1975-07-06 0800 GMT. However, the first 15 points of the series have NA for the sealevel value, so that the first non-missing point is at 2300 GMT, or 1600 MST as in Foreman's Appendix 7.2.
data(sealevel.tuk)
T_TIDE
dataset (which
Pawlowicz et al. (2002) seems to have based on Appendix 7.2 of Foreman
(1977), to build a test case), edited a bit to get the dates in a
format that Rcould scan, and put into a data file. This file was
then used to create sealevel
object follows.)# Note: the data file is not supplied with the package, # and the 7h offset is to convert from Mountain Standard Time. tuk <- read.table("tuk_time_elev.dat", header=FALSE, as.is=TRUE) t <- as.POSIXlt(strptime(tuk$V1, "%d-%b-%Y %H:%M:%S", tz="GMT")+7*3600, tz="GMT") eta <- tuk$V2 sealevel.tuk <- as.sealevel(eta=eta, t=t, station.name="Tuktoyaktuk", region="NWT", station.number=6485, longitude=133.0292, latitude=69.43889, year=1975, GMT.offset=0) #save(sealevel.tuk, file="oce/data/sealevel.tuk.rda")
Pawlowicz, Rich, Bob Beardsley, and Steve Lentz, 2002.
Classical tidal harmonic analysis including error estimates in MATLAB using T_TIDE
.
Computers and Geosciences, 28, 929-937.