Impossible to modify mass rate with additional state derivatives
When using AdditionalStateDerivatives
in propagation, it is possible via combinedDerivatives
to update the rates for the dependent variables of integration e.g. equinoctial elements. However there is no way to change the one for mass. This is really a drawback for optimal control which is an application mentioned in the doc.