updates examples to accommodate changes due to nummaterial_velocity_f…

…ile usage
danielpeter · Feb 8, 2024 · 7d45cfd · 7d45cfd
1 parent 784966c
commit 7d45cfd
Show file tree

Hide file tree

Showing 170 changed files with 1,090 additions and 577 deletions.
diff --git a/DATA/Par_file b/DATA/Par_file
@@ -279,15 +279,16 @@ read_external_mesh              = .false.
 
 # data concerning mesh, when generated using third-party app (more info in README)
 # (see also absorbing_conditions above)
-mesh_file                       = ./DATA/mesh_file          # file containing the mesh
-nodes_coords_file               = ./DATA/nodes_coords_file  # file containing the nodes coordinates
-materials_file                  = ./DATA/materials_file     # file containing the material number for each element
-free_surface_file               = ./DATA/free_surface_file  # file containing the free surface
-axial_elements_file             = ./DATA/axial_elements_file   # file containing the axial elements if AXISYM is true
-absorbing_surface_file          = ./DATA/absorbing_surface_file   # file containing the absorbing surface
-acoustic_forcing_surface_file   = ./DATA/MSH/Surf_acforcing_Bottom_enforcing_mesh   # file containing the acoustic forcing surface
-absorbing_cpml_file             = ./DATA/absorbing_cpml_file   # file containing the CPML element numbers
-tangential_detection_curve_file = ./DATA/courbe_eros_nodes  # file containing the curve delimiting the velocity model
+mesh_file                       = dummy                     # file containing the mesh
+nodes_coords_file               = dummy                     # file containing the nodes coordinates
+materials_file                  = dummy                     # file containing the material number for each element
+nummaterial_velocity_file       = dummy                     # file containing the material properties for each material
+free_surface_file               = dummy                     # file containing the free surface
+axial_elements_file             = dummy                     # file containing the axial elements if AXISYM is true
+absorbing_surface_file          = dummy                     # file containing the absorbing surface
+acoustic_forcing_surface_file   = dummy                     # file containing the acoustic forcing surface
+absorbing_cpml_file             = dummy                     # file containing the CPML element numbers
+tangential_detection_curve_file = dummy                     # file containing the curve delimiting the velocity model
 
 #-----------------------------------------------------------
 #

diff --git a/EXAMPLES/BENCHIE_water_PVC/DATA/Par_file b/EXAMPLES/BENCHIE_water_PVC/DATA/Par_file
@@ -278,15 +278,16 @@ read_external_mesh              = .true.
 
 # data concerning mesh, when generated using third-party app (more info in README)
 # (see also absorbing_conditions above)
-mesh_file                       = ./DATA/mesh_file          # file containing the mesh
-nodes_coords_file               = ./DATA/nodes_coords_file  # file containing the nodes coordinates
-materials_file                  = ./DATA/materials_file     # file containing the material number for each element
-free_surface_file               = ./DATA/free_surface_file  # file containing the free surface
-axial_elements_file             = ./DATA/axial_elements_file   # file containing the axial elements if AXISYM is true
-absorbing_surface_file          = ./DATA/absorbing_surface_file   # file containing the absorbing surface
-acoustic_forcing_surface_file   = ./DATA/MSH/Surf_acforcing_Bottom_enforcing_mesh   # file containing the acoustic forcing surface
-absorbing_cpml_file             = ./DATA/absorbing_cpml_file   # file containing the CPML element numbers
-tangential_detection_curve_file = ./DATA/courbe_eros_nodes  # file containing the curve delimiting the velocity model
+mesh_file                       = ./MESH/mesh_file                             # file containing the mesh
+nodes_coords_file               = ./MESH/nodes_coords_file                     # file containing the nodes coordinates
+materials_file                  = ./MESH/materials_file                        # file containing the material number for each element
+nummaterial_velocity_file       = ./MESH/nummaterial_velocity_file             # file containing the material properties for each material
+free_surface_file               = ./MESH/free_surface_file                     # file containing the free surface
+axial_elements_file             = ./MESH/axial_elements_file                   # file containing the axial elements if AXISYM is true
+absorbing_surface_file          = ./MESH/absorbing_surface_file                # file containing the absorbing surface
+acoustic_forcing_surface_file   = ./MESH/Surf_acforcing_Bottom_enforcing_mesh  # file containing the acoustic forcing surface
+absorbing_cpml_file             = ./MESH/absorbing_cpml_file                   # file containing the CPML element numbers
+tangential_detection_curve_file = ./MESH/courbe_eros_nodes                     # file containing the curve delimiting the velocity model
 
 #-----------------------------------------------------------
 #

