[cig-commits] [commit] devel, master: Paul fixed problems in Par_file when using PML with external mesh (58325a0)

cig_noreply at geodynamics.org cig_noreply at geodynamics.org
Wed Jun 18 15:21:32 PDT 2014


Repository : https://github.com/geodynamics/specfem2d

On branches: devel,master
Link       : https://github.com/geodynamics/specfem2d/compare/fc67e6fd7ad890705b2b72b4b3c509accb22249e...e9ca46c40131588d89d7b0883250bc6584ce6b4c

>---------------------------------------------------------------

commit 58325a08ca778271540da64e775e5d8fd5a940ac
Author: Xie Zhinan <xiezhinan1984 at gmail.com>
Date:   Thu Aug 9 13:46:38 2012 +0000

    Paul fixed problems in Par_file when using PML with external mesh


>---------------------------------------------------------------

58325a08ca778271540da64e775e5d8fd5a940ac
 .../Par_file_Abel_Balanche_bathy_source_solid      |  1 +
 Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in       |  1 +
 Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in    |  1 +
 INDUSTRIAL_FORMAT/Par_file                         |  1 +
 M2_UPPA/Par_file_M2_UPPA                           |  1 +
 Rayleigh_wave_no_crack/Par_file_Rayleigh_2D        |  1 +
 Rayleigh_wave_with_crack/Par_file_Rayleigh_2D      |  1 +
 Tape2007/Par_file_Tape2007_132rec_checker          |  1 +
 Tape2007/Par_file_Tape2007_onerec                  |  1 +
 Tape2007_kernel/Par_file_Tape2007_onerec           |  1 +
 Tromp2005/Par_file_Tromp2005                       |  1 +
 Tromp2005/Par_file_Tromp2005_s100                  |  1 +
 Tromp2005_kernel/Par_file_Tromp2005                |  1 +
 acoustic_poroelastic/Par_file_acoustic_poroelastic |  1 +
 attenuation/Par_file_attenuation_2D                |  1 +
 canyon/Par_file_canyon                             |  1 +
 .../fluid_solid_external_mesh/Par_file_fluid_solid |  1 +
 .../Par_file_fluid_solid                           |  1 +
 .../Par_file_fluid_solid                           |  1 +
 init_plane/Par_file_Slave                          |  1 +
 init_plane/Par_file_Slave_for                      |  1 +
 init_plane/Par_file_Slave_kernel                   |  1 +
 noise_layered/model_0/Par_file_fair                |  1 +
 noise_layered/model_0/Par_file_good                |  1 +
 noise_layered/model_1/Par_file_best                |  1 +
 noise_layered/model_1/Par_file_fair                |  1 +
 noise_layered/model_1/Par_file_good                |  1 +
 noise_layered/model_2/Par_file_fair                |  1 +
 noise_layered/model_2/Par_file_good                |  1 +
 noise_uniform/Par_file_noise_1                     |  1 +
 noise_uniform/Par_file_noise_2                     |  1 +
 noise_uniform/Par_file_noise_3                     |  1 +
 ...e_their_format_when_new_parameters_are_added.py |  9 +++--
 semi_infinite_homo/Par_file_elastic_2D             | 11 +++---
 .../plot_points_per_wavelength_histogram.gnu       |  8 ++---
 .../points_per_wavelength_histogram_S_in_solid.txt | 40 +++++++++++-----------
 36 files changed, 66 insertions(+), 34 deletions(-)

diff --git a/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid b/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid
index 4d610fe..bb9a460 100644
--- a/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid
+++ b/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid
@@ -139,6 +139,7 @@ mesh_file                       = ./DATA/Mesh_canyon/canyon_mesh_file   # file c
 nodes_coords_file               = ./DATA/Mesh_canyon/canyon_nodes_coords_file   # file containing the nodes coordinates
 materials_file                  = ./DATA/Mesh_canyon/canyon_materials_file   # file containing the material number for each element
 free_surface_file               = ./DATA/Mesh_canyon/canyon_free_surface_file   # file containing the free surface
+absorbing_surface_file          = ./DATA/Mesh_canyon/canyon_absorbing_surface_file   # file containing the absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in b/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in
index 6573c30..c010dc9 100644
--- a/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in
+++ b/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in
@@ -135,6 +135,7 @@ mesh_file                       = Mesh_SqrCirc   # file containing the mesh
 nodes_coords_file               = Nodes_SqrCirc    # file containing the nodes coordinates
 materials_file                  = Material_SqrCirc  # file containing the material number for each element
 free_surface_file               = Surf_free_SqrCirc   # file containing the free surface
