def METHOD_NAME(ensemble):
    mean_ensemble = np.zeros(len(ensemble[0]))
    for iEns in range(len(ensemble)):
        mean_ensemble += ensemble[iEns]
    # Compute mean of ensemble
    mean_ensemble /= float(len(ensemble))
    return mean_ensemble