compute analysis of mean
ua = ub + K @ (w-ObsOp(ub))
xbi = np.zeros([n,N]) #ensemble of forecast anomalies
xai = np.zeros([n,N]) #ensemble of analysis anomalies
for i in range(N):