        msSat = [msSat ; msSol.s(c2,1)];
        refSat = [refSat ; refSol.s(c2,1)];
        
        % Increase time and continue if we do not want to plot saturation
        if doPlotSim
            clf;
            % Saturation plot
            subplot(2,2,1)
            plotGrid_DFM(g_fine, 'FaceColor', 'None', 'EdgeAlpha', 0)
            plotCellData_DFM(g_fine, refSol.s(:,1), refSol.s(:,1) > 1e-4);
            plotFractures(g_fine,find(refSol.s(:,1) > 1e-4),refSol.s(:,1)); axis tight; colorbar;
      