%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Open and drilling
figure(fignum)
clf

% recovery
subplot(1,3,1)
ig=find(finite(Dcata(:,4)));igf=ig(length(ig));
fill([0;Dcata(ig,4);0],[Dcata(1,1);Dcata(ig,1);Dcata(igf,1)],[0 0 1])
hold on
plot(Dcata(:,4),Dcata(:,1),'b-')
axis('ij'); axis([0 100 0 1420])
title('Recovery');xlabel('%');ylabel('Depth, mbsf');

% OPEN
subplot(1,3,2)
for i=1:length(doi(:,1)),
    plot([0 doi(i,[3 3]) 0],doi(i,[1 1 2 2]),'r-');
    hold on
end
plot(Dopen(:,2)+Dopen(:,3)/2,Dopen(:,1),'g-')
plot(Dopen(:,2)-Dopen(:,3)/2,Dopen(:,1),'g-')
plot(Dopen(:,2),Dopen(:,1),'b-')
axis('ij'); axis([0 3 0 1420])
title('Open fractures');xlabel('Intensity');


% DRILL
subplot(1,3,3)
for i=1:length(ddi(:,1)),
    plot([0 ddi(i,[3 3]) 0],ddi(i,[1 1 2 2]),'r-');
    hold on
end
plot(Ddril(:,2)+Ddril(:,3)/2,Ddril(:,1),'g-')
plot(Ddril(:,2)-Ddril(:,3)/2,Ddril(:,1),'g-')
plot(Ddril(:,2),Ddril(:,1),'b-')
axis('ij'); axis([0 3 0 1420])
title('Drilling fractures');xlabel('Intensity')
