function nbrhood = cosmo_spherical_neighborhood(ds, varargin)
cosmo sphere offsets hdr
cosmo split hdr