Merge pull request #81 from keiyamamo/save_mesh

Save mesh

Author: keiyamamo

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)