 mug * (rhog / rhol);
        <span class="keyword">end</span>

        <span class="comment">%------------------------------------------------------------------</span>
        <span class="comment">% Compute the compressibility coefficients</span>
        <span class="keyword">function</span> [cll, clg, cgl, cgg] = compressibilityCoeffs(this,Sl,pl,pg)
            <span class="comment">% Get porous media parameters</span>
            biot = this.porousMedia.biot;
            phi  = this.porousMedia.phi;
   