[cig-commits] [commit] devel, master: initial version of titled PML commited, please becarful to use it, since the current pml_init.F90 code using to define the pml parameter do not support titled mesh. (2d68b73)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Wed Jun 18 15:22:38 PDT 2014
Repository : https://github.com/geodynamics/specfem2d
On branches: devel,master
Link : https://github.com/geodynamics/specfem2d/compare/fc67e6fd7ad890705b2b72b4b3c509accb22249e...e9ca46c40131588d89d7b0883250bc6584ce6b4c
>---------------------------------------------------------------
commit 2d68b7391231f0d52b64f4f40a8840de524fdca0
Author: Xie Zhinan <xiezhinan1984 at gmail.com>
Date: Mon Jan 7 15:44:21 2013 +0000
initial version of titled PML commited, please becarful to use it, since the current pml_init.F90 code using to define the pml parameter do not support titled mesh.
>---------------------------------------------------------------
2d68b7391231f0d52b64f4f40a8840de524fdca0
.../Par_file_Abel_Balanche_bathy_source_solid | 2 ++
Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in | 2 ++
Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in | 2 ++
INDUSTRIAL_FORMAT/Par_file | 2 ++
M2_UPPA/Par_file_M2_UPPA | 2 ++
Rayleigh_wave_no_crack/Par_file_Rayleigh_2D | 2 ++
Rayleigh_wave_with_crack/Par_file_Rayleigh_2D | 2 ++
Tape2007/Par_file_Tape2007_132rec_checker | 2 ++
Tape2007/Par_file_Tape2007_onerec | 2 ++
Tape2007_kernel/Par_file_Tape2007_onerec | 2 ++
Tromp2005/Par_file_Tromp2005 | 2 ++
Tromp2005/Par_file_Tromp2005_s100 | 2 ++
Tromp2005_kernel/Par_file_Tromp2005 | 2 ++
acoustic_poroelastic/Par_file_acoustic_poroelastic | 2 ++
attenuation/Par_file_attenuation_2D | 2 ++
canyon/Par_file_canyon | 2 ++
.../fluid_solid_external_mesh/Par_file_fluid_solid | 2 ++
.../Par_file_fluid_solid | 2 ++
.../Par_file_fluid_solid | 2 ++
global_Earth_ak135f/Par_file | 2 ++
init_plane/Par_file_Slave | 2 ++
init_plane/Par_file_Slave_for | 2 ++
init_plane/Par_file_Slave_kernel | 2 ++
noise_layered/model_0/Par_file_fair | 2 ++
noise_layered/model_0/Par_file_good | 2 ++
noise_layered/model_1/Par_file_best | 2 ++
noise_layered/model_1/Par_file_fair | 2 ++
noise_layered/model_1/Par_file_good | 2 ++
noise_layered/model_2/Par_file_fair | 2 ++
noise_layered/model_2/Par_file_good | 2 ++
noise_uniform/Par_file_noise_1 | 2 ++
noise_uniform/Par_file_noise_2 | 2 ++
noise_uniform/Par_file_noise_3 | 2 ++
salt_dome/CPML_homogeneous/Par_file | 2 ++
salt_dome/CPML_normal/Par_file | 2 ++
salt_dome/Stacey_homogeneous/Par_file | 4 ++-
salt_dome/Stacey_normal/Par_file | 2 ++
semi_infinite_homo/Par_file_elastic_2D | 4 ++-
.../plot_points_per_wavelength_histogram.gnu | 8 ++---
.../points_per_wavelength_histogram_S_in_solid.txt | 40 +++++++++++-----------
40 files changed, 102 insertions(+), 26 deletions(-)
diff --git a/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid b/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid
index 715264c..3eed3aa 100644
--- a/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid
+++ b/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid
@@ -118,6 +118,8 @@ read_external_mesh = .false.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .true.
diff --git a/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in b/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in
index 4c30457..6997c49 100644
--- a/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in
+++ b/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in
@@ -114,6 +114,8 @@ read_external_mesh = .true.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .true.
ADD_SPRING_TO_STACEY = .true.
diff --git a/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in b/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in
index 92322d4..3e7a06f 100644
--- a/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in
+++ b/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in
@@ -114,6 +114,8 @@ read_external_mesh = .true.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .true.
ADD_SPRING_TO_STACEY = .true.
diff --git a/INDUSTRIAL_FORMAT/Par_file b/INDUSTRIAL_FORMAT/Par_file
index 901fc72..f1a5dbb 100644
--- a/INDUSTRIAL_FORMAT/Par_file
+++ b/INDUSTRIAL_FORMAT/Par_file
@@ -112,6 +112,8 @@ read_external_mesh = .false.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .true.
diff --git a/M2_UPPA/Par_file_M2_UPPA b/M2_UPPA/Par_file_M2_UPPA
index 1da1715..caf3bf2 100644
--- a/M2_UPPA/Par_file_M2_UPPA
+++ b/M2_UPPA/Par_file_M2_UPPA
@@ -112,6 +112,8 @@ read_external_mesh = .false.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .true.
diff --git a/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D b/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D
index 725174c..f7ed58c 100644
--- a/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D
+++ b/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D
@@ -109,6 +109,8 @@ read_external_mesh = .false.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .true.
diff --git a/Rayleigh_wave_with_crack/Par_file_Rayleigh_2D b/Rayleigh_wave_with_crack/Par_file_Rayleigh_2D
index 03d62f4..52a1664 100644
--- a/Rayleigh_wave_with_crack/Par_file_Rayleigh_2D
+++ b/Rayleigh_wave_with_crack/Par_file_Rayleigh_2D
@@ -113,6 +113,8 @@ read_external_mesh = .false.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .true.
diff --git a/Tape2007/Par_file_Tape2007_132rec_checker b/Tape2007/Par_file_Tape2007_132rec_checker
index 873fbee..ed53ff2 100644
--- a/Tape2007/Par_file_Tape2007_132rec_checker
+++ b/Tape2007/Par_file_Tape2007_132rec_checker
@@ -895,6 +895,8 @@ read_external_mesh = .false.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .true.
diff --git a/Tape2007/Par_file_Tape2007_onerec b/Tape2007/Par_file_Tape2007_onerec
index 21d3515..1233a62 100644
--- a/Tape2007/Par_file_Tape2007_onerec
+++ b/Tape2007/Par_file_Tape2007_onerec
@@ -110,6 +110,8 @@ read_external_mesh = .false.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .true.
diff --git a/Tape2007_kernel/Par_file_Tape2007_onerec b/Tape2007_kernel/Par_file_Tape2007_onerec
index d6e3af1..7482ea8 100644
--- a/Tape2007_kernel/Par_file_Tape2007_onerec
+++ b/Tape2007_kernel/Par_file_Tape2007_onerec
@@ -110,6 +110,8 @@ read_external_mesh = .false.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .true.
diff --git a/Tromp2005/Par_file_Tromp2005 b/Tromp2005/Par_file_Tromp2005
index 8c5070d..321c977 100644
--- a/Tromp2005/Par_file_Tromp2005
+++ b/Tromp2005/Par_file_Tromp2005
@@ -109,6 +109,8 @@ read_external_mesh = .false.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .true.
diff --git a/Tromp2005/Par_file_Tromp2005_s100 b/Tromp2005/Par_file_Tromp2005_s100
index 7573579..da3453c 100644
--- a/Tromp2005/Par_file_Tromp2005_s100
+++ b/Tromp2005/Par_file_Tromp2005_s100
@@ -109,6 +109,8 @@ read_external_mesh = .false.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .true.
diff --git a/Tromp2005_kernel/Par_file_Tromp2005 b/Tromp2005_kernel/Par_file_Tromp2005
index f21ea92..55c23a2 100644
--- a/Tromp2005_kernel/Par_file_Tromp2005
+++ b/Tromp2005_kernel/Par_file_Tromp2005
@@ -109,6 +109,8 @@ read_external_mesh = .false.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .true.
diff --git a/acoustic_poroelastic/Par_file_acoustic_poroelastic b/acoustic_poroelastic/Par_file_acoustic_poroelastic
index 47ac9e1..44e290e 100644
--- a/acoustic_poroelastic/Par_file_acoustic_poroelastic
+++ b/acoustic_poroelastic/Par_file_acoustic_poroelastic
@@ -118,6 +118,8 @@ read_external_mesh = .false.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .true.
diff --git a/attenuation/Par_file_attenuation_2D b/attenuation/Par_file_attenuation_2D
index 867fe7b..6763db1 100644
--- a/attenuation/Par_file_attenuation_2D
+++ b/attenuation/Par_file_attenuation_2D
@@ -109,6 +109,8 @@ read_external_mesh = .false.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .true.
diff --git a/canyon/Par_file_canyon b/canyon/Par_file_canyon
index 7ad2660..548a2d6 100644
--- a/canyon/Par_file_canyon
+++ b/canyon/Par_file_canyon
@@ -109,6 +109,8 @@ read_external_mesh = .true.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .true.
ADD_SPRING_TO_STACEY = .true.
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 49a6e77..6de5538 100644
--- a/fluid_solid/fluid_solid_external_mesh/Par_file_fluid_solid
+++ b/fluid_solid/fluid_solid_external_mesh/Par_file_fluid_solid
@@ -114,6 +114,8 @@ read_external_mesh = .true.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .true.
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 4a8f892..bef6692 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
@@ -111,6 +111,8 @@ read_external_mesh = .false.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .true.
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 02b4edf..db8ed49 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
@@ -111,6 +111,8 @@ read_external_mesh = .false.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .true.
diff --git a/global_Earth_ak135f/Par_file b/global_Earth_ak135f/Par_file
index e99973f..db3d02f 100644
--- a/global_Earth_ak135f/Par_file
+++ b/global_Earth_ak135f/Par_file
@@ -115,6 +115,8 @@ read_external_mesh = .true.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .false.
diff --git a/init_plane/Par_file_Slave b/init_plane/Par_file_Slave
index e9bdc38..07600a3 100644
--- a/init_plane/Par_file_Slave
+++ b/init_plane/Par_file_Slave
@@ -110,6 +110,8 @@ read_external_mesh = .false.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .true.
diff --git a/init_plane/Par_file_Slave_for b/init_plane/Par_file_Slave_for
index a26860e..12754bb 100644
--- a/init_plane/Par_file_Slave_for
+++ b/init_plane/Par_file_Slave_for
@@ -110,6 +110,8 @@ read_external_mesh = .false.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .true.
diff --git a/init_plane/Par_file_Slave_kernel b/init_plane/Par_file_Slave_kernel
index 55c664b..2327fe3 100644
--- a/init_plane/Par_file_Slave_kernel
+++ b/init_plane/Par_file_Slave_kernel
@@ -110,6 +110,8 @@ read_external_mesh = .false.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .true.
diff --git a/noise_layered/model_0/Par_file_fair b/noise_layered/model_0/Par_file_fair
index 10eb253..5e15ea8 100644
--- a/noise_layered/model_0/Par_file_fair
+++ b/noise_layered/model_0/Par_file_fair
@@ -109,6 +109,8 @@ read_external_mesh = .false.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .true.
diff --git a/noise_layered/model_0/Par_file_good b/noise_layered/model_0/Par_file_good
index 17431d9..1cdbc22 100644
--- a/noise_layered/model_0/Par_file_good
+++ b/noise_layered/model_0/Par_file_good
@@ -109,6 +109,8 @@ read_external_mesh = .false.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .true.
diff --git a/noise_layered/model_1/Par_file_best b/noise_layered/model_1/Par_file_best
index 7cd300c..7187799 100644
--- a/noise_layered/model_1/Par_file_best
+++ b/noise_layered/model_1/Par_file_best
@@ -110,6 +110,8 @@ read_external_mesh = .false.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .true.
diff --git a/noise_layered/model_1/Par_file_fair b/noise_layered/model_1/Par_file_fair
index b7ed9c0..a49b785 100644
--- a/noise_layered/model_1/Par_file_fair
+++ b/noise_layered/model_1/Par_file_fair
@@ -110,6 +110,8 @@ read_external_mesh = .false.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .true.
diff --git a/noise_layered/model_1/Par_file_good b/noise_layered/model_1/Par_file_good
index 9dec092..e4076e1 100644
--- a/noise_layered/model_1/Par_file_good
+++ b/noise_layered/model_1/Par_file_good
@@ -110,6 +110,8 @@ read_external_mesh = .false.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .true.
diff --git a/noise_layered/model_2/Par_file_fair b/noise_layered/model_2/Par_file_fair
index 31cb4b1..85c2c25 100644
--- a/noise_layered/model_2/Par_file_fair
+++ b/noise_layered/model_2/Par_file_fair
@@ -111,6 +111,8 @@ read_external_mesh = .false.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .true.
diff --git a/noise_layered/model_2/Par_file_good b/noise_layered/model_2/Par_file_good
index cdbea27..e0bd998 100644
--- a/noise_layered/model_2/Par_file_good
+++ b/noise_layered/model_2/Par_file_good
@@ -111,6 +111,8 @@ read_external_mesh = .false.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .true.
diff --git a/noise_uniform/Par_file_noise_1 b/noise_uniform/Par_file_noise_1
index 06e287e..231dc86 100644
--- a/noise_uniform/Par_file_noise_1
+++ b/noise_uniform/Par_file_noise_1
@@ -109,6 +109,8 @@ read_external_mesh = .false.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .true.
diff --git a/noise_uniform/Par_file_noise_2 b/noise_uniform/Par_file_noise_2
index a9f8fb0..0b90559 100644
--- a/noise_uniform/Par_file_noise_2
+++ b/noise_uniform/Par_file_noise_2
@@ -109,6 +109,8 @@ read_external_mesh = .false.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .true.
diff --git a/noise_uniform/Par_file_noise_3 b/noise_uniform/Par_file_noise_3
index 0fd6a5e..db186ee 100644
--- a/noise_uniform/Par_file_noise_3
+++ b/noise_uniform/Par_file_noise_3
@@ -109,6 +109,8 @@ read_external_mesh = .false.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .true.
diff --git a/salt_dome/CPML_homogeneous/Par_file b/salt_dome/CPML_homogeneous/Par_file
index 362c2b2..479039a 100644
--- a/salt_dome/CPML_homogeneous/Par_file
+++ b/salt_dome/CPML_homogeneous/Par_file
@@ -116,6 +116,8 @@ read_external_mesh = .true.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .true.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .false.
diff --git a/salt_dome/CPML_normal/Par_file b/salt_dome/CPML_normal/Par_file
index dc8195f..e9cd45f 100644
--- a/salt_dome/CPML_normal/Par_file
+++ b/salt_dome/CPML_normal/Par_file
@@ -123,6 +123,8 @@ read_external_mesh = .true.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .true.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .false.
diff --git a/salt_dome/Stacey_homogeneous/Par_file b/salt_dome/Stacey_homogeneous/Par_file
index a58dfd7..8120728 100644
--- a/salt_dome/Stacey_homogeneous/Par_file
+++ b/salt_dome/Stacey_homogeneous/Par_file
@@ -24,7 +24,7 @@ p_sv = .true. # set the type of calculation (
# time step parameters
nt = 10000 # total number of time steps
-deltat = 4.25e-4 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
+deltat = 4.25e-4 # duration of a time step (see section "How to choose the time step" of the manual for how to do this) (for the choice of deltat please refer to section 4.5 of the user manual)
USER_T0 = 0.0d0 # 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 RK4 4th-order 4-stage Runge-Kutta
@@ -116,6 +116,8 @@ read_external_mesh = .true.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .true.
ADD_SPRING_TO_STACEY = .true.
diff --git a/salt_dome/Stacey_normal/Par_file b/salt_dome/Stacey_normal/Par_file
index 7cd33b9..9ead0ff 100644
--- a/salt_dome/Stacey_normal/Par_file
+++ b/salt_dome/Stacey_normal/Par_file
@@ -123,6 +123,8 @@ read_external_mesh = .true.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .true.
ADD_SPRING_TO_STACEY = .false.
diff --git a/semi_infinite_homo/Par_file_elastic_2D b/semi_infinite_homo/Par_file_elastic_2D
index d8193ca..019b76e 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 = 15000 # total number of time steps
+nt = 1500 # total number of time steps
deltat = 0.8e-3 # duration of a time step (see section "How to choose the time step" of the manual for how to do this)
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
@@ -109,6 +109,8 @@ read_external_mesh = .false.
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .true.
NELEM_PML_THICKNESS = 3
+ROTATE_PML_ACTIVATE = .false.
+ROTATE_PML_ANGLE = 30.
STACEY_ABSORBING_CONDITIONS = .false.
ADD_SPRING_TO_STACEY = .false.
diff --git a/semi_infinite_homo/plot_points_per_wavelength_histogram.gnu b/semi_infinite_homo/plot_points_per_wavelength_histogram.gnu
index a108347..a382333 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_P_in_fluid.gif"
+ #set output "points_per_wavelength_histogram_S_in_solid.gif"
- set boxwidth 1.5000000E-02
- set xlabel "Range of min number of points per P wavelength in fluid"
+ set boxwidth 8.66025034E-03
+ set xlabel "Range of min number of points per S wavelength in solid"
set ylabel "Percentage of elements (%)"
- plot "points_per_wavelength_histogram_P_in_fluid.txt" with boxes
+ plot "points_per_wavelength_histogram_S_in_solid.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 2e24674..6906e54 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.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
+ 8.5779781 0.0000000
+ 8.5866375 0.0000000
+ 8.5952988 0.0000000
+ 8.6039581 0.0000000
+ 8.6126194 0.0000000
+ 8.6212788 0.0000000
+ 8.6299400 0.0000000
+ 8.6385994 0.0000000
+ 8.6472588 0.0000000
+ 8.6559200 14.000000
+ 8.6645794 86.000000
+ 8.6732407 0.0000000
+ 8.6819000 0.0000000
+ 8.6905613 0.0000000
+ 8.6992207 0.0000000
+ 8.7078819 0.0000000
+ 8.7165413 0.0000000
+ 8.7252026 0.0000000
+ 8.7338619 0.0000000
+ 8.7425232 0.0000000
More information about the CIG-COMMITS
mailing list