init_verlet_sweep Subroutine

    subroutine init_verlet_sweep()
        ! here, we initiate the first dpsi_cache from a runge-kutta calculation
        character(*), parameter :: this_routine = "init_verlet_sweep"

        write(stdout, *) "Prepared initial delta_psi, starting verlet calculation"
        call build_initial_delta_psi()
        ! rescale the timestep
        call assign_value_to_tau(iterInit * tau, this_routine)
        tau_imag = iterInit * tau_imag
        tau_real = iterInit * tau_real
        ! There is only one step now (we might log the second spawns as quasi-second
        ! step later on
        runge_kutta_step = 1

    end subroutine init_verlet_sweep