Finding facet index in imported mesh

You can use distribute_entity_data

as done in dolfinx/gmshio.py at 513d933eee37062b7847f51c48cbe1a6ee4083da · FEniCS/dolfinx · GitHub