I just noticed that the example notebook uses:
csm = {'omega_M': 0.3,
'omega_B': 0.3,
'omega_L': 0.7,
'h': 0.7,
'w': -1.0,
'ns': 0.96,
'sigma_8': 0.8}
leading to omega_C = omega_M - omega_B = 0. Is that intentional? The resulting linear matter power spectrum is quite weird:
