[cig-commits] [commit] devel, master: used Paul Cristini's script to add the new PML parameters to all the Par_files (3616a41)

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


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

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

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

commit 3616a417bd0c22f9391deba68d88b978e6baf796
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date:   Wed Jun 13 15:43:33 2012 +0000

    used Paul Cristini's script to add the new PML parameters to all the Par_files


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

3616a417bd0c22f9391deba68d88b978e6baf796
 .../Par_file_Abel_Balanche_bathy_source_solid      |  4 ++-
 Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in       |  4 ++-
 Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in    |  4 ++-
 INDUSTRIAL_FORMAT/Par_file                         |  4 ++-
 M2_UPPA/Par_file_M2_UPPA                           |  4 ++-
 Rayleigh_wave_no_crack/Par_file_Rayleigh_2D        |  4 ++-
 Rayleigh_wave_with_crack/Par_file_Rayleigh_2D      |  4 ++-
 Tape2007/Par_file_Tape2007_132rec_checker          |  4 ++-
 Tape2007/Par_file_Tape2007_onerec                  |  4 ++-
 Tape2007_kernel/Par_file_Tape2007_onerec           |  4 ++-
 Tromp2005/Par_file_Tromp2005                       |  4 ++-
 Tromp2005/Par_file_Tromp2005_s100                  |  4 ++-
 Tromp2005_kernel/Par_file_Tromp2005                |  4 ++-
 acoustic_poroelastic/Par_file_acoustic_poroelastic |  4 ++-
 attenuation/Par_file_attenuation_2D                |  4 ++-
 canyon/Par_file_canyon                             |  4 ++-
 .../fluid_solid_external_mesh/Par_file_fluid_solid |  4 ++-
 .../Par_file_fluid_solid                           |  4 ++-
 .../Par_file_fluid_solid                           |  4 ++-
 init_plane/Par_file_Slave                          |  4 ++-
 init_plane/Par_file_Slave_for                      |  4 ++-
 init_plane/Par_file_Slave_kernel                   |  4 ++-
 noise_layered/model_0/Par_file_fair                |  4 ++-
 noise_layered/model_0/Par_file_good                |  4 ++-
 noise_layered/model_1/Par_file_best                |  4 ++-
 noise_layered/model_1/Par_file_fair                |  4 ++-
 noise_layered/model_1/Par_file_good                |  4 ++-
 noise_layered/model_2/Par_file_fair                |  4 ++-
 noise_layered/model_2/Par_file_good                |  4 ++-
 noise_uniform/Par_file_noise_1                     |  4 ++-
 noise_uniform/Par_file_noise_2                     |  4 ++-
 noise_uniform/Par_file_noise_3                     |  4 ++-
 ...e_their_format_when_new_parameters_are_added.py | 31 ++++++++++++++++++++++
 semi_infinite_homo/Par_file_elastic_2D             |  4 ++-
 34 files changed, 130 insertions(+), 33 deletions(-)

diff --git a/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid b/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid
index c8c409e..4dddd55 100644
--- a/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid
+++ b/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid
@@ -116,7 +116,9 @@ nbmodels                        = 2              # nb of different models
 read_external_mesh              = .false.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
diff --git a/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in b/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in
index 90184cd..6139b94 100644
--- a/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in
+++ b/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in
@@ -112,7 +112,9 @@ nbmodels                        = 2              # nb of different models
 read_external_mesh              = .true.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
diff --git a/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in b/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in
index 31653c1..2aff56f 100644
--- a/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in
+++ b/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in
@@ -112,7 +112,9 @@ nbmodels                        = 2              # nb of different models
 read_external_mesh              = .true.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
diff --git a/INDUSTRIAL_FORMAT/Par_file b/INDUSTRIAL_FORMAT/Par_file
index 357a232..f3b4b17 100644
--- a/INDUSTRIAL_FORMAT/Par_file
+++ b/INDUSTRIAL_FORMAT/Par_file
@@ -110,7 +110,9 @@ nbmodels                        = 4              # nb of different models
 read_external_mesh              = .false.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
diff --git a/M2_UPPA/Par_file_M2_UPPA b/M2_UPPA/Par_file_M2_UPPA
index bd92106..9c92dba 100644
--- a/M2_UPPA/Par_file_M2_UPPA
+++ b/M2_UPPA/Par_file_M2_UPPA
@@ -110,7 +110,9 @@ nbmodels                        = 4              # nb of different models
 read_external_mesh              = .false.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
