Modern Fortran implementation of DOP853 ODE solver
