anchor.diff:     real(kind = xp) :: E                 ! Total energy
      integer         :: iter      = 0     ! Number of iterations ran
      integer         :: dump_freq = 10    ! Frequency at which the system is sampled
  -   integer         :: maxiter   = 500   ! Maximum number of iterations
  +   integer         :: maxtime = npoints/100 ! Maximum time (ad hoc)
   
      !---------------------------------------------
      ! Read initial positions
  
positive.code:    !---------------------------------------------
     ! Loop over time
     !---------------------------------------------
     do while (iter < maxiter)
  
        call integrate(v, a, dt/2)  ! Compute v(t+dt/2)
        call integrate(r, v, dt)    ! Compute r(t+dt)
  
negative.code:       iter = iter + 1
  
        if (mod(iter, dump_freq) == 0) then
           print *, 'Dump', iter, t
           !call compute_energy(m, r, v, Ec, Ep, E) ! Compute Ec, Ep, E at t+dt
           call compute_energy_omp(m, r, v, Ec, Ep, E) ! Compute Ec, Ep, E at t+dt
           call write_dump(una, un, iter, Ec, Ep, E, t, r, v)
  