diff --git a/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D b/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D
index 57ff706..be01ce5 100644
--- a/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D
+++ b/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D
@@ -107,7 +107,9 @@ nbmodels                        = 1              # nb of different models
 read_external_mesh              = .false.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
diff --git a/Rayleigh_wave_with_crack/Par_file_Rayleigh_2D b/Rayleigh_wave_with_crack/Par_file_Rayleigh_2D
index 7910f18..ea5f4fe 100644
--- a/Rayleigh_wave_with_crack/Par_file_Rayleigh_2D
+++ b/Rayleigh_wave_with_crack/Par_file_Rayleigh_2D
@@ -111,7 +111,9 @@ nbmodels                        = 3              # nb of different models
 read_external_mesh              = .false.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
diff --git a/Tape2007/Par_file_Tape2007_132rec_checker b/Tape2007/Par_file_Tape2007_132rec_checker
index 5a2090a..bd63107 100644
--- a/Tape2007/Par_file_Tape2007_132rec_checker
+++ b/Tape2007/Par_file_Tape2007_132rec_checker
@@ -893,7 +893,9 @@ nbmodels                        = 1              # nb of different models
 read_external_mesh              = .false.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
diff --git a/Tape2007/Par_file_Tape2007_onerec b/Tape2007/Par_file_Tape2007_onerec
index 99e25a5..92eb79f 100644
--- a/Tape2007/Par_file_Tape2007_onerec
+++ b/Tape2007/Par_file_Tape2007_onerec
@@ -108,7 +108,9 @@ nbmodels                        = 1              # nb of different models
 read_external_mesh              = .false.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
diff --git a/Tape2007_kernel/Par_file_Tape2007_onerec b/Tape2007_kernel/Par_file_Tape2007_onerec
index 2fc2b21..8f2bc56 100644
--- a/Tape2007_kernel/Par_file_Tape2007_onerec
+++ b/Tape2007_kernel/Par_file_Tape2007_onerec
@@ -108,7 +108,9 @@ nbmodels                        = 1              # nb of different models
 read_external_mesh              = .false.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
diff --git a/Tromp2005/Par_file_Tromp2005 b/Tromp2005/Par_file_Tromp2005
index d677f14..5e4c4f0 100644
--- a/Tromp2005/Par_file_Tromp2005
+++ b/Tromp2005/Par_file_Tromp2005
@@ -107,7 +107,9 @@ nbmodels                        = 1              # nb of different models
 read_external_mesh              = .false.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
diff --git a/Tromp2005/Par_file_Tromp2005_s100 b/Tromp2005/Par_file_Tromp2005_s100
index b685496..72dd97a 100644
--- a/Tromp2005/Par_file_Tromp2005_s100
+++ b/Tromp2005/Par_file_Tromp2005_s100
@@ -107,7 +107,9 @@ nbmodels                        = 1              # nb of different models
 read_external_mesh              = .false.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
diff --git a/Tromp2005_kernel/Par_file_Tromp2005 b/Tromp2005_kernel/Par_file_Tromp2005
index 88e524a..550127b 100644
--- a/Tromp2005_kernel/Par_file_Tromp2005
+++ b/Tromp2005_kernel/Par_file_Tromp2005
@@ -107,7 +107,9 @@ nbmodels                        = 1              # nb of different models
 read_external_mesh              = .false.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
diff --git a/acoustic_poroelastic/Par_file_acoustic_poroelastic b/acoustic_poroelastic/Par_file_acoustic_poroelastic
index 381bf4e..5645523 100644
--- a/acoustic_poroelastic/Par_file_acoustic_poroelastic
+++ b/acoustic_poroelastic/Par_file_acoustic_poroelastic
@@ -116,7 +116,9 @@ nbmodels                        = 2              # nb of different models
 read_external_mesh              = .false.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
diff --git a/attenuation/Par_file_attenuation_2D b/attenuation/Par_file_attenuation_2D
index 8927950..1bce0c4 100644
--- a/attenuation/Par_file_attenuation_2D
+++ b/attenuation/Par_file_attenuation_2D
@@ -107,7 +107,9 @@ nbmodels                        = 1              # nb of different models
 read_external_mesh              = .false.
 
 # absorbing boundary active or not
