Converting mesh from XML format to hd5 format

[dokken]

As I said previously, there are plenty of examples of usecases of meshio on the forum, for instance:

Transitioning from mesh.xml to mesh.xdmf, from dolfin-convert to meshio

The problem is that meshio defines multiple cell block for facets. The following works: import meshio import numpy as np msh = meshio.read("juanfe360.msh") line_cells = [] for cell in msh.cells: if cell.type == "triangle": triangle_cells = cell.data elif cell.type == "line": if len(line_cells) == 0: line_cells = cell.data else: line_cells = np.vstack([line_cells, cell.data]) line_data = [] for key in msh.cell_data_dict["gmsh:physical"].…

Note that most of these problems reads the meshes from msh. However, a similar procedure can be done with xml meshes.