Get predation mortality rate for plankton

Calculates the predation mortality rate \(\mu_p(w)\) on the plankton spectrum by plankton size (in units 1/year).

getPlanktonMort(
  params,
  n = params@initial_n,
  n_pp = params@initial_n_pp,
  n_other = params@initial_n_other,
  pred_rate = getPredRate(params, n = n, n_pp = n_pp, n_other = n_other)
)

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 of abundances for other dynamical components of the ecosystem
pred_rate	An array of predation rates of dimension no. sp x no. community size bins x no. of size bins in whole spectra (i.e. community + plankton, the w_full slot). The array is optional. If it is not provided it is calculated by the getPredRate function.

Value

A vector of mortality rate by plankton size.

Details

Used by the project function for running size based simulations.

See also

Other rate functions: getEGrowth(), getEReproAndGrowth(), getERepro(), getEncounter(), getFMortGear(), getFMort(), getFeedingLevel(), getMort(), getPredMort(), getPredRate(), getRDD(), getRDI(), getRates(), getStarvMort()

Examples

if (FALSE) {
data(NS_species_params_gears)
data(inter)
params <- newMultispeciesParams(NS_species_params_gears, inter)
# With constant fishing effort for all gears for 20 time steps
sim <- project(params, t_max = 20, effort = 0.5)
# Get plankton mortality at one time step
n <- sim@n[21,,]
n_pp <- sim@n_pp[21,]
getPlanktonMort(params,n,n_pp)
}