ant to plot saturation</span>
        <span class="keyword">if</span> doPlotSim
            clf;
            <span class="comment">% Saturation plot</span>
            subplot(2,2,1)
            plotGrid_DFM(g_fine, <span class="string">'FaceColor'</span>, <span class="string">'None'</span>, <span class="string">'EdgeAlpha'</span>, 0)
            plotCellData_DFM(g_fine, refSol.s(:,1), refSol.s(:,1) &gt; 1e-4);
            plotFractures(g_fine,find(refSol.s(:,1) &gt; 1e-4),refSol.s(:,1)); axis <span class="