Turbulence spectrum \(phi\), as a function of wavenumber \(k\) for a run starring at 0905h on October 3, 1959, with data in their Appendix 1 and graph in their Figure 9 and Appendix 1. The units are c.g.s., that is \(k\) is in \(cm^{-2}\) and \(phi\) is in \(cm^3/s^2\).
data(turbulence, package="ocedata")
H. L. Grant and R. W. Stewart and A. Moilliet, 1962. Turbulence spectra from a tidal channel. Journal of Fluid Mechanics, 12(2): 241-268.
# Recreate Fig 9 of Grant et al. (1962), with intercept
# chosen to mimic the published figure
data(turbulence, package="ocedata")
plot(log10(turbulence$k), log10(turbulence$phi),
xlab=expression(log[10]*k), ylab=expression(log[10]*phi), pch=20)
grid()
power <- -5/3
intercept <- -0.3
abline(intercept, power)