The cbc or power spectrum is not working properly due to an issue with the interpolate function

[DQuanR2]

Hi, I found an issue when running the code from the examples e.g.,
python example.py -n 32 34 32 -l 0.565 0.6 0.565, this would call the cbc-spectrum as a default.
This crashes with a
ValueError: A value in x_new is below the interpolation range. given by the "scipy" function "interpolate".

This can be solved by adding the flag fill_value = "extrapolate" in the interp1d function:
self.especf = interpolate.interp1d(kcbc, ecbc,'cubic', fill_value="extrapolate")

Note:
I am unaware of the physical implication this might have.