[cig-commits] [commit] devel, master: added the missing C Delta_t / 2 term in the mass matrix for Clayton-Engquist absorbing elements (mentioned by Jean-Paul Ampuero in 2003, see the book of Hughes 1987 chapter 9). Zhinan Xie also added a spring to further stabilize the Clayton-Engquist absorbing conditions. (bae684c)

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


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

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

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

commit bae684c7b645dd62c9ab7d019b473b7fdc574095
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date:   Mon Mar 19 02:12:09 2012 +0000

    added the missing C Delta_t / 2 term in the mass matrix for Clayton-Engquist absorbing elements (mentioned by Jean-Paul Ampuero in 2003, see the book of Hughes 1987 chapter 9).
    Zhinan Xie also added a spring to further stabilize the Clayton-Engquist absorbing conditions.


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

bae684c7b645dd62c9ab7d019b473b7fdc574095
 Abel_Brest/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 +
 change_all_2D_Par_files.csh                                            | 3 ++-
 fluid_solid/fluid_solid_external_mesh/Par_file_fluid_solid             | 1 +
 .../from_2000_Geophysics_paper_flat_ocean_bottom/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 +
 semi_infinite_homo/Par_file_elastic_2D                                 | 1 +
 34 files changed, 35 insertions(+), 1 deletion(-)

diff --git a/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid b/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid
index 8af3730..a18c766 100644
--- a/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid
+++ b/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid
@@ -99,6 +99,7 @@ read_external_mesh              = .false.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
diff --git a/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in b/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in
index 98c55e9..cf1ae34 100644
--- a/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in
+++ b/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in
@@ -95,6 +95,7 @@ read_external_mesh              = .true.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
diff --git a/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in b/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in
index b88dedb..121327f 100644
--- a/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in
+++ b/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in
@@ -95,6 +95,7 @@ read_external_mesh              = .true.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
diff --git a/INDUSTRIAL_FORMAT/Par_file b/INDUSTRIAL_FORMAT/Par_file
index 059de1d..e84e157 100644
--- a/INDUSTRIAL_FORMAT/Par_file
+++ b/INDUSTRIAL_FORMAT/Par_file
@@ -93,6 +93,7 @@ read_external_mesh              = .false.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
diff --git a/M2_UPPA/Par_file_M2_UPPA b/M2_UPPA/Par_file_M2_UPPA
index 3065cf7..b2e2735 100644
--- a/M2_UPPA/Par_file_M2_UPPA
+++ b/M2_UPPA/Par_file_M2_UPPA
@@ -93,6 +93,7 @@ read_external_mesh              = .false.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
diff --git a/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D b/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D
index a68e20c..a8e2fe4 100644
--- a/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D
+++ b/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D
@@ -90,6 +90,7 @@ read_external_mesh              = .false.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
diff --git a/Rayleigh_wave_with_crack/Par_file_Rayleigh_2D b/Rayleigh_wave_with_crack/Par_file_Rayleigh_2D
index ed10434..691e2f4 100644
--- a/Rayleigh_wave_with_crack/Par_file_Rayleigh_2D
+++ b/Rayleigh_wave_with_crack/Par_file_Rayleigh_2D
@@ -94,6 +94,7 @@ read_external_mesh              = .false.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
diff --git a/Tape2007/Par_file_Tape2007_132rec_checker b/Tape2007/Par_file_Tape2007_132rec_checker
index 7887b4d..ea927a7 100644
--- a/Tape2007/Par_file_Tape2007_132rec_checker
+++ b/Tape2007/Par_file_Tape2007_132rec_checker
@@ -876,6 +876,7 @@ read_external_mesh              = .false.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
diff --git a/Tape2007/Par_file_Tape2007_onerec b/Tape2007/Par_file_Tape2007_onerec
index 6526740..b39ecdc 100644
--- a/Tape2007/Par_file_Tape2007_onerec
+++ b/Tape2007/Par_file_Tape2007_onerec
@@ -91,6 +91,7 @@ read_external_mesh              = .false.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
diff --git a/Tape2007_kernel/Par_file_Tape2007_onerec b/Tape2007_kernel/Par_file_Tape2007_onerec
index 9c9f962..9220eee 100644
--- a/Tape2007_kernel/Par_file_Tape2007_onerec
+++ b/Tape2007_kernel/Par_file_Tape2007_onerec
@@ -91,6 +91,7 @@ read_external_mesh              = .false.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
diff --git a/Tromp2005/Par_file_Tromp2005 b/Tromp2005/Par_file_Tromp2005
index e71fbae..b027f25 100644
--- a/Tromp2005/Par_file_Tromp2005
+++ b/Tromp2005/Par_file_Tromp2005
@@ -90,6 +90,7 @@ read_external_mesh              = .false.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
diff --git a/Tromp2005/Par_file_Tromp2005_s100 b/Tromp2005/Par_file_Tromp2005_s100
index 83cdb55..c6ea35f 100644
--- a/Tromp2005/Par_file_Tromp2005_s100
+++ b/Tromp2005/Par_file_Tromp2005_s100
@@ -90,6 +90,7 @@ read_external_mesh              = .false.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
diff --git a/Tromp2005_kernel/Par_file_Tromp2005 b/Tromp2005_kernel/Par_file_Tromp2005
index 5e1b3ab..9ff3e45 100644
--- a/Tromp2005_kernel/Par_file_Tromp2005
+++ b/Tromp2005_kernel/Par_file_Tromp2005
@@ -90,6 +90,7 @@ read_external_mesh              = .false.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
diff --git a/acoustic_poroelastic/Par_file_acoustic_poroelastic b/acoustic_poroelastic/Par_file_acoustic_poroelastic
index 5567df1..55c3791 100644
--- a/acoustic_poroelastic/Par_file_acoustic_poroelastic
+++ b/acoustic_poroelastic/Par_file_acoustic_poroelastic
@@ -99,6 +99,7 @@ read_external_mesh              = .false.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
diff --git a/attenuation/Par_file_attenuation_2D b/attenuation/Par_file_attenuation_2D
index dad2450..4834b8c 100644
--- a/attenuation/Par_file_attenuation_2D
+++ b/attenuation/Par_file_attenuation_2D
@@ -90,6 +90,7 @@ read_external_mesh              = .false.
 
 # absorbing boundary active or not
 absorbing_conditions            = .false.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
