       title(<span class="string">'Saturation inexact'</span>)
            <span class="comment">% Align colorbars</span>
            caxis([0 1])


            <span class="comment">% Pressure plot</span>
            subplot(2,2,3)
            plotCellData_DFM(g_fine, refSol.pressure); axis <span class="string">tight</span>; colorbar;
            title(<span class="string">'Pressure Ref'</span>)
            cbar = caxis();
            subplot(2,2,4)
            hs = plotCellData_DFM(g_fine, msSol.pressure)