Integration measures for coupled bulk-surface equations

Here you are trying to do a mixed variational form with a codim-1 mesh, which is currently under development in the following branches:

With these two PRs you should be able to create a mixed variational form with variables living in cells, and variables that only live in a function space on a subset of facets.