matthiasm@8: function bc13 = basschroma13(bc)
matthiasm@8: bc(isnan(sum(bc,2)),:) = 1;
matthiasm@8: bc = qnormalise(bc,1,2);
matthiasm@8: bc13 = [bc,(1./12./max(bc,[],2)).^2];
matthiasm@8: bc13 = qnormalise(bc13, inf, 2);