Please read through my post carefully.
What I stated is:
- 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
- You can now use your mesh function (
mf
In the link) in boundary conditions such as integral measuresds_C = Measure(“ds”, domain=mesh, subdomain_marker=mf)
or in Dirichlet conditionsbc = DirichletBC(V, u, mf, 1)