I have a problem with a minus sign in the integration formula, but I don't know where it comes from.
Consider the domain $Q=\{(x,y)\in \mathbb R^2|x^2+y^2<1,y<0\}$, ie the lower half disk. Consider as well a $C^\infty$ function $g$ such that $g\equiv 0$ inside a neighborhood of the circumference $S=\{(x,y)|x^2+y^2=1\}$.
I want to compute the integral $\int_Q g_{y}$. With the usual orientation of $\mathbb R^2$ and using Fubini we have:
$\int_Q g_{y}=\int_{-1}^{1}\int_{-\sqrt{1-x^2}}^0g_{y}(x,y)dydx=\int_{-1}^1g(x,0)dx$
However we can also compute this integral integrating by parts and we get: $\int_Q g_{y}=-\int_Q 0\cdot g+\int_{\partial Q} 1 g\nu_2$
where $\nu_2$ is an outer normal to the boundary. On the boundary, $g$ is zero outside the axis so we only integrate on it, there the outer normal is $\nu=(0,1)$ so $\nu_2=1$ and the orientation of the border is from right to left since it has the induced orientation. So we get:
$\int_Q g_{y}=\int_{1}^{-1}g(x,0)dx$
Which has clearly the opposite sign. Where is the problem?