debugging optimizer

da/tools/pipeline.py

@@ -175,6 +175,7 @@ def SampleState(DaCycle,Samples,StateVector, ObservationOperator):
 def SampleOneCycle(DaCycle,Samples,StateVector, ObservationOperator,lag):
     """ Perform all actions needed to sample one cycle """
     from da.tools.general import AdvanceTime
+    import copy
 
     if lag>0:print "1", lag, StateVector.EnsembleMembers[0][0].ModelSample.Data.getvalues('obs').sum()
 
@@ -239,7 +240,7 @@ def SampleOneCycle(DaCycle,Samples,StateVector, ObservationOperator,lag):
     for Member in members:
         filename               = os.path.join(DaCycle['dir.output'],'samples.%03d.nc'%Member.membernumber)
         dummy                  = Samples.AddSimulations(filename,silent)
-        Member.ModelSample     = copy(Samples)
+        Member.ModelSample     = copy.copy(Samples)
         silent=True
 
     StateVector.nobs       += Samples.getlength()