anchor.diff:        if (mod(iter, dump_freq) == 0) then
            print*, 'Dump!'
            call compute_energy(m, r, v, Ec, Ep, E) ! Compute Ec, Ep, E at t+dt
  -         call write_dump(11, 12, iter, Ec, Ep, E, t, r, v)
  +         call write_dump(una, un, iter, Ec, Ep, E, t, r, v)
         end if
   
      end do
   
  -   close(unit=11)
  -   close(unit=12)
  -
      !---------------------------------------------
  -   ! Deallocations
  +   ! Close files
      !---------------------------------------------
  -   deallocate(r)
  -   deallocate(v)
  -   deallocate(a)
  +   close(un)
  +   close(una)
   
      !---------------------------------------------
      ! Elapsed time
  
positive.code:    !---------------------------------------------
     call cpu_time(start_time)
  
     !---------------------------------------------
     ! Read number of points
     !---------------------------------------------
     open(unit=10, file='initial_conditions.dat')
     call read_npoints(10)
  
     !---------------------------------------------
     ! Allocations
     !---------------------------------------------
     allocate(m(npoints))
     allocate(r(npoints, 3))
     allocate(v(npoints, 3))
     allocate(a(npoints, 3))
  
     !---------------------------------------------
     ! Read initial positions
     !---------------------------------------------
     call read_mpos(10, m, r)
     close(unit=10)
  
     !---------------------------------------------
     ! Compute initial speeds, accelerations, variables
  
negative.code:    call compute_initial_variables()
  
     !---------------------------------------------
     ! Compute time step as a fraction of dynamic time
     !---------------------------------------------
     dt = 1.e-3_xp
     print *, dt
  
     !---------------------------------------------
     ! Open files for output, add headers
     !---------------------------------------------
     open(newunit=una, file='output.dat', status="replace")
     open(newunit=un, file='output_int.dat', status="replace")
  
     call write_dump_headers(un)
  
     !---------------------------------------------
     ! Loop over time
     !---------------------------------------------
     do while (iter < maxiter)
  
        call integrate(v, a, dt/2)  ! Compute v(t+dt/2)
  