Phytoplankton
P(i), i=1:np
Light/temperature
Advection/mixing
River discharges
Boundary flux
Y(i), i=1: ny    (Auxiliary state variables)
Nutrients
N(i), i=1:nn
Detritus
D(i), i=1:nd
Zooplankton
Z(i), i=1:nz
Bacteria
B(i), i=1:nb
DOM(i), i=1:nm