+absorbing_surface_file          = Surf_abs_SqrCirc   # file containing the absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in b/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in
index 21a5ea5..72ae513 100644
--- a/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in
+++ b/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in
@@ -135,6 +135,7 @@ mesh_file                       = Mesh_SqrCirc   # file containing the mesh
 nodes_coords_file               = Nodes_SqrCirc    # file containing the nodes coordinates
 materials_file                  = Material_SqrCirc  # file containing the material number for each element
 free_surface_file               = Surf_free_SqrCirc   # file containing the free surface
+absorbing_surface_file          = Surf_abs_SqrCirc   # file containing the absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/INDUSTRIAL_FORMAT/Par_file b/INDUSTRIAL_FORMAT/Par_file
index 6b7fae0..bc5bb94 100644
--- a/INDUSTRIAL_FORMAT/Par_file
+++ b/INDUSTRIAL_FORMAT/Par_file
@@ -133,6 +133,7 @@ mesh_file                       = ./DATA/Mesh_canyon/canyon_mesh_file   # file c
 nodes_coords_file               = ./DATA/Mesh_canyon/canyon_nodes_coords_file   # file containing the nodes coordinates
 materials_file                  = ./DATA/Mesh_canyon/canyon_materials_file   # file containing the material number for each element
 free_surface_file               = ./DATA/Mesh_canyon/canyon_free_surface_file   # file containing the free surface
+absorbing_surface_file          = ./DATA/Mesh_canyon/canyon_absorbing_surface_file   # file containing the absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/M2_UPPA/Par_file_M2_UPPA b/M2_UPPA/Par_file_M2_UPPA
index 3e5b9b4..cf64618 100644
--- a/M2_UPPA/Par_file_M2_UPPA
+++ b/M2_UPPA/Par_file_M2_UPPA
@@ -133,6 +133,7 @@ mesh_file                       = ./DATA/Mesh_canyon/canyon_mesh_file   # file c
 nodes_coords_file               = ./DATA/Mesh_canyon/canyon_nodes_coords_file   # file containing the nodes coordinates
 materials_file                  = ./DATA/Mesh_canyon/canyon_materials_file   # file containing the material number for each element
 free_surface_file               = ./DATA/Mesh_canyon/canyon_free_surface_file   # file containing the free surface
+absorbing_surface_file          = ./DATA/Mesh_canyon/canyon_absorbing_surface_file   # file containing the absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D b/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D
index bc3f0a4..b3251a9 100644
--- a/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D
+++ b/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D
@@ -130,6 +130,7 @@ mesh_file                       = ./DATA/mesh/canyon_mesh_file   # file containi
 nodes_coords_file               = ./DATA/mesh/canyon_nodes_coords_file   # file containing the nodes coordinates
 materials_file                  = ./DATA/mesh/canyon_materials_file   # file containing the material number for each element
 free_surface_file               = ./DATA/mesh/canyon_free_surface_file   # file containing the free surface
+absorbing_surface_file          = ./DATA/mesh/canyon_absorbing_surface_file   # file containing the absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/Rayleigh_wave_with_crack/Par_file_Rayleigh_2D b/Rayleigh_wave_with_crack/Par_file_Rayleigh_2D
index f106dce..6ecd0c6 100644
--- a/Rayleigh_wave_with_crack/Par_file_Rayleigh_2D
+++ b/Rayleigh_wave_with_crack/Par_file_Rayleigh_2D
@@ -134,6 +134,7 @@ mesh_file                       = ./DATA/mesh/canyon_mesh_file   # file containi
 nodes_coords_file               = ./DATA/mesh/canyon_nodes_coords_file   # file containing the nodes coordinates
 materials_file                  = ./DATA/mesh/canyon_materials_file   # file containing the material number for each element
 free_surface_file               = ./DATA/mesh/canyon_free_surface_file   # file containing the free surface
+absorbing_surface_file          = ./DATA/mesh/canyon_absorbing_surface_file   # file containing the absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/Tape2007/Par_file_Tape2007_132rec_checker b/Tape2007/Par_file_Tape2007_132rec_checker
index 0db92d0..07531d8 100644
--- a/Tape2007/Par_file_Tape2007_132rec_checker
+++ b/Tape2007/Par_file_Tape2007_132rec_checker
@@ -916,6 +916,7 @@ mesh_file                       = ./DATA/ice_water_rock_1D/ice_water_rock_1D_ele
 nodes_coords_file               = ./DATA/ice_water_rock_1D/ice_water_rock_1D_nodes          # file containing node coordinates
 materials_file                  = ./DATA/ice_water_rock_1D/ice_water_rock_1D_material       # file containing material index for each element
 free_surface_file               = ./DATA/ice_water_rock_1D/ice_water_rock_1D_surface_free   # file containing free surface
