[cig-commits] [commit] devel, master: adds C-PML boundary conditions in CPU code for elastic and acoustic domains. BEWARE!!! runs ok but no fully tested -> desactivated by default in all examples. (a324b48)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Wed Jun 18 16:50:51 PDT 2014
Repository : https://github.com/geodynamics/specfem3d
On branches: devel,master
Link : https://github.com/geodynamics/specfem3d/compare/6026e367984905ab133865f62fa6293b343759b9...47f703851338234f96397e7da9fbff63d8178b8a
>---------------------------------------------------------------
commit a324b488fb65c1673da2fb7e441effb2d98febe3
Author: Joseph Charles <jcharles at lma.cnrs-mrs.fr>
Date: Sat Jan 12 21:24:23 2013 +0000
adds C-PML boundary conditions in CPU code for elastic and acoustic domains. BEWARE!!! runs ok but no fully tested -> desactivated by default in all examples.
>---------------------------------------------------------------
a324b488fb65c1673da2fb7e441effb2d98febe3
BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC/DATA/Par_file | 25 ++++++++++++++++++----
BENCHMARK_CLAERBOUT_ADJOINT/ELASTIC/DATA/Par_file | 25 ++++++++++++++++++----
Mount_StHelens/DATA/Par_file | 23 +++++++++++++++++---
homogeneous_halfspace_HEX27/DATA/Par_file | 20 +++++++++++++++--
homogeneous_halfspace_HEX8/DATA/Par_file | 20 +++++++++++++++--
homogeneous_poroelastic/DATA/Par_file | 24 +++++++++++++++++----
layered_halfspace/DATA/Par_file | 23 +++++++++++++++++---
meshfem3D_examples/many_interfaces/DATA/Par_file | 23 +++++++++++++++++---
meshfem3D_examples/simple_model/DATA/Par_file | 23 +++++++++++++++++---
meshfem3D_examples/socal1D/DATA/Par_file | 23 +++++++++++++++++---
splay_faults/DATA/Par_file | 20 +++++++++++++++--
tomographic_model/DATA/Par_file | 23 +++++++++++++++++---
tpv102/DATA/Par_file | 20 +++++++++++++++--
tpv103/DATA/Par_file | 20 +++++++++++++++--
tpv15/DATA/Par_file | 20 +++++++++++++++--
tpv16/DATA/Par_file | 20 +++++++++++++++--
tpv5/DATA/Par_file | 20 +++++++++++++++--
waterlayered_halfspace/DATA/Par_file | 23 +++++++++++++++++---
18 files changed, 346 insertions(+), 49 deletions(-)
diff --git a/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC/DATA/Par_file b/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC/DATA/Par_file
index 3c91e05..444956a 100644
--- a/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC/DATA/Par_file
+++ b/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC/DATA/Par_file
@@ -54,6 +54,23 @@ ABSORBING_CONDITIONS = .true.
# absorbing top surface (defined in mesh as 'free_surface_file')
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
+# C-PML boundary conditions for a regional simulation
+PML_CONDITIONS = .false.
+
+# C-PML top surface
+PML_INSTEAD_OF_FREE_SURFACE = .false.
+
+# C-PML thickness min/max
+PML_WIDTH_MIN = 10000.0d0
+PML_WIDTH_MAX = 25000.0d0
+
+# C-PML dominant frequency
+f0_FOR_PML = 12.7
+
+# parameters used to rotate C-PML boundary conditions by a given angle (not completed yet)
+# ROTATE_PML_ACTIVATE = .false.
+# ROTATE_PML_ANGLE = 0.
+
# save AVS or OpenDX movies
# MOVIE_TYPE = 1 to show the top surface
# MOVIE_TYPE = 2 to show all the external faces of the mesh
@@ -70,7 +87,7 @@ HDUR_MOVIE = 0.0
SAVE_MESH_FILES = .true.
# path to store the local database file on each node
-LOCAL_PATH = ../OUTPUT_FILES/DATABASES_MPI/
+LOCAL_PATH = ../OUTPUT_FILES/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 100
@@ -79,10 +96,10 @@ NTSTEP_BETWEEN_OUTPUT_INFO = 100
NTSTEP_BETWEEN_OUTPUT_SEISMOS = 10000
# interval in time steps for reading adjoint traces
-NTSTEP_BETWEEN_READ_ADJSRC = 0 # 0 = read the whole adjoint sources at the same time
+# 0 = read the whole adjoint sources at the same time
+NTSTEP_BETWEEN_READ_ADJSRC = 0
-# use a (tilted) FORCESOLUTION force point source (or several) located exactly at a grid point
-# instead of a CMTSOLUTION moment-tensor source.
+# use a (tilted) FORCESOLUTION force point source (or several) instead of a CMTSOLUTION moment-tensor source.
# This can be useful e.g. for oil industry foothills simulations or asteroid simulations
# in which the source is a vertical force, normal force, inclined force, impact etc.
# If this flag is turned on, the FORCESOLUTION file must be edited by precising:
diff --git a/BENCHMARK_CLAERBOUT_ADJOINT/ELASTIC/DATA/Par_file b/BENCHMARK_CLAERBOUT_ADJOINT/ELASTIC/DATA/Par_file
index 3b90548..0c0a648 100644
--- a/BENCHMARK_CLAERBOUT_ADJOINT/ELASTIC/DATA/Par_file
+++ b/BENCHMARK_CLAERBOUT_ADJOINT/ELASTIC/DATA/Par_file
@@ -54,6 +54,23 @@ ABSORBING_CONDITIONS = .true.
# absorbing top surface (defined in mesh as 'free_surface_file')
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
+# C-PML boundary conditions for a regional simulation
+PML_CONDITIONS = .false.
+
+# C-PML top surface
+PML_INSTEAD_OF_FREE_SURFACE = .false.
+
+# C-PML thickness min/max
+PML_WIDTH_MIN = 10000.0d0
+PML_WIDTH_MAX = 25000.0d0
+
+# C-PML dominant frequency
+f0_FOR_PML = 12.7
+
+# parameters used to rotate C-PML boundary conditions by a given angle (not completed yet)
+# ROTATE_PML_ACTIVATE = .false.
+# ROTATE_PML_ANGLE = 0.
+
# save AVS or OpenDX movies
# MOVIE_TYPE = 1 to show the top surface
# MOVIE_TYPE = 2 to show all the external faces of the mesh
@@ -70,7 +87,7 @@ HDUR_MOVIE = 0.0
SAVE_MESH_FILES = .true.
# path to store the local database file on each node
-LOCAL_PATH = ../OUTPUT_FILES/DATABASES_MPI/
+LOCAL_PATH = ../OUTPUT_FILES/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 100
@@ -79,10 +96,10 @@ NTSTEP_BETWEEN_OUTPUT_INFO = 100
NTSTEP_BETWEEN_OUTPUT_SEISMOS = 10000
# interval in time steps for reading adjoint traces
-NTSTEP_BETWEEN_READ_ADJSRC = 0 # 0 = read the whole adjoint sources at the same time
+# 0 = read the whole adjoint sources at the same time
+NTSTEP_BETWEEN_READ_ADJSRC = 0
-# use a (tilted) FORCESOLUTION force point source (or several) located exactly at a grid point
-# instead of a CMTSOLUTION moment-tensor source.
+# use a (tilted) FORCESOLUTION force point source (or several) instead of a CMTSOLUTION moment-tensor source.
# This can be useful e.g. for oil industry foothills simulations or asteroid simulations
# in which the source is a vertical force, normal force, inclined force, impact etc.
# If this flag is turned on, the FORCESOLUTION file must be edited by precising:
diff --git a/Mount_StHelens/DATA/Par_file b/Mount_StHelens/DATA/Par_file
index 456c998..7be25fa 100644
--- a/Mount_StHelens/DATA/Par_file
+++ b/Mount_StHelens/DATA/Par_file
@@ -54,6 +54,23 @@ ABSORBING_CONDITIONS = .true.
# absorbing top surface (defined in mesh as 'free_surface_file')
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
+# C-PML boundary conditions for a regional simulation
+PML_CONDITIONS = .false.
+
+# C-PML top surface
+PML_INSTEAD_OF_FREE_SURFACE = .false.
+
+# C-PML thickness min/max
+PML_WIDTH_MIN = 10000.0d0
+PML_WIDTH_MAX = 25000.0d0
+
+# C-PML dominant frequency
+f0_FOR_PML = 12.7
+
+# parameters used to rotate C-PML boundary conditions by a given angle (not completed yet)
+# ROTATE_PML_ACTIVATE = .false.
+# ROTATE_PML_ANGLE = 0.
+
# save AVS or OpenDX movies
# MOVIE_TYPE = 1 to show the top surface
# MOVIE_TYPE = 2 to show all the external faces of the mesh
@@ -79,10 +96,10 @@ NTSTEP_BETWEEN_OUTPUT_INFO = 500
NTSTEP_BETWEEN_OUTPUT_SEISMOS = 10000
# interval in time steps for reading adjoint traces
-NTSTEP_BETWEEN_READ_ADJSRC = 0 # 0 = read the whole adjoint sources at the same time
+# 0 = read the whole adjoint sources at the same time
+NTSTEP_BETWEEN_READ_ADJSRC = 0
-# use a (tilted) FORCESOLUTION force point source (or several) located exactly at a grid point
-# instead of a CMTSOLUTION moment-tensor source.
+# use a (tilted) FORCESOLUTION force point source (or several) instead of a CMTSOLUTION moment-tensor source.
# This can be useful e.g. for oil industry foothills simulations or asteroid simulations
# in which the source is a vertical force, normal force, inclined force, impact etc.
# If this flag is turned on, the FORCESOLUTION file must be edited by precising:
diff --git a/homogeneous_halfspace_HEX27/DATA/Par_file b/homogeneous_halfspace_HEX27/DATA/Par_file
index aa5442e..a2ff73e 100644
--- a/homogeneous_halfspace_HEX27/DATA/Par_file
+++ b/homogeneous_halfspace_HEX27/DATA/Par_file
@@ -54,6 +54,23 @@ ABSORBING_CONDITIONS = .false.
# absorbing top surface (defined in mesh as 'free_surface_file')
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
+# C-PML boundary conditions for a regional simulation
+PML_CONDITIONS = .false.
+
+# C-PML top surface
+PML_INSTEAD_OF_FREE_SURFACE = .false.
+
+# C-PML thickness min/max
+PML_WIDTH_MIN = 10000.0d0
+PML_WIDTH_MAX = 25000.0d0
+
+# C-PML dominant frequency
+f0_FOR_PML = 12.7
+
+# parameters used to rotate C-PML boundary conditions by a given angle (not completed yet)
+# ROTATE_PML_ACTIVATE = .false.
+# ROTATE_PML_ANGLE = 0.
+
# save AVS or OpenDX movies
# MOVIE_TYPE = 1 to show the top surface
# MOVIE_TYPE = 2 to show all the external faces of the mesh
@@ -82,8 +99,7 @@ NTSTEP_BETWEEN_OUTPUT_SEISMOS = 10000
# 0 = read the whole adjoint sources at the same time
NTSTEP_BETWEEN_READ_ADJSRC = 0
-# use a (tilted) FORCESOLUTION force point source (or several) located exactly at a grid point
-# instead of a CMTSOLUTION moment-tensor source.
+# use a (tilted) FORCESOLUTION force point source (or several) instead of a CMTSOLUTION moment-tensor source.
# This can be useful e.g. for oil industry foothills simulations or asteroid simulations
# in which the source is a vertical force, normal force, inclined force, impact etc.
# If this flag is turned on, the FORCESOLUTION file must be edited by precising:
diff --git a/homogeneous_halfspace_HEX8/DATA/Par_file b/homogeneous_halfspace_HEX8/DATA/Par_file
index d2478ea..943b0cd 100644
--- a/homogeneous_halfspace_HEX8/DATA/Par_file
+++ b/homogeneous_halfspace_HEX8/DATA/Par_file
@@ -54,6 +54,23 @@ ABSORBING_CONDITIONS = .false.
# absorbing top surface (defined in mesh as 'free_surface_file')
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
+# C-PML boundary conditions for a regional simulation
+PML_CONDITIONS = .false.
+
+# C-PML top surface
+PML_INSTEAD_OF_FREE_SURFACE = .false.
+
+# C-PML thickness min/max
+PML_WIDTH_MIN = 10000.0d0
+PML_WIDTH_MAX = 25000.0d0
+
+# C-PML dominant frequency
+f0_FOR_PML = 12.7
+
+# parameters used to rotate C-PML boundary conditions by a given angle (not completed yet)
+# ROTATE_PML_ACTIVATE = .false.
+# ROTATE_PML_ANGLE = 0.
+
# save AVS or OpenDX movies
# MOVIE_TYPE = 1 to show the top surface
# MOVIE_TYPE = 2 to show all the external faces of the mesh
@@ -82,8 +99,7 @@ NTSTEP_BETWEEN_OUTPUT_SEISMOS = 10000
# 0 = read the whole adjoint sources at the same time
NTSTEP_BETWEEN_READ_ADJSRC = 0
-# use a (tilted) FORCESOLUTION force point source (or several) located exactly at a grid point
-# instead of a CMTSOLUTION moment-tensor source.
+# use a (tilted) FORCESOLUTION force point source (or several) instead of a CMTSOLUTION moment-tensor source.
# This can be useful e.g. for oil industry foothills simulations or asteroid simulations
# in which the source is a vertical force, normal force, inclined force, impact etc.
# If this flag is turned on, the FORCESOLUTION file must be edited by precising:
diff --git a/homogeneous_poroelastic/DATA/Par_file b/homogeneous_poroelastic/DATA/Par_file
index 59c6a98..22b663c 100644
--- a/homogeneous_poroelastic/DATA/Par_file
+++ b/homogeneous_poroelastic/DATA/Par_file
@@ -1,4 +1,4 @@
-*# simulation input parameters
+# simulation input parameters
#
# forward or adjoint simulation
# 1 = forward, 2 = adjoint, 3 = both simultaneously
@@ -54,6 +54,23 @@ ABSORBING_CONDITIONS = .false.
# absorbing top surface (defined in mesh as 'free_surface_file')
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
+# C-PML boundary conditions for a regional simulation
+PML_CONDITIONS = .false.
+
+# C-PML top surface
+PML_INSTEAD_OF_FREE_SURFACE = .false.
+
+# C-PML thickness min/max
+PML_WIDTH_MIN = 10000.0d0
+PML_WIDTH_MAX = 25000.0d0
+
+# C-PML dominant frequency
+f0_FOR_PML = 12.7
+
+# parameters used to rotate C-PML boundary conditions by a given angle (not completed yet)
+# ROTATE_PML_ACTIVATE = .false.
+# ROTATE_PML_ANGLE = 0.
+
# save AVS or OpenDX movies
# MOVIE_TYPE = 1 to show the top surface
# MOVIE_TYPE = 2 to show all the external faces of the mesh
@@ -79,11 +96,10 @@ NTSTEP_BETWEEN_OUTPUT_INFO = 250
NTSTEP_BETWEEN_OUTPUT_SEISMOS = 10000
# interval in time steps for reading adjoint traces
-# # 0 = read the whole adjoint sources at the same time
+# 0 = read the whole adjoint sources at the same time
NTSTEP_BETWEEN_READ_ADJSRC = 0
-# use a (tilted) FORCESOLUTION force point source (or several) located exactly at a grid point
-# instead of a CMTSOLUTION moment-tensor source.
+# use a (tilted) FORCESOLUTION force point source (or several) instead of a CMTSOLUTION moment-tensor source.
# This can be useful e.g. for oil industry foothills simulations or asteroid simulations
# in which the source is a vertical force, normal force, inclined force, impact etc.
# If this flag is turned on, the FORCESOLUTION file must be edited by precising:
diff --git a/layered_halfspace/DATA/Par_file b/layered_halfspace/DATA/Par_file
index 935428d..2a882f5 100644
--- a/layered_halfspace/DATA/Par_file
+++ b/layered_halfspace/DATA/Par_file
@@ -54,6 +54,23 @@ ABSORBING_CONDITIONS = .true.
# absorbing top surface (defined in mesh as 'free_surface_file')
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
+# C-PML boundary conditions for a regional simulation
+PML_CONDITIONS = .false.
+
+# C-PML top surface
+PML_INSTEAD_OF_FREE_SURFACE = .false.
+
+# C-PML thickness min/max
+PML_WIDTH_MIN = 10000.0d0
+PML_WIDTH_MAX = 25000.0d0
+
+# C-PML dominant frequency
+f0_FOR_PML = 12.7
+
+# parameters used to rotate C-PML boundary conditions by a given angle (not completed yet)
+# ROTATE_PML_ACTIVATE = .false.
+# ROTATE_PML_ANGLE = 0.
+
# save AVS or OpenDX movies
# MOVIE_TYPE = 1 to show the top surface
# MOVIE_TYPE = 2 to show all the external faces of the mesh
@@ -79,10 +96,10 @@ NTSTEP_BETWEEN_OUTPUT_INFO = 500
NTSTEP_BETWEEN_OUTPUT_SEISMOS = 10000
# interval in time steps for reading adjoint traces
-NTSTEP_BETWEEN_READ_ADJSRC = 0 # 0 = read the whole adjoint sources at the same time
+# 0 = read the whole adjoint sources at the same time
+NTSTEP_BETWEEN_READ_ADJSRC = 0
-# use a (tilted) FORCESOLUTION force point source (or several) located exactly at a grid point
-# instead of a CMTSOLUTION moment-tensor source.
+# use a (tilted) FORCESOLUTION force point source (or several) instead of a CMTSOLUTION moment-tensor source.
# This can be useful e.g. for oil industry foothills simulations or asteroid simulations
# in which the source is a vertical force, normal force, inclined force, impact etc.
# If this flag is turned on, the FORCESOLUTION file must be edited by precising:
diff --git a/meshfem3D_examples/many_interfaces/DATA/Par_file b/meshfem3D_examples/many_interfaces/DATA/Par_file
index 367fc34..f2cf55b 100644
--- a/meshfem3D_examples/many_interfaces/DATA/Par_file
+++ b/meshfem3D_examples/many_interfaces/DATA/Par_file
@@ -54,6 +54,23 @@ ABSORBING_CONDITIONS = .true.
# absorbing top surface (defined in mesh as 'free_surface_file')
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
+# C-PML boundary conditions for a regional simulation
+PML_CONDITIONS = .false.
+
+# C-PML top surface
+PML_INSTEAD_OF_FREE_SURFACE = .false.
+
+# C-PML thickness min/max
+PML_WIDTH_MIN = 10000.0d0
+PML_WIDTH_MAX = 25000.0d0
+
+# C-PML dominant frequency
+f0_FOR_PML = 12.7
+
+# parameters used to rotate C-PML boundary conditions by a given angle (not completed yet)
+# ROTATE_PML_ACTIVATE = .false.
+# ROTATE_PML_ANGLE = 0.
+
# save AVS or OpenDX movies
# MOVIE_TYPE = 1 to show the top surface
# MOVIE_TYPE = 2 to show all the external faces of the mesh
@@ -79,10 +96,10 @@ NTSTEP_BETWEEN_OUTPUT_INFO = 500
NTSTEP_BETWEEN_OUTPUT_SEISMOS = 10000
# interval in time steps for reading adjoint traces
-NTSTEP_BETWEEN_READ_ADJSRC = 0 # 0 = read the whole adjoint sources at the same time
+# 0 = read the whole adjoint sources at the same time
+NTSTEP_BETWEEN_READ_ADJSRC = 0
-# use a (tilted) FORCESOLUTION force point source (or several) located exactly at a grid point
-# instead of a CMTSOLUTION moment-tensor source.
+# use a (tilted) FORCESOLUTION force point source (or several) instead of a CMTSOLUTION moment-tensor source.
# This can be useful e.g. for oil industry foothills simulations or asteroid simulations
# in which the source is a vertical force, normal force, inclined force, impact etc.
# If this flag is turned on, the FORCESOLUTION file must be edited by precising:
diff --git a/meshfem3D_examples/simple_model/DATA/Par_file b/meshfem3D_examples/simple_model/DATA/Par_file
index c2bcb8f..e14db13 100644
--- a/meshfem3D_examples/simple_model/DATA/Par_file
+++ b/meshfem3D_examples/simple_model/DATA/Par_file
@@ -54,6 +54,23 @@ ABSORBING_CONDITIONS = .false.
# absorbing top surface (defined in mesh as 'free_surface_file')
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
+# C-PML boundary conditions for a regional simulation
+PML_CONDITIONS = .false.
+
+# C-PML top surface
+PML_INSTEAD_OF_FREE_SURFACE = .false.
+
+# C-PML thickness min/max
+PML_WIDTH_MIN = 10000.0d0
+PML_WIDTH_MAX = 25000.0d0
+
+# C-PML dominant frequency
+f0_FOR_PML = 12.7
+
+# parameters used to rotate C-PML boundary conditions by a given angle (not completed yet)
+# ROTATE_PML_ACTIVATE = .false.
+# ROTATE_PML_ANGLE = 0.
+
# save AVS or OpenDX movies
# MOVIE_TYPE = 1 to show the top surface
# MOVIE_TYPE = 2 to show all the external faces of the mesh
@@ -79,10 +96,10 @@ NTSTEP_BETWEEN_OUTPUT_INFO = 100
NTSTEP_BETWEEN_OUTPUT_SEISMOS = 10000
# interval in time steps for reading adjoint traces
-NTSTEP_BETWEEN_READ_ADJSRC = 0 # 0 = read the whole adjoint sources at the same time
+# 0 = read the whole adjoint sources at the same time
+NTSTEP_BETWEEN_READ_ADJSRC = 0
-# use a (tilted) FORCESOLUTION force point source (or several) located exactly at a grid point
-# instead of a CMTSOLUTION moment-tensor source.
+# use a (tilted) FORCESOLUTION force point source (or several) instead of a CMTSOLUTION moment-tensor source.
# This can be useful e.g. for oil industry foothills simulations or asteroid simulations
# in which the source is a vertical force, normal force, inclined force, impact etc.
# If this flag is turned on, the FORCESOLUTION file must be edited by precising:
diff --git a/meshfem3D_examples/socal1D/DATA/Par_file b/meshfem3D_examples/socal1D/DATA/Par_file
index 196796e..581fe54 100644
--- a/meshfem3D_examples/socal1D/DATA/Par_file
+++ b/meshfem3D_examples/socal1D/DATA/Par_file
@@ -54,6 +54,23 @@ ABSORBING_CONDITIONS = .true.
# absorbing top surface (defined in mesh as 'free_surface_file')
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
+# C-PML boundary conditions for a regional simulation
+PML_CONDITIONS = .false.
+
+# C-PML top surface
+PML_INSTEAD_OF_FREE_SURFACE = .false.
+
+# C-PML thickness min/max
+PML_WIDTH_MIN = 10000.0d0
+PML_WIDTH_MAX = 25000.0d0
+
+# C-PML dominant frequency
+f0_FOR_PML = 12.7
+
+# parameters used to rotate C-PML boundary conditions by a given angle (not completed yet)
+# ROTATE_PML_ACTIVATE = .false.
+# ROTATE_PML_ANGLE = 0.
+
# save AVS or OpenDX movies
# MOVIE_TYPE = 1 to show the top surface
# MOVIE_TYPE = 2 to show all the external faces of the mesh
@@ -79,10 +96,10 @@ NTSTEP_BETWEEN_OUTPUT_INFO = 100
NTSTEP_BETWEEN_OUTPUT_SEISMOS = 10000
# interval in time steps for reading adjoint traces
-NTSTEP_BETWEEN_READ_ADJSRC = 0 # 0 = read the whole adjoint sources at the same time
+# 0 = read the whole adjoint sources at the same time
+NTSTEP_BETWEEN_READ_ADJSRC = 0
-# use a (tilted) FORCESOLUTION force point source (or several) located exactly at a grid point
-# instead of a CMTSOLUTION moment-tensor source.
+# use a (tilted) FORCESOLUTION force point source (or several) instead of a CMTSOLUTION moment-tensor source.
# This can be useful e.g. for oil industry foothills simulations or asteroid simulations
# in which the source is a vertical force, normal force, inclined force, impact etc.
# If this flag is turned on, the FORCESOLUTION file must be edited by precising:
diff --git a/splay_faults/DATA/Par_file b/splay_faults/DATA/Par_file
index 4c4c946..6e3117a 100644
--- a/splay_faults/DATA/Par_file
+++ b/splay_faults/DATA/Par_file
@@ -54,6 +54,23 @@ ABSORBING_CONDITIONS = .true.
# absorbing top surface (defined in mesh as 'free_surface_file')
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
+# C-PML boundary conditions for a regional simulation
+PML_CONDITIONS = .false.
+
+# C-PML top surface
+PML_INSTEAD_OF_FREE_SURFACE = .false.
+
+# C-PML thickness min/max
+PML_WIDTH_MIN = 10000.0d0
+PML_WIDTH_MAX = 25000.0d0
+
+# C-PML dominant frequency
+f0_FOR_PML = 12.7
+
+# parameters used to rotate C-PML boundary conditions by a given angle (not completed yet)
+# ROTATE_PML_ACTIVATE = .false.
+# ROTATE_PML_ANGLE = 0.
+
# save AVS or OpenDX movies
# MOVIE_TYPE = 1 to show the top surface
# MOVIE_TYPE = 2 to show all the external faces of the mesh
@@ -82,8 +99,7 @@ NTSTEP_BETWEEN_OUTPUT_SEISMOS = 10000
# 0 = read the whole adjoint sources at the same time
NTSTEP_BETWEEN_READ_ADJSRC = 0
-# use a (tilted) FORCESOLUTION force point source (or several) located exactly at a grid point
-# instead of a CMTSOLUTION moment-tensor source.
+# use a (tilted) FORCESOLUTION force point source (or several) instead of a CMTSOLUTION moment-tensor source.
# This can be useful e.g. for oil industry foothills simulations or asteroid simulations
# in which the source is a vertical force, normal force, inclined force, impact etc.
# If this flag is turned on, the FORCESOLUTION file must be edited by precising:
diff --git a/tomographic_model/DATA/Par_file b/tomographic_model/DATA/Par_file
index 29acdfb..2d27ec6 100644
--- a/tomographic_model/DATA/Par_file
+++ b/tomographic_model/DATA/Par_file
@@ -54,6 +54,23 @@ ABSORBING_CONDITIONS = .true.
# absorbing top surface (defined in mesh as 'free_surface_file')
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
+# C-PML boundary conditions for a regional simulation
+PML_CONDITIONS = .false.
+
+# C-PML top surface
+PML_INSTEAD_OF_FREE_SURFACE = .false.
+
+# C-PML thickness min/max
+PML_WIDTH_MIN = 10000.0d0
+PML_WIDTH_MAX = 25000.0d0
+
+# C-PML dominant frequency
+f0_FOR_PML = 12.7
+
+# parameters used to rotate C-PML boundary conditions by a given angle (not completed yet)
+# ROTATE_PML_ACTIVATE = .false.
+# ROTATE_PML_ANGLE = 0.
+
# save AVS or OpenDX movies
# MOVIE_TYPE = 1 to show the top surface
# MOVIE_TYPE = 2 to show all the external faces of the mesh
@@ -79,10 +96,10 @@ NTSTEP_BETWEEN_OUTPUT_INFO = 500
NTSTEP_BETWEEN_OUTPUT_SEISMOS = 10000
# interval in time steps for reading adjoint traces
-NTSTEP_BETWEEN_READ_ADJSRC = 0 # 0 = read the whole adjoint sources at the same time
+# 0 = read the whole adjoint sources at the same time
+NTSTEP_BETWEEN_READ_ADJSRC = 0
-# use a (tilted) FORCESOLUTION force point source (or several) located exactly at a grid point
-# instead of a CMTSOLUTION moment-tensor source.
+# use a (tilted) FORCESOLUTION force point source (or several) instead of a CMTSOLUTION moment-tensor source.
# This can be useful e.g. for oil industry foothills simulations or asteroid simulations
# in which the source is a vertical force, normal force, inclined force, impact etc.
# If this flag is turned on, the FORCESOLUTION file must be edited by precising:
diff --git a/tpv102/DATA/Par_file b/tpv102/DATA/Par_file
index 64ca619..7fe6985 100644
--- a/tpv102/DATA/Par_file
+++ b/tpv102/DATA/Par_file
@@ -54,6 +54,23 @@ ABSORBING_CONDITIONS = .true.
# absorbing top surface (defined in mesh as 'free_surface_file')
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
+# C-PML boundary conditions for a regional simulation
+PML_CONDITIONS = .false.
+
+# C-PML top surface
+PML_INSTEAD_OF_FREE_SURFACE = .false.
+
+# C-PML thickness min/max
+PML_WIDTH_MIN = 10000.0d0
+PML_WIDTH_MAX = 25000.0d0
+
+# C-PML dominant frequency
+f0_FOR_PML = 12.7
+
+# parameters used to rotate C-PML boundary conditions by a given angle (not completed yet)
+# ROTATE_PML_ACTIVATE = .false.
+# ROTATE_PML_ANGLE = 0.
+
# save AVS or OpenDX movies
# MOVIE_TYPE = 1 to show the top surface
# MOVIE_TYPE = 2 to show all the external faces of the mesh
@@ -82,8 +99,7 @@ NTSTEP_BETWEEN_OUTPUT_SEISMOS = 10000
# 0 = read the whole adjoint sources at the same time
NTSTEP_BETWEEN_READ_ADJSRC = 0
-# use a (tilted) FORCESOLUTION force point source (or several) located exactly at a grid point
-# instead of a CMTSOLUTION moment-tensor source.
+# use a (tilted) FORCESOLUTION force point source (or several) instead of a CMTSOLUTION moment-tensor source.
# This can be useful e.g. for oil industry foothills simulations or asteroid simulations
# in which the source is a vertical force, normal force, inclined force, impact etc.
# If this flag is turned on, the FORCESOLUTION file must be edited by precising:
diff --git a/tpv103/DATA/Par_file b/tpv103/DATA/Par_file
index 64ca619..7fe6985 100644
--- a/tpv103/DATA/Par_file
+++ b/tpv103/DATA/Par_file
@@ -54,6 +54,23 @@ ABSORBING_CONDITIONS = .true.
# absorbing top surface (defined in mesh as 'free_surface_file')
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
+# C-PML boundary conditions for a regional simulation
+PML_CONDITIONS = .false.
+
+# C-PML top surface
+PML_INSTEAD_OF_FREE_SURFACE = .false.
+
+# C-PML thickness min/max
+PML_WIDTH_MIN = 10000.0d0
+PML_WIDTH_MAX = 25000.0d0
+
+# C-PML dominant frequency
+f0_FOR_PML = 12.7
+
+# parameters used to rotate C-PML boundary conditions by a given angle (not completed yet)
+# ROTATE_PML_ACTIVATE = .false.
+# ROTATE_PML_ANGLE = 0.
+
# save AVS or OpenDX movies
# MOVIE_TYPE = 1 to show the top surface
# MOVIE_TYPE = 2 to show all the external faces of the mesh
@@ -82,8 +99,7 @@ NTSTEP_BETWEEN_OUTPUT_SEISMOS = 10000
# 0 = read the whole adjoint sources at the same time
NTSTEP_BETWEEN_READ_ADJSRC = 0
-# use a (tilted) FORCESOLUTION force point source (or several) located exactly at a grid point
-# instead of a CMTSOLUTION moment-tensor source.
+# use a (tilted) FORCESOLUTION force point source (or several) instead of a CMTSOLUTION moment-tensor source.
# This can be useful e.g. for oil industry foothills simulations or asteroid simulations
# in which the source is a vertical force, normal force, inclined force, impact etc.
# If this flag is turned on, the FORCESOLUTION file must be edited by precising:
diff --git a/tpv15/DATA/Par_file b/tpv15/DATA/Par_file
index 0182a3f..8cb3837 100644
--- a/tpv15/DATA/Par_file
+++ b/tpv15/DATA/Par_file
@@ -54,6 +54,23 @@ ABSORBING_CONDITIONS = .true.
# absorbing top surface (defined in mesh as 'free_surface_file')
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
+# C-PML boundary conditions for a regional simulation
+PML_CONDITIONS = .false.
+
+# C-PML top surface
+PML_INSTEAD_OF_FREE_SURFACE = .false.
+
+# C-PML thickness min/max
+PML_WIDTH_MIN = 10000.0d0
+PML_WIDTH_MAX = 25000.0d0
+
+# C-PML dominant frequency
+f0_FOR_PML = 12.7
+
+# parameters used to rotate C-PML boundary conditions by a given angle (not completed yet)
+# ROTATE_PML_ACTIVATE = .false.
+# ROTATE_PML_ANGLE = 0.
+
# save AVS or OpenDX movies
# MOVIE_TYPE = 1 to show the top surface
# MOVIE_TYPE = 2 to show all the external faces of the mesh
@@ -82,8 +99,7 @@ NTSTEP_BETWEEN_OUTPUT_SEISMOS = 10000
# 0 = read the whole adjoint sources at the same time
NTSTEP_BETWEEN_READ_ADJSRC = 0
-# use a (tilted) FORCESOLUTION force point source (or several) located exactly at a grid point
-# instead of a CMTSOLUTION moment-tensor source.
+# use a (tilted) FORCESOLUTION force point source (or several) instead of a CMTSOLUTION moment-tensor source.
# This can be useful e.g. for oil industry foothills simulations or asteroid simulations
# in which the source is a vertical force, normal force, inclined force, impact etc.
# If this flag is turned on, the FORCESOLUTION file must be edited by precising:
diff --git a/tpv16/DATA/Par_file b/tpv16/DATA/Par_file
index d576134..50daba6 100644
--- a/tpv16/DATA/Par_file
+++ b/tpv16/DATA/Par_file
@@ -54,6 +54,23 @@ ABSORBING_CONDITIONS = .true.
# absorbing top surface (defined in mesh as 'free_surface_file')
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
+# C-PML boundary conditions for a regional simulation
+PML_CONDITIONS = .false.
+
+# C-PML top surface
+PML_INSTEAD_OF_FREE_SURFACE = .false.
+
+# C-PML thickness min/max
+PML_WIDTH_MIN = 10000.d0
+PML_WIDTH_MAX = 25000.d0
+
+# C-PML dominant frequency
+f0_FOR_PML = 12.7
+
+# parameters used to rotate C-PML boundary conditions by a given angle (not completed yet)
+# ROTATE_PML_ACTIVATE = .false.
+# ROTATE_PML_ANGLE = 0.
+
# save AVS or OpenDX movies
# MOVIE_TYPE = 1 to show the top surface
# MOVIE_TYPE = 2 to show all the external faces of the mesh
@@ -82,8 +99,7 @@ NTSTEP_BETWEEN_OUTPUT_SEISMOS = 10000
# 0 = read the whole adjoint sources at the same time
NTSTEP_BETWEEN_READ_ADJSRC = 0
-# use a (tilted) FORCESOLUTION force point source (or several) located exactly at a grid point
-# instead of a CMTSOLUTION moment-tensor source.
+# use a (tilted) FORCESOLUTION force point source (or several) instead of a CMTSOLUTION moment-tensor source.
# This can be useful e.g. for oil industry foothills simulations or asteroid simulations
# in which the source is a vertical force, normal force, inclined force, impact etc.
# If this flag is turned on, the FORCESOLUTION file must be edited by precising:
diff --git a/tpv5/DATA/Par_file b/tpv5/DATA/Par_file
index 04cfdb3..7722794 100644
--- a/tpv5/DATA/Par_file
+++ b/tpv5/DATA/Par_file
@@ -54,6 +54,23 @@ ABSORBING_CONDITIONS = .true.
# absorbing top surface (defined in mesh as 'free_surface_file')
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
+# C-PML boundary conditions for a regional simulation
+PML_CONDITIONS = .false.
+
+# C-PML top surface
+PML_INSTEAD_OF_FREE_SURFACE = .false.
+
+# C-PML thickness min/max
+PML_WIDTH_MIN = 10000.0d0
+PML_WIDTH_MAX = 25000.0d0
+
+# C-PML dominant frequency
+f0_FOR_PML = 12.7
+
+# parameters used to rotate C-PML boundary conditions by a given angle (not completed yet)
+# ROTATE_PML_ACTIVATE = .false.
+# ROTATE_PML_ANGLE = 0.
+
# save AVS or OpenDX movies
# MOVIE_TYPE = 1 to show the top surface
# MOVIE_TYPE = 2 to show all the external faces of the mesh
@@ -82,8 +99,7 @@ NTSTEP_BETWEEN_OUTPUT_SEISMOS = 10000
# 0 = read the whole adjoint sources at the same time
NTSTEP_BETWEEN_READ_ADJSRC = 0
-# use a (tilted) FORCESOLUTION force point source (or several) located exactly at a grid point
-# instead of a CMTSOLUTION moment-tensor source.
+# use a (tilted) FORCESOLUTION force point source (or several) instead of a CMTSOLUTION moment-tensor source.
# This can be useful e.g. for oil industry foothills simulations or asteroid simulations
# in which the source is a vertical force, normal force, inclined force, impact etc.
# If this flag is turned on, the FORCESOLUTION file must be edited by precising:
diff --git a/waterlayered_halfspace/DATA/Par_file b/waterlayered_halfspace/DATA/Par_file
index 66d15a9..1531f40 100644
--- a/waterlayered_halfspace/DATA/Par_file
+++ b/waterlayered_halfspace/DATA/Par_file
@@ -54,6 +54,23 @@ ABSORBING_CONDITIONS = .true.
# absorbing top surface (defined in mesh as 'free_surface_file')
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
+# C-PML boundary conditions for a regional simulation
+PML_CONDITIONS = .false.
+
+# C-PML top surface
+PML_INSTEAD_OF_FREE_SURFACE = .false.
+
+# C-PML thickness min/max
+PML_WIDTH_MIN = 10000.0d0
+PML_WIDTH_MAX = 25000.0d0
+
+# C-PML dominant frequency
+f0_FOR_PML = 12.7
+
+# parameters used to rotate C-PML boundary conditions by a given angle (not completed yet)
+# ROTATE_PML_ACTIVATE = .false.
+# ROTATE_PML_ANGLE = 0.
+
# save AVS or OpenDX movies
# MOVIE_TYPE = 1 to show the top surface
# MOVIE_TYPE = 2 to show all the external faces of the mesh
@@ -79,10 +96,10 @@ NTSTEP_BETWEEN_OUTPUT_INFO = 500
NTSTEP_BETWEEN_OUTPUT_SEISMOS = 10000
# interval in time steps for reading adjoint traces
-NTSTEP_BETWEEN_READ_ADJSRC = 0 # 0 = read the whole adjoint sources at the same time
+# 0 = read the whole adjoint sources at the same time
+NTSTEP_BETWEEN_READ_ADJSRC = 0
-# use a (tilted) FORCESOLUTION force point source (or several) located exactly at a grid point
-# instead of a CMTSOLUTION moment-tensor source.
+# use a (tilted) FORCESOLUTION force point source (or several) instead of a CMTSOLUTION moment-tensor source.
# This can be useful e.g. for oil industry foothills simulations or asteroid simulations
# in which the source is a vertical force, normal force, inclined force, impact etc.
# If this flag is turned on, the FORCESOLUTION file must be edited by precising:
More information about the CIG-COMMITS
mailing list