           plotCellData_DFM(g_fine, refSol.pressure); axis tight; colorbar;
            title('Pressure Ref')
            cbar = caxis();
            subplot(2,2,4)
            hs = plotCellData_DFM(g_fine, msSol.pressure); axis tight; colorbar;
            title('Pressure MSFV')
            caxis(cbar)
            
            
            drawnow
        end
        t = t + dT;

    end

   
     caseData.refSat{casenr} = refSat;
     caseData.msSat{casenr} = msSat;
      
end

%% Plot results
if doPlotRe