Bug fixes and saturated zone improvements

[waternumbers]

Fix a number of small bugs and change the saturated zone code to keep mass balance when root is only solved approximately