# NOT RUN {
data(NL_6min_2013)
### an elevation of interest: 0.9117 m
MHW <- 0.9117
wave.dur(NL_6min_2013[, 2], MHW)
# }
# NOT RUN {
### function also accepts and returns vectors
elevs <- data.frame(elevs = c(-1000:1500) / 10^3)
elevs$wave.dur <- wave.dur(NL_6min_2013[, 2], elevs$elevs)
plot(elevs ~ wave.dur, data = elevs)
# }
Run the code above in your browser using DataLab