diff --git a/...ENCHIE_water_PVC/DATA/absorbing_cpml_file → ...ENCHIE_water_PVC/MESH/absorbing_cpml_file b/...ENCHIE_water_PVC/DATA/absorbing_cpml_file → ...ENCHIE_water_PVC/MESH/absorbing_cpml_file
diff --git a/...HIE_water_PVC/DATA/absorbing_surface_file → ...HIE_water_PVC/MESH/absorbing_surface_file b/...HIE_water_PVC/DATA/absorbing_surface_file → ...HIE_water_PVC/MESH/absorbing_surface_file
diff --git a/.../BENCHIE_water_PVC/DATA/free_surface_file → .../BENCHIE_water_PVC/MESH/free_surface_file b/.../BENCHIE_water_PVC/DATA/free_surface_file → .../BENCHIE_water_PVC/MESH/free_surface_file
diff --git a/...LES/BENCHIE_water_PVC/DATA/materials_file → ...LES/BENCHIE_water_PVC/MESH/materials_file b/...LES/BENCHIE_water_PVC/DATA/materials_file → ...LES/BENCHIE_water_PVC/MESH/materials_file
diff --git a/EXAMPLES/BENCHIE_water_PVC/DATA/mesh_file → EXAMPLES/BENCHIE_water_PVC/MESH/mesh_file b/EXAMPLES/BENCHIE_water_PVC/DATA/mesh_file → EXAMPLES/BENCHIE_water_PVC/MESH/mesh_file
diff --git a/.../BENCHIE_water_PVC/DATA/nodes_coords_file → .../BENCHIE_water_PVC/MESH/nodes_coords_file b/.../BENCHIE_water_PVC/DATA/nodes_coords_file → .../BENCHIE_water_PVC/MESH/nodes_coords_file
diff --git a/EXAMPLES/BENCHIE_water_PVC/MESH/nummaterial_velocity_file b/EXAMPLES/BENCHIE_water_PVC/MESH/nummaterial_velocity_file
@@ -0,0 +1,33 @@
+# nummaterial_velocity_file
+# format:
+#(1)domain_id #(2)material_id #(3)rho #(4)vp #(5)vs #(6)Q_k #(7)Q_mu #(8)ani
+#
+#  where
+#     domain_id          : 1=acoustic / 2=elastic / 3=poroelastic
+#     material_id        : POSITIVE integer identifier of material block
+#     rho                : density
+#     vp                 : P-velocity
+#     vs                 : S-velocity
+#     Q_k                : 9999 = no Q_kappa attenuation
+#     Q_mu               : 9999 = no Q_mu attenuation
+#     ani                : 0=no anisotropy/ 1,2,.. check with aniso_model.f90
+#
+# example:
+# 2   1 2300 2800 1500 9999.0 9999.0 0
+#
+# or
+#
+#(1)domain_id #(2)material_id  tomography elastic  #(3)filename #(4)positive
+#
+#  where
+#     domain_id : 1=acoustic / 2=elastic / 3=poroelastic
+#     material_id        : NEGATIVE integer identifier of material block
+#     filename           : filename of the tomography file
+#     positive           : a positive unique identifier
+#
+# example:
+# 2  -1 tomography elastic tomo.xyz 1
+#
+# materials
+2 1 1412.d0 2260.d0 1050.d0 9999.0 9999.0 0
+1 2 1000.d0 1477.d0    0.d0 9999.0 9999.0 0  
diff --git a/EXAMPLES/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC/DATA/Par_file b/EXAMPLES/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC/DATA/Par_file
@@ -263,6 +263,7 @@ read_external_mesh              = .false.
 mesh_file                       = dummy                     # file containing the mesh
 nodes_coords_file               = dummy                     # file containing the nodes coordinates
 materials_file                  = dummy                     # file containing the material number for each element
+nummaterial_velocity_file       = dummy                     # file containing the material properties for each material
 free_surface_file               = dummy                     # file containing the free surface
 axial_elements_file             = dummy                     # file containing the axial elements if AXISYM is true
 absorbing_surface_file          = dummy                     # file containing the absorbing surface

diff --git a/EXAMPLES/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC_ELASTIC/DATA/Par_file b/EXAMPLES/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC_ELASTIC/DATA/Par_file
@@ -264,6 +264,7 @@ read_external_mesh              = .false.
 mesh_file                       = dummy                     # file containing the mesh
 nodes_coords_file               = dummy                     # file containing the nodes coordinates
 materials_file                  = dummy                     # file containing the material number for each element
+nummaterial_velocity_file       = dummy                     # file containing the material properties for each material
 free_surface_file               = dummy                     # file containing the free surface
 axial_elements_file             = dummy                     # file containing the axial elements if AXISYM is true
 absorbing_surface_file          = dummy                     # file containing the absorbing surface

diff --git a/EXAMPLES/BENCHMARK_CLAERBOUT_ADJOINT/ELASTIC/DATA/Par_file b/EXAMPLES/BENCHMARK_CLAERBOUT_ADJOINT/ELASTIC/DATA/Par_file
@@ -263,6 +263,7 @@ read_external_mesh              = .false.
 mesh_file                       = dummy                     # file containing the mesh
 nodes_coords_file               = dummy                     # file containing the nodes coordinates
 materials_file                  = dummy                     # file containing the material number for each element
+nummaterial_velocity_file       = dummy                     # file containing the material properties for each material
 free_surface_file               = dummy                     # file containing the free surface
 axial_elements_file             = dummy                     # file containing the axial elements if AXISYM is true
 absorbing_surface_file          = dummy                     # file containing the absorbing surface

diff --git a/EXAMPLES/Gmsh_example_Stacey_MPI/DATA/Par_file b/EXAMPLES/Gmsh_example_Stacey_MPI/DATA/Par_file
@@ -269,15 +269,16 @@ read_external_mesh              = .true.
 
 # data concerning mesh, when generated using third-party app (more info in README)
 # (see also absorbing_conditions above)
-mesh_file                       = MESH/Mesh_SqrCirc         # file containing the mesh
-nodes_coords_file               = MESH/Nodes_SqrCirc        # file containing the nodes coordinates
-materials_file                  = MESH/Material_SqrCirc     # file containing the material number for each element
-free_surface_file               = MESH/Surf_free_SqrCirc    # file containing the free surface
-axial_elements_file             = dummy                     # file containing the axial elements if AXISYM is true
-absorbing_surface_file          = MESH/Surf_abs_SqrCirc     # file containing the absorbing surface
-acoustic_forcing_surface_file   = dummy                     # file containing the acoustic forcing surface
-absorbing_cpml_file             = dummy                     # file containing the CPML element numbers
-tangential_detection_curve_file = dummy                     # file containing the curve delimiting the velocity model
+mesh_file                       = ./MESH/Mesh_SqrCirc               # file containing the mesh
+nodes_coords_file               = ./MESH/Nodes_SqrCirc              # file containing the nodes coordinates
+materials_file                  = ./MESH/Material_SqrCirc           # file containing the material number for each element
+nummaterial_velocity_file       = ./MESH/nummaterial_velocity_file  # file containing the material properties for each material
+free_surface_file               = ./MESH/Surf_free_SqrCirc          # file containing the free surface
+axial_elements_file             = dummy                             # file containing the axial elements if AXISYM is true
+absorbing_surface_file          = ./MESH/Surf_abs_SqrCirc           # file containing the absorbing surface
+acoustic_forcing_surface_file   = dummy                             # file containing the acoustic forcing surface
+absorbing_cpml_file             = dummy                             # file containing the CPML element numbers
+tangential_detection_curve_file = dummy                             # file containing the curve delimiting the velocity model
 
 #-----------------------------------------------------------
 #

