eSaturationDegree(pc, this.porousMedia);
            % Derivative of the gas density wrt to pg
            drhogdpg = this.derivativeGasDensityWrtGasPressure(rhog,pg);
            % Compressibility coefficients
            ccc =  phi * dSldpc;
            ccg = 0.0;
            cgc = -phi * dSldpc * (rhog / rhol);
            cgg =  phi * (1.0 - Sl) * drhogdpg / rhol;
        end
        %REPLACE_WITH_DASH_DASHREPLACE_WITH_DASH_DASHREPLACE_WITH_DASH_DASHREPLACE_WITH_DASH_DASHREPLACE_WITH_DASH_DASHREPLACE_WI