Fixing a bug to find radial interpolation points in assemble_sph