-absorbing_conditions            = .false.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
diff --git a/canyon/Par_file_canyon b/canyon/Par_file_canyon
index 5a8d55b..b80e5ba 100644
--- a/canyon/Par_file_canyon
+++ b/canyon/Par_file_canyon
@@ -107,7 +107,9 @@ nbmodels                        = 1              # nb of different models
 read_external_mesh              = .true.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
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 2a9610d..ee4573c 100644
--- a/fluid_solid/fluid_solid_external_mesh/Par_file_fluid_solid
+++ b/fluid_solid/fluid_solid_external_mesh/Par_file_fluid_solid
@@ -112,7 +112,9 @@ nbmodels                        = 2              # nb of different models
 read_external_mesh              = .true.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
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 fc510af..423f002 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
@@ -109,7 +109,9 @@ nbmodels                        = 2              # nb of different models
 read_external_mesh              = .false.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
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 194ef75..21aee31 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
@@ -109,7 +109,9 @@ nbmodels                        = 2              # nb of different models
 read_external_mesh              = .false.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
diff --git a/init_plane/Par_file_Slave b/init_plane/Par_file_Slave
index 0c8d704..d856db4 100644
--- a/init_plane/Par_file_Slave
+++ b/init_plane/Par_file_Slave
@@ -108,7 +108,9 @@ nbmodels                        = 2              # nb of different models
 read_external_mesh              = .false.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
diff --git a/init_plane/Par_file_Slave_for b/init_plane/Par_file_Slave_for
index e4f66c1..6d37f1f 100644
--- a/init_plane/Par_file_Slave_for
+++ b/init_plane/Par_file_Slave_for
@@ -108,7 +108,9 @@ nbmodels                        = 2              # nb of different models
 read_external_mesh              = .false.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
diff --git a/init_plane/Par_file_Slave_kernel b/init_plane/Par_file_Slave_kernel
index d7fd79f..69f2f86 100644
--- a/init_plane/Par_file_Slave_kernel
+++ b/init_plane/Par_file_Slave_kernel
@@ -108,7 +108,9 @@ nbmodels                        = 2              # nb of different models
 read_external_mesh              = .false.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
diff --git a/noise_layered/model_0/Par_file_fair b/noise_layered/model_0/Par_file_fair
index baaed09..995b334 100644
--- a/noise_layered/model_0/Par_file_fair
+++ b/noise_layered/model_0/Par_file_fair
@@ -107,7 +107,9 @@ nbmodels                        = 1              # nb of different models
 read_external_mesh              = .false.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
diff --git a/noise_layered/model_0/Par_file_good b/noise_layered/model_0/Par_file_good
index 34906b2..e0f65f3 100644
--- a/noise_layered/model_0/Par_file_good
+++ b/noise_layered/model_0/Par_file_good
@@ -107,7 +107,9 @@ nbmodels                        = 1              # nb of different models
 read_external_mesh              = .false.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
diff --git a/noise_layered/model_1/Par_file_best b/noise_layered/model_1/Par_file_best
index a0361cc..fbf553a 100644
--- a/noise_layered/model_1/Par_file_best
+++ b/noise_layered/model_1/Par_file_best
@@ -108,7 +108,9 @@ nbmodels                        = 2              # nb of different models
 read_external_mesh              = .false.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
diff --git a/noise_layered/model_1/Par_file_fair b/noise_layered/model_1/Par_file_fair
index 95296bb..6e0506a 100644
--- a/noise_layered/model_1/Par_file_fair
+++ b/noise_layered/model_1/Par_file_fair
@@ -108,7 +108,9 @@ nbmodels                        = 2              # nb of different models
 read_external_mesh              = .false.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
diff --git a/noise_layered/model_1/Par_file_good b/noise_layered/model_1/Par_file_good
index 7b0a484..e56ca64 100644
--- a/noise_layered/model_1/Par_file_good
+++ b/noise_layered/model_1/Par_file_good
@@ -108,7 +108,9 @@ nbmodels                        = 2              # nb of different models
 read_external_mesh              = .false.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
