Reduced order upper ocean biogeochemical flux model with 17 state variables