Read Tensor Mesh

Read a UBC tensor mesh file

import pooch

# sphinx_gallery_thumbnail_number = 1
import PVGeo

Download sample data files and keep track of names:

url = "https://github.com/OpenGeoVis/PVGeo/raw/main/tests/data/Craig-Chile/craig_chile.msh"
mesh_file = pooch.retrieve(url=url, known_hash=None)
url = "https://github.com/OpenGeoVis/PVGeo/raw/main/tests/data/Craig-Chile/Lpout.mod"
model_file = pooch.retrieve(url=url, known_hash=None)

Read the mesh and model

reader = PVGeo.ubc.TensorMeshReader()
reader.set_mesh_filename(mesh_file)
reader.add_model_file_name(model_file)
mesh = reader.apply()
mesh

Header	Data Arrays
RectilinearGrid Information N Cells 190440 N Points 200900 X Bounds 3.550e+05, 3.722e+05 Y Bounds 5.999e+06, 6.016e+06 Z Bounds -5.250e+03, 3.000e+03 Dimensions 70, 70, 41 N Arrays 1	Name Field Type N Comp Min Max 93c3bf0f879573c07f76c0a0cbc55705-Lpout.mod Cells float64 1 -6.000e-01 3.000e-01

Use a PyVista threshold filter to remove NaN data values

mesh.threshold().plot()

Static Scene

Interactive Scene

Or inspect slices of the model

mesh.slice_orthogonal().plot()

Static Scene

Interactive Scene

Or threshold a data range

mesh.threshold([-0.6, -0.3]).plot(clim=[-0.6, 0.3])

Static Scene

Interactive Scene