anchor.diff:          v_com = v_com / tmp_mass
           kin_e = kin_e / 2
   
  -        call thermo_data%append(hfile, temperature, e1+e2+e3, kin_e, e1+e2+e3+kin_e, v_com)
  +        call thermo_data%append(hfile, temperature, e1+e2+e3+e_wall, kin_e, e1+e2+e3+e_wall+kin_e, v_com)
           call axial_cf%add(i-equilibration_loops, com_pos, unit_r)
           call varia%tac()
        end if
  
positive.code:   else
       e3 = 0
    end if
    solvent% force_old = solvent% force
    colloids% force_old = colloids% force
    if (do_quaternion) call rigid_janus%compute_force_torque(colloids)
  
negative.code:   type(neighbor_list_t) :: neigh
    type(lj_params_t) :: solvent_colloid_lj
    type(lj_params_t) :: colloid_lj
    type(lj_params_t) :: walls_colloid_lj
  
    integer, parameter :: N_species = 2
  
  