On the parameterization of phytoplankton primary production in water ecosystem models