This is NOT a question but the explanation of how the loss terms for phytoplankton are calculated.
I hope this will help someone who want to understand this.
There are three terms of respiration, mortality and excretion for the phytoplankton loss.
Each loss term is basically done by [rate] * [concentration of phytoplankton] * [delta t (time step)].
For respiration,
1) the initial value is calculated by the [power_temperature] and [salinity_freshwater] functions =[a]
2) the above value [a] is multiplied by [respiration_coefficient] = [b]
3) the final value can be added with a fraction of growth defined by the parameter of [respiration_fraction_growth] =[c]
4) the loss term by respiration is calculated by the [c] * [concentration of phytoplankton] * [delta t (time step)].
For mortality and excretion,
1) the default value is [b] for respiration.
2) the final value is multiplied by the fraction of respiration defined by the parameter of [mortality_fraction_respiration] or [excretion_fraction_respiration] [d]
3) the loss term by mortality or excretion is calculated by the [d] * [concentration of phytoplankton] * [delta t (time step)].
[mortality_coefficient] and [excretion_coefficient] are not used.
Hiroshi
