Skip to content
Snippets Groups Projects
Commit 0e975e47 authored by Pooria Danaeifar's avatar Pooria Danaeifar
Browse files

Fix get_dimensions for multiblocks

parent 7e84d58d
No related branches found
No related tags found
1 merge request!85Add pyvista plot
......@@ -24,6 +24,8 @@
params.control_dict.end_time = 0.002
sim = Simul(params)
sim.make.exec("run")
sim.output.fields.plot_mesh(color="w", show=True)
sim.make.exec("clean")
......
......@@ -73,9 +73,18 @@
def get_dimensions(mesh):
interior_mesh = mesh[0]
centers = interior_mesh.cell_centers()
points = centers.points
n_blocks = mesh.n_blocks
if n_blocks > 2:
data = mesh[0][0]
points = data.cell_centers().points
for block in range(1, n_blocks):
data0 = mesh[block][0]
points0 = data0.cell_centers().points
points = np.concatenate((points, points0), axis=0)
else:
interior_mesh = mesh[0]
centers = interior_mesh.cell_centers()
points = centers.points
cell_coords = x, y, z = points[:, 0], points[:, 1], points[:, 2]
dimensions = ""
for letter, coord in zip("xyz", cell_coords):
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment