debugging optimizer

da/tools/pipeline.py

@@ -161,6 +161,7 @@ def SampleState(DaCycle,Samples,StateVector, ObservationOperator):
         msg     = header+".....Ensemble Kalman Filter at lag %d"% (int(lag)+1,)           ; logging.info(msg) 
 
         SampleOneCycle(DaCycle,Samples,StateVector, ObservationOperator,lag)
+        print StateVector.EnsembleMembers[lag][0].ModelSample.Data.getvalues('obs')
 
     # Optionally, post-processing of the model output can be added that deals for instance with
     # sub-sampling of time series, vertical averaging, etc.