Index: wflow-py/wflow/reservoir_Si.py
===================================================================
diff -u -rfacc3a7f729f6910dac37d4761992c70d238fcbd -r5304872d008eedd1549905bf187ecd067a92a73a
--- wflow-py/wflow/reservoir_Si.py	(.../reservoir_Si.py)	(revision facc3a7f729f6910dac37d4761992c70d238fcbd)
+++ wflow-py/wflow/reservoir_Si.py	(.../reservoir_Si.py)	(revision 5304872d008eedd1549905bf187ecd067a92a73a)
@@ -39,7 +39,7 @@
     self.Pe_[k] = max(self.Precipitation, 0)
     self.Ei_[k] = 0.
     self.Si_[k] = 0.
-    self.wbSi_[k] = self.Precipitation - self.Ei - self.Pe - self.Si[k] + self.Si_t[k]
+    self.wbSi_[k] = self.Precipitation - self.Ei_[k] - self.Pe_[k] - self.Si[k] + self.Si_t[k]
 
 def interception_overflow2(self, k):
     """
@@ -55,8 +55,8 @@
     
     self.wbSi_[k] = self.Precipitation - self.Ei - self.Pe - self.Si[k] + self.Si_t[k]        
     
-    self.Pe = self.Pe + self.Qw         #added on 3-11-2015 for snow module
-    self.Ei = self.Ei + self.Ew
+    self.Pe = self.Pe + self.Qw_[k]         #added on 3-11-2015 for snow module
+    self.Ei = self.Ei + self.Ew_[k]
     
     self.Ei_[k]=self.Ei
     self.Pe_[k]=self.Pe
@@ -87,8 +87,8 @@
     
     self.wbSi_[k] = self.Precipitation - self.Ei - self.Pe - self.Si[k] + self.Si_t[k]        
     
-    self.Pe = self.Pe + self.Qw         #added on 3-11-2015 for snow module
-    self.Ei = self.Ei + self.Ew
+    self.Pe = self.Pe + self.Qw_[k]         #added on 3-11-2015 for snow module
+    self.Ei = self.Ei + self.Ew_[k]
     
     self.Ei_[k]=self.Ei
     self.Pe_[k]=self.Pe