Index: wflow-py/wflow/wflow_fit.py
===================================================================
diff -u -rc1ec494ff17c68d792fe1cfd1bc74854b420cbc5 -rcca16d276b7a138c6857735b88102f24cab28bf2
--- wflow-py/wflow/wflow_fit.py	(.../wflow_fit.py)	(revision c1ec494ff17c68d792fe1cfd1bc74854b420cbc5)
+++ wflow-py/wflow/wflow_fit.py	(.../wflow_fit.py)	(revision cca16d276b7a138c6857735b88102f24cab28bf2)
@@ -129,6 +129,7 @@
             mod = __import__(modeltofit, globals(), locals(), [], -1)
             self.WF = mod
 
+       
         self.results= []
         self.runId = runId
         self.caseName = casename
@@ -144,6 +145,8 @@
         self.dynModelFw = self.WF.wf_DynamicFramework(self.myModel, self.stopTime,self.startTime)
         # Load model config from files and check directory structure
         self.dynModelFw.createRunId(NoOverWrite=False,level=30)
+        
+    
         #self.dynModelFw.logger.setLevel(20)
         self.log = self.dynModelFw._userModel().logger
         self.conf  = self.dynModelFw._userModel().config
@@ -229,11 +232,17 @@
         
         for ts in range(self.startTime,self.stopTime):
             self.dynModelFw._runDynamic(ts,ts) # runs for all timesteps
+        
+        
 
-        results, head = pcrut.readtss(os.path.join(self.caseName,self.runId,"run.tss"))
+      
         # save output state
         self.dynModelFw._runSuspend()
         self.dynModelFw._wf_shutdown()
+
+        tssfile = os.path.join(self.caseName,self.runId,"run.tss")
+        
+        results, head = pcrut.readtss(tssfile)
         return results[self.WarmUpSteps:,self.ColSim].astype(np.float64)
       
     def savemaps(self,pars,savetoinput=False):