[cig-commits] [commit] devel, master: adds point source flags in Par_file related to the use of an inclined force (2d6496b)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Wed Jun 18 16:49:44 PDT 2014
Repository : https://github.com/geodynamics/specfem3d
On branches: devel,master
Link : https://github.com/geodynamics/specfem3d/compare/6026e367984905ab133865f62fa6293b343759b9...47f703851338234f96397e7da9fbff63d8178b8a
>---------------------------------------------------------------
commit 2d6496b6b5cd47eb3f26a4cca0835e7e529b0cf6
Author: Joseph Charles <jcharles at lma.cnrs-mrs.fr>
Date: Wed Sep 19 18:47:12 2012 +0000
adds point source flags in Par_file related to the use of an inclined force
>---------------------------------------------------------------
2d6496b6b5cd47eb3f26a4cca0835e7e529b0cf6
.../ACOUSTIC/in_data_files/Par_file | 8 +++++---
.../ELASTIC/in_data_files/Par_file | 8 +++++---
Mount_StHelens/in_data_files/Par_file | 8 +++++---
homogeneous_halfspace/in_data_files/Par_file | 8 +++++---
homogeneous_poroelastic/in_data_files/Par_file | 8 +++++---
layered_halfspace/README | 8 --------
layered_halfspace/in_data_files/Par_file | 8 +++++---
layered_halfspace/process.sh | 16 ----------------
meshfem3D_examples/many_interfaces/Par_file | 8 +++++---
meshfem3D_examples/simple_model/in_data_files/Par_file | 8 +++++---
meshfem3D_examples/socal1D/in_data_files/Par_file | 8 +++++---
tomographic_model/in_data_files/Par_file | 8 +++++---
waterlayered_halfspace/in_data_files/Par_file | 8 +++++---
13 files changed, 55 insertions(+), 57 deletions(-)
diff --git a/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC/in_data_files/Par_file b/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC/in_data_files/Par_file
index faa4fae..f327661 100644
--- a/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC/in_data_files/Par_file
+++ b/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC/in_data_files/Par_file
@@ -63,11 +63,13 @@ NTSTEP_BETWEEN_READ_ADJSRC = 0 # 0 = read the whole adjoint sources at
# use a force source located exactly at a grid point instead of a CMTSOLUTION 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, impact etc.
+# in which the source is a vertical force, normal force, inclined force, impact etc.
USE_FORCE_POINT_SOURCE = .false.
FACTOR_FORCE_SOURCE = 1.d15
-# direction in comp E/N/Z = 1/2/3
-COMPONENT_FORCE_SOURCE = 3
+# components of a (non-unitary) direction vector for the force source on the E/N/Z_UP basis
+COMPONENT_DIR_VECT_SOURCE_E = 1.d0
+COMPONENT_DIR_VECT_SOURCE_N = -2.d0
+COMPONENT_DIR_VECT_SOURCE_Z_UP = -1.d0
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .true.
diff --git a/BENCHMARK_CLAERBOUT_ADJOINT/ELASTIC/in_data_files/Par_file b/BENCHMARK_CLAERBOUT_ADJOINT/ELASTIC/in_data_files/Par_file
index b365d86..2b1bfc7 100644
--- a/BENCHMARK_CLAERBOUT_ADJOINT/ELASTIC/in_data_files/Par_file
+++ b/BENCHMARK_CLAERBOUT_ADJOINT/ELASTIC/in_data_files/Par_file
@@ -63,11 +63,13 @@ NTSTEP_BETWEEN_READ_ADJSRC = 0 # 0 = read the whole adjoint sources at
# use a force source located exactly at a grid point instead of a CMTSOLUTION 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, impact etc.
+# in which the source is a vertical force, normal force, inclined force, impact etc.
USE_FORCE_POINT_SOURCE = .false.
FACTOR_FORCE_SOURCE = 1.d15
-# direction in comp E/N/Z = 1/2/3
-COMPONENT_FORCE_SOURCE = 3
+# components of a (non-unitary) direction vector for the force source on the E/N/Z_UP basis
+COMPONENT_DIR_VECT_SOURCE_E = 1.d0
+COMPONENT_DIR_VECT_SOURCE_N = -2.d0
+COMPONENT_DIR_VECT_SOURCE_Z_UP = -1.d0
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .true.
diff --git a/Mount_StHelens/in_data_files/Par_file b/Mount_StHelens/in_data_files/Par_file
index 081bb56..dc758e4 100644
--- a/Mount_StHelens/in_data_files/Par_file
+++ b/Mount_StHelens/in_data_files/Par_file
@@ -63,11 +63,13 @@ NTSTEP_BETWEEN_READ_ADJSRC = 0 # 0 = read the whole adjoint sources at
# use a force source located exactly at a grid point instead of a CMTSOLUTION 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, impact etc.
+# in which the source is a vertical force, normal force, inclined force, impact etc.
USE_FORCE_POINT_SOURCE = .false.
FACTOR_FORCE_SOURCE = 1.d15
-# direction in comp E/N/Z = 1/2/3
-COMPONENT_FORCE_SOURCE = 3
+# components of a (non-unitary) direction vector for the force source on the E/N/Z_UP basis
+COMPONENT_DIR_VECT_SOURCE_E = 1.d0
+COMPONENT_DIR_VECT_SOURCE_N = -2.d0
+COMPONENT_DIR_VECT_SOURCE_Z_UP = -1.d0
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .false.
diff --git a/homogeneous_halfspace/in_data_files/Par_file b/homogeneous_halfspace/in_data_files/Par_file
index 1f9a162..3a380ef 100644
--- a/homogeneous_halfspace/in_data_files/Par_file
+++ b/homogeneous_halfspace/in_data_files/Par_file
@@ -65,11 +65,13 @@ NTSTEP_BETWEEN_READ_ADJSRC = 0
# use a force source located exactly at a grid point instead of a CMTSOLUTION 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, impact etc.
+# in which the source is a vertical force, normal force, inclined force, impact etc.
USE_FORCE_POINT_SOURCE = .false.
FACTOR_FORCE_SOURCE = 1.d15
-# direction in comp E/N/Z = 1/2/3
-COMPONENT_FORCE_SOURCE = 3
+# components of a (non-unitary) direction vector for the force source on the E/N/Z_UP basis
+COMPONENT_DIR_VECT_SOURCE_E = 1.d0
+COMPONENT_DIR_VECT_SOURCE_N = -2.d0
+COMPONENT_DIR_VECT_SOURCE_Z_UP = -1.d0
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .false.
diff --git a/homogeneous_poroelastic/in_data_files/Par_file b/homogeneous_poroelastic/in_data_files/Par_file
index d821556..e15fa33 100644
--- a/homogeneous_poroelastic/in_data_files/Par_file
+++ b/homogeneous_poroelastic/in_data_files/Par_file
@@ -65,11 +65,13 @@ NTSTEP_BETWEEN_READ_ADJSRC = 0
# use a force source located exactly at a grid point instead of a CMTSOLUTION 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, impact etc.
+# in which the source is a vertical force, normal force, inclined force, impact etc.
USE_FORCE_POINT_SOURCE = .false.
FACTOR_FORCE_SOURCE = 1.d15
-# direction in comp E/N/Z = 1/2/3
-COMPONENT_FORCE_SOURCE = 3
+# components of a (non-unitary) direction vector for the force source on the E/N/Z_UP basis
+COMPONENT_DIR_VECT_SOURCE_E = 1.d0
+COMPONENT_DIR_VECT_SOURCE_N = -2.d0
+COMPONENT_DIR_VECT_SOURCE_Z_UP = -1.d0
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .true.
diff --git a/layered_halfspace/README b/layered_halfspace/README
index e4743ba..2548487 100644
--- a/layered_halfspace/README
+++ b/layered_halfspace/README
@@ -107,14 +107,6 @@ step-by-step tutorial:
5. run simulation:
- - modify src/shared/constants.h:
- to use a vertical force source, with a Ricker wavelet source time function,
- turn flag on for parameter USE_FORCE_POINT_SOURCE:
- ...
- logical, parameter :: USE_FORCE_POINT_SOURCE = .true.
- double precision, parameter :: FACTOR_FORCE_SOURCE = -1.d15
- ...
-
- compile specfem3D:
> make xspecfem3D
diff --git a/layered_halfspace/in_data_files/Par_file b/layered_halfspace/in_data_files/Par_file
index fe893cd..314b71d 100644
--- a/layered_halfspace/in_data_files/Par_file
+++ b/layered_halfspace/in_data_files/Par_file
@@ -64,11 +64,13 @@ NTSTEP_BETWEEN_READ_ADJSRC = 0 # 0 = read the whole adjoint sources at
# use a force source located exactly at a grid point instead of a CMTSOLUTION 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, impact etc.
+# in which the source is a vertical force, normal force, inclined force, impact etc.
USE_FORCE_POINT_SOURCE = .false.
FACTOR_FORCE_SOURCE = 1.d15
-# direction in comp E/N/Z = 1/2/3
-COMPONENT_FORCE_SOURCE = 3
+# components of a (non-unitary) direction vector for the force source on the E/N/Z_UP basis
+COMPONENT_DIR_VECT_SOURCE_E = 1.d0
+COMPONENT_DIR_VECT_SOURCE_N = -2.d0
+COMPONENT_DIR_VECT_SOURCE_Z_UP = -1.d0
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .false.
diff --git a/layered_halfspace/process.sh b/layered_halfspace/process.sh
index 6b26de3..7b36e96 100755
--- a/layered_halfspace/process.sh
+++ b/layered_halfspace/process.sh
@@ -34,22 +34,6 @@ mkdir -p in_out_files/DATABASES_MPI
rm -rf in_out_files/OUTPUT_FILES/*
rm -rf in_out_files/DATABASES_MPI/*
-# compiles executables in root directory
-cd ../../
-
-# compiles with flag for a point force (with a ricker source time function)
-cd src/shared/
-cp constants.h constants.h.org
-sed -e "s:USE_FORCE_POINT_SOURCE.*:USE_FORCE_POINT_SOURCE = .true. :" constants.h.org > constants.h.tmp
-sed -e "s:FACTOR_FORCE_SOURCE.*:FACTOR_FORCE_SOURCE = -1.d15 :" constants.h.tmp > constants.h
-rm -f constants.h.tmp
-cd ../../
-make > tmp.log
-
-# backup & restores original file again
-cp src/shared/constants.h $currentdir/in_out_files/OUTPUT_FILES/constants.h.bak
-mv src/shared/constants.h.org src/shared/constants.h
-
cd $currentdir
# links executables
diff --git a/meshfem3D_examples/many_interfaces/Par_file b/meshfem3D_examples/many_interfaces/Par_file
index 8ac6d6f..1babb06 100644
--- a/meshfem3D_examples/many_interfaces/Par_file
+++ b/meshfem3D_examples/many_interfaces/Par_file
@@ -64,11 +64,13 @@ NTSTEP_BETWEEN_READ_ADJSRC = 0 # 0 = read the whole adjoint sources at
# use a force source located exactly at a grid point instead of a CMTSOLUTION 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, impact etc.
+# in which the source is a vertical force, normal force, inclined force, impact etc.
USE_FORCE_POINT_SOURCE = .false.
FACTOR_FORCE_SOURCE = 1.d15
-# direction in comp E/N/Z = 1/2/3
-COMPONENT_FORCE_SOURCE = 3
+# components of a (non-unitary) direction vector for the force source on the E/N/Z_UP basis
+COMPONENT_DIR_VECT_SOURCE_E = 1.d0
+COMPONENT_DIR_VECT_SOURCE_N = -2.d0
+COMPONENT_DIR_VECT_SOURCE_Z_UP = -1.d0
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .false.
diff --git a/meshfem3D_examples/simple_model/in_data_files/Par_file b/meshfem3D_examples/simple_model/in_data_files/Par_file
index 8e3fff1..a0e7114 100644
--- a/meshfem3D_examples/simple_model/in_data_files/Par_file
+++ b/meshfem3D_examples/simple_model/in_data_files/Par_file
@@ -64,11 +64,13 @@ NTSTEP_BETWEEN_READ_ADJSRC = 0 # 0 = read the whole adjoint sources at
# use a force source located exactly at a grid point instead of a CMTSOLUTION 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, impact etc.
+# in which the source is a vertical force, normal force, inclined force, impact etc.
USE_FORCE_POINT_SOURCE = .false.
FACTOR_FORCE_SOURCE = 1.d15
-# direction in comp E/N/Z = 1/2/3
-COMPONENT_FORCE_SOURCE = 3
+# components of a (non-unitary) direction vector for the force source on the E/N/Z_UP basis
+COMPONENT_DIR_VECT_SOURCE_E = 1.d0
+COMPONENT_DIR_VECT_SOURCE_N = -2.d0
+COMPONENT_DIR_VECT_SOURCE_Z_UP = -1.d0
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .false.
diff --git a/meshfem3D_examples/socal1D/in_data_files/Par_file b/meshfem3D_examples/socal1D/in_data_files/Par_file
index 053a064..0c646a9 100644
--- a/meshfem3D_examples/socal1D/in_data_files/Par_file
+++ b/meshfem3D_examples/socal1D/in_data_files/Par_file
@@ -64,11 +64,13 @@ NTSTEP_BETWEEN_READ_ADJSRC = 0 # 0 = read the whole adjoint sources at
# use a force source located exactly at a grid point instead of a CMTSOLUTION 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, impact etc.
+# in which the source is a vertical force, normal force, inclined force, impact etc.
USE_FORCE_POINT_SOURCE = .false.
FACTOR_FORCE_SOURCE = 1.d15
-# direction in comp E/N/Z = 1/2/3
-COMPONENT_FORCE_SOURCE = 3
+# components of a (non-unitary) direction vector for the force source on the E/N/Z_UP basis
+COMPONENT_DIR_VECT_SOURCE_E = 1.d0
+COMPONENT_DIR_VECT_SOURCE_N = -2.d0
+COMPONENT_DIR_VECT_SOURCE_Z_UP = -1.d0
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .true.
diff --git a/tomographic_model/in_data_files/Par_file b/tomographic_model/in_data_files/Par_file
index e8ac208..f47900b 100644
--- a/tomographic_model/in_data_files/Par_file
+++ b/tomographic_model/in_data_files/Par_file
@@ -64,11 +64,13 @@ NTSTEP_BETWEEN_READ_ADJSRC = 0 # 0 = read the whole adjoint sources at
# use a force source located exactly at a grid point instead of a CMTSOLUTION 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, impact etc.
+# in which the source is a vertical force, normal force, inclined force, impact etc.
USE_FORCE_POINT_SOURCE = .false.
FACTOR_FORCE_SOURCE = 1.d15
-# direction in comp E/N/Z = 1/2/3
-COMPONENT_FORCE_SOURCE = 3
+# components of a (non-unitary) direction vector for the force source on the E/N/Z_UP basis
+COMPONENT_DIR_VECT_SOURCE_E = 1.d0
+COMPONENT_DIR_VECT_SOURCE_N = -2.d0
+COMPONENT_DIR_VECT_SOURCE_Z_UP = -1.d0
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .false.
diff --git a/waterlayered_halfspace/in_data_files/Par_file b/waterlayered_halfspace/in_data_files/Par_file
index 2469e2b..9189e75 100644
--- a/waterlayered_halfspace/in_data_files/Par_file
+++ b/waterlayered_halfspace/in_data_files/Par_file
@@ -64,11 +64,13 @@ NTSTEP_BETWEEN_READ_ADJSRC = 0 # 0 = read the whole adjoint sources at
# use a force source located exactly at a grid point instead of a CMTSOLUTION 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, impact etc.
+# in which the source is a vertical force, normal force, inclined force, impact etc.
USE_FORCE_POINT_SOURCE = .false.
FACTOR_FORCE_SOURCE = 1.d15
-# direction in comp E/N/Z = 1/2/3
-COMPONENT_FORCE_SOURCE = 3
+# components of a (non-unitary) direction vector for the force source on the E/N/Z_UP basis
+COMPONENT_DIR_VECT_SOURCE_E = 1.d0
+COMPONENT_DIR_VECT_SOURCE_N = -2.d0
+COMPONENT_DIR_VECT_SOURCE_Z_UP = -1.d0
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .false.
More information about the CIG-COMMITS
mailing list