uses second derivative information to find range for saturation at front for BL solution