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:   call h5gclose_f(timers_group, error)
  
    call dimer_io%close()
    call solvent_io%close()
    call n_solvent_el%close()
    call link_angle_el%close()
    call hfile%close()
  
negative.code:         end do
          call varia%tac()
  
       end do md_loop
  
       call solvent_cells%random_shift(state(1))
  