I am trying to use OrdinaryDiffEq for the time integration of a hyperbolic PDE system in a finite-volume context. I want to control the size of the time-step through a CFL condition, rather than by some error norm. The calculation of this desired time-step is simple for a given state U, but I can’t see a way to specify that time-step to the solver.