anchor.diff:    double precision :: current_time
     integer :: j, k
     type(timer_t), target :: varia, main, time_flag, time_change
  +  type(timer_t), target :: time_select, time_release, time_reset
     double precision :: total_time
     type(timer_list_t) :: timer_list
     integer(HID_T) :: timers_group
  
positive.code: 
      integer :: enzyme_i, enz_2
      integer :: s_found, p_found
  
      select_substrate_loop: do enzyme_i = 1, N_enzymes
  
  
negative.code:            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)
  