Solves the Mussel population bioenergetic balance with a 4th order Runge Kutta method
Mussel_pop_RKsolver(
Param,
times,
IC,
Tint,
Phyint,
DTint,
POCint,
Ccont,
Ncont,
Pcont,
POMint,
TSSint,
N
)a list containing the weights of the mussel, the excreted CNP, the mussel CNP, temperature limitation functions, metabolic rates, oxygen consumption
a vector containing model parameters
integration extremes and integration timestep
initial condition
the interpolated water temperature time series
the interpolated phytoplankton time series
the interpolated detritus time series
the interpolated POC time series
the C/C content of the POC
the N/C content of POC
the P/C content of POC
the interpolated POM time series
the interpolated TSS time series
the number of indivduals time series