search

nschloe / meshio

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

Storing a scalar property of the whole grid in a VTK file #1446

Open maneesh29s opened 6 months ago

maneesh29s commented 6 months ago

I want to store 100x100 grid cells in a VTK file. Along with the cells and points, I want to store 2 scalar values, X = 100 Y = 100, which represent the X and Y dimensions of the grid

The above values can not be stored as point_data or cell_data as X and Y are properties of the whole grid.

https://examples.vtk.org/site/Cxx/PolyData/FieldData/ The above link mentions that "For data that does not align with either points or cells, FieldData should be used" The Mesh class has a field_data attribute

But the vtk_sections list present in __vtk51.py does not have any section related to field data.

So what is the best way to store this?

As a workaround, I tried adding a FIELD_DATA section in the __vtk51.py and made some changes to the write() and _read_section() methods. This approach solves my issue, but I am not sure if this is valid from the domain point of view.

chrisflesher commented 5 months ago

I am interested in something similar, storing a single scalar metadata value in a PLY file.