diff --git a/EXAMPLES/Gmsh_example_Stacey_MPI/DATA/Par_file.serial b/EXAMPLES/Gmsh_example_Stacey_MPI/DATA/Par_file.serial
@@ -269,15 +269,16 @@ read_external_mesh              = .true.
 
 # data concerning mesh, when generated using third-party app (more info in README)
 # (see also absorbing_conditions above)
-mesh_file                       = MESH/Mesh_SqrCirc         # file containing the mesh
-nodes_coords_file               = MESH/Nodes_SqrCirc        # file containing the nodes coordinates
-materials_file                  = MESH/Material_SqrCirc     # file containing the material number for each element
-free_surface_file               = MESH/Surf_free_SqrCirc    # file containing the free surface
-axial_elements_file             = dummy                     # file containing the axial elements if AXISYM is true
-absorbing_surface_file          = MESH/Surf_abs_SqrCirc     # file containing the absorbing surface
-acoustic_forcing_surface_file   = dummy                     # file containing the acoustic forcing surface
-absorbing_cpml_file             = dummy                     # file containing the CPML element numbers
-tangential_detection_curve_file = dummy                     # file containing the curve delimiting the velocity model
+mesh_file                       = ./MESH/Mesh_SqrCirc               # file containing the mesh
+nodes_coords_file               = ./MESH/Nodes_SqrCirc              # file containing the nodes coordinates
+materials_file                  = ./MESH/Material_SqrCirc           # file containing the material number for each element
+nummaterial_velocity_file       = ./MESH/nummaterial_velocity_file  # file containing the material properties for each material
+free_surface_file               = ./MESH/Surf_free_SqrCirc          # file containing the free surface
+axial_elements_file             = dummy                             # file containing the axial elements if AXISYM is true
+absorbing_surface_file          = ./MESH/Surf_abs_SqrCirc           # file containing the absorbing surface
+acoustic_forcing_surface_file   = dummy                             # file containing the acoustic forcing surface
+absorbing_cpml_file             = dummy                             # file containing the CPML element numbers
+tangential_detection_curve_file = dummy                             # file containing the curve delimiting the velocity model
 
 #-----------------------------------------------------------
 #

diff --git a/EXAMPLES/Gmsh_example_Stacey_MPI/MESH/nummaterial_velocity_file b/EXAMPLES/Gmsh_example_Stacey_MPI/MESH/nummaterial_velocity_file
@@ -0,0 +1,33 @@
+# nummaterial_velocity_file
+# format:
+#(1)domain_id #(2)material_id #(3)rho #(4)vp #(5)vs #(6)Q_k #(7)Q_mu #(8)ani
+#
+#  where
+#     domain_id          : 1=acoustic / 2=elastic / 3=poroelastic
+#     material_id        : POSITIVE integer identifier of material block
+#     rho                : density
+#     vp                 : P-velocity
+#     vs                 : S-velocity
+#     Q_k                : 9999 = no Q_kappa attenuation
+#     Q_mu               : 9999 = no Q_mu attenuation
+#     ani                : 0=no anisotropy/ 1,2,.. check with aniso_model.f90
+#
+# example:
+# 2   1 2300 2800 1500 9999.0 9999.0 0
+#
+# or
+#
+#(1)domain_id #(2)material_id  tomography elastic  #(3)filename #(4)positive
+#
+#  where
+#     domain_id : 1=acoustic / 2=elastic / 3=poroelastic
+#     material_id        : NEGATIVE integer identifier of material block
+#     filename           : filename of the tomography file
+#     positive           : a positive unique identifier
+#
+# example:
+# 2  -1 tomography elastic tomo.xyz 1
+#
+# materials
+1 1 2000.d0 1680.d0 0.0d0 10.d0 10.d0 0
+1 2 1000.d0 1477.d0 0.0d0 10.d0 10.d0 0
diff --git a/EXAMPLES/Industrial_Format_SEP/DATA/Par_file b/EXAMPLES/Industrial_Format_SEP/DATA/Par_file
@@ -268,15 +268,16 @@ read_external_mesh              = .false.
 
 # data concerning mesh, when generated using third-party app (more info in README)
 # (see also absorbing_conditions above)
