anchor.diff:    end if
     solvent% force_old = solvent% force
     colloids% force_old = colloids% force
  +  if (do_quaternion) call rigid_janus%compute_force_torque(colloids, solvent_cells%edges)
   
     write(*,*) 'Running for', equilibration_loops, '+', N_loop, 'loops'
     call main%tic()
  
positive.code:         end do
       end do
    end if
  
    write(*,*) 'number of links:', i_link
  
  
negative.code:   double precision :: unit_r(3)
    double precision :: com_pos(3), head_pos(3)
    integer :: head
    type(rigid_body_t) :: rigid_janus
    double precision :: quaternion_treshold
  
    integer, parameter :: block_length = 8
    type(axial_correlator_t) :: axial_cf
  