+absorbing_surface_file          = ./DATA/ice_water_rock_1D/ice_water_rock_1D_surface_absorb # file containing absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/Tape2007/Par_file_Tape2007_onerec b/Tape2007/Par_file_Tape2007_onerec
index bac2027..d94cc74 100644
--- a/Tape2007/Par_file_Tape2007_onerec
+++ b/Tape2007/Par_file_Tape2007_onerec
@@ -131,6 +131,7 @@ mesh_file                       = ./DATA/ice_water_rock_1D/ice_water_rock_1D_ele
 nodes_coords_file               = ./DATA/ice_water_rock_1D/ice_water_rock_1D_nodes          # file containing node coordinates
 materials_file                  = ./DATA/ice_water_rock_1D/ice_water_rock_1D_material       # file containing material index for each element
 free_surface_file               = ./DATA/ice_water_rock_1D/ice_water_rock_1D_surface_free   # file containing free surface
+absorbing_surface_file          = ./DATA/ice_water_rock_1D/ice_water_rock_1D_surface_absorb # file containing absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/Tape2007_kernel/Par_file_Tape2007_onerec b/Tape2007_kernel/Par_file_Tape2007_onerec
index b907cb7..4979d59 100644
--- a/Tape2007_kernel/Par_file_Tape2007_onerec
+++ b/Tape2007_kernel/Par_file_Tape2007_onerec
@@ -131,6 +131,7 @@ mesh_file                       = ./DATA/ice_water_rock_1D/ice_water_rock_1D_ele
 nodes_coords_file               = ./DATA/ice_water_rock_1D/ice_water_rock_1D_nodes          # file containing node coordinates
 materials_file                  = ./DATA/ice_water_rock_1D/ice_water_rock_1D_material       # file containing material index for each element
 free_surface_file               = ./DATA/ice_water_rock_1D/ice_water_rock_1D_surface_free   # file containing free surface
+absorbing_surface_file          = ./DATA/ice_water_rock_1D/ice_water_rock_1D_surface_absorb # file containing absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/Tromp2005/Par_file_Tromp2005 b/Tromp2005/Par_file_Tromp2005
index 0a0e2da..adeef85 100644
--- a/Tromp2005/Par_file_Tromp2005
+++ b/Tromp2005/Par_file_Tromp2005
@@ -130,6 +130,7 @@ mesh_file                       = ./DATA/ice_water_rock_1D/ice_water_rock_1D_ele
 nodes_coords_file               = ./DATA/ice_water_rock_1D/ice_water_rock_1D_nodes          # file containing node coordinates
 materials_file                  = ./DATA/ice_water_rock_1D/ice_water_rock_1D_material       # file containing material index for each element
 free_surface_file               = ./DATA/ice_water_rock_1D/ice_water_rock_1D_surface_free   # file containing free surface
+absorbing_surface_file          = ./DATA/ice_water_rock_1D/ice_water_rock_1D_surface_absorb # file containing absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/Tromp2005/Par_file_Tromp2005_s100 b/Tromp2005/Par_file_Tromp2005_s100
index efe35ed..d906477 100644
--- a/Tromp2005/Par_file_Tromp2005_s100
+++ b/Tromp2005/Par_file_Tromp2005_s100
@@ -130,6 +130,7 @@ mesh_file                       = ./DATA/ice_water_rock_1D/ice_water_rock_1D_ele
 nodes_coords_file               = ./DATA/ice_water_rock_1D/ice_water_rock_1D_nodes          # file containing node coordinates
 materials_file                  = ./DATA/ice_water_rock_1D/ice_water_rock_1D_material       # file containing material index for each element
 free_surface_file               = ./DATA/ice_water_rock_1D/ice_water_rock_1D_surface_free   # file containing free surface
+absorbing_surface_file          = ./DATA/ice_water_rock_1D/ice_water_rock_1D_surface_absorb # file containing absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/Tromp2005_kernel/Par_file_Tromp2005 b/Tromp2005_kernel/Par_file_Tromp2005
index d54005c..617825a 100644
--- a/Tromp2005_kernel/Par_file_Tromp2005
+++ b/Tromp2005_kernel/Par_file_Tromp2005
@@ -130,6 +130,7 @@ mesh_file                       = ./DATA/ice_water_rock_1D/ice_water_rock_1D_ele
 nodes_coords_file               = ./DATA/ice_water_rock_1D/ice_water_rock_1D_nodes          # file containing node coordinates
 materials_file                  = ./DATA/ice_water_rock_1D/ice_water_rock_1D_material       # file containing material index for each element
 free_surface_file               = ./DATA/ice_water_rock_1D/ice_water_rock_1D_surface_free   # file containing free surface
