load amo_indices

%
% remove trailing NaNs
%
amo_smooth_date = amo_date;

II = find(isnan(amo));
amo(II) = [];
amo_date(II) = [];

amo_smooth_date = amo_date;

II = find(isnan(amo_smooth));
amo_smooth(II) = [];


[amo_butter_ann, amo_butter_bia, amo_butter_mos, amo_lp_ann, amo_lp_mos, amo_prolate_ann, amo_prolate_bia, amo_prolate_mos, amo_rmn_ann, amo_rmn_bia, amo_rmn_mos] = make_filtered_monthly_series (amo, amo_date);

[amo_smooth_butter_ann, amo_smooth_butter_bia, amo_smooth_butter_mos, amo_smooth_lp_ann, amo_smooth_lp_mos, amo_smooth_prolate_ann, amo_smooth_prolate_bia, amo_smooth_prolate_mos, amo_smooth_rmn_ann, amo_smooth_rmn_bia, amo_smooth_rmn_mos] = make_filtered_monthly_series (amo_smooth, amo_smooth_date);


save amo_smooth_filtered.mat amo_date amo*

