anchor.diff:    call varia%init('varia')
     call time_flag%init('flag')
     call time_change%init('change')
  -
  -  call timer_list%init(17)
  +  call time_select%init('select_substrate')
  +  call time_release%init('release_substrate')
  +  call time_reset%init('reset_bit')
  +
  +  call timer_list%init(20)
     call timer_list%append(varia)
     call timer_list%append(time_flag)
     call timer_list%append(time_change)
  +  call timer_list%append(time_select)
  +  call timer_list%append(time_release)
  +  call timer_list%append(time_reset)
   
     call h5open_f(error)
     call hfile%create(args%output_file, 'RMPCDMD::three_bead_enzyme', &
  
positive.code: 
      end do select_substrate_loop
  
    end subroutine select_substrate
  
  
  
negative.code:      end do md_loop
  
  
       call time_release%tic()
       ! Check if unbinding should occur
       do enzyme_i = 1, N_enzymes
          if (current_time >= next_reaction_time(enzyme_i)) then
  