anchor.diff:     print *, '# Simulation parameters'
      print *, '# dt', dt
      print *, '# npoints', npoints
  -   print *, '# maxiter', maxiter
  +   print *, '# maxtime', maxtime
   
      !---------------------------------------------
      ! Open files for output, add headers
  
positive.code:       iter = iter + 1
  
        if (mod(iter, dump_freq) == 0) then
           print *, 'Dump', iter
           !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)
  
negative.code:    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         :: maxtime = npoints/100 ! Maximum time (ad hoc)
  
     !---------------------------------------------
     ! Read initial positions
  