diff --git a/canyon/Par_file_canyon b/canyon/Par_file_canyon
index 33426f5..a22e9d7 100644
--- a/canyon/Par_file_canyon
+++ b/canyon/Par_file_canyon
@@ -90,6 +90,7 @@ read_external_mesh              = .true.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
diff --git a/change_all_2D_Par_files.csh b/change_all_2D_Par_files.csh
index 0c69a42..ad211f3 100755
--- a/change_all_2D_Par_files.csh
+++ b/change_all_2D_Par_files.csh
@@ -11,11 +11,12 @@ foreach file ( DATA/Par_file EXAMPLES/Abel_Brest/Par_file_Abel_Balanche_bathy_so
 set newfile = "../"$file
 
 # DK DK leave the white spaces, in order to have the right alignment with other variables
+###sed -e '1,$s/PERFORM_CUTHILL_MCKEE           = .true.         # perform inverse Cuthill-McKee (1969) optimization\/permutation for mesh numbering/PERFORM_CUTHILL_MCKEE           = .false.        # perform inverse Cuthill-McKee (1969) optimization\/permutation for mesh numbering (can be very costly and not very useful)/g' < $newfile > __________temp_27_zzzyyyyxxxx__________
 sed -e '1,$s/PERFORM_CUTHILL_MCKEE           = .true.         # perform inverse Cuthill-McKee (1969) optimization\/permutation for mesh numbering/PERFORM_CUTHILL_MCKEE           = .false.        # perform inverse Cuthill-McKee (1969) optimization\/permutation for mesh numbering (can be very costly and not very useful)/g' < $newfile > __________temp_27_zzzyyyyxxxx__________
 
 clean_a_Par_file.py __________temp_27_zzzyyyyxxxx__________ $file
 
-mv __________temp_27_zzzyyyyxxxx__________ $newfile
+#mv __________temp_27_zzzyyyyxxxx__________ $newfile
 
 end
 
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 eedd0cb..7944191 100644
--- a/fluid_solid/fluid_solid_external_mesh/Par_file_fluid_solid
+++ b/fluid_solid/fluid_solid_external_mesh/Par_file_fluid_solid
@@ -95,6 +95,7 @@ read_external_mesh              = .true.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
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 b0a3889..faa4d66 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
@@ -92,6 +92,7 @@ read_external_mesh              = .false.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
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 08866d9..6388a4b 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
@@ -92,6 +92,7 @@ read_external_mesh              = .false.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
diff --git a/init_plane/Par_file_Slave b/init_plane/Par_file_Slave
index f60abd9..6f9deac 100644
--- a/init_plane/Par_file_Slave
+++ b/init_plane/Par_file_Slave
@@ -91,6 +91,7 @@ read_external_mesh              = .false.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
diff --git a/init_plane/Par_file_Slave_for b/init_plane/Par_file_Slave_for
index 1c7125f..1c90e77 100644
--- a/init_plane/Par_file_Slave_for
+++ b/init_plane/Par_file_Slave_for
@@ -91,6 +91,7 @@ read_external_mesh              = .false.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
diff --git a/init_plane/Par_file_Slave_kernel b/init_plane/Par_file_Slave_kernel
index 464dae2..b07e30b 100644
--- a/init_plane/Par_file_Slave_kernel
+++ b/init_plane/Par_file_Slave_kernel
@@ -91,6 +91,7 @@ read_external_mesh              = .false.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
diff --git a/noise_layered/model_0/Par_file_fair b/noise_layered/model_0/Par_file_fair
index ac2c3d1..74d514c 100644
--- a/noise_layered/model_0/Par_file_fair
+++ b/noise_layered/model_0/Par_file_fair
@@ -90,6 +90,7 @@ read_external_mesh              = .false.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
diff --git a/noise_layered/model_0/Par_file_good b/noise_layered/model_0/Par_file_good
index 9e7042d..b7e1e98 100644
--- a/noise_layered/model_0/Par_file_good
+++ b/noise_layered/model_0/Par_file_good
@@ -90,6 +90,7 @@ read_external_mesh              = .false.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
diff --git a/noise_layered/model_1/Par_file_best b/noise_layered/model_1/Par_file_best
index f3a08cb..a2d4203 100644
--- a/noise_layered/model_1/Par_file_best
+++ b/noise_layered/model_1/Par_file_best
@@ -91,6 +91,7 @@ read_external_mesh              = .false.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
diff --git a/noise_layered/model_1/Par_file_fair b/noise_layered/model_1/Par_file_fair
index ebae388..60c3a58 100644
--- a/noise_layered/model_1/Par_file_fair
+++ b/noise_layered/model_1/Par_file_fair
@@ -91,6 +91,7 @@ read_external_mesh              = .false.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
diff --git a/noise_layered/model_1/Par_file_good b/noise_layered/model_1/Par_file_good
index ce3f70c..fe38332 100644
--- a/noise_layered/model_1/Par_file_good
+++ b/noise_layered/model_1/Par_file_good
@@ -91,6 +91,7 @@ read_external_mesh              = .false.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
diff --git a/noise_layered/model_2/Par_file_fair b/noise_layered/model_2/Par_file_fair
index d459a65..db75d7a 100644
--- a/noise_layered/model_2/Par_file_fair
+++ b/noise_layered/model_2/Par_file_fair
@@ -92,6 +92,7 @@ read_external_mesh              = .false.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
diff --git a/noise_layered/model_2/Par_file_good b/noise_layered/model_2/Par_file_good
index 61259f1..2cea418 100644
--- a/noise_layered/model_2/Par_file_good
+++ b/noise_layered/model_2/Par_file_good
@@ -92,6 +92,7 @@ read_external_mesh              = .false.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
diff --git a/noise_uniform/Par_file_noise_1 b/noise_uniform/Par_file_noise_1
index 318e88c..36d8df4 100644
--- a/noise_uniform/Par_file_noise_1
+++ b/noise_uniform/Par_file_noise_1
@@ -90,6 +90,7 @@ read_external_mesh              = .false.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
diff --git a/noise_uniform/Par_file_noise_2 b/noise_uniform/Par_file_noise_2
index 32c5163..7b12ef9 100644
--- a/noise_uniform/Par_file_noise_2
+++ b/noise_uniform/Par_file_noise_2
@@ -90,6 +90,7 @@ read_external_mesh              = .false.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
diff --git a/noise_uniform/Par_file_noise_3 b/noise_uniform/Par_file_noise_3
index 8ea307e..7514754 100644
--- a/noise_uniform/Par_file_noise_3
+++ b/noise_uniform/Par_file_noise_3
@@ -90,6 +90,7 @@ read_external_mesh              = .false.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.
diff --git a/semi_infinite_homo/Par_file_elastic_2D b/semi_infinite_homo/Par_file_elastic_2D
index 016ffd8..3339e6e 100644
--- a/semi_infinite_homo/Par_file_elastic_2D
+++ b/semi_infinite_homo/Par_file_elastic_2D
@@ -90,6 +90,7 @@ read_external_mesh              = .false.
 
 # absorbing boundary active or not
 absorbing_conditions            = .true.
+ADD_SPRING_TO_STACEY            = .true.
 
 # for horizontal periodic conditions: detect common points between left and right edges
 ADD_PERIODIC_CONDITIONS         = .false.



More information about the CIG-COMMITS mailing list