IRP=36 Faust for two-phase Buckley-Leverett problem

IRP = 36 Functions used by Faust (1985) for two-phase Buckley-Leverett problem

$k_{rl}$ = ( $S_l$ - 0.16)2/0.64

$k_{rg}$ = 0

$k_{rn}$ = (0.8 - $S_l$ )2/0.64

where $k_{rl}$ and $k_{rn}$ are limited to values between 0 and 1, no parameters.