Get tags from gmsh
I made changes to the GmshMeshReceiver class in meshmode/mesh/io.py
Basically, I recorded a map from nodes to the elements that contain them, then populate boundary tags from "lower-dimensional" elements (which get thrown away) up to the higher dimensional elements that contain them.
Only physical tags of boundary dimension (i.e. dimension of the mesh - 1) are marked
Edited by Ben Sepanski