diff --git a/noise_layered/model_2/Par_file_fair b/noise_layered/model_2/Par_file_fair
index ad175bf..faac520 100644
--- a/noise_layered/model_2/Par_file_fair
+++ b/noise_layered/model_2/Par_file_fair
@@ -109,7 +109,9 @@ nbmodels                        = 3              # nb of different models
 read_external_mesh              = .false.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
diff --git a/noise_layered/model_2/Par_file_good b/noise_layered/model_2/Par_file_good
index 916ea77..c669aee 100644
--- a/noise_layered/model_2/Par_file_good
+++ b/noise_layered/model_2/Par_file_good
@@ -109,7 +109,9 @@ nbmodels                        = 3              # nb of different models
 read_external_mesh              = .false.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
diff --git a/noise_uniform/Par_file_noise_1 b/noise_uniform/Par_file_noise_1
index 6a226bf..2b8f341 100644
--- a/noise_uniform/Par_file_noise_1
+++ b/noise_uniform/Par_file_noise_1
@@ -107,7 +107,9 @@ nbmodels                        = 1              # nb of different models
 read_external_mesh              = .false.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
diff --git a/noise_uniform/Par_file_noise_2 b/noise_uniform/Par_file_noise_2
index 14ad00a..265c11b 100644
--- a/noise_uniform/Par_file_noise_2
+++ b/noise_uniform/Par_file_noise_2
@@ -107,7 +107,9 @@ nbmodels                        = 1              # nb of different models
 read_external_mesh              = .false.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
diff --git a/noise_uniform/Par_file_noise_3 b/noise_uniform/Par_file_noise_3
index 8f9ef7d..4b98758 100644
--- a/noise_uniform/Par_file_noise_3
+++ b/noise_uniform/Par_file_noise_3
@@ -107,7 +107,9 @@ nbmodels                        = 1              # nb of different models
 read_external_mesh              = .false.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
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 f820e33..d3554d0 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
@@ -385,6 +385,36 @@ def ProcessParfile_r20307(fic):
     print 'xxxxx------> '+fic+' processed to r20307'
     return
 #------------------------------------------------------------------------------
+def ProcessParfile_r20359(fic):
+    # Open the file and get all lines from Par_file
+    ligs= LoadLig(fic)
+
+    # Test if already processed
+    for lig in ligs:
+        if 'NELEM_PML_THICKNESS' in lig:
+            print '----> '+fic+' already processed to r20359'            
+            return
+    #
+    a1='PML_BOUNDARY_CONDITIONS         = .true.\nNELEM_PML_THICKNESS  ' + \
+    '           = 3\nSTACEY_ABSORBING_CONDITIONS     = .false.\n'
+
+    #--------------------------------------------------------------------------
+    # Add new parameters
+    # 
+    for ilg, lig in enumerate(ligs):
+        if lig.startswith('absorbing_conditions'):
+            ligs.pop(ilg)
+            ligs.insert(ilg,a1)
+    #
+    move(fic,fic+'.before_update_to_r20359')
+    #
+    fm = open(fic,'w')
+    fm.writelines(ligs)
+    fm.close()
+    #
+    print 'xxxxx------> '+fic+' processed to r20359'
+    return
+#------------------------------------------------------------------------------
 if __name__=='__main__':
     ## List of all files of current directory
     Fichiers=[]
@@ -407,6 +437,7 @@ if __name__=='__main__':
                     ProcessParfile_r19521(fic)
                     ProcessParfile_r19804(fic)
                     ProcessParfile_r20307(fic)
+                    ProcessParfile_r20359(fic)
                 print '~'*80
     #                
     print 'Number of Par_file analysed : ', Ct_Par_file   
diff --git a/semi_infinite_homo/Par_file_elastic_2D b/semi_infinite_homo/Par_file_elastic_2D
index 052f93d..5eba7b9 100644
--- a/semi_infinite_homo/Par_file_elastic_2D
+++ b/semi_infinite_homo/Par_file_elastic_2D
@@ -107,7 +107,9 @@ nbmodels                        = 1              # nb of different models
 read_external_mesh              = .false.
 
 # absorbing boundary active or not
-absorbing_conditions            = .true.
+PML_BOUNDARY_CONDITIONS         = .true.
+NELEM_PML_THICKNESS             = 3
+STACEY_ABSORBING_CONDITIONS     = .false.
 ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges



More information about the CIG-COMMITS mailing list