Mussel_pop_RKsolver {RAC}

Solves the Mussel population bioenergetic balance with a 4th order Runge Kutta method

Description

Solves the Mussel population bioenergetic balance with a 4th order Runge Kutta method

Usage

Mussel_pop_RKsolver(
  Param,
  times,
  IC,
  Tint,
  Phyint,
  DTint,
  POCint,
  Ccont,
  Ncont,
  Pcont,
  POMint,
  TSSint,
  N
)

Arguments

Param

a vector containing model parameters

times

integration extremes and integration timestep

IC

initial condition

Tint

the interpolated water temperature time series

Phyint

the interpolated phytoplankton time series

DTint

the interpolated detritus time series

POCint

the interpolated POC time series

Ccont

the C/C content of the POC

Ncont

the N/C content of POC

Pcont

the P/C content of POC

POMint

the interpolated POM time series

TSSint

the interpolated TSS time series

N

the number of indivduals time series

Value

a list containing the weights of the mussel, the excreted CNP, the mussel CNP, temperature limitation functions, metabolic rates, oxygen consumption