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