function [offsets, distances] = cosmo_sphere_offsets(radius)
cosmo slice hdr
cosmo spherical neighborhood hdr