ensemble size
N = 10
#initialize ensemble
uai = np.zeros([3,N])
for i in range(N):
uai[:,i] = u0b + np.random.multivariate_normal(np.zeros(n), B)
km = 0
for k in range(nt):