function h = krrg_gw(sw, SW, KRO, varargin)
n = numel(varargin);%表示没有输入裂缝参数
if n == 0
    h = interptable(SW, KRO, sw);
else
    SWF = varargin{1};
    KROF = varargin{2};
    rpt = varargin{3};
    hm = interptable(SW, KRO, sw);
    hf = interptable(SWF, KROF, sw);
    h = rpt .* hm + ~rpt .* hf;%分别为裂缝与基质
end
% h=sw+(1-sw);
end