j=1;
x=1;

for i=1:(length(spmr)-1)
    
    if spmr(1,i)==0|i==length((spmr)-1)
        j=1;
        if ~exist('SKu','var')
            Ku(1:6,x)=0;
            Kd(1:6,x)=0;
            x=x+1;
        else
        for ik=1:6
        Ku(ik,x)=mean(SKu(ik,:));
        Kd(ik,x)=mean(SKd(ik,:));
        end;
        x=x+1;
        clear SKu;
        clear SKd;
        end;
    end;
    
        
    
    if spmr(1,i)~=0&spmr(1,i+1)~=0
    
        for ij=1:6
            
            Depth=spmr(1,i)- spmr(1,i+1);
            LuS=spmr((ij+1),i);
            LuD=spmr((ij+1),i+1);
            EdS=spmr((ij+13),i);
            EdD=spmr((ij+13),i+1);
            
            SKu(ij,j)=calc_k(LuS,LuD,Depth);
            SKd(ij,j)=calc_k(EdS,EdD,Depth);
            
            
        end;
    end;
    j=j+1;
end;
    