%% 

figure
    [bbb,bbb1] = maxflat(20,5,1/1.7/128); 
load data_BARCAN/run85
ix=find(xwg==18.5); 
etaLF = filtfilt(bbb,bbb1,etawg(:,ix));  
ax(1)=subplot(511); 
plot(twg,etawg(:,ix),'color',[1 0 0]), hold on 
plot(twg,etaLF,'color',[1 0 0],'linewidth',2) 
ylabel('elevation (m)'), legend('C1','location','northwest')
load data_BARCAN/run100
ix=find(xwg==18.5); 
etaLF = filtfilt(bbb,bbb1,etawg(:,ix));  
ax(2)=subplot(512); 
plot(twg,etawg(:,ix),'color',[0 0 1]), hold on 
plot(twg,etaLF,'color',[0 0 1],'linewidth',2)  
ylabel('elevation (m)'), legend('C2','location','northwest')
load data_BARCAN/run128
ix=find(xwg==18.5); 
etaLF = filtfilt(bbb,bbb1,etawg(:,ix));  
ax(3)=subplot(513); 
plot(twg,etawg(:,ix),'color',[0 0.7 0]), hold on 
plot(twg,etaLF,'color',[0 0.7 0],'linewidth',2)  
ylabel('elevation (m)'), legend('C3','location','northwest')
load data_BARCAN/run111
ix=find(xwg==18.5); 
etaLF = filtfilt(bbb,bbb1,etawg(:,ix));  
ax(4)=subplot(514); 
plot(twg,etawg(:,ix),'color',[1 0. 1]), hold on 
plot(twg,etaLF,'color',[1 0. 1],'linewidth',2)  
ylabel('elevation (m)'), legend('C0','location','northwest')
load data_BARCAN/run149
ix=find(xwg==20); 
etaLF = filtfilt(bbb,bbb1,etawg(:,ix));  
ax(5)=subplot(515); 
plot(twg,etawg(:,ix),'color',[1 0.7 0]), hold on 
plot(twg,etaLF,'color',[1 0.7 0],'linewidth',2)  
ylabel('elevation (m)'), legend('C5','location','northwest')
xlabel('time (s)') 
linkaxes(ax,'xy'), xlim([0 50]), ylim([-.06 .08])

print Fig_A_0 -dpdf