anchor.diff:                     end if
                      list_p(n_p) = m
                   end if
  -             end if
  -          else
  -             if (dist > enzyme_capture_radius) then
  -                solvent%flags(m) = ibclr(solvent%flags(m), ENZYME_REGION_BIT)
  +             else
  +                if (dist > enzyme_capture_radius) then
  +                   solvent%flags(m) = ibclr(solvent%flags(m), ENZYME_REGION_BIT)
  +                end if
                end if
             end if
          end do select_loop
  
positive.code:   skin = 1
    n_extra_sorting = 0
  
    call neigh% make_stencil(solvent_cells, max_cut+skin)
  
    call neigh% update_list(colloids, solvent, max_cut+skin, solvent_cells, solvent_colloid_lj)
  
    e1 = compute_force(colloids, solvent, neigh, solvent_cells% edges, solvent_colloid_lj)
    e2 = compute_force_n2(colloids, solvent_cells% edges, colloid_lj)
  
negative.code:            call apply_pbc(solvent, solvent_cells% edges)
             call apply_pbc(colloids, solvent_cells% edges)
             call solvent% sort(solvent_cells)
             call neigh% update_list(colloids, solvent, enzyme_capture_radius + skin, solvent_cells, solvent_colloid_lj)
             call varia%tic()
             !$omp parallel do
             do k = 1, solvent%Nmax
  