Real
PorousFlowCapillaryPressureBC::d2CapillaryPressureCurve(Real saturation, unsigned /*qp*/) const
{
  Real seff = effectiveSaturationFromSaturation(saturation);
  return PorousFlowBrooksCorey::d2CapillaryPressure(seff, _pe, _lambda, _pc_max) * _dseff_ds *
         _dseff_ds;
}
Real
PorousFlowCapillaryPressureBC::capillaryPressureCurve(Real saturation) const
{
  Real seff = effectiveSaturationFromSaturation(saturation);
  return PorousFlowBrooksCorey::capillaryPressure(seff, _pe, _lambda, _pc_max);
}