correction_function {soilfoodwebs}

A function to correct stoichiometry dynamically

Description

A function to correct stoichiometry dynamically

Usage

correction_function(
  biomasses,
  cij,
  CN,
  p,
  a,
  canIMM,
  dietlimits,
  diet_correct = TRUE,
  Conly = FALSE,
  Immobilizationlimit = Inf
)

Arguments

biomasses

A vector of biomasses.

cij

The consumption matrix from Cijfcn

CN

A vector of C:N ratios.

p

A vector of production efficiency.

a

A vector of assimilation efficiency.

canIMM

A Boolean vector of whether the nodes can immobilize nitrogen.

dietlimits

The diet limits matrix for the stoichiometry correction (proportion of diet)?

diet_correct

Boolean: Does the organism correct it's diet?

Conly

Boolean: Is the model meant for carbon only?

Immobilizationlimit

This is the limit of the amount of nitrogen the food web can immobilize nitrogen (NOT PLANTS). This will impact the calculations of inorganic nitrogen dynamics.

Details

This function takes inputs from the ODE and outputs corrected consumption rates. The key difference from 'corrstoich' is that the prey DO NOT correct their feeding rates to compensate for higher consumption from the predators, so the system can leave equilibrium if a diet shift occurs.

Value

Returns the consumption rates (FMAT) and production efficiencies (p).