# \donttest{
library(oce)
data(adv)
# Velocity time-series
plot(adv)
# Spectrum of upward component of velocity, with ``turbulent'' reference line
s <- spectrum(adv[["v"]][, 3], plot = FALSE)
plot(log10(s$freq), log10(s$spec), type = "l")
for (a in seq(-20, 20, by = 1)) {
abline(a = a, b = -5 / 3, col = "gray", lty = "dotted")
}
# }
Run the code above in your browser using DataLab