fluid-matrix interactions: add a tabulated two-phase pc/kr law which uses splines

compared to the variant which uses linear interpolation, this makes
the resulting curves smoother and thus hopefully improves the
convergence rates of the simulations.