To the GMSH authors, have you incorporated an option where it's possible to define an internal surface boundary inside a 3D domain (meshed volume)? I noticed that it's possible to have this feature in 2D with 1D lines using Line … in Surface. By the way great program! Gianluca