load MBHsmooths1.txt % http://www.climateaudit.info/wp-content/uploads/2009/11/mbhsmooths1.txt % MBHsmooths1 = [ Year AnnualRecon Instrumental MB98SmoothTrick MBH98SmoothNoTrick MBH99SmoothTrick MBH99SmoothNoTrick ] [B98,A98]=butter(10,2/50); % ’50 year lowpass’ in98=MBHsmooths1(401:981,2); % Annual Recon MBH98, %source: % http://www.nature.com/nature/journal/v430/n6995/extref/FigureData/nhmean.txt , % col 2 Trick=MBHsmooths1(982:996,3); % pad with instrumental (1981..1995) .. % source: % http://www.nature.com/nature/journal/v430/n6995/extref/FigureData/nhmean.txt % col 3 out98t=flipud(filter(B98,A98,flipud([filter(B98,A98,[in98;Trick;zeros(100,1)])]))); %..then smooth % Note that filter.m initializes … Continue reading Code: The Trick