ts.buoyancy.freq: Calculate the buoyancy (Brunt-Vaisala) frequency for a...

In clairervh/GLEON: Lake Physics Tools

Description

Function for simplifying the calculation of buoyancy frequency. Can usually be called directly on data loaded directly using load.ts and load.bathy.

Usage

ts.buoyancy.freq(wtr, at.thermo=TRUE, na.rm=FALSE, ...)

Arguments

wtr	A data frame of water temperatures (in Celsius). Loaded using load.ts
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 thermo.depth function when extracting buoyancy frequency at only the thermocline. Common parameters to supply would be seasonal and slope.

Value

Returns a data frame with the timeseries of buoyancy frequency in units sec^-2. Includes a ‘datetime’ column.

Author(s)

Luke Winslow

References

Imberger, J., Patterson, J.C., 1990. Physical limnology. Advances in Applied Mechanics 27, 353-370.

See Also

buoyancy.freq

Examples

	#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')
	

clairervh/GLEON documentation built on May 12, 2019, 2:04 p.m.