The bug was the following: the reordering transport solvers accounts for sources and sinks, and modifies saturation accordingly. This saturation is the one we must use for computing sink outflow (= fracflow(s)*outflux). If we wait until after segregation, saturation is changed, and we compute fracflow(s')*outflux instead.