% reproduce Fig 2a from the WDCA submission:
load BEB14.dat
t=[1995:2009];
subplot(311);boxplot(BEB14(:,2:13)',t);grid on
set(gca,'Fontname','Times')
axis([0 16 26 34]);
set(gca,'YTick',[26:2:34],'YAxisLocation','left');grid on
ylabel({'Nov-Apr';'\delta^{18}O (permil, SMOW)'})
set(gca,'XTick',[0:16])
set(gca,'XTickLabel',{'    ';'95/96';'96/97';'97/98';'98/99';'99/00';'00/01';'01/02';'02/03';'03/04';'04/05';'05/06';'06/07';'07/08';'08/09';'09/10';'     '})
text(0.3,32.5,'(a)')
% add the simulated o18 to this plot:
hold on
plot(1:15,BEB14(:,14),'ko-','MarkerFaceColor','k')
hold off

