fit_epsilon: fit epsilon from block of ADV data

Description Usage Arguments Details Value References Examples

Description

solves for epsilon from ADV data in the inertial subrange assuming Zappa et al. 2003

Usage

1
2
fit_epsilon(chunk.adv, freq = 32, lower = 20, upper = 80,
  diagnostic = FALSE)

Arguments

chunk.adv

a data.frame created with load.ADV, with the window.idx column

freq

a numeric, in Hz, that represents the input data

lower

a numeric, in Hz, that represents the lower bounds for range of epsilon calculation

upper

a numeric, in Hz, that represents the upper bounds for range of epsilon calculation

diagnostic

a boolean for diagnostic plot

Details

a GDopp function for splitting sensor data into temporal segments.

Value

a single value of epsilon estimating from the temporal block given.

References

Zappa, Christopher J., Peter A. Raymond, Eugene A. Terray, and Wade R. McGillis. "Variation in surface turbulence and the gas transfer velocity over a tidal cycle in a macro-tidal estuary." Estuaries 26, no. 6 (2003): 1401-1415.

Examples

1
2
3
4
5
6
7
8
9
## Not run: 
folder.nm <- system.file('extdata', package = 'GDopp')
file.nm <- "ALQ102.dat"
data.adv <- load_adv(file.nm=file.nm, folder.nm =folder.nm)
window.adv <- window_adv(data.adv,freq=32,window.mins=10)
fit_epsilon(window.adv[window.adv$window.idx==7, ],freq=32)
fit_epsilon(window.adv[window.adv$window.idx==7, ],freq=32, diagnostic = TRUE)

## End(Not run)

USGS-R/GDopp documentation built on May 9, 2019, 6:09 p.m.