The integral of Flux density (dFluxTh) should be computed on face regions in 2D and on volume regions in 3D.