h = np.vstack([np.atleast_1d(r).transpose(axes).ravel() for r in steps])