Index: wflow-py/Sandbox/wflow_topoflex/wflow_topoflex.py
===================================================================
diff -u -r2b8f6107e313a4df0cc695b144b6249d4afea7f1 -r878b2936e851fa3890ca517ade8a9effe1c6c63a
--- wflow-py/Sandbox/wflow_topoflex/wflow_topoflex.py	(.../wflow_topoflex.py)	(revision 2b8f6107e313a4df0cc695b144b6249d4afea7f1)
+++ wflow-py/Sandbox/wflow_topoflex/wflow_topoflex.py	(.../wflow_topoflex.py)	(revision 878b2936e851fa3890ca517ade8a9effe1c6c63a)
@@ -316,8 +316,8 @@
                                       "model", "wflow_surfaceArea", "staticmaps/wflow_surfaceArea.map")
         wflow_transit = configget(self.config,
                                   "model", "wflow_transit", "staticmaps/wflow_transit.map")
-        wflow_velocity = configget(self.config,
-                                  "model", "wflow_velocity", "staticmaps/wflow_velocity.map")
+        #wflow_velocity = configget(self.config,
+                                  #"model", "wflow_velocity", "staticmaps/wflow_velocity.map")
         wflow_percent = [configget(self.config,
                                    "model", "wflow_percent_" + str(self.Classes[i]),
                                    "staticmaps/wflow_percent" + str(self.Classes[i]) + ".map") for i in self.Classes]
@@ -340,7 +340,7 @@
         self.totalArea = areatotal(self.surfaceArea, nominal(self.TopoId))
         self.percentArea = self.surfaceArea / self.totalArea
         self.Transit = scalar(readmap(os.path.join(self.Dir, wflow_transit)))  #: Map with surface area per cell
-        self.velocity = scalar(readmap(os.path.join(self.Dir, wflow_velocity)))  #: Map with surface area per cell
+        #self.velocity = scalar(readmap(os.path.join(self.Dir, wflow_velocity)))  #: Map with surface area per cell
         self.gaugesR = nominal(readmap(os.path.join(self.Dir, wflow_gauges)))
         self.percent = []
         for i in self.Classes:
@@ -607,14 +607,14 @@
         eval(eval_str)
 
         # ROUTING
-        self.Qtot = self.Qftotal + self.Qs_  # total local discharge in mm/hour
-        self.Qtotal = self.Qtot / 1000 * self.surfaceArea / self.timestepsecs  # total local discharge in m3/s
-        self.Qstate_t = self.Qstate
-        self.Qrout = accutraveltimeflux(self.TopoLdd, self.Qstate + self.Qtotal, self.velocity)
-        self.Qstate = accutraveltimestate(self.TopoLdd, self.Qstate + self.Qtotal, self.velocity)
-        # water balance of flux routing
-        self.dSdt = self.Qstate-self.Qstate_t
-        self.WB_rout = (accuflux(self.TopoLdd, self.Qtotal - self.dSdt)-self.Qrout)/accuflux(self.TopoLdd, self.Qtotal)
+        #self.Qtot = self.Qftotal + self.Qs_  # total local discharge in mm/hour
+        #self.Qtotal = self.Qtot / 1000 * self.surfaceArea / self.timestepsecs  # total local discharge in m3/s
+        #self.Qstate_t = self.Qstate
+        #self.Qrout = accutraveltimeflux(self.TopoLdd, self.Qstate + self.Qtotal, self.velocity)
+        #self.Qstate = accutraveltimestate(self.TopoLdd, self.Qstate + self.Qtotal, self.velocity)
+        ## water balance of flux routing
+        #self.dSdt = self.Qstate-self.Qstate_t
+        #self.WB_rout = (accuflux(self.TopoLdd, self.Qtotal - self.dSdt)-self.Qrout)/accuflux(self.TopoLdd, self.Qtotal)
         self.Qtlag = self.Qflag_ / self.timestepsecs + self.Qs_ / 1000 * self.surfaceArea / self.timestepsecs
         self.QLagTot = areatotal(self.Qtlag, nominal(self.TopoId))  # catchment total runoff with looptijd
 
@@ -641,8 +641,8 @@
         self.Ei = areatotal(sum(multiply(self.Ei_, self.percent)) / 1000 * self.surfaceArea, nominal(self.TopoId))
         self.Eu = areatotal(sum(multiply(self.Eu_, self.percent)) / 1000 * self.surfaceArea, nominal(self.TopoId))
         self.Er = areatotal(sum(multiply(self.Er_, self.percent)) / 1000 * self.surfaceArea, nominal(self.TopoId))
-        #self.Qtot = self.QLagTot * self.timestepsecs
-        self.QtotnoRout = areatotal(self.Qtotal, nominal(self.TopoId))
+        self.Qtot = self.QLagTot * self.timestepsecs
+        #self.QtotnoRout = areatotal(self.Qtotal, nominal(self.TopoId))
         self.SiWB = areatotal(sum(multiply(self.Si, self.percent)) / 1000 * self.surfaceArea, nominal(self.TopoId))
         self.Si_WB = areatotal(sum(multiply(self.Si_t, self.percent)) / 1000 * self.surfaceArea, nominal(self.TopoId))
         self.SuWB = areatotal(sum(multiply(self.Su, self.percent)) / 1000 * self.surfaceArea, nominal(self.TopoId))