When the model is running under a NBP constraint, it has a residual between the constraint value and the calculated NBP due to NPP + RH + LUC. This PR changes the model behavior so that the residual is sent to the deep ocean (following the behavior of the CO2 constraint) rather than added to or debited from the atmosphere, because the old behavior could result in runaway atmospheric CO2 under certain conditions. This was fully documented in #659 .
This does not change any behavior except when running under this constraint.
When the model is running under a NBP constraint, it has a residual between the constraint value and the calculated NBP due to NPP + RH + LUC. This PR changes the model behavior so that the residual is sent to the deep ocean (following the behavior of the CO2 constraint) rather than added to or debited from the atmosphere, because the old behavior could result in runaway atmospheric CO2 under certain conditions. This was fully documented in #659 .
This does not change any behavior except when running under this constraint.
Closes #659