anchor.diff:    double precision :: d,prob
     double precision :: skin, co_max, so_max
     integer :: N_MD_steps, N_loop
  +  integer :: colloid_sampling
     integer :: n_extra_sorting
     double precision :: kin_e, temperature
     integer, dimension(N_species) :: n_solvent, catalytic_change, bulk_change
  
positive.code:            end if
          end if
  
       end do md_loop
  
       call random_number(solvent_cells% origin)
  
negative.code:      call catalytic_change_el%append(catalytic_change)
       call bulk_change_el%append(bulk_change)
  
       call h5fflush_f(hfile%id, H5F_SCOPE_GLOBAL_F, error)
  
       if (fixed) then
  