Bugfix, correcting incomplete merge.
Guard against "too small" max saturation.

This parameter might well be user configurable ...
Fixing some syntax.
Update of hysteresis model is now performed prior

to each time step, thus also capturing the initial state.