function ds_splits = cosmo_split(ds, split_by, dim, check)
cosmo spherical neighborhood hdr
cosmo squareform hdr