ber of x-dimension
% points (i.e. length of x input), ny corresponds to the number of
% different datasets, and nens is the number of ensemble members in each
% dataset. (See 'dims' parameter to permute these dimensions)  
%
% |[A, h] = ensemble2bnd(x,y, 'plot', plottype)| plots the datasets as
% either lines with shaded patches indicating percentiles ('boundedline'),
% lines with errorbars ('errorbar'), unstacked bar plots with errorbars
% ('bar'), or boxplots ('boxplot').  The default, 'none', indicates n