Ashwin V. Mohanan
Ashwin V. Mohanan
Support for `open_dataset` with non-box meshes is not complete. See #13 for an attempt at this. Hints, suggestions or contributions welcome.
This [comment](https://github.com/pydata/xarray/issues/4470#issuecomment-700670814) advises loop indices for curvilinear structured meshes.
For most cases this seems to be the reason, but it appears to happen (rarely) for certain runs even with box simulations as seen in #56. In v1.5.0 we will...
@jcanton Let me know if you can do this. You were the original author, afterall :smile:
`open_dataset` works only for cubes or cuboids. Is your mesh like that? This is because I have not found a nice way to extract the 1D coordinates from the 3D...
I could help you more if you can generate a minimal example to reproduce.
See also #13 if you are interested in non-cubical geometries.
We need to check if this was a duplicate of #56 which was fixed
Some source code to reproduce this exists here, but the field files need to be generated https://github.com/eX-Mech/pymech-test-data/blob/gmsh-bug
This feature should perhaps wait for explicit indexing support from xarray https://github.com/pydata/xarray/projects/1.