Add discontinuity time list to `odesolve`

Let a user specify a list of times for `odesolve` that are matched precisely during time stepping. Spitfire makes sure to nail the exact final time specified, so a user can already do this themselves by restarting with several `odesolve` calls, but it would be nice to provide the feature for them and automatically handle output concatenation, etc.