AMReX/BoxLib Grid Reader

[hgopalan]

Describe the feature you would like to be added.

I use this plugin in Paraview to read the output data from AMR-Wind and ERF codes, being developed for energy applications. There is also a native python reader: https://yt-project.org which can handle the output files. However, I like using pyvista. I have been using it within PyQt5 interface for the last few years for urban climate modeling. I would request the ability to use the amrex data format within PyVista to avoid using paraview.

Links to VTK Documentation, Examples, or Class Definitions.

https://kitware.github.io/paraview-docs/latest/python/paraview.simple.AMReXBoxLibGridReader.html

Pseudocode or Screenshots

Sample code in yt:

ds=yt.load(sys.argv[1])
print("ds:",ds.field_list)
left_edge=ds.domain_left_edge
right_edge=ds.domain_right_edge
domain_dimension=[ds.domain_dimensions[0],ds.domain_dimensions[1],ds.domain_dimensions[2]]
nsize=ds.domain_dimensions[0]
vsize=ds.domain_dimensions[2]
all_data_level_0 = ds.arbitrary_grid(left_edge=left_edge,right_edge=right_edge,dims=domain_dimension)
velocityx = all_data_level_0["boxlib", "x_velocity"]
velocityy = all_data_level_0["boxlib", "y_velocity"]
velocityz = all_data_level_0["boxlib", "z_velocity"]
terrainBlank=all_data_level_0["boxlib", "z_phys"]

Outputs are numpy arrays and cane be plotted using matplotlib. I would prefer to natively plot the results using pyvista UI.

[MatthewFlamm]

I'm not sure whether this vtk reader is relevant here https://vtk.org/doc/nightly/html/classvtkAMReXGridReader.html

[tkoyama010]

I'm not sure whether this vtk reader is relevant here https://vtk.org/doc/nightly/html/classvtkAMReXGridReader.html

As a reference, I found the Sample code and data of vtk.

• https://gitlab.kitware.com/vtk/vtk/-/blob/v9.3.0/Filters/AMR/Testing/Python/TestAMRResampleFilter.py
• https://github.com/open-cv/VTKData/tree/master/Data/AMR/Enzo/DD0010