-mesh_file                       = ./mesh/mesh               # file containing the mesh
-nodes_coords_file               = ./mesh/nodes_coords       # file containing the nodes coordinates
-materials_file                  = ./mesh/mat                # file containing the material number for each element
-free_surface_file               = ./mesh/surface_free       # file containing the free surface
-axial_elements_file             = ./mesh/axial_elements_file   # file containing the axial elements if AXISYM is true
-absorbing_surface_file          = ./mesh/surface_abs        # file containing the absorbing surface
-acoustic_forcing_surface_file   = ./mesh/Surf_acforcing_Bottom_enforcing_mesh   # file containing the acoustic forcing surface
-absorbing_cpml_file             = ./mesh/Elements_CPML_list # file containing the CPML element numbers
-tangential_detection_curve_file = ./mesh/courbe_eros_nodes  # file containing the curve delimiting the velocity model
+mesh_file                       = dummy                     # file containing the mesh
+nodes_coords_file               = dummy                     # file containing the nodes coordinates
+materials_file                  = dummy                     # file containing the material number for each element
+nummaterial_velocity_file       = dummy                     # file containing the material properties for each material
+free_surface_file               = dummy                     # file containing the free surface
+axial_elements_file             = dummy                     # file containing the axial elements if AXISYM is true
+absorbing_surface_file          = dummy                     # file containing the absorbing surface
+acoustic_forcing_surface_file   = dummy                     # file containing the acoustic forcing surface
+absorbing_cpml_file             = dummy                     # file containing the CPML element numbers
+tangential_detection_curve_file = dummy                     # file containing the curve delimiting the velocity model
 
 #-----------------------------------------------------------
 #

diff --git a/EXAMPLES/Industrial_Format_SEP/REF_SEIS/Par_file.step_1 b/EXAMPLES/Industrial_Format_SEP/REF_SEIS/Par_file.step_1
@@ -268,15 +268,16 @@ read_external_mesh              = .false.
 
 # data concerning mesh, when generated using third-party app (more info in README)
 # (see also absorbing_conditions above)
-mesh_file                       = ./mesh/mesh               # file containing the mesh
-nodes_coords_file               = ./mesh/nodes_coords       # file containing the nodes coordinates
-materials_file                  = ./mesh/mat                # file containing the material number for each element
-free_surface_file               = ./mesh/surface_free       # file containing the free surface
-axial_elements_file             = ./mesh/axial_elements_file   # file containing the axial elements if AXISYM is true
-absorbing_surface_file          = ./mesh/surface_abs        # file containing the absorbing surface
-acoustic_forcing_surface_file   = ./mesh/Surf_acforcing_Bottom_enforcing_mesh   # file containing the acoustic forcing surface
-absorbing_cpml_file             = ./mesh/Elements_CPML_list # file containing the CPML element numbers
-tangential_detection_curve_file = ./mesh/courbe_eros_nodes  # file containing the curve delimiting the velocity model
+mesh_file                       = dummy                     # file containing the mesh
+nodes_coords_file               = dummy                     # file containing the nodes coordinates
+materials_file                  = dummy                     # file containing the material number for each element
+nummaterial_velocity_file       = dummy                     # file containing the material properties for each material
+free_surface_file               = dummy                     # file containing the free surface
+axial_elements_file             = dummy                     # file containing the axial elements if AXISYM is true
+absorbing_surface_file          = dummy                     # file containing the absorbing surface
+acoustic_forcing_surface_file   = dummy                     # file containing the acoustic forcing surface
+absorbing_cpml_file             = dummy                     # file containing the CPML element numbers
+tangential_detection_curve_file = dummy                     # file containing the curve delimiting the velocity model
 
 #-----------------------------------------------------------
 #