- d = EnsembleY[:,iens] - Hm * Xb[iens]
- Xa = Xb[iens] + K*d
-
- self.StoredVariables["CurrentState"].store( Xa.A1 )
- self.StoredVariables["Innovation"].store( d.A1 )
+ HXb = numpy.ravel(numpy.dot(Hm, Xb[iens]))
+ if self._toStore("SimulatedObservationAtBackground"):
+ self.StoredVariables["SimulatedObservationAtBackground"].store( HXb )
+ d = numpy.ravel(EnsembleY[:,iens]) - HXb
+ if self._toStore("Innovation"):
+ self.StoredVariables["Innovation"].store( d )
+ Xa = numpy.ravel(Xb[iens]) + numpy.dot(K, d)
+ self.StoredVariables["CurrentState"].store( Xa )
+ if self._toStore("SimulatedObservationAtCurrentState"):
+ self.StoredVariables["SimulatedObservationAtCurrentState"].store( numpy.dot(Hm, Xa) )