|
4 | 4 | from mpi4py import MPI
|
5 | 5 | import meshio
|
6 | 6 | import dolfinx
|
| 7 | +from os import path |
7 | 8 | from petsc4py import PETSc
|
8 | 9 | import ufl
|
9 | 10 | import sys
|
@@ -167,7 +168,6 @@ def boundary(x):
|
167 | 168 | T_per_vertex = [ self.uh.x.array[i] for i in range(self.mesh.geometry.x.shape[0]) if i in p_to_keep ]
|
168 | 169 | age_per_vertex = [ self.mesh_vertices_age[i] for i in range(self.mesh.geometry.x.shape[0]) if i in p_to_keep ]
|
169 | 170 |
|
170 |
| - from os import path |
171 | 171 | filename = path.join(out_path, self.modelName+'_'+str(self.tti)+'.epc')
|
172 | 172 | write_tetra_grid_with_properties(filename, np.array(points_cached), tet_renumbered, "tetramesh",
|
173 | 173 | np.array(T_per_vertex), np.array(age_per_vertex), poro0_per_cell, decay_per_cell, density_per_cell,
|
@@ -290,7 +290,6 @@ def write_hexa_mesh_resqml( self, out_path, tti):
|
290 | 290 | T_per_vertex_keep = [ T_per_vertex[i] for i in range(nv) if i in p_to_keep ]
|
291 | 291 | age_per_vertex_keep = [ self.age_per_vertex[i] for i in range(nv) if i in p_to_keep ]
|
292 | 292 |
|
293 |
| - from os import path |
294 | 293 | filename_hex = path.join(out_path, self.modelName+'_hexa_'+str(tti)+'.epc')
|
295 | 294 | points_cached=np.array(points_cached)
|
296 | 295 | write_hexa_grid_with_properties(filename_hex, points_cached, hexa_renumbered, "hexamesh",
|
@@ -355,13 +354,11 @@ def write_hexa_mesh_timeseries( self, out_path):
|
355 | 354 | for i in range(Temp_per_vertex_series.shape[1]):
|
356 | 355 | ts = Temp_per_vertex_series[:,i]
|
357 | 356 | ro = VR.easyRoDL(ts)
|
358 |
| - if (i%100==0): |
359 |
| - print(f"index {i} RO {ro[0:10]} {ro[-10:]}", flush=True) |
360 | 357 | Ro_per_vertex_series[:,i] = ro.flatten()
|
361 | 358 |
|
362 | 359 | hexa_renumbered = [ [point_original_to_cached[i] for i in hexa] for hexa in hexa_to_keep ]
|
363 | 360 |
|
364 |
| - from os import path |
| 361 | + |
365 | 362 |
|
366 | 363 | filename_hex = path.join(out_path, self.modelName+'_hexa_ts_'+str(self.tti)+'.epc')
|
367 | 364 | write_hexa_grid_with_timeseries(filename_hex, points_cached_series, hexa_renumbered, "hexamesh",
|
|
0 commit comments