nschloe / meshio

:spider_web: input/output for many mesh formats
MIT License
1.96k stars 402 forks source link

IndexError with msh2vtu #1360

Open emartincb opened 2 years ago

emartincb commented 2 years ago

Hello all,

I was trying to convert .msh file to vtu files with msh2vtu.py but when I run it, the terminal return this error:

--> python3 msh2vtu.py mesh.msh

MeshIO 5.3.4 found, MSH2VTU was tested with MeshIO 4.4.8. Newer version of MeshIO than supported. Backward compatibility may be missing! MeshIO version > 5.0 brings relevant changes. MSH2VTU will catch up as PyVista does so.

Original mesh (read) 1086 points in 3 dimensions; cells: 87 line, 2017 triangle; point_data=['gmsh:dim_tags']; cell_data=['gmsh:physical', 'gmsh:geometrical']; cell_sets=['null flux boundary', 'hydraulic head', 'tailings', 'starter_dam', 'erosion', 'gmsh:bounding_entities']

Trying to save original mesh as vtu-file (possibly not all features may be saved) Info: VTU format cannot write cell_sets. Converting them to cell_data... Traceback (most recent call last): File "/home/eloy/Documentos/Merriespruit_mesh/msh2vtu.py", line 210, in meshio.write(output_basename + "_original.vtu", mesh, binary=not args.ascii) File "/home/eloy/.local/lib/python3.10/site-packages/meshio/_helpers.py", line 188, in write return writer(filename, mesh, **kwargs) File "/home/eloy/.local/lib/python3.10/site-packages/meshio/vtu/_vtu.py", line 644, in write mesh.cell_sets_to_data(key) File "/home/eloy/.local/lib/python3.10/site-packages/meshio/_mesh.py", line 328, in cell_sets_to_data arr[cc] = i IndexError: index 14 is out of bounds for axis 0 with size 13

Is it a geometry issue?

Thanks a lot.

Regards, MC

FrankFrank9 commented 1 year ago

I'm having the exact same issue, same error message. Have you solved? Any idea? Actualy I was using meshio convert test.msh test.vtu