);
iceS = reshapetimeseries(t, extent_S, <span class="string">'bin'</span>, 52);

ice = cat(3, iceN, iceS);
</pre><p>By default, the function calculates the mean and upper/lower bounds of the data at each x-location in the data.  We can plot either <a href="boundedline_documentation.html">bounded lines</a>:</p><pre class="codeinput">[A, h] = ensemble2bnd(1:52,ice, <span class="string">'dims'</span>, <span class="string">'xey'</span>, <span class="string">'plot'</span>, <span class="string">'boundedline'</sp