anchor.diff:    type(PTo) :: config
   
     integer :: i, L(3)
  +  double precision :: current_time
     integer :: j, k
     type(timer_t), target :: varia, main, time_flag, time_change
     double precision :: total_time
  
positive.code: 
       if (i>equilibration_loops) call select_substrate
  
       call varia%tic()
  
       if (sampling) then
  
negative.code: 
    proba_s = PTread_d(config, 'proba_s', loc=params_group)
    proba_p = PTread_d(config, 'proba_p', loc=params_group)
    rate_release_s = PTread_d(config, 'rate_release_s', loc=params_group)
    rate_release_p = PTread_d(config, 'rate_release_p', loc=params_group)
    total_rate = rate_release_s + rate_release_p
    enzyme_capture_radius = PTread_d(config, 'enzyme_capture_radius', loc=params_group)
    substrate_fraction = PTread_d(config, 'substrate_fraction', loc=params_group)
  
  