cosmo cluster neighborhood hdr

function nbrhood = cosmo_cluster_neighborhood(ds, varargin)