subroutine end_verlet_sweep() ! switch back to runge-kutta after adjusting alpha to get a new delta_psi write(stdout, *) "Switching to runge-kutta for update of alpha" tVerletSweep = .false. call assign_value_to_tau(tau / iterInit, 'end_verlet_sweep') end subroutine end_verlet_sweep