[karam_j]

Hello, I am implementing a diffusion term using Gauss Theorem and adding all the face fluxes in the cell to get a source term. This is not working properly as i have to know when to add the flux and when to subtract based on the face normal direction, which is sometimes pointing in a direction and sometimes on another one for internal faces.

Is there a clean way to implement this without using F_C0 and F_C1, which are pretty confusing for me given the assignment of C0 and C1 does not seem consistent throughout the domain.