function nbrhood = cosmo_interval_neighborhood(ds, label, varargin)
cosmo index unique hdr
cosmo isequaln hdr