mcyst_batch: Microcystin Production Model (Batch Version)

In simecolModels: Model collection for the simecol package

Description

This is a batch version of the MCYST model with logistic growth.

Usage

mcyst_batch()

Details

To see all details, please have a look into the implementation and the original publications.

Value

S4 object according to the odeModel specification. The object contains the following slots:

main	The differential equations for cell abundance and Microcystin concentration: cells abundance of cells (1/L), mcyst microcystin concentration (fg/L).

parms	Vector with the named parameters of the model: mu Growth rate of cells (1/d), K Carrying Capacity (cells /L), p MCYST production coefficient (fg/cell), dM MCYST depletion rate (1/d), s sedimentation rate (1/d).
times	Simulation time and integration interval.
init	Vector with start values for cells and mcyst.
solver	Character string with the integration method.

References

Jähnichen, S., Petzoldt, T. and Benndorf, J. (2001). Evidence for Control of Microcystin Dynamics in Bautzen Reservoir (Germany) by Cyanobacterial Population Growth Rates and Dissolved Inorganic Carbon. Archiv für Hydrobiologie, 2, 150, 177-196.

Jähnichen, S., Ihle, T. and Petzoldt, T. (2008). Variability of microcystin cell quota: A small model explains dynamics and equilibria. Limnologica - Ecology and Management of Inland Waters, 38, 339-349. URL http://dx.doi.org/10.1016/j.limno.2008.05.003

See Also

sim, parms, init, times.

Examples

## create model object
m <- mcyst_batch()
## inspect parameters, start values, time steps
parms(m)
init(m)
times(m)
## simulate the model
m <- sim(m)
o <- out(m)

par(mfrow=c(3,1))
plot(o$time, o$cells, col="green", type="l",
  xlab="", ylab="Microcystis", ylim=c(0, max(o$cells)))
plot(o$time, o$mcyst/1000, col="cyan",  type="l",
  xlab="time (d)", ylab="Microcystin mug/l")
plot(o$time, o$mcyst/o$cells, col="red",  type="l",
  xlab="time", ylab="spec. MCYST")

simecolModels documentation built on May 2, 2019, 4:59 p.m.