next, and averages Feb 28-29 data in leap years.  Monthly binning can also be achieved by setting the 'bin' input appropriately. (Note that the first year in this dataset started mid-October, which leads to the outlier-like point in the monthly averages).</p><pre class="codeinput">[icem, yr, tmidm] = reshapetimeseries(t, extent_N, <span class="string">'bin'</span>, <span class="string">'month'</span>);

plot(tmidm, icem, <span class="string">'r'</span>)
</pre><img vspace="5" hspace="5" src="reshapetimeserie