Dampen inputs of rpmodel — dampen

Applies an exponential dampening input time series with specified time scale.

dampen_vec(vec, tau)

Arguments

vec	A numeric vector for the time series of a daily meteorological variable used as input for `rpmodel` (temperature, vapour pressure deficit, CO2, or atmospheric pressure). The length of `x` must be at least 365, i.e., corresponding to one year.
tau	The time scale of dampening (e-folding time scale of a perturbation). Must be smaller or equal to 365 d.

Value

A numeric vector of equal length as x with damped variation. The dampening is calculated as: $$ S(t+1) - S(t) = (X(t+1) - S(t)) / \tau $$ Where $X$ is the daily varying time series given by argument x, $S$ is the dampened time returned by this function, and $\tau$ is the decay time scale of a perturbation, given by argument tau.

Examples

if (FALSE) {
dampen_vec(
 vec = 20 * (sin(doy*pi/(365)))^2 + rnorm(365, mean = 0, sd = 5),
 tau = 40 
 )
}