es 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).  

[icem, yr, tmidm] = reshapetimeseries(t, extent_N, 'bin', 'month');

plot(tmidm, icem, 'r')

%%
% While Jan 1 is usually a good place to place the year-to-year wrap,
% sometimes your dataset may call for a different placement.  For example,
% we can place the w