bility = RelativePermeabilityPolynomialLiquid();
            <span class="keyword">elseif</span> strcmp(<span class="string">'UMAT'</span>,matData.porousMedia.liqRelPermeability)
                curve = matData.porousMedia.klr_umat;
                this.liqRelativePermeability = RelativePermeabilityUMAT(curve(:,1),curve(:,2));
            <span class="keyword">end</span>
            <span class="comment">% Gas phase relative permeability function</span>
            <span class="keyword">if</span> strcmp(<sp