    <span class="comment">%------------------------------------------------------------------</span>
        <span class="comment">% Compute the compressibility coefficients</span>
        <span class="keyword">function</span> [ccc, ccg, cgc,cgg] = compressibilityCoeffsPgPc(this,Sl,pl,pg)
            <span class="comment">% Get porous media parameters</span>
            phi  = this.porousMedia.phi;
            <span class="comment">% Get fluid densities</span>
            rhol = this.liquidFluid.getDensity(