def get_ind_sphere(mesh, ind_active, origin, radius):
    """Retreives the indices of a sphere object coordintes in a mesh."""

    return (
            (mesh.gridCC[ind_active, 0] <= origin[0] + radius)
            & (mesh.gridCC[ind_active, 0] >= origin[0] - radius)
            & (mesh.gridCC[ind_active, 1] <= origin[1] + radius)
            & (mesh.gridCC[ind_active, 1] >= origin[1] - radius)
            & (mesh.gridCC[ind_active, 2] <= origin[2] + radius)
            & (mesh.gridCC[ind_active, 2] >= origin[2] - radius)
    )