https://github.com/AquaticEcoDynamics/libaed2/blob/b7ac2df190dddccabf5e21264f3510b0ce07812a/src/aed2_carbon.F90#L147