anchor.diff:     ! Compute time step as a fraction of dynamic time
      !---------------------------------------------
      dt = 1.e-3_xp
  -   print *, dt
  +   print *, '# Simulation parameters'
  +   print *, '# dt', dt
  +   print *, '# npoints', npoints
  +   print *, '# maxiter', maxiter
   
      !---------------------------------------------
      ! Open files for output, add headers
  
positive.code:    close(un)
     close(una)
  
     !---------------------------------------------
     ! Elapsed time
     !---------------------------------------------
     call cpu_time(total_time)
     total_time = total_time - start_time
     print*, 'Elapsed time:', total_time
  
  end program n_body
  
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         :: maxiter   = 500   ! Maximum number of iterations
  
     !---------------------------------------------
     ! Read initial positions
  