+absorbing_surface_file          = ./DATA/ice_water_rock_1D/ice_water_rock_1D_surface_absorb # file containing absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/acoustic_poroelastic/Par_file_acoustic_poroelastic b/acoustic_poroelastic/Par_file_acoustic_poroelastic
index be9b91d..463b040 100644
--- a/acoustic_poroelastic/Par_file_acoustic_poroelastic
+++ b/acoustic_poroelastic/Par_file_acoustic_poroelastic
@@ -139,6 +139,7 @@ mesh_file                       = ./DATA/Mesh_canyon/canyon_mesh_file   # file c
 nodes_coords_file               = ./DATA/Mesh_canyon/canyon_nodes_coords_file   # file containing the nodes coordinates
 materials_file                  = ./DATA/Mesh_canyon/canyon_materials_file   # file containing the material number for each element
 free_surface_file               = ./DATA/Mesh_canyon/canyon_free_surface_file   # file containing the free surface
+absorbing_surface_file          = ./DATA/Mesh_canyon/canyon_absorbing_surface_file   # file containing the absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/attenuation/Par_file_attenuation_2D b/attenuation/Par_file_attenuation_2D
index 81c83e9..f35e3ae 100644
--- a/attenuation/Par_file_attenuation_2D
+++ b/attenuation/Par_file_attenuation_2D
@@ -130,6 +130,7 @@ mesh_file                       = ./DATA/Mesh_canyon/canyon_mesh_file   # file c
 nodes_coords_file               = ./DATA/Mesh_canyon/canyon_nodes_coords_file   # file containing the nodes coordinates
 materials_file                  = ./DATA/Mesh_canyon/canyon_materials_file   # file containing the material number for each element
 free_surface_file               = ./DATA/Mesh_canyon/canyon_free_surface_file   # file containing the free surface
+absorbing_surface_file          = ./DATA/Mesh_canyon/canyon_absorbing_surface_file   # file containing the absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/canyon/Par_file_canyon b/canyon/Par_file_canyon
index 90f762f..9a1fcd1 100644
--- a/canyon/Par_file_canyon
+++ b/canyon/Par_file_canyon
@@ -130,6 +130,7 @@ mesh_file                       = ./DATA/mesh/canyon_mesh_file   # file containi
 nodes_coords_file               = ./DATA/mesh/canyon_nodes_coords_file   # file containing the nodes coordinates
 materials_file                  = ./DATA/mesh/canyon_materials_file   # file containing the material number for each element
 free_surface_file               = ./DATA/mesh/canyon_free_surface_file   # file containing the free surface
+absorbing_surface_file          = ./DATA/mesh/canyon_absorbing_surface_file   # file containing the absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/fluid_solid/fluid_solid_external_mesh/Par_file_fluid_solid b/fluid_solid/fluid_solid_external_mesh/Par_file_fluid_solid
index 34df167..058dddd 100644
--- a/fluid_solid/fluid_solid_external_mesh/Par_file_fluid_solid
+++ b/fluid_solid/fluid_solid_external_mesh/Par_file_fluid_solid
@@ -135,6 +135,7 @@ mesh_file                       = ./mesh/mesh          # file containing the mes
 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
+absorbing_surface_file          = ./mesh/surface_abs   # file containing the absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/Par_file_fluid_solid b/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/Par_file_fluid_solid
index db11f5d..f62ca48 100644
--- a/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/Par_file_fluid_solid
+++ b/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/Par_file_fluid_solid
@@ -132,6 +132,7 @@ mesh_file                       = ./DATA/Mesh_canyon/canyon_mesh_file   # file c
 nodes_coords_file               = ./DATA/Mesh_canyon/canyon_nodes_coords_file   # file containing the nodes coordinates
 materials_file                  = ./DATA/Mesh_canyon/canyon_materials_file   # file containing the material number for each element
 free_surface_file               = ./DATA/Mesh_canyon/canyon_free_surface_file   # file containing the free surface
+absorbing_surface_file          = ./DATA/Mesh_canyon/canyon_absorbing_surface_file   # file containing the absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/Par_file_fluid_solid b/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/Par_file_fluid_solid
index b8dd0d5..3c4123d 100644
--- a/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/Par_file_fluid_solid
+++ b/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/Par_file_fluid_solid
@@ -132,6 +132,7 @@ mesh_file                       = ./DATA/Mesh_canyon/canyon_mesh_file   # file c
 nodes_coords_file               = ./DATA/Mesh_canyon/canyon_nodes_coords_file   # file containing the nodes coordinates
 materials_file                  = ./DATA/Mesh_canyon/canyon_materials_file   # file containing the material number for each element
 free_surface_file               = ./DATA/Mesh_canyon/canyon_free_surface_file   # file containing the free surface
