uidFluid.getDensity(pl);
            rhog = this.gasFluid.getDensity(pg);
            % Compute relative permeability coefficients
            klr = this.liqRelativePermeability.calculate(Sl, this.porousMedia);
            kgr = this.gasRelativePermeability.calculate(Sl, this.porousMedia);
            % Permeability matrices
            Kll = - K * klr / mul;
            Klg = K * klr / mul;
            Kgl = zeros(2);
            Kgg = K * kgr / mug * (rhog / rhol);
        end

        %REPLACE_WITH_DASH_