anchor.diff:          end do
           call n_solvent_el%append(n_solvent)
           call link_angle_el%append(link_angle)
  +        if (modulo(i, colloid_sampling)==0) then
  +           call dimer_io%position%append(colloids%pos)
  +           call dimer_io%velocity%append(colloids%vel)
  +           call dimer_io%image%append(colloids%image)
  +        end if
  +
        end if
   
     end do
  
positive.code:   type(thermo_t) :: thermo_data
    double precision :: temperature, kin_e
    double precision :: v_com(3)
    double precision :: com_pos(3)
    double precision :: unit_r(3)
    double precision :: tmp_vec(3), normal_vec(3), tmp_q(4)
    type(particle_system_io_t) :: dimer_io
    type(particle_system_io_t) :: solvent_io
    double precision :: bulk_rate(2)
    logical :: bulk_rmpcd
  
  
negative.code:   !! Histogram variable
    type(histogram_t), allocatable :: radial_hist(:)
    double precision, allocatable :: dummy_hist(:,:,:)
  
    integer :: equilibration_loops
    integer :: colloid_sampling
  