 ensemble2bnd(1:52,ice, 'dims', 'xey', 'plot', 'errorbar')

%%
% or bar plots with errorbars:

cla
[A,h] = ensemble2bnd(1:52, ice, 'dims', 'xey', 'plot', 'bar')
set(h.bar, 'edgecolor', 'none');

%%
% Note that the field names in the handle output structures vary
% based upon the plotting option you choose.

%%
% We can also increase the number of intervals shown.  For example, perhaps
% you'd like to see interquartile ranges as well as the minimum and
% maximum.  In this case, it may make more sense to use 