      title('Saturation Ref')
            caxis([0 1]);
            
            subplot(2,2,2)
            plotGrid_DFM(g_fine, 'FaceColor', 'None', 'EdgeAlpha', 0)
            plotCellData_DFM(g_fine, msSol.s(:,1), msSol.s(:,1) > 1e-4);
            plotFractures(g_fine,find(msSol.s(:,1) > 1e-4),msSol.s(:,1)); axis tight; colorbar;
            title('Saturation inexact')
            % Align colorbars
            caxis([0 1])
            
            
            % Pressure plot
            subplot(2,2,3)
 