anchor.diff:          co_max = colloids% maximum_displacement()
   
           if ( (co_max >= skin/2) .or. (so_max >= skin/2) ) then
  -           call varia%tic()
              call apply_pbc(solvent, solvent_cells% edges)
              call apply_pbc(colloids, solvent_cells% edges)
  -           call varia%tac()
              call solvent% sort(solvent_cells)
              call neigh% update_list(colloids, solvent, max_cut + skin, solvent_cells, solvent_colloid_lj)
              call varia%tic()
  
positive.code: 
       end do md_loop
  
       call varia%tic()
  
       if (sampling) then
          temperature = compute_temperature(solvent, solvent_cells)
          kin_e = sum(solvent% vel**2)
          v_com = sum(solvent%vel, dim=2)
          tmp_mass = solvent%Nmax
  
negative.code: 
          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 apply_pbc(solvent, solvent_cells% edges)
       call apply_pbc(colloids, solvent_cells% edges)
  