reading modeldatamismatch right after reading obs now

da/tools/pipeline.py

@@ -207,6 +207,10 @@ def SampleOneCycle(DaCycle,Samples,StateVector, ObservationOperator,lag):
 
     status   = Samples.AddObs() 
 
+    # Add model-data mismatch to all samples, this *might* use output from the ensemble in the future??
+
+    status   = Samples.AddModelDataMismatch()
+
     filename   = Samples.WriteSampleInfo() 
 
     # Write filename to DaCycle, and to output collection list
@@ -219,9 +223,9 @@ def SampleOneCycle(DaCycle,Samples,StateVector, ObservationOperator,lag):
 
     # Add model-data mismatch to all samples, this *might* use output from the ensemble in the future??
 
-    status   = Samples.AddModelDataMismatch()
+    #status   = Samples.AddModelDataMismatch()
 
-    msg = "Added Model Data Mismatch to all samples " ; logging.debug(msg)
+    #msg = "Added Model Data Mismatch to all samples " ; logging.debug(msg)
 
     # Read forecast model samples that were written to NetCDF files by each member. Add them to the exisiting
     # Observation object for each sample loop. This data fill be written to file in the output folder for each sample cycle.