anchor.diff:  
     double precision :: total_time
     type(timer_t), target :: varia, main, time_flag, time_refuel, time_change
  +  type(timer_t), target :: bulk_reac_timer
     type(timer_list_t) :: timer_list
     integer(HID_T) :: timers_group
   
  
positive.code: 
       call simple_mpcd_step(solvent, solvent_cells, state)
  
       call bulk_reaction(solvent, solvent_cells, 2, 1, bulk_rate, tau, state)
  
       if (sampling) then
          n_solvent = 0
  
negative.code:   call timer_list%append(solvent%time_md_vel)
    call timer_list%append(solvent%time_max_disp)
    call timer_list%append(colloids%time_self_force)
    call timer_list%append(colloids%time_rattle_pos)
    call timer_list%append(colloids%time_rattle_vel)
    call timer_list%append(neigh%time_update)
    call timer_list%append(neigh%time_force)
    call timer_list%append(colloids%time_max_disp)
  