## 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')
# ## End(Not run)
Run the code above in your browser using DataLab