cosmo singleton neighborhood hdr

function nh = cosmo_singleton_neighborhood(ds)