anchor.diff:                 (dist < solvent_colloid_lj%cut(to_species, colloids%species(i)))
          end do
       end do
  -    solvent%pos(:,bound_molecule_idx) = x_new
  +
  +    i = solvent%id_to_idx(bound_molecule_id)
  +
  +    solvent%species(i) = to_species
  +    solvent%pos(:,i) = x_new
   
       ! Use c.o.m. velocity so that no kinetic energy exchange must take place
  -    solvent%vel(:,bound_molecule_idx) = colloids%vel(:,2)
  +    solvent%vel(:,i) = colloids%vel(:,2)
   
       enzyme_bound = .false.
   
  
positive.code:     ! transfer mass
      colloids%mass(1) = colloids%mass(1) + 1
  
      bound_molecule_idx = idx
      solvent%species(idx) = 0
  
      call add_energy_to_cell(cell_idx, excess_kinetic_energy)
  
      enzyme_bound = .true.
  
    end subroutine bind_molecule
  
  
negative.code:   type(PTo) :: config
  
    integer :: i, L(3)
    double precision :: current_time
    integer :: j, k
    type(timer_t), target :: varia, main, time_flag, time_change
    double precision :: total_time
  