Identifying 3D domain boundaries

Check out this post as well,

You can also found boundary labels more conveniently with built-in function “labels” as
int[int] labs = labels(Th);

if

1 Like