anchor.diff:  
     call h5gclose_f(timers_group, error)
   
  +  call radial_hist_el%close()
     call colloids_io%close()
     call hfile%close()
     call h5close_f(error)
  
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 thermo_data%append(hfile, temperature, e1+e2, kin_e, e1+e2+kin_e, v_com)
  
       radial_hist%data = 0
       call compute_radial_histogram(radial_hist, colloids%pos(:,1), &
            solvent_cells%edges, solvent)
       call radial_hist_el%append(radial_hist%data)
  
    end do
  
    call thermo_data%append(hfile, temperature, e1+e2, kin_e, e1+e2+kin_e, &
  