There was an issue in grouping.spherical_average: if keyword little_h is False in function, then the input k-bins are automatically adjusted to be in h Mpc-1 units:

# transform kgrid to little_h units
 if not little_h:
      kbins = kbins / uvp.cosmo.h
      bin_widths = bin_widths / uvp.cosmo.h

However this change occurred after defining kbin_left and kbin_right that are later used to bin. I changed this and moved the above code block up.