%% 

figure
    [bbb,bbb1] = maxflat(20,5,1/1.7/128); 
    [ccc,ccc1] = maxflat(20,5,1/1.7/50); 
load data_BARCAN/run92
etawgLF = filtfilt(bbb,bbb1,etawgG);  
etanumG(:,xnum>25.3)=NaN;
etanumLF = filtfilt(ccc,ccc1,etanumG);  
ax(1)=subplot(211); 
plot(xwg,400*std(etawgG-etawgLF),'^','color',[1 0 0],'markersize',3), hold on 
plot(xwg,400*std(etawgLF),'^','color',[1 0 0],'markersize',6)
plot(xnum,400*std(etanumG-etanumLF),'-','color',[1 0 0])
plot(xnum,400*std(etanumLF),'--','color',[1 0 0])
ylabel('wave height (m)'), xlabel('cross-shore distance (m)')
ax(2)=subplot(212); 
plot(xwg,skewness(etawgG),'^','color',[1 0 0],'markerfacecolor',[1 0 0],'markersize',3), hold on 
plot(xnum,skewness(etanumG),'-','color',[1 0 0])
ylabel('wave skewness'), xlabel('cross-shore distance (m)')
linkaxes(ax,'x'), xlim([5 25.5])

load data_BARCAN/run104
etawgLF = filtfilt(bbb,bbb1,etawgG);  
etanumG(:,xnum>25.2)=NaN;
etanumLF = filtfilt(ccc,ccc1,etanumG);  
ax(1)=subplot(211); 
plot(xwg,400*std(etawgG-etawgLF),'o','color',[0 0 1],'markersize',3), hold on 
plot(xwg,400*std(etawgLF),'o','color',[0 0 1],'markersize',6)
plot(xnum,400*std(etanumG-etanumLF),'-','color',[0 0 1])
plot(xnum,400*std(etanumLF),'--','color',[0 0 1])
ax(2)=subplot(212); 
plot(xwg,skewness(etawgG),'o','color',[0 0 1],'markerfacecolor',[0 0 1],'markersize',3), hold on 
plot(xnum,skewness(etanumG),'-','color',[0 0 1])

load data_BARCAN/run146
etawgLF = filtfilt(bbb,bbb1,etawgG);  
etanumG(:,xnum>25.0)=NaN;
etanumLF = filtfilt(ccc,ccc1,etanumG);  
ax(1)=subplot(211); 
plot(xwg,400*std(etawgG-etawgLF),'s','color',[0 0.7 0],'markersize',3), hold on 
plot(xwg,400*std(etawgLF),'s','color',[0 0.7 0],'markersize',6)
plot(xnum,400*std(etanumG-etanumLF),'-','color',[0 0.7 0])
plot(xnum,400*std(etanumLF),'--','color',[0 0.7 0])
ax(2)=subplot(212); 
plot(xwg,skewness(etawgG),'s','color',[0 0.7 0],'markerfacecolor',[0 0.7 0],'markersize',3), hold on 
plot(xnum,skewness(etanumG),'-','color',[0 0.7 0])

load data_BARCAN/run149
etawgLF = filtfilt(bbb,bbb1,etawgG);  
etanumG(:,xnum>25.0)=NaN;
etanumLF = filtfilt(ccc,ccc1,etanumG);  
ax(1)=subplot(211); 
plot(xwg,400*std(etawgG-etawgLF),'d','color',[1 0.7 0],'markersize',3), hold on 
plot(xwg,400*std(etawgLF),'d','color',[1 0.7 0],'markersize',6)
plot(xnum,400*std(etanumG-etanumLF),'-','color',[1 0.7 0])
plot(xnum,400*std(etanumLF),'--','color',[1 0.7 0])
ax(2)=subplot(212); 
plot(xwg,skewness(etawgG),'d','color',[1 0.7 0],'markerfacecolor',[1 0.7 0],'markersize',3), hold on 
plot(xnum,skewness(etanumG),'-','color',[1 0.7 0])

ylim([-.6 1.6])
linkaxes(ax,'x'), xlim([15 25.5])

print Fig_A_00 -dpdf