anchor.diff:             call apply_pbc(colloids, solvent_cells% edges)
              call solvent% sort(solvent_cells)
              call neigh% update_list(colloids, solvent, max_cut + skin, solvent_cells, solvent_colloid_lj)
  -           call varia%tic()
  +           call so_timer%tic()
              !$omp parallel do
              do k = 1, solvent%Nmax
                 solvent% pos_old(:,k) = solvent% pos(:,k)
              end do
  +           call so_timer%tac()
  +           call varia%tic()
              colloids% pos_old = colloids% pos
              n_extra_sorting = n_extra_sorting + 1
              call varia%tac()
  
positive.code: 
       if (sampling) then
          temperature = compute_temperature(solvent, solvent_cells)
          call varia%tic()
          kin_e = 0
          v_com = 0
          !$omp parallel do private(k) reduction(+:kin_e) reduction(+:v_com)
          do k = 1, solvent%Nmax
             kin_e = kin_e + sum(solvent%vel(:,k)**2)
             v_com = v_com + solvent%vel(:,k)
          end do
          tmp_mass = solvent%Nmax
          do k = 1, colloids%Nmax
             v_com = v_com + colloids%mass(colloids%species(k)) * colloids%vel(:,k)
  
negative.code: 
    call main%init('main')
    call varia%init('varia')
    call so_timer%init('so_timer')
    call time_flag%init('flag')
    call time_refuel%init('refuel')
    call time_change%init('change')
    call bulk_reac_timer%init('bulk_reac')
    call q_timer%init('quaternion_vv')
  
    call timer_list%init(24)
    call timer_list%append(varia)
    call timer_list%append(so_timer)
    call timer_list%append(time_flag)
    call timer_list%append(time_refuel)
    call timer_list%append(time_change)
  