----------------------------</span>
        <span class="comment">% Compute the permeability matrices</span>
        <span class="keyword">function</span> [Kll, Klg, Kgl, Kgg] = permeabilityMtrcsPgPc(this,Sl,pl,pg)
            K = this.porousMedia.intrinsicPermeabilityMatrix();
            <span class="comment">% Get fluids dynamic viscosity</span>
            mul = this.liquidFluid.mu;
            mug = this.gasFluid.mu;
            <span class="comment">% Get fluid densities</span>
            rhol = this