You cant expect the “+” restriction to be consistent, see: Integrating over an interior surface - #4 by MiroK
I guess you want to integrate over every cell, and then for each facet of the cell compute the dS integrals. This would be that you would integrate over each facet twice (one from each side). Is that what you want?