For surfaces, instead of printing "1" as the number of physical tags the surface entity belongs to, loop through the values available in the tesselation data (Tess.DomTessFaceNb) for coincidences with the current entity ID.
If this ID is found in the tesselation data, the surface is part of the domain boundary, and so two physical tags should be printed: the ID itself and the tag of the domain face.
Refers to Issue #791
For surfaces, instead of printing "1" as the number of physical tags the surface entity belongs to, loop through the values available in the tesselation data (Tess.DomTessFaceNb) for coincidences with the current entity ID.
If this ID is found in the tesselation data, the surface is part of the domain boundary, and so two physical tags should be printed: the ID itself and the tag of the domain face.