<div dir="ltr"><div>Hi all,</div><div><br></div><div><div class="gmail-post-text">

<p>I'm having some troubles setting the ID for the physical group of the boundary and then extract them later in Python. Here's my <code>main.geo</code> file:</p>
<pre><code>// Gmsh project created on Sat Jun 20 11:28:19 2020
Mesh.Algorithm = 6;
Mesh.CharacteristicLengthMin = 0.4;
Mesh.CharacteristicLengthMax = 0.4;

//+
Box(1) = {-2, -2, -2, 4, 4, 4};
// Surface for Box
box_bnd[] = Boundary{Volume{1};};
Physical Surface("boundary", 1) = box_bnd[];
//+
Sphere(3) = {0, 0, 0, 1.0, -Pi/2, Pi/2, 2*Pi};
sphere_bnd[] = Boundary{Volume{3};};
Physical Surface("boundary", 1) += sphere_bnd[];
// Cut
BooleanDifference(5) = {Volume{1}; Delete;}{Volume{3}; Delete;};
Physical Volume("internal", 2) = {5};
</code></pre>
<p>Having run this inside Gmsh and generate <code>main.msh</code>, I try to explore this mesh inside Python and extract out the coordinates on the boundary for checking. Here is my <code>extract.py</code>:</p>
<pre><code>import matplotlib.pyplot as plt
import gmsh

# Mesh characteristic length
char_length = 0.4
# Some reference dimension
max_coord, min_coord = -2.0, 2.0
a = 1.0
# Tolerance allowed in the mesh point
tol = 1e-12

gmsh.initialize()
gmsh.option.setNumber("General.Terminal", 1)
gmsh.open("main.msh".format(char_length))

# all mesh node coordinates
nodeTags, nodeCoords, _ = gmsh.model.mesh.getNodes()

# just the boundary coordinates
bTags, bCoords, = gmsh.model.mesh.getNodesForPhysicalGroup(3, 1)

gmsh.finalize()
</code></pre>
<p>which seems to give me empty arrays of <code>bTags</code> and <code>bCoords</code>, i.e.:</p>
<pre><code>In : run extract.py
Warning : Gmsh has aleady been initialized
Info    : 33 entities
Info    : 1260 nodes
Info    : 6524 elements