Why the integral computed on face or on 2D grid gives different results?

Yann_Le_Floch

Such a difference could be explained by the difference of the discretization on 2D grid and on face. In Flux, 2D grids have their own discretization (number of discretization elements) imposed by the user which could be different from the mesh discretization on the face.