function intadd(kf,kc)

fprintf('interpolate and add correction on level %d to solution on level %d \n',...
    kc, kf)
