anchor.diff:             next_reaction_time(enzyme_i) = huge(next_reaction_time(enzyme_i))
           end if
        end do
  +     call time_release%tac()
   
        call solvent_cells%random_shift(state(1))
        call apply_pbc(solvent, solvent_cells% edges)
  
positive.code:   double precision :: current_time
    integer :: j, k
    type(timer_t), target :: varia, main, time_flag, time_change
    double precision :: total_time
    type(timer_list_t) :: timer_list
    integer(HID_T) :: timers_group
  
negative.code:   call varia%init('varia')
    call time_flag%init('flag')
    call time_change%init('change')
    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', &
  