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