debugging optimizer

da/tools/pipeline.py

@@ -204,6 +204,8 @@ def SampleOneCycle(DaCycle,Samples,StateVector, ObservationOperator,lag):
 
     dummy   = WriteEnsembleData(DaCycle, StateVector, lag)
 
+    if lag>0:print "2", lag, StateVector.EnsembleMembers[0][0].ModelSample.Data.getvalues('obs').sum()
+
     dummy   = Samples.Initialize(DaCycle)
 
     dummy   = Samples.Validate() 
@@ -214,7 +216,6 @@ def SampleOneCycle(DaCycle,Samples,StateVector, ObservationOperator,lag):
 
     DaCycle['ObsOperator.inputfile'] = filename
 
-    if lag>0:print "2", lag, StateVector.EnsembleMembers[0][0].ModelSample.Data.getvalues('obs').sum()
     # Run the observation operator
 
     dummy   = RunForecastModel(DaCycle,ObservationOperator)