%
%
load a01.dat
load a02.dat
load a03.dat
load a04.dat
load a05.dat
load a06.dat
load a07.dat
load a08.dat
load a09.dat
load a10.dat
load a11.dat
load a12.dat
load a13.dat
load a14.dat
load a15.dat
load b01.dat
load b02.dat
load b03.dat
load b04.dat
load b05.dat
load b06.dat
load b07.dat
load b08.dat
load b09.dat
load b10.dat
load b11.dat
load b12.dat
load b13.dat
load b14.dat
load b15.dat
load c01.dat
load c02.dat
load c03.dat
load c04.dat
load c05.dat
load c06.dat
load c07.dat
load c08.dat
load c09.dat
load c10.dat
load c11.dat
load c12.dat
load c13.dat
load c14.dat
load c15.dat
load d01.dat
load d02.dat
load d03.dat
load d04.dat
load d05.dat
load d06.dat
load d07.dat
load d08.dat
load d09.dat
load d10.dat
load d11.dat
load d12.dat
load d13.dat
load d14.dat
load d15.dat
%
%  make the strain profiles for Imperial and Carrizo
%
V=40;
d=6;
x=c01(:,1);
xd2=1.+(x/d).^2;
str=1.e3*V./(pi*d*xd2);
%
%   plot the Imperial and Carrizo profiles
%
orient tall
subplot(2,1,2),plot(x,str,'--',c01(:,1),c01(:,2),':',c02(:,1),c02(:,2),':',c03(:,1),c03(:,2),':',c04(:,1),c04(:,2),':',c05(:,1),c05(:,2),':',c06(:,1),c06(:,2),':',c07(:,1),c07(:,2),c08(:,1),c08(:,2),c09(:,1),c09(:,2),c10(:,1),c10(:,2),c11(:,1),c11(:,2),c12(:,1),c12(:,2),c13(:,1),c13(:,2),c14(:,1),c14(:,2),c15(:,1),c15(:,2));
xlabel('distance (km)')
legend('GPS-2D-model','mccaffrey','tape','becker','shen','parsons','bormann\_hammond','gpsgridder','zeng','hackl','holt','loveless\_meade','kreemer','smith\_konter','tong','bird')
ylabel('strain (nanostrain/yr)')
title('B - B` Imperial')
%
%  make the strain profiles for Imperial and Carrizo
%
V=35;
d=17;
x=d01(:,1);
xd2=1.+(x/d).^2;
str=1.e3*V./(pi*d*xd2);
subplot(2,1,1),plot(x,str,'--',d01(:,1),d01(:,2),':',d02(:,1),d02(:,2),':',d03(:,1),d03(:,2),':',d04(:,1),d04(:,2),':',d05(:,1),d05(:,2),':',d06(:,1),d06(:,2),':',d07(:,1),d07(:,2),d08(:,1),d08(:,2),d09(:,1),d09(:,2),d10(:,1),d10(:,2),d11(:,1),d11(:,2),d12(:,1),d12(:,2),d13(:,1),d13(:,2),d14(:,1),d14(:,2),d15(:,1),d15(:,2));
axis([-80,80,0,1000])
ylabel('strain (nanostrain/yr)')
xlabel('distance (km)')
legend('GPS-2D-model','mccaffrey','tape','becker','shen','parsons','bormann\_hammond','gpsgridder','zeng','hackl','holt','loveless\_meade','kreemer','smith\_konter','tong','bird')
title('A - A` Carrizo')
