with virtual observations
n,N = ubi.shape # n is the state dimension and N is the size of ensemble
m = w.shape[0] # m is the size of measurement vector