Parameter interaction between `theta_sat` and `theta_mp` may inadvertently turn off macropore flow

[wknoben]

If theta_sat < theta_mp, macropore flow will never activate. This is effectively a model decision and should (probably) at least throw a warning if not a graceful exit.