So my first guess is that you didn’t prune the z-coordinate from your mesh, see: Converting simple 2D mesh from gmsh to dolfin I would also suggest to use XDMF, and not VTK to save your mesh, see: Transitioning from mesh.xml to mesh.xdmf, from dolfin-convert to meshio for usecases and syntax