Keep plankton abundance constant

This function can be used instead of the standard plankton_semichemostat in order to keep the Plankton spectrum constant over time.

plankton_constant(params, n, n_pp, n_other, rates, t, dt, ...)

Arguments

params	A MizerParams object
n	A matrix of species abundances (species x size)
n_pp	A vector of the plankton abundance by size
n_other	A list with the abundances of other components
rates	A list of rates as returned by getRates()
t	The current time
dt	Time step
...	Unused

Value

Vector containing plankton spectrum at next timestep

Examples

if (FALSE) {
data(NS_species_params_gears)
data(inter)
params <- newMultispeciesParams(NS_species_params_gears, inter,
                                plankton_dynamics = "plankton_constant")
}