%%%%%%%%%%%%%%%%%%%%%%%% Serpentine foliations
figure(fignum)
clf

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


subplot(1,3,2)
load U1309D_serpdip.dat
sd=U1309D_serpdip;clear U1309D_serpdip
inan=find(isnan(sd(:,3)));
plot(sd(:,3),sd(:,1),'ro')
hold on
plot(sd(inan,1)*0,sd(inan,1),'b>')
axis('ij');axis([0 90 0 1420])
xlabel('Dip')
title('Foliation dip')
