Seasonal semichemostat resource dynamics — seasonal_resource

This implements the standard semichemostat dynamics for the resource (see mizer::resource_semichemostat()) but with a time-dependent carrying capacity. The carrying capacity is given by $c_{R} (w, t) = (1 + maxR \cdot vonMises (t, μ, κ)) c_{R} (w)$ where $c_{R} (w)$ is the standard carrying capacity, $μ$ is the mean of the von-Mises distribution, and $κ$ is the concentration parameter of the von-Mises distribution. These parameters must be supplied in the rp$maxR, rp$mu and rp$kappa slots of the resource_params slot of the params object.

Usage

seasonal_resource_semichemostat(
  params,
  n,
  n_pp,
  n_other,
  rates,
  t,
  dt,
  resource_rate,
  resource_capacity,
  ...
)

Arguments

params: A MizerParams object
n: A matrix of species abundances (species x size)
n_pp: A vector of the resource abundance by size
n_other: A list with the abundances of other components
rates: A list of rates as returned by mizerRates()
t: The current time
dt: Time step
resource_rate: Resource replenishment rate
resource_capacity: Resource carrying capacity
...: Unused

Value

Vector containing resource spectrum at next timestep