We need to add a check in the python layer in nu_masses which throws an error if you pass a combination of Omeganuh2 and the mode of neutrino splitting which is not physical. Right now it returns nan or a meaningless error.

This is because this test used to be in the C layer and now relies on checks at a python level which haven't ever been added in this function.