%%%%%%%%%%% figure 6
figure(fignum)
clf
subplot(1,3,2)
load faults_zsdrp.dat
fd=faults_zsdrp;clear faults_zsdrp;
ip=find(fd(:,4)==2);id=find(fd(:,4)==1);
fp=fd(ip,:);fd=fd(id,:);
plot(fd(:,3),fd(:,1),'b.')
hold on

% separate data by plunge intervals
ia=find(fd(:,6)<=30);
ib=find(fd(:,6)>30 & fd(:,6)<=60);
ic=find(fd(:,6)>60);[length(fd(:,1)) length(ia) length(ic) sum(fd(:,6)<45) sum(finite(fd(:,6)))]
plot(fd(ia,3),fd(ia,1),'r*')
plot(fd(ib,3),fd(ib,1),'b>')
plot(fd(ic,3),fd(ic,1),'kx')

fall=[fp;fd];
axis('ij');axis([0 90 0 1420])
title('Green veins')
xlabel('Dip')

subplot(1,3,3)
plot(fp(:,3),fp(:,1),'b.')
hold on
ia=find(fp(:,6)<=30);
ib=find(fp(:,6)>30 & fp(:,6)<=60);
ic=find(fp(:,6)>60);[length(fp(:,1)) length(ia) length(ic) sum(fp(:,6)<45) sum(finite(fp(:,6)))]
plot(fp(ia,3),fp(ia,1),'r*')
plot(fp(ib,3),fp(ib,1),'b>')
plot(fp(ic,3),fp(ic,1),'kx')

% fault zone
load fault_zone.dat;fz=fault_zone;clear fault_zone
plot(fz(:,3),fz(:,1),'ko')
plot(fz(:,3),fz(:,1),'k*')

axis('ij');axis([0 90 0 1420])
title('Pale green veins')
xlabel('Dip')

subplot(1,3,1)
plot(Dcata(:,2),Dcata(:,1),'r-')
hold on
plot(Dvein(:,2),Dvein(:,1),'b-')
for i=1:length(fz(:,1));
    for j=1:2;
        subplot(1,3,j)
        plot([0;5],fz(i,[1 1]),'k:')
    end
end
subplot(1,3,1)

axis('ij');axis([0 4 0 1420])
title('Cataclastic deformation')
xlabel('Intensity')