save(outfile,'recon','-ascii')

% make smooth versions too

mt=size(recon,2);
recon_sm=recon;

%recon_sm(:,2:mt)=lowpassmin(recon(:,2:mt),0.05);  % 20 yr smooth
%save(strcat(outfile,'_20yrsmooth'),'recon_sm','-ascii')

% for comparisons only
recon_sm(:,2:mt)=lowpassmin(recon(:,2:mt),0.1);  % 10 yr smooth
save(strcat(outfile,'_10yrsmooth'),'recon_sm','-ascii')
