the median as the
% center line rather than the mean:

cla
ensemble2bnd(1:52,ice, 'dims', 'xey', 'plot', 'boundedline', ...
    'cent', 'median', 'prc', [0 25 75 100]);

%%
% If you add a large number of evenly-spaced percentiles, you can mimic a
% shaded distribution plot.

cla
[~,h] = ensemble2bnd(1:52,ice, 'dims', 'xey', 'plot', 'boundedline', ...
    'cent', 'median', 'prc', [0:2:48 52:2:100], 'tlim', [0.1 0.9]);

%%
% The plotting order of the patch objects moves from outer bounds to inner
% bounds, wh