Hi. I'm trying to read a mesh of higher element order, in this example a mesh consisting of 10-node tetrahedral elements, from gmsh, into PETSC. But It looks like the mesh is not properly being loaded and converted into a DMPlex. gmsh tells me it has generated a mesh with 7087 nodes, but when I view my dm object it tells me it has 1081 0-cells. This is the printout I get ... Info : Done meshing order 2 (Wall 0.0169823s, CPU 0.016662s) Info : 7087 nodes 5838 elements ... DM Object: DM_0x84000000_0 1 MPI process type: plex DM_0x84000000_0 in 3 dimensions: Number of 0-cells per rank: 1081 Number of 1-cells per rank: 6006 Number of 2-cells per rank: 9104 Number of 3-cells per rank: 4178 Labels: celltype: 4 strata with value/size (0 (1081), 6 (4178), 3 (9104), 1 (6006)) depth: 4 strata with value/size (0 (1081), 1 (6006), 2 (9104), 3 (4178)) Cell Sets: 1 strata with value/size (2 (4178)) Face Sets: 6 strata with value/size (12 (190), 21 (242), 20 (242), 11 (192), 22 (242), 10 (188)) Field P2: adjacency FEM ... To replicate the error try generating a mesh according to https://gmsh.info/doc/texinfo/gmsh.html#t5 <https://gmsh.info/doc/texinfo/gmsh.html#t5> setting the element order to 2, and then loading the mesh using DMPlexCreateGmshFromFile I don't have any issues when i set the element order to 1. Thanks in advance, Vilmer