anchor.diff:    type(h5md_element_t) :: dummy_element
     integer(HID_T) :: fields_group, params_group
     integer(HID_T) :: box_group
  +  integer(HID_T) :: dummy_id
     type(thermo_t) :: thermo_data
     double precision :: temperature, kin_e
     double precision :: v_com(3)
  
positive.code:         com_pos = colloids%pos(:,2)
          call axial_cf%add(i-equilibration_loops, com_pos, unit_r)
  
       end if
  
       call varia%tac()
  
negative.code:   call solvent_io%velocity%append(solvent%vel)
    call solvent_io%image%append(solvent%image)
    call solvent_io%species%append(solvent%species)
  
    call h5gcreate_f(hfile%id, 'fields', fields_group, error)
  
    ! copy radial histogram data
    allocate(dummy_hist(size(radial_hist(1)%data, 1), &
         size(radial_hist(1)%data, 2), size(radial_hist)))
    do i = 1, size(radial_hist)
       call correct_radial_histogram(radial_hist(i))
       dummy_hist(:,:,i) = radial_hist(i)%data / dble(N_loop)
    end do
  
    call dummy_element%create_fixed(fields_group, 'radial_histogram', dummy_hist)
  
    call h5oopen_f(fields_group, 'radial_histogram', dummy_id, error)
    call h5md_write_attribute(dummy_id, 'xmin', radial_hist(1)%xmin)
    call h5md_write_attribute(dummy_id, 'dx', radial_hist(1)%dx)
    call h5oclose_f(dummy_id, error)
  
    call h5gclose_f(fields_group, error)
  
  
    ! Store timing data
    call timer_list%append(solvent%time_stream)
  