Calculates the energy rate (grams/year) available for reproduction after growth and metabolism have been accounted for.
getERepro( params, n = params@initial_n, n_pp = params@initial_n_pp, n_other = params@initial_n_other, e = getEReproAndGrowth(params, n = n, n_pp = n_pp, n_other = n_other) )
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 |
e | The energy available for reproduction and growth (optional). A
matrix of size no. species x no. size bins. If not supplied, is calculated
internally using |
A two dimensional array (prey species x prey size) holding
$$\psi_i(w)E_{r.i}(w)$$
where \(E_{r.i}(w)\) is the rate at which energy becomes available for
growth and reproduction, calculated with getEReproAndGrowth
,
and \(\psi_i(w)\) is the proportion of this energy that is used for
reproduction. This proportion is taken from the params
object and is
set with setReproduction
.
Other rate functions:
getEGrowth()
,
getEReproAndGrowth()
,
getEncounter()
,
getFMortGear()
,
getFMort()
,
getFeedingLevel()
,
getMort()
,
getPlanktonMort()
,
getPredMort()
,
getPredRate()
,
getRDD()
,
getRDI()
,
getRates()
,
getStarvMort()
if (FALSE) { data(NS_species_params_gears) data(inter) params <- newMultispeciesParams(NS_species_params_gears, inter) # Project with constant fishing effort for all gears for 20 time steps sim <- project(params, t_max = 20, effort = 0.5) # Get the energy at a particular time step getERepro(params,sim@n[21,,],sim@n_pp[21,]) }