Description Usage Arguments Value References See Also Examples
Function for simplifying the calculation of buoyancy frequency. Can usually
be called directly on data loaded directly using load.ts
and
load.bathy
.
1 | ts.buoyancy.freq(wtr, at.thermo = TRUE, na.rm = FALSE, ...)
|
wtr |
A data frame of water temperatures (in Celsius). Loaded using
|
at.thermo |
Boolean indicating if only buoyancy frequency at the thermocline should be returned. If false, full profile is returned. |
na.rm |
Boolean indicated if step-by-step removal of NA's should be tried. If false, a timestep with any NA values will likely return an NA value. If true, best effort will be made to calculate indices despite NA values. |
... |
Additional parameters will be passed on to |
Returns a data frame with the timeseries of buoyancy frequency in
units sec^-2
. Includes a ‘datetime’ column.
Imberger, J., Patterson, J.C., 1990. Physical limnology. Advances in Applied Mechanics 27, 353-370.
1 2 3 4 5 6 7 8 9 10 11 12 | #Get the path for the package example file included
wtr.path <- system.file('extdata', 'Sparkling.daily.wtr', package="rLakeAnalyzer")
#Load data for example lake, Sparkilng Lake, Wisconsin.
sp.wtr = load.ts(wtr.path)
N2 = ts.buoyancy.freq(sp.wtr, seasonal=FALSE)
SN2 = ts.buoyancy.freq(sp.wtr, seasonal=TRUE)
plot(N2, type='l', ylab='Buoyancy Frequency', xlab='Date')
lines(SN2, col='red')
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.