# NOT RUN {
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")
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab