anchor.diff:      integer :: m
       double precision :: x(3)
   
  +    call change_timer%tic()
       catalytic_change = 0
       !$omp parallel do private(x, dist_to_C_sq, dist_to_N_sq) reduction(+:catalytic_change)
       do m = 1, solvent% Nmax
  
positive.code:           end if
         end if
      end do
  
    end subroutine change_species
  
  
negative.code:     integer :: r, s
      double precision :: x(3)
  
      call flag_timer%tic()
      do s = 1,neigh% n(1)
         r = neigh%list(s,1)
         if (solvent% species(r) == 1) then
  