subplot(141)
l=length(c(:,1))+3;
z(2:l-2)=c(:,1);
c1(2:l-2)=c(:,2);
c4(2:l-2)=c(:,2)+c(:,5);
c2(2:l-2)=c(:,2)+c(:,5)+c(:,3);
c5(2:l-2)=c(:,2)+c(:,5)+c(:,3)+c(:,6);
c3(2:l-2)=c(:,2)+c(:,5)+c(:,3)+c(:,6)+c(:,4);
c6(2:l-2)=c(:,2)+c(:,5)+c(:,3)+c(:,6)+c(:,4)+c(:,7);

zmn=max(z);
zmx=min(z);
z(1) = zmn; z(l-1)=zmx; z(l)=zmn;
c1(1)=0;c1(l-1)=0;c1(l)=0;
c4(1)=0;c4(l-1)=0;c4(l)=0;
c3(1)=0;c3(l-1)=0;c3(l)=0;
c2(1)=0;c2(l-1)=0;c2(l)=0;
c5(1)=0;c5(l-1)=0;c5(l)=0;
c6(1)=0;c6(l-1)=0;c6(l)=0;

cm1=[255 202 153]/255;
cm2=[204 254 255]/255;
cm4=[0 0 0];
cm3=[0 0 0.5625];
cm5=[43 184 239]/255;
cm6=[161 204 41]/255;

fill(c6,z,cm6),hold on
fill(c3,z,cm3)
fill(c5,z,cm5)
fill(c2,z,cm2)
fill(c4,z,cm4)
fill(c1,z,cm1)
axis([0 1 -1400 0])
title('U1309D')
ylabel('depth (mbsf), 20 m running average, 5 m step')
xlabel('fraction of cored interval')
