anchor.diff:    ! Print out initial particle positions
     frame = 0
     call particle_write
  +
  +  ! Initialize timers
  +  neigh_t = 0.0_WP
  +  dens_t = 0.0_WP
  +  force_t = 0.0_WP
  +  ef_t = 0.0_WP
  +  ps_t = 0.0_WP
  +  tot_t = 0.0_WP
   
     return
   end subroutine init
  
positive.code:   use posvel
    use neighbors
    use forces
    implicit none
    character(30) :: init_file, proc_string
    integer :: i, unit
  
negative.code:                 cbinz = binz + stz
                  do j = 1, part_count(cbinx,cbiny, cbinz)
                    nb = binpart(j,cbinx,cbiny,cbinz)
                    if(nb .le. p) cycle
                    distx = px(p) - px(nb)
                    disty = py(p) - py(nb)
                    distz = pz(p) - pz(nb)
  