anchor.diff:    integer(HID_T) :: box_group
     type(h5md_element_t) :: dummy_element
   
  +  integer(HID_T) :: fields_group
  +  type(histogram_t) :: radial_hist
  +  type(h5md_element_t) :: radial_hist_el
  +
     double precision :: skin
     double precision :: rsq
     logical :: tooclose
  
positive.code:   call dummy_element%create_fixed(box_group, 'edges', solvent_cells%edges)
    call h5gclose_f(box_group, error)
  
    i = 1
    place_colloids: do while (i<=N_colloids)
       colloids%pos(1,i) = threefry_double(state(1))*L(1)
  
negative.code: 
    call h5gclose_f(timers_group, error)
  
    call radial_hist_el%close()
    call colloids_io%close()
    call hfile%close()
    call h5close_f(error)
  