%
%
%
load a1.dat
load a2.dat
load a3.dat
load a4.dat
load a5.dat
load a6.dat
load a7.dat
load a8.dat
load a9.dat
load a10.dat
load b1.dat
load b2.dat
load b3.dat
load b4.dat
load b5.dat
load b6.dat
load b7.dat
load b8.dat
load b9.dat
load b10.dat
load c1.dat
load c2.dat
load c3.dat
load c4.dat
load c5.dat
load c6.dat
load c7.dat
load c8.dat
load c9.dat
load c10.dat
%
%   plot the three profiles
%
orient tall
subplot(3,1,1),plot(a1(:,3),a1(:,4),'k',a2(:,3),a2(:,4),'r',a3(:,3),a3(:,4),'c',a4(:,3),a4(:,4),'b',a5(:,3),a5(:,4),'--k',a6(:,3),a6(:,4),'--r',a7(:,3),a7(:,4),'--c',a8(:,3),a8(:,4),'--b',a9(:,3),a9(:,4),':k',a10(:,3),a10(:,4),'g');
xlabel('distance (km)')
legend('becker','fialko','freed','harvard','holt','kreemer','mccaffrey','pollitz','smith-konter','zheng')
ylabel('strain (nanostrain/yr)')
title('A - A`')
subplot(3,1,2),plot(b1(:,3),b1(:,4),'k',b2(:,3),b2(:,4),'r',b3(:,3),b3(:,4),'c',b4(:,3),b4(:,4),'b',b5(:,3),b5(:,4),'--k',b6(:,3),b6(:,4),'--r',b7(:,3),b7(:,4),'--c',b8(:,3),b8(:,4),'--b',b9(:,3),b9(:,4),':k',b10(:,3),b10(:,4),'g');
xlabel('distance (km)')
legend('becker','fialko','freed','harvard','holt','kreemer','mccaffrey','pollitz','smith-konter','zheng')
ylabel('strain (nanostrain/yr)')
title('B - B`')
subplot(3,1,3),plot(c1(:,3),c1(:,4),'k',c2(:,3),c2(:,4),'r',c3(:,3),c3(:,4),'c',c4(:,3),c4(:,4),'b',c5(:,3),c5(:,4),'--k',c6(:,3),c6(:,4),'--r',c7(:,3),c7(:,4),'--c',c8(:,3),c8(:,4),'--b',c9(:,3),c9(:,4),':k',c10(:,3),c10(:,4),'g');
xlabel('distance (km)')
legend('becker','fialko','freed','harvard','holt','kreemer','mccaffrey','pollitz','smith-konter','zheng')
ylabel('strain (nanostrain/yr)')
title('C - C`')
