MeshSeries can't read xdmf if the mesh contains quadratic order elements

Seems to be a meshio issue, so we will have to work aorund this.

Happens form quadratic triangles and quadratic quadrilaterals. I would guess also fr 3D quadratic elements but haven't tested.

File ~/ogstools/ogstools/meshlib/mesh_series.py:90, in MeshSeries._read_xdmf(self, timestep)
     89 def _read_xdmf(self, timestep: int) -> pv.UnstructuredGrid:
---> 90     points, cells = self._xdmf_reader.read_points_cells()
     91     _, point_data, cell_data, field_data = self._xdmf_reader.read_data(
     92         timestep
     93     )
     94     meshio_mesh = meshio.Mesh(
     95         points, cells, point_data, cell_data
     96     )

File ~/ogstools/.venv/lib/python3.10/site-packages/meshio/xdmf/time_series.py:116, in TimeSeriesReader.read_points_cells(self)
    113     cell_type = c.attrib["TopologyType"]
    115 if cell_type == "Mixed":
--> 116     cells = translate_mixed_cells(data)
    117 else:
    118     cells.append(CellBlock(xdmf_to_meshio_type[cell_type], data))

File ~/ogstools/.venv/lib/python3.10/site-packages/meshio/xdmf/common.py:143, in translate_mixed_cells(data)
    141         r += 1
    142     r += 1
--> 143     r += xdmf_idx_to_num_nodes[xdmf_type]
    145 types = np.array(types)
    146 offsets = np.array(offsets)

KeyError: 37