%%%%%%%%%%%%% Figure 5 cataclasitc and veins
figure(fignum)
clf

subplot(1,3,1)
for i=1:length(dvi(:,1)),
    plot([0 dvi(i,[3 3]) 0],dvi(i,[1 1 2 2]),'r-');
    hold on
end
hold on
plot(Dvein(:,2)+Dvein(:,3)/2,Dvein(:,1),'g-')
plot(Dvein(:,2)-Dvein(:,3)/2,Dvein(:,1),'g-')
plot(Dvein(:,2),Dvein(:,1),'b-')
axis('ij');axis([0 5 0 1420])
ylabel('Depth, mbsf');xlabel('Intensity')
title('Vein intensity')

subplot(1,3,2)
for i=1:length(dci(:,1)),
    plot([0 dci(i,[3 3]) 0],dci(i,[1 1 2 2]),'r-');
    hold on
end
hold on
plot(Dcata(:,2)+Dcata(:,3)/2,Dcata(:,1),'g-')
plot(Dcata(:,2)-Dcata(:,3)/2,Dcata(:,1),'g-')
plot(Dcata(:,2),Dcata(:,1),'b-')

% fault zones
load fault_zone.dat;fz=fault_zone;clear fault_zone

for i=1:length(fz(:,1));
    for j=1:2;
        subplot(1,3,j)
        plot([0;5],fz(i,[1 1]),'k:')
    end
end

axis('ij');axis([0 5 0 1420])
xlabel('Intensity')
title('Cataclastic intensity')