anchor.diff:    double precision :: link_treshold
     integer :: i_link
     double precision :: dist, rattle_pos_tolerance, rattle_vel_tolerance
  +  double precision :: unit_r(3)
  +  double precision :: com_pos(3)
  +  integer :: head
  +
  +  integer, parameter :: block_length = 8
  +  type(axial_correlator_t) :: axial_cf
  +  integer(HID_T) :: correlator_group
  +
  +  integer :: equilibration_loops
  +  integer :: colloid_sampling
  +  logical :: sampling
   
     type(args_t) :: args
     character(len=144) :: data_filename
  
positive.code:     end do
    end subroutine refuel
  
  end program single_janus_pbc
  
negative.code:   solvent% force_old = solvent% force
    colloids% force_old = colloids% force
  
    write(*,*) 'Running for', equilibration_loops, '+', N_loop, 'loops'
    call main%tic()
    sampling = .false.
    do i = 0, N_loop+equilibration_loops
       if (i==equilibration_loops) sampling = .true.
       if (modulo(i,32) == 0) write(*,'(i05)',advance='no') i
       md_loop: do j = 1, N_MD_steps
          call md_pos(solvent, dt)
  
  