compute the mean of forecast ensemble
ub = np.mean(ubi,1)