Rev 591

* Upstream closes #7
* Proper initializations in readsub.f
* Fix for reservoir outflow