anchor.diff:          n_solvent = 0
           do k = 1, solvent%Nmax
              j = solvent%species(k)
  -           if (j <= 0) continue
  +           if (j <= 0) cycle
              n_solvent(j) = n_solvent(j) + 1
           end do
           call n_solvent_el%append(n_solvent)
  
positive.code:   type(PTo) :: config
  
    integer :: i, L(3)
    integer :: j, k
    type(timer_t), target :: varia, main, time_flag, time_change
    double precision :: total_time
  
negative.code:     ! transfer mass
      colloids%mass(1) = colloids%mass(1) + 1
  
      bound_molecule_id = solvent%id(idx)
      solvent%species(idx) = 0
  
      call add_energy_to_cell(cell_idx, excess_kinetic_energy)
  
      enzyme_bound = .true.
  
      next_reaction_time = current_time - log(threefry_double(state(1)))/total_rate ! sample from Poisson process
  
    end subroutine bind_molecule
  
  