s

    % compare concentration at outlet
    clf
    plot(caseData.refSat{1},'*')
    hold all
    plot(caseData.msSat{1})
    plot(caseData.msSat{2})
    legend({'Reference solution','Conforming coarse grid','Non-conforming coarse grid'},'Location','SouthEast')
    xlabel('PVI')
    ylabel('Concentration')
    
    % compare the gmres residuals
    figure
    for i = 1:2
        semilogy(0:caseData.msSol{i}.iter,caseData.msSol{i}.residuals(1:caseData.msSol{i}.iter+1))
        hold all
    end
    legend({'