+absorbing_surface_file          = ./DATA/Mesh_canyon/canyon_absorbing_surface_file   # file containing the absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/init_plane/Par_file_Slave b/init_plane/Par_file_Slave
index 988f243..8a62491 100644
--- a/init_plane/Par_file_Slave
+++ b/init_plane/Par_file_Slave
@@ -131,6 +131,7 @@ mesh_file                       = ./DATA/ice_water_rock_1D/ice_water_rock_1D_ele
 nodes_coords_file               = ./DATA/ice_water_rock_1D/ice_water_rock_1D_nodes          # file containing node coordinates
 materials_file                  = ./DATA/ice_water_rock_1D/ice_water_rock_1D_material       # file containing material index for each element
 free_surface_file               = ./DATA/ice_water_rock_1D/ice_water_rock_1D_surface_free   # file containing free surface
+absorbing_surface_file          = ./DATA/ice_water_rock_1D/ice_water_rock_1D_surface_absorb # file containing absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/init_plane/Par_file_Slave_for b/init_plane/Par_file_Slave_for
index 73f4c1e..dde9f7b 100644
--- a/init_plane/Par_file_Slave_for
+++ b/init_plane/Par_file_Slave_for
@@ -131,6 +131,7 @@ mesh_file                       = ./DATA/ice_water_rock_1D/ice_water_rock_1D_ele
 nodes_coords_file               = ./DATA/ice_water_rock_1D/ice_water_rock_1D_nodes          # file containing node coordinates
 materials_file                  = ./DATA/ice_water_rock_1D/ice_water_rock_1D_material       # file containing material index for each element
 free_surface_file               = ./DATA/ice_water_rock_1D/ice_water_rock_1D_surface_free   # file containing free surface
+absorbing_surface_file          = ./DATA/ice_water_rock_1D/ice_water_rock_1D_surface_absorb # file containing absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/init_plane/Par_file_Slave_kernel b/init_plane/Par_file_Slave_kernel
index b565bef..5adf56e 100644
--- a/init_plane/Par_file_Slave_kernel
+++ b/init_plane/Par_file_Slave_kernel
@@ -131,6 +131,7 @@ mesh_file                       = ./DATA/ice_water_rock_1D/ice_water_rock_1D_ele
 nodes_coords_file               = ./DATA/ice_water_rock_1D/ice_water_rock_1D_nodes          # file containing node coordinates
 materials_file                  = ./DATA/ice_water_rock_1D/ice_water_rock_1D_material       # file containing material index for each element
 free_surface_file               = ./DATA/ice_water_rock_1D/ice_water_rock_1D_surface_free   # file containing free surface
+absorbing_surface_file          = ./DATA/ice_water_rock_1D/ice_water_rock_1D_surface_absorb # file containing absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/noise_layered/model_0/Par_file_fair b/noise_layered/model_0/Par_file_fair
index 68626b1..4c8742d 100644
--- a/noise_layered/model_0/Par_file_fair
+++ b/noise_layered/model_0/Par_file_fair
@@ -130,6 +130,7 @@ mesh_file                       = ./DATA/Mesh_canyon/canyon_mesh_file   # file c
 nodes_coords_file               = ./DATA/Mesh_canyon/canyon_nodes_coords_file   # file containing the nodes coordinates
 materials_file                  = ./DATA/Mesh_canyon/canyon_materials_file   # file containing the material number for each element
 free_surface_file               = ./DATA/Mesh_canyon/canyon_free_surface_file   # file containing the free surface
+absorbing_surface_file          = ./DATA/Mesh_canyon/canyon_absorbing_surface_file   # file containing the absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/noise_layered/model_0/Par_file_good b/noise_layered/model_0/Par_file_good
index 061e3d6..254df64 100644
--- a/noise_layered/model_0/Par_file_good
+++ b/noise_layered/model_0/Par_file_good
@@ -130,6 +130,7 @@ mesh_file                       = ./DATA/Mesh_canyon/canyon_mesh_file   # file c
 nodes_coords_file               = ./DATA/Mesh_canyon/canyon_nodes_coords_file   # file containing the nodes coordinates
 materials_file                  = ./DATA/Mesh_canyon/canyon_materials_file   # file containing the material number for each element
 free_surface_file               = ./DATA/Mesh_canyon/canyon_free_surface_file   # file containing the free surface
+absorbing_surface_file          = ./DATA/Mesh_canyon/canyon_absorbing_surface_file   # file containing the absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/noise_layered/model_1/Par_file_best b/noise_layered/model_1/Par_file_best
index 3c0a148..09b2e98 100644
--- a/noise_layered/model_1/Par_file_best
+++ b/noise_layered/model_1/Par_file_best
@@ -131,6 +131,7 @@ mesh_file                       = ./DATA/Mesh_canyon/canyon_mesh_file   # file c
 nodes_coords_file               = ./DATA/Mesh_canyon/canyon_nodes_coords_file   # file containing the nodes coordinates
 materials_file                  = ./DATA/Mesh_canyon/canyon_materials_file   # file containing the material number for each element
 free_surface_file               = ./DATA/Mesh_canyon/canyon_free_surface_file   # file containing the free surface
