Added Free surface

PrincetonUniversity · Feb 21, 2025 · fee63b1 · fee63b1
1 parent e1e98bc
commit fee63b1
Show file tree

Hide file tree

Showing 2 changed files with 12 additions and 0 deletions.
diff --git a/include/IO/fortranio/fortran_io.tpp b/include/IO/fortranio/fortran_io.tpp
@@ -54,6 +54,7 @@ void fortran_read_line(std::ifstream &stream, Args... values) {
   } catch (const std::exception &e) {
     std::ostringstream error_message;
     error_message << "Error reading fortran line with buffer length: "
+                  << buffer_length << "\n"
                   << e.what() << "(" << __FILE__ << ":" << __LINE__ << ")";
     throw std::runtime_error(error_message.str());
   }

diff --git a/src/IO/mesh/impl/fortran/dim3/mesh.cpp b/src/IO/mesh/impl/fortran/dim3/mesh.cpp
@@ -409,6 +409,12 @@ specfem::IO::read_3d_mesh(const std::string mesh_parameters_file,
 #endif
   }
 
+  // Read the number of free surface faces
+  int num_free_surface_faces;
+  try_read_line("num_free_surface_faces", stream, &num_free_surface_faces);
+  check_values("num_free_surface_faces", num_free_surface_faces,
+               mesh.parameters.num_free_surface_faces);
+
   // Create the free surface object
   mesh.free_surface =
       specfem::mesh::free_surface<specfem::dimension::type::dim3>(
@@ -421,6 +427,11 @@ specfem::IO::read_3d_mesh(const std::string mesh_parameters_file,
                  mesh.free_surface.jacobian2Dw);
   try_read_array("free_surface_normal", stream, mesh.free_surface.normal);
 
+#ifndef NDEBUG
+  // Print the free surface
+  mesh.free_surface.print();
+#endif
+
   stream.close();
 
   return mesh;