anchor.diff:    call time_flag%init('flag')
     call time_refuel%init('refuel')
     call time_change%init('change')
  -
  -  call timer_list%init(16)
  +  call bulk_reac_timer%init('bulk_reac')
  +
  +  call timer_list%init(19)
     call timer_list%append(varia)
     call timer_list%append(time_flag)
     call timer_list%append(time_refuel)
     call timer_list%append(time_change)
  +  call timer_list%append(bulk_reac_timer)
   
     call h5open_f(error)
     call hfile%create(args%output_file, 'RMPCDMD::single_janus_pbc', &
  
positive.code: 
    double precision :: total_time
    type(timer_t), target :: varia, main, time_flag, time_refuel, time_change
    type(timer_list_t) :: timer_list
    integer(HID_T) :: timers_group
  
  
negative.code: 
       call simple_mpcd_step(solvent, solvent_cells, state)
  
       call bulk_reac_timer%tic()
       call bulk_reaction(solvent, solvent_cells, 2, 1, bulk_rate, tau, state)
       call bulk_reac_timer%tac()
  
       if (sampling) then
          n_solvent = 0
  