* Use l in krw/krn instead of hard-coded l=0.5
* Correct the non-wetting relperm to comply with the Mualem model
* Correct the non-wetting relperm derivative