p_c p_nw-p_w Capillary pressure
Hc p_c/(rho_water*g_const) Capillary pressure head
Se_w (1+(alpha*rm1(Hc[1/m])[m])^N)^(-M) Effective saturation, wetting phase
Se_nw 1-Se_w Effective saturation, non-wetting phase
theta_w theta_r+Se_w*(theta_s-theta_r) Volume fraction, wetting phase
theta_nw theta_s-theta_w Volume fraction, non-wetting phase
kr_w (Se_w^L*(1-(1-Se_w^(1/M))^M)^2)+eps Relative permeability, wetting phase
kr_nw ((1-Se_w)^L*(1-Se_w^(1/M))^(2*M))+eps Relative permeability, non-wetting phase
Cp ((alpha*M/(1-M)*(theta_s-theta_r)*Se_w^(1/M)*(1-Se_w^(1/M))^M))/(rho_water*g_const) Specific capacity