The H-exchange BC is Neumann by default. We need to add a Dirichlet for M and generalize this for the boundaries.