Closes #1002 (closed).
Allow using Richards model with the TwoPImmiscble fluid system. Uses more general indices and asserts that the phases are as assumed before.
Theoretically not backwards compatible if someone had a custom fluid system where phases are exchanged. But the error message is clear due to the static assert and it should be easy to fix then.