%%%%%%%%%%%%%%%%%%%%%% Plastic/magmatic foliations
figure(fignum)
clf

subplot(1,4,1)
plot(Dmagm(:,2),Dmagm(:,1),'r-')
hold on
plot(Dplas(:,2),Dplas(:,1),'b-')
axis('ij');axis([0 2 0 1420])
xlabel('Intensity')
title('Plastic & magmatic foliation')
ylabel ('Depth, mbsf')

subplot(1,4,2)
load fol_mag.dat;fm=fol_mag;
load fol_plas.dat;fp=fol_plas;
clear  fol_mag fol_plas
plot(fm(:,4),fm(:,1),'ro')
hold on
plot(fp(:,4),fp(:,1),'pb')
%%%% 304

% for shift correction
load U1309D_shift.dat; s= U1309D_shift;clear U1309D_shift
load U1309D_cores.dat; d=U1309D_cores;clear U1309D_cores

load fol_mag304.dat;fmb=fol_mag304;clear fol_mag304
nz=dodepth(fmb(:,[1 2 3]),d,s);nz=nz(:,4);
plot(fmb(:,4),nz,'ro')
load fol_plas304.dat;fpb=fol_plas304;
nz=dodepth(fpb(:,[1 2 3]),d,s);nz=nz(:,4);
plot(fpb(:,4),nz,'bp')

axis('ij');axis([0 90 0 1420])
xlabel('Dip')
title('Foliation dip')


%%%%%%%%% Layering
load layering_comp.dat;lc=layering_comp;clear layering_comp;
load layering_grain.dat;lg=layering_grain;clear layering_grain;
load fol_ox.dat;fo=fol_ox;clear fol_ox;
subplot(1,4,3)
plot(lc(:,3),lc(:,1),'ro')
hold on
plot(lg(:,3),lg(:,1),'bp')
plot(fo(:,4),fo(:,1),'k*')

xlabel('Dip')
title('Layering')
axis('ij');axis([0 90 0 1420])

%%%%%%%%%%%% other alteration magmatic vein igneous contact
load dip_altfront.dat; daf=dip_altfront;clear dip_altfront;
load dip_igcont.dat; dic=dip_igcont;clear dip_igcont;
load dip_magvein.dat; dmv=dip_magvein;clear dip_magvein;

subplot(1,4,4)
plot(daf(:,3),daf(:,1),'bo')
hold on
plot(dic(:,3),dic(:,1),'kx')
plot(dmv(:,3),dmv(:,1),'r*')
xlabel('Dip')
title('Veins, contacts & alteration')
axis('ij');axis([0 90 0 1420])