
% Read in the QBO time series and make a plot.
%
% Todd Mitchell, October 2005

a = load( 'qbo19532001.asci' );
whos a
[ yr1, yr2 ] = range( a(:,1) )
nyr = yr2 - yr1 + 1
nt  = nyr * 12;
ts = a(:,3);

clf
fillbar3( ts, yr1, yr2 );

ts2 = running_ave( ts, 6 );
ts2 = ts2(8:12:nt);
clf
fillbar3( ts2, yr1, yr2, 1, 1950 );

b = time_shift( ts2, yr1, yr2, 1900, 2004, 1 );
save qbojjason.dat b -ascii
