anchor.diff:  
     call colloids% init(2,2, mass) !there will be 2 species of colloids
   
  -  !call PTkill(config)
  -  
  +  call PTkill(config)
  +
     open(17,file ='dimerdata_FullExp_1.txt')
     open(18,file ='dimerdata_FullExp_2.txt')
     open(19,file ='dimerdata_vx_flow_wall.txt')  
  
positive.code: 
    call h5open_f(error)
  
    g = [0.001d0, 0.d0, 0.d0] !PTread_ivec(config, 'g', 3) !place in inputfile
    bufferlength = 20 !PTread_i(config, 'bufferlength')
    prob = 1.d0 !PTread_d(config,'probability')
  
    L = [50,50,15] !PTread_ivec(config, 'L', 3)
    L(1) = L(1)+ bufferlength
    
    rho = 10 !PTread_i(config, 'rho')
    N = rho *L(1)*L(2)*L(3)
  
    T = 5.d0 !PTread_d(config, 'T')
    d = 4.5d0 !PTread_d(config, 'd')
  
    wall_v = 0
    wall_t = [T, T]
    
    tau =0.2d0 !PTread_d(config, 'tau')
    N_MD_steps = 100 !PTread_i(config, 'N_MD')
    dt = tau / N_MD_steps
    N_loop = 20000 !PTread_i(config, 'N_loop')
  
    
  
    sigma_C = 2.d0 !PTread_d(config, 'sigma_C')
    sigma_N = 2.d0 !PTread_d(config, 'sigma_N')
    
    epsilon(:,1) = 1 !PTread_dvec(config, 'epsilon_C', 2)
    epsilon(:,2) = [1.d0, 0.1d0, 1.d0] !PTread_dvec(config, 'epsilon_N', 2)
  
    sigma(:,1) = sigma_C
    sigma(:,2) = sigma_N
  
negative.code: 
    call colloid_lj% init(epsilon(1:2,:), sigma(1:2,:), sigma_cut(1:2,:))
  
    epsilon(1,1) = 1.d0
    epsilon(1,2) = 1.d0
    epsilon(2,1) = 1.d0
    epsilon(2,2) = 1.d0
    sigma(1,:) = [sigma_C, sigma_N]
    sigma(2,:) = [sigma_C, sigma_N]
    sigma_cut = sigma*2**(1.d0/6.d0)
  