+absorbing_surface_file          = ./DATA/Mesh_canyon/canyon_absorbing_surface_file   # file containing the absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/noise_layered/model_1/Par_file_fair b/noise_layered/model_1/Par_file_fair
index 913c28d..2cc307f 100644
--- a/noise_layered/model_1/Par_file_fair
+++ b/noise_layered/model_1/Par_file_fair
@@ -131,6 +131,7 @@ mesh_file                       = ./DATA/Mesh_canyon/canyon_mesh_file   # file c
 nodes_coords_file               = ./DATA/Mesh_canyon/canyon_nodes_coords_file   # file containing the nodes coordinates
 materials_file                  = ./DATA/Mesh_canyon/canyon_materials_file   # file containing the material number for each element
 free_surface_file               = ./DATA/Mesh_canyon/canyon_free_surface_file   # file containing the free surface
+absorbing_surface_file          = ./DATA/Mesh_canyon/canyon_absorbing_surface_file   # file containing the absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/noise_layered/model_1/Par_file_good b/noise_layered/model_1/Par_file_good
index ba241d6..d7f63ed 100644
--- a/noise_layered/model_1/Par_file_good
+++ b/noise_layered/model_1/Par_file_good
@@ -131,6 +131,7 @@ mesh_file                       = ./DATA/Mesh_canyon/canyon_mesh_file   # file c
 nodes_coords_file               = ./DATA/Mesh_canyon/canyon_nodes_coords_file   # file containing the nodes coordinates
 materials_file                  = ./DATA/Mesh_canyon/canyon_materials_file   # file containing the material number for each element
 free_surface_file               = ./DATA/Mesh_canyon/canyon_free_surface_file   # file containing the free surface
+absorbing_surface_file          = ./DATA/Mesh_canyon/canyon_absorbing_surface_file   # file containing the absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/noise_layered/model_2/Par_file_fair b/noise_layered/model_2/Par_file_fair
index 71b07c0..2c0ae75 100644
--- a/noise_layered/model_2/Par_file_fair
+++ b/noise_layered/model_2/Par_file_fair
@@ -132,6 +132,7 @@ mesh_file                       = ./DATA/Mesh_canyon/canyon_mesh_file   # file c
 nodes_coords_file               = ./DATA/Mesh_canyon/canyon_nodes_coords_file   # file containing the nodes coordinates
 materials_file                  = ./DATA/Mesh_canyon/canyon_materials_file   # file containing the material number for each element
 free_surface_file               = ./DATA/Mesh_canyon/canyon_free_surface_file   # file containing the free surface
+absorbing_surface_file          = ./DATA/Mesh_canyon/canyon_absorbing_surface_file   # file containing the absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/noise_layered/model_2/Par_file_good b/noise_layered/model_2/Par_file_good
index 4becc4c..c85fb7f 100644
--- a/noise_layered/model_2/Par_file_good
+++ b/noise_layered/model_2/Par_file_good
@@ -132,6 +132,7 @@ mesh_file                       = ./DATA/Mesh_canyon/canyon_mesh_file   # file c
 nodes_coords_file               = ./DATA/Mesh_canyon/canyon_nodes_coords_file   # file containing the nodes coordinates
 materials_file                  = ./DATA/Mesh_canyon/canyon_materials_file   # file containing the material number for each element
 free_surface_file               = ./DATA/Mesh_canyon/canyon_free_surface_file   # file containing the free surface
+absorbing_surface_file          = ./DATA/Mesh_canyon/canyon_absorbing_surface_file   # file containing the absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/noise_uniform/Par_file_noise_1 b/noise_uniform/Par_file_noise_1
index 1fba4af..19a538b 100644
--- a/noise_uniform/Par_file_noise_1
+++ b/noise_uniform/Par_file_noise_1
@@ -130,6 +130,7 @@ mesh_file                       = ./DATA/Mesh_canyon/canyon_mesh_file   # file c
 nodes_coords_file               = ./DATA/Mesh_canyon/canyon_nodes_coords_file   # file containing the nodes coordinates
 materials_file                  = ./DATA/Mesh_canyon/canyon_materials_file   # file containing the material number for each element
 free_surface_file               = ./DATA/Mesh_canyon/canyon_free_surface_file   # file containing the free surface
