Index: examples/wflow_maas_topoflex/wflow_topoflex.ini
===================================================================
diff -u -r9cda4f050ab505ec75d1cf9064980e9182f83aca -rebc36a6b97f3a09586461f97bf0ef371f9fb669c
--- examples/wflow_maas_topoflex/wflow_topoflex.ini	(.../wflow_topoflex.ini)	(revision 9cda4f050ab505ec75d1cf9064980e9182f83aca)
+++ examples/wflow_maas_topoflex/wflow_topoflex.ini	(.../wflow_topoflex.ini)	(revision ebc36a6b97f3a09586461f97bf0ef371f9fb669c)
@@ -25,6 +25,8 @@
 spinUp_time = 5
 NSEbreak = 0
 
+#intbl = intbl_tanja
+
 #W=wetland(0) H=hillslope(1)  P=plateau(2)  WD = drained wetland(3) 
 classes = ['W','H','P']
 timestepsecs = 3600
@@ -67,13 +69,16 @@
 
 [run]
 runinfofile = inmaps/runinfo.xml
+#reinit = 1
 
 [layout]
 # if set to zero the cell-size is given in lat/long (the default)
 sizeinmetres = 1
 
 [outputmaps]
 
+self.Ss = Ss
+self.Si[0] = Si
 # List all timeseries in tss format to save in this section. Timeseries are
 # produced per subcatchment. 
 
Index: wflow-py/wflow/wflow_topoflex.py
===================================================================
diff -u -r9cda4f050ab505ec75d1cf9064980e9182f83aca -rebc36a6b97f3a09586461f97bf0ef371f9fb669c
--- wflow-py/wflow/wflow_topoflex.py	(.../wflow_topoflex.py)	(revision 9cda4f050ab505ec75d1cf9064980e9182f83aca)
+++ wflow-py/wflow/wflow_topoflex.py	(.../wflow_topoflex.py)	(revision ebc36a6b97f3a09586461f97bf0ef371f9fb669c)
@@ -816,8 +816,13 @@
             usage()
             return
 
-    opts, args = getopt.getopt(argv, 'C:S:T:Ic:s:R:F:fl:L:P:p:')
+    try:
+        opts, args = getopt.getopt(argv, 'C:S:T:Ic:s:R:F:fl:L:P:p:i:') #'XF:L:hC:Ii:v:S:T:WR:u:s:EP:p:Xx:U:fOc:l:')
+    except getopt.error, msg:
+        pcrut.usage(msg)
 
+    # opts, args = getopt.getopt(argv, 'C:S:T:Ic:s:R:F:fl:L:P:p:')
+
     for o, a in opts:
         if o == '-F': 
             runinfoFile = a
@@ -865,7 +870,7 @@
             right = a.split('=')[1]
             configset(myModel.config,'variable_change_timestep',left,right,overwrite=True)
         if o == '-X': configset(myModel.config, 'model', 'OverWriteInit', '1', overwrite=True)
-        if o == '-I': configset(myModel.config, 'model', 'reinit', '1', overwrite=True)
+        if o == '-I': configset(myModel.config, 'run', 'reinit', '1', overwrite=True)
         if o == '-i': configset(myModel.config, 'model', 'intbl', a, overwrite=True)
         if o == '-s': configset(myModel.config, 'model', 'timestepsecs', a, overwrite=True)
         if o == '-x': configset(myModel.config, 'model', 'sCatch', a, overwrite=True)