function [kf, kn] = divide_k(k,qfit)
kf = k(logical(qfit));
kn=k(~logical(qfit));
