Consider: https://jorgensd.github.io/dolfinx-tutorial/chapter3/subdomains.html?highlight=meshio#subdomains-defined-from-external-mesh-data