anchor.diff:    call init_genrand(mt, int(clock, c_long))
   
     call h5open_f(error)
  +
  +  prob = PTread_d(config,'probability')
   
     L = PTread_ivec(config, 'L', 3)
     rho = PTread_i(config, 'rho')
  
positive.code: 
    subroutine flag_particles
    double precision :: dist_to_C_sq
    real :: prob 
    real :: rndnumbers(solvent% Nmax)
    integer :: r
    prob = 1.0
    
    call random_number(prob)
    
    do  r = 1,solvent% Nmax
       if (solvent% species(r) == 1) then
         dist_to_C_sq = dot_product(colloids% pos(:,1) - solvent% pos(:,r),colloids% pos(:,1) - solvent% pos(:,r))
         if (dist_to_C_sq < 3**2) then
           if (rndnumbers(r) <= prob) then
             solvent% flag(r) = 1 
           end if
  
negative.code: 
    double precision :: e1, e2
    double precision :: tau, dt , T
    double precision :: d,prob
    double precision :: skin, co_max, so_max
    integer :: N_MD_steps, N_loop
    integer :: n_extra_sorting
  