Merge pull request #81 from keiyamamo/save_mesh

Save mesh
KVSlab · Sep 4, 2023 · 803b318 · 803b318
2 parents a209f92 + be09bdd
commit 803b318
Showing 1 changed file with 8 additions and 0 deletions.
diff --git a/turtleFSI/monolithic.py b/turtleFSI/monolithic.py
@@ -62,6 +62,14 @@
 # Get mesh information
 mesh, domains, boundaries = get_mesh_domain_and_boundaries(**vars())
 
+# Save mesh, domains, and boundaries for post-processing
+if restart_folder is None:
+    mesh_path = results_folder.joinpath("Mesh", "mesh.h5")
+    with HDF5File(mesh.mpi_comm(), mesh_path.__str__(), "w") as hdf:
+        hdf.write(mesh, "/mesh")
+        hdf.write(boundaries, "/boundaries")
+        hdf.write(domains, "/domains")
+
 # Control FEniCS output
 set_log_level(loglevel)