Marking surfaces of a 3D mesh (mshr) to apply b.c

If your are using CG 1 elements the solution is continuous, and you can chose either side. Note that the ("+") and ("-") restriction is arbitrary, and if you want to Get the integration to be on a specific side, a hack has to be made, see Integrating over an interior surface