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

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

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.

Examples

Run this code

# NOT RUN {

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