added diagnostic output from optimizer

da/baseclasses/optimizer.py

@@ -85,8 +85,9 @@ class Optimizer(object):
             members                                             = StateVector.EnsembleMembers[n]
             self.x[n*self.nparams:(n+1)*self.nparams]           = members[0].ParameterValues[n]
             self.X_prime[n*self.nparams:(n+1)*self.nparams,:]   = np.transpose(np.array([m.ParameterValues for m in members]))
-            allobs.extend(                                        StateVector.EnsembleMembers[n][0].ModelSample.Data.getvalues('obs') )
-            allsamples.extend(                                    StateVector.EnsembleMembers[n][0].ModelSample.Data.getvalues('simulated') )
+            allobs.extend(                                        members[0].ModelSample.Data.getvalues('obs') )
+            allsamples.extend(                                    members[0].ModelSample.Data.getvalues('simulated') )
+            print n,len(allobs)
 
         print self.obs.shape
         print len(allobs)

da/tools/pipeline.py

@@ -238,8 +238,6 @@ def SampleOneCycle(DaCycle,Samples,StateVector, ObservationOperator,lag):
 def Invert(DaCycle, StateVector, Optimizer ):
     """ Perform the inverse calculation """
 
-    nobs      = len(StateVector.EnsembleMembers[-1][0].ModelSample.Data.getvalues('obs'))
-
     dims      = ( int(DaCycle['time.nlag']),
                   int(DaCycle['forecast.nmembers']),
                   int(DaCycle.DaSystem['nparameters']),