For non-zero neumann boundary conditions we did a scaling to recover
DgradU at the boundary from the prescribed flux. We don't do this anymore
and assume that when the UseTpfaBoundary property is set to false, that the
user provides -DgradU at the boundary instead of, e.g., mass fluxes.