Establish subsets for the boundary marks.
Cause external facet integrals to occur only over the relevant subset.
Minor bugfix to cause the correct subset index to be found. Introduce test of facet integrals.