ts.wedderburn.number: Calculate Wedderburn number for a timeseries.

In rLakeAnalyzer: Lake Physics Tools

Description

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

Usage

ts.wedderburn.number(wtr, wnd, wnd.height, bathy, Ao, seasonal = TRUE)

Arguments

wtr	A data frame of water temperatures (in Celsius). Loaded using load.ts
wnd	A data frame of wind speeds (in m/s). Loaded using load.ts
wnd.height	Height of the anemometer above the lake surface in meters
bathy	A data frame containing hypsometric data. Loaded using load.bathy
Ao	Numeric value for the water body surface area (m^2) at zero meters depth
seasonal	Boolean indicating if seasonal thermocline should be used in calculation.

Value

Returns a data frame with the timeseries of Wedderburn number. Includes a ‘datetime’ column.

References

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

See Also

wedderburn.number,ts.lake.number

Examples

	#Get the path for the package example file included
	wtr.path <- system.file('extdata', 'Sparkling.daily.wtr', package="rLakeAnalyzer")
	wnd.path <- system.file('extdata', 'Sparkling.daily.wnd', package="rLakeAnalyzer")
	bathy.path <- system.file('extdata', 'Sparkling.bth', package="rLakeAnalyzer")
	
	#Load data for example lake, Sparkilng lake, in Wisconsin.
	sp.wtr = load.ts(wtr.path)
	sp.wnd = load.ts(wnd.path)
	sp.bathy = load.bathy(bathy.path)
	
	sp.area = 64e4  #Area of Sparkling lake in m^2
	wnd.height = 2  #Height of Sparkling lake anemometer
	
	w.n = ts.wedderburn.number(sp.wtr, sp.wnd, wnd.height, sp.bathy, sp.area)
	plot(w.n$datetime, w.n$wedderburn.number, type='l', ylab='Wedderburn Number', xlab='Date')
	

rLakeAnalyzer documentation built on June 10, 2019, 1:02 a.m.