+absorbing_surface_file          = ./DATA/Mesh_canyon/canyon_absorbing_surface_file   # file containing the absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/noise_uniform/Par_file_noise_2 b/noise_uniform/Par_file_noise_2
index bf8c473..d501b6d 100644
--- a/noise_uniform/Par_file_noise_2
+++ b/noise_uniform/Par_file_noise_2
@@ -130,6 +130,7 @@ mesh_file                       = ./DATA/Mesh_canyon/canyon_mesh_file   # file c
 nodes_coords_file               = ./DATA/Mesh_canyon/canyon_nodes_coords_file   # file containing the nodes coordinates
 materials_file                  = ./DATA/Mesh_canyon/canyon_materials_file   # file containing the material number for each element
 free_surface_file               = ./DATA/Mesh_canyon/canyon_free_surface_file   # file containing the free surface
+absorbing_surface_file          = ./DATA/Mesh_canyon/canyon_absorbing_surface_file   # file containing the absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/noise_uniform/Par_file_noise_3 b/noise_uniform/Par_file_noise_3
index c5fa1c4..317c44e 100644
--- a/noise_uniform/Par_file_noise_3
+++ b/noise_uniform/Par_file_noise_3
@@ -130,6 +130,7 @@ mesh_file                       = ./DATA/Mesh_canyon/canyon_mesh_file   # file c
 nodes_coords_file               = ./DATA/Mesh_canyon/canyon_nodes_coords_file   # file containing the nodes coordinates
 materials_file                  = ./DATA/Mesh_canyon/canyon_materials_file   # file containing the material number for each element
 free_surface_file               = ./DATA/Mesh_canyon/canyon_free_surface_file   # file containing the free surface
+absorbing_surface_file          = ./DATA/Mesh_canyon/canyon_absorbing_surface_file   # file containing the absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/process_DATA_Par_files_to_update_their_format_when_new_parameters_are_added.py b/process_DATA_Par_files_to_update_their_format_when_new_parameters_are_added.py
index ae1ec3a..1d0917a 100755
--- a/process_DATA_Par_files_to_update_their_format_when_new_parameters_are_added.py
+++ b/process_DATA_Par_files_to_update_their_format_when_new_parameters_are_added.py
@@ -38,7 +38,7 @@ def OuvreParFile(fic):
         return
     #
     for var in VariableNum:
-        print var
+        #print var
         for ligne in lignes:
             lsplit=string.split(ligne)
             if lsplit!= []:
@@ -432,9 +432,8 @@ def ProcessParfile_r20561(fic):
     # Add new parameters
     # 
     for ilg, lig in enumerate(ligs):
-        if lig.startswith('absorbing_surface_file'):
-            ligs.pop(ilg)
-            ligs.insert(ilg,a1)
+        if lig.startswith('absorbing_surface'):
+            ligs.insert(ilg+1,a1)
     #
     move(fic,fic+'.before_update_to_r20561')
     #
@@ -469,7 +468,7 @@ if __name__=='__main__':
                     ProcessParfile_r20307(fic)
                     ProcessParfile_r20359(fic)
                     ProcessParfile_r20561(fic)
-                print '~'*80
+                    print '~'*80
     #                
     print 'Number of Par_file analysed : ', Ct_Par_file   
     print 'END OF Par_file PROCESSING'
diff --git a/semi_infinite_homo/Par_file_elastic_2D b/semi_infinite_homo/Par_file_elastic_2D
index 4a22f97..2944b9d 100644
--- a/semi_infinite_homo/Par_file_elastic_2D
+++ b/semi_infinite_homo/Par_file_elastic_2D
@@ -23,7 +23,7 @@ freq0                           =  10            # frequency for viscous attenua
 p_sv                            = .true.         # set the type of calculation (P-SV or SH/membrane waves)
 
 # time step parameters
-nt                              = 1500           # total number of time steps
+nt                              = 15000           # total number of time steps
 deltat                          = 0.8e-3         # duration of a time step
 USER_T0                         = 0.15d0         # use this t0 as earliest starting time rather than the automatically calculated one
 time_stepping_scheme            = 1              # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical 4th-order 4-stage Runge-Kutta
@@ -91,7 +91,7 @@ use_binary_for_wavefield_dumps  = .false.        # use ASCII or single-precision
 ####
 output_grid_Gnuplot              = .false.        # save the grid in a text file or not
 output_grid_ASCII               = .false.        # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
-output_energy                   = .false.        # compute and output acoustic and elastic energy (slows down the code significantly)
+output_energy                   = .true.        # compute and output acoustic and elastic energy (slows down the code significantly)
 
 # velocity and density models
 nbmodels                        = 1              # nb of different models
