How to specify Boundary conditions when the mesh is imported from GMSH?

Please read through my post carefully.
What I stated is:

  1. If you have marked your boundary in gmsh, you can read in this data to a mesh function as described in: Transitioning from mesh.xml to mesh.xdmf, from dolfin-convert to meshio - #3 by dokken
  2. You can now use your mesh function (mf In the link) in boundary conditions such as integral measures ds_C = Measure(“ds”, domain=mesh, subdomain_marker=mf) or in Dirichlet conditions bc = DirichletBC(V, u, mf, 1)
1 Like