Shifted van-Genuchten curve added.  Fixes #3350

This is for 2-phase flow.  It may make the Jacobian nonsingular in
2phase models when one phase disappears and residual saturation is
zero.  Time will tell whether this is useful.

Included are doco and tests