anchor.diff:  
           call thermo_data%append(hfile, temperature, e1+e2+e3, kin_e, e1+e2+e3+kin_e, v_com)
           call axial_cf%add(i-equilibration_loops, com_pos, unit_r)
  -
  +        call varia%tac()
        end if
   
        call solvent_cells%random_shift(state(1))
  -
  -     call varia%tac()
   
        call apply_pbc(solvent, solvent_cells% edges)
        call apply_pbc(colloids, solvent_cells% edges)
  
positive.code:   call time_change%init('change')
    call bulk_reac_timer%init('bulk_reac')
  
    call timer_list%init(20)
    call timer_list%append(varia)
    call timer_list%append(time_flag)
    call timer_list%append(time_refuel)
  
negative.code: 
       end do md_loop
  
       if (sampling) then
          temperature = compute_temperature(solvent, solvent_cells)
          call varia%tic()
          kin_e = sum(solvent% vel**2)
          v_com = sum(solvent%vel, dim=2)
          tmp_mass = solvent%Nmax
  