No commit message

No commit message

da/baseclasses/optimizer.py

@@ -77,7 +77,10 @@ class Optimizer(object):
         # species type
   #      self.species       = np.zeros(self.nobs,str)
         self.species       = np.zeros(self.nobs,dtype='<U6')
-	
+
+        # species mask
+        self.speciesmask    = {}	
+
         # Total covariance of fluxes and obs in units of obs [H P H^t + R]
         self.HPHR            =  np.zeros( (self.nobs,self.nobs,), float)
         # Kalman Gain matrix
@@ -124,9 +127,12 @@ class Optimizer(object):
 
                 simulatedensemble =                               Samples.Data.getvalues('simulated')
 
+                print 'simulatedensemble',np.array(simulatedensemble).shape
+
                 if allsimulated == None :
                     allsimulated = np.array(simulatedensemble)
                 else:
+                    print 'simulatedensemble else',np.array(simulatedensemble).shape,allsimulated.shape
                     allsimulated = np.concatenate((allsimulated,np.array(simulatedensemble)),axis=0)
 
 

da/baseclasses/statevector.py

@@ -64,7 +64,6 @@ class EnsembleMember(object):
 
         self.membernumber           = membernumber   # the member number
         self.ParameterValues        = None           # Parameter values of this member
-        self.ModelSample            = None           # Model Sampled Parameter values of this member
 
     def __str__(self):
         return "%03d"%self.membernumber