@@ -101,15 +101,15 @@ nbmodels                        = 1              # nb of different models
 # III: (model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu).
 # For istropic elastic/acoustic material use I and set Vs to zero to make a given model acoustic, for anisotropic elastic use II,
 # and for isotropic poroelastic material use III. The mesh can contain acoustic, elastic, and poroelastic models simultaneously.
-1 1 2700.d0 3000.d0 1732.051d0 0 0 20. 10. 0 0 0 0 0 0
+1 1 2700.d0 3000.d0 0.d0 0 0 20. 10. 0 0 0 0 0 0
 
 # external mesh or not
 read_external_mesh              = .false.
 
 # absorbing boundary active or not
-PML_BOUNDARY_CONDITIONS         = .false.
+PML_BOUNDARY_CONDITIONS         = .true.
 NELEM_PML_THICKNESS             = 3
-STACEY_ABSORBING_CONDITIONS     = .true.
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .false.
 
 # for horizontal periodic conditions: detect common points between left and right edges
@@ -130,6 +130,7 @@ mesh_file                       = ./DATA/Mesh_canyon/canyon_mesh_file   # file c
 nodes_coords_file               = ./DATA/Mesh_canyon/canyon_nodes_coords_file   # file containing the nodes coordinates
 materials_file                  = ./DATA/Mesh_canyon/canyon_materials_file   # file containing the material number for each element
 free_surface_file               = ./DATA/Mesh_canyon/canyon_free_surface_file   # file containing the free surface
+absorbing_surface_file          = ./DATA/Mesh_canyon/canyon_absorbing_surface_file   # file containing the absorbing surface
 CPML_element_file               = EltPML_xxxxxx  # file containing the CPML element numbers
 tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
 
diff --git a/semi_infinite_homo/plot_points_per_wavelength_histogram.gnu b/semi_infinite_homo/plot_points_per_wavelength_histogram.gnu
index 0f77db4..a108347 100644
--- a/semi_infinite_homo/plot_points_per_wavelength_histogram.gnu
+++ b/semi_infinite_homo/plot_points_per_wavelength_histogram.gnu
@@ -1,9 +1,9 @@
  set term x11
  #set term gif
- #set output "points_per_wavelength_histogram_S_in_solid.gif"
+ #set output "points_per_wavelength_histogram_P_in_fluid.gif"
  
- set boxwidth   8.6602550E-03
- set xlabel "Range of min number of points per S wavelength in solid"
+ set boxwidth   1.5000000E-02
+ set xlabel "Range of min number of points per P wavelength in fluid"
  set ylabel "Percentage of elements (%)"
- plot "points_per_wavelength_histogram_S_in_solid.txt" with boxes
+ plot "points_per_wavelength_histogram_P_in_fluid.txt" with boxes
  pause -1 "hit any key..."
diff --git a/semi_infinite_homo/points_per_wavelength_histogram_S_in_solid.txt b/semi_infinite_homo/points_per_wavelength_histogram_S_in_solid.txt
index dd6e32f..2e24674 100644
--- a/semi_infinite_homo/points_per_wavelength_histogram_S_in_solid.txt
+++ b/semi_infinite_homo/points_per_wavelength_histogram_S_in_solid.txt
@@ -1,20 +1,20 @@
-   8.577982        0.0000000E+00
-   8.586643        0.0000000E+00
-   8.595303        0.0000000E+00
-   8.603964        0.0000000E+00
-   8.612623        0.0000000E+00
-   8.621284        0.0000000E+00
-   8.629944        0.0000000E+00
-   8.638605        0.0000000E+00
-   8.647264        0.0000000E+00
-   8.655926         2.000000    
-   8.664585         98.00000    
-   8.673244        0.0000000E+00
-   8.681906        0.0000000E+00
-   8.690565        0.0000000E+00
-   8.699226        0.0000000E+00
-   8.707886        0.0000000E+00
-   8.716547        0.0000000E+00
-   8.725206        0.0000000E+00
-   8.733868        0.0000000E+00
-   8.742527        0.0000000E+00
+   8.5779819         0.0000000    
+   8.5866432         0.0000000    
+   8.5953026         0.0000000    
+   8.6039639         0.0000000    
+   8.6126232         0.0000000    
+   8.6212845         0.0000000    
+   8.6299438         0.0000000    
+   8.6386051         0.0000000    
+   8.6472645         0.0000000    
+   8.6559258         14.000000    
+   8.6645851         86.000000    
+   8.6732445         0.0000000    
+   8.6819057         0.0000000    
+   8.6905651         0.0000000    
+   8.6992264         0.0000000    
+   8.7078857         0.0000000    
+   8.7165470         0.0000000    
+   8.7252064         0.0000000    
+   8.7338676         0.0000000    
+   8.7425270         0.0000000    



More information about the CIG-COMMITS mailing list