Solves the Mussel individual bioenergetic balance with a 4th order Runge Kutta method for spatialized model
Mussel_spatial_RKsolver(
Param,
times,
IC,
Tint,
Phyint,
DTint,
POCint,
Ccont,
Ncont,
Pcont,
POMint,
TSSint
)
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