DifferentialEquations - Derivatives in ODE function/ nesting AD

Even with the new AD, this will be needed because the missing component cannot be calculated without a derivative rule to catch the discontinuity. The issue to follow is:

An implementation of Enright’s corrections is needed regardless of what AD is used, so it would be much appreciated!