[cig-commits] [commit] devel, master: moves parameter 'NGNOD' from constants.h to Par_file, updates the 'simple_model' example (0849d5d)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Wed Jun 18 16:50:16 PDT 2014
Repository : https://github.com/geodynamics/specfem3d
On branches: devel,master
Link : https://github.com/geodynamics/specfem3d/compare/6026e367984905ab133865f62fa6293b343759b9...47f703851338234f96397e7da9fbff63d8178b8a
>---------------------------------------------------------------
commit 0849d5dbcdc0643d429dbd0012191a3c9b06cf20
Author: Joseph Charles <jcharles at lma.cnrs-mrs.fr>
Date: Thu Oct 25 00:01:02 2012 +0000
moves parameter 'NGNOD' from constants.h to Par_file, updates the 'simple_model' example
>---------------------------------------------------------------
0849d5dbcdc0643d429dbd0012191a3c9b06cf20
BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC/in_data_files/Par_file | 9 ++++++++-
BENCHMARK_CLAERBOUT_ADJOINT/ELASTIC/in_data_files/Par_file | 9 ++++++++-
Mount_StHelens/in_data_files/Par_file | 9 ++++++++-
homogeneous_halfspace_HEX27/in_data_files/Par_file | 9 ++++++++-
homogeneous_halfspace_HEX8/DATA/Par_file | 9 ++++++++-
homogeneous_poroelastic/in_data_files/Par_file | 9 ++++++++-
layered_halfspace/in_data_files/Par_file | 9 ++++++++-
meshfem3D_examples/many_interfaces/Par_file | 9 ++++++++-
.../simple_model/{in_data_files => DATA}/CMTSOLUTION | 0
.../simple_model/{in_data_files => DATA}/FORCESOLUTION | 0
.../simple_model/{in_data_files => DATA}/Par_file | 9 ++++++++-
.../simple_model/{in_data_files => DATA}/STATIONS | 0
.../{in_data_files => DATA}/meshfem3D_files/Mesh_Par_file | 0
.../{in_data_files => DATA}/meshfem3D_files/interface1.dat | 0
.../{in_data_files => DATA}/meshfem3D_files/interface2.dat | 0
.../{in_data_files => DATA}/meshfem3D_files/interfaces.dat | 0
meshfem3D_examples/simple_model/process.sh | 10 ++++++----
meshfem3D_examples/socal1D/in_data_files/Par_file | 9 ++++++++-
tomographic_model/in_data_files/Par_file | 9 ++++++++-
waterlayered_halfspace/in_data_files/Par_file | 9 ++++++++-
20 files changed, 102 insertions(+), 16 deletions(-)
diff --git a/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC/in_data_files/Par_file b/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC/in_data_files/Par_file
index 4a179f7..9bd8a54 100644
--- a/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC/in_data_files/Par_file
+++ b/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC/in_data_files/Par_file
@@ -18,6 +18,13 @@ NPROC = 4
NSTEP = 10000
DT = 0.001
+# number of nodes for 2D and 3D shape functions for hexahedra
+# we use either 8-node mesh elements (bricks) or 27-node elements.
+# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
+# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
+# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
+NGNOD = 8
+
# models:
# available options are:
# default (model parameters described by mesh properties)
@@ -67,7 +74,6 @@ NTSTEP_BETWEEN_READ_ADJSRC = 0 # 0 = read the whole adjoint sources at
# use a force source (or several) 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, inclined force, impact etc.
-USE_FORCE_POINT_SOURCE = .false.
# If this flag is turned on, the FORCESOLUTION file must be edited by precising:
# - the corresponding time-shift parameter,
# - the half duration parameter of the source,
@@ -76,6 +82,7 @@ USE_FORCE_POINT_SOURCE = .false.
# - the components of a (non-unitary) direction vector for the force source in the E/N/Z_UP basis.
# The direction vector is made unitary internally in the code and thus only its direction matters here;
# its norm is ignored and the norm of the force used is the factor force source times the source time function.
+USE_FORCE_POINT_SOURCE = .false.
# set to use a Ricker source time function instead of a gaussian (when using CMTSOLUTION files)
# or a quasi-Heaviside (when using FORCESOLUTION files)
diff --git a/BENCHMARK_CLAERBOUT_ADJOINT/ELASTIC/in_data_files/Par_file b/BENCHMARK_CLAERBOUT_ADJOINT/ELASTIC/in_data_files/Par_file
index bf4999a..d9f6ee8 100644
--- a/BENCHMARK_CLAERBOUT_ADJOINT/ELASTIC/in_data_files/Par_file
+++ b/BENCHMARK_CLAERBOUT_ADJOINT/ELASTIC/in_data_files/Par_file
@@ -18,6 +18,13 @@ NPROC = 4
NSTEP = 4000
DT = 0.0006
+# number of nodes for 2D and 3D shape functions for hexahedra
+# we use either 8-node mesh elements (bricks) or 27-node elements.
+# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
+# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
+# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
+NGNOD = 8
+
# models:
# available options are:
# default (model parameters described by mesh properties)
@@ -67,7 +74,6 @@ NTSTEP_BETWEEN_READ_ADJSRC = 0 # 0 = read the whole adjoint sources at
# use a force source (or several) 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, inclined force, impact etc.
-USE_FORCE_POINT_SOURCE = .false.
# If this flag is turned on, the FORCESOLUTION file must be edited by precising:
# - the corresponding time-shift parameter,
# - the half duration parameter of the source,
@@ -76,6 +82,7 @@ USE_FORCE_POINT_SOURCE = .false.
# - the components of a (non-unitary) direction vector for the force source in the E/N/Z_UP basis.
# The direction vector is made unitary internally in the code and thus only its direction matters here;
# its norm is ignored and the norm of the force used is the factor force source times the source time function.
+USE_FORCE_POINT_SOURCE = .false.
# set to use a Ricker source time function instead of a gaussian (when using CMTSOLUTION files)
# or a quasi-Heaviside (when using FORCESOLUTION files)
diff --git a/Mount_StHelens/in_data_files/Par_file b/Mount_StHelens/in_data_files/Par_file
index 5cd4343..0ecdceb 100644
--- a/Mount_StHelens/in_data_files/Par_file
+++ b/Mount_StHelens/in_data_files/Par_file
@@ -18,6 +18,13 @@ NPROC = 4
NSTEP = 2500
DT = 0.005d0
+# number of nodes for 2D and 3D shape functions for hexahedra
+# we use either 8-node mesh elements (bricks) or 27-node elements.
+# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
+# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
+# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
+NGNOD = 8
+
# models:
# available options are:
# default (model parameters described by mesh properties)
@@ -67,7 +74,6 @@ NTSTEP_BETWEEN_READ_ADJSRC = 0 # 0 = read the whole adjoint sources at
# use a force source (or several) 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, inclined force, impact etc.
-USE_FORCE_POINT_SOURCE = .false.
# If this flag is turned on, the FORCESOLUTION file must be edited by precising:
# - the corresponding time-shift parameter,
# - the half duration parameter of the source,
@@ -76,6 +82,7 @@ USE_FORCE_POINT_SOURCE = .false.
# - the components of a (non-unitary) direction vector for the force source in the E/N/Z_UP basis.
# The direction vector is made unitary internally in the code and thus only its direction matters here;
# its norm is ignored and the norm of the force used is the factor force source times the source time function.
+USE_FORCE_POINT_SOURCE = .false.
# set to use a Ricker source time function instead of a gaussian (when using CMTSOLUTION files)
# or a quasi-Heaviside (when using FORCESOLUTION files)
diff --git a/homogeneous_halfspace_HEX27/in_data_files/Par_file b/homogeneous_halfspace_HEX27/in_data_files/Par_file
index f09b00d..6bc5b42 100644
--- a/homogeneous_halfspace_HEX27/in_data_files/Par_file
+++ b/homogeneous_halfspace_HEX27/in_data_files/Par_file
@@ -18,6 +18,13 @@ NPROC = 4
NSTEP = 1000
DT = 0.05
+# number of nodes for 2D and 3D shape functions for hexahedra
+# we use either 8-node mesh elements (bricks) or 27-node elements.
+# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
+# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
+# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
+NGNOD = 8
+
# models:
# available options are:
# default (model parameters described by mesh properties)
@@ -69,7 +76,6 @@ NTSTEP_BETWEEN_READ_ADJSRC = 0
# use a force source (or several) 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, inclined force, impact etc.
-USE_FORCE_POINT_SOURCE = .false.
# If this flag is turned on, the FORCESOLUTION file must be edited by precising:
# - the corresponding time-shift parameter,
# - the half duration parameter of the source,
@@ -78,6 +84,7 @@ USE_FORCE_POINT_SOURCE = .false.
# - the components of a (non-unitary) direction vector for the force source in the E/N/Z_UP basis.
# The direction vector is made unitary internally in the code and thus only its direction matters here;
# its norm is ignored and the norm of the force used is the factor force source times the source time function.
+USE_FORCE_POINT_SOURCE = .false.
# set to use a Ricker source time function instead of a gaussian (when using CMTSOLUTION files)
# or a quasi-Heaviside (when using FORCESOLUTION files)
diff --git a/homogeneous_halfspace_HEX8/DATA/Par_file b/homogeneous_halfspace_HEX8/DATA/Par_file
index f09b00d..6bc5b42 100644
--- a/homogeneous_halfspace_HEX8/DATA/Par_file
+++ b/homogeneous_halfspace_HEX8/DATA/Par_file
@@ -18,6 +18,13 @@ NPROC = 4
NSTEP = 1000
DT = 0.05
+# number of nodes for 2D and 3D shape functions for hexahedra
+# we use either 8-node mesh elements (bricks) or 27-node elements.
+# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
+# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
+# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
+NGNOD = 8
+
# models:
# available options are:
# default (model parameters described by mesh properties)
@@ -69,7 +76,6 @@ NTSTEP_BETWEEN_READ_ADJSRC = 0
# use a force source (or several) 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, inclined force, impact etc.
-USE_FORCE_POINT_SOURCE = .false.
# If this flag is turned on, the FORCESOLUTION file must be edited by precising:
# - the corresponding time-shift parameter,
# - the half duration parameter of the source,
@@ -78,6 +84,7 @@ USE_FORCE_POINT_SOURCE = .false.
# - the components of a (non-unitary) direction vector for the force source in the E/N/Z_UP basis.
# The direction vector is made unitary internally in the code and thus only its direction matters here;
# its norm is ignored and the norm of the force used is the factor force source times the source time function.
+USE_FORCE_POINT_SOURCE = .false.
# set to use a Ricker source time function instead of a gaussian (when using CMTSOLUTION files)
# or a quasi-Heaviside (when using FORCESOLUTION files)
diff --git a/homogeneous_poroelastic/in_data_files/Par_file b/homogeneous_poroelastic/in_data_files/Par_file
index 57b52fe..2e9c886 100644
--- a/homogeneous_poroelastic/in_data_files/Par_file
+++ b/homogeneous_poroelastic/in_data_files/Par_file
@@ -18,6 +18,13 @@ NPROC = 14
NSTEP = 4000
DT = 5.d-5 #0.5d-3
+# number of nodes for 2D and 3D shape functions for hexahedra
+# we use either 8-node mesh elements (bricks) or 27-node elements.
+# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
+# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
+# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
+NGNOD = 8
+
# models:
# available options are:
# default (model parameters described by mesh properties)
@@ -69,7 +76,6 @@ NTSTEP_BETWEEN_READ_ADJSRC = 0
# use a force source (or several) 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, inclined force, impact etc.
-USE_FORCE_POINT_SOURCE = .false.
# If this flag is turned on, the FORCESOLUTION file must be edited by precising:
# - the corresponding time-shift parameter,
# - the half duration parameter of the source,
@@ -78,6 +84,7 @@ USE_FORCE_POINT_SOURCE = .false.
# - the components of a (non-unitary) direction vector for the force source in the E/N/Z_UP basis.
# The direction vector is made unitary internally in the code and thus only its direction matters here;
# its norm is ignored and the norm of the force used is the factor force source times the source time function.
+USE_FORCE_POINT_SOURCE = .false.
# set to use a Ricker source time function instead of a gaussian (when using CMTSOLUTION files)
# or a quasi-Heaviside (when using FORCESOLUTION files)
diff --git a/layered_halfspace/in_data_files/Par_file b/layered_halfspace/in_data_files/Par_file
index 2a7e18c..1182b7c 100644
--- a/layered_halfspace/in_data_files/Par_file
+++ b/layered_halfspace/in_data_files/Par_file
@@ -18,6 +18,13 @@ NPROC = 4
NSTEP = 3000
DT = 0.0075
+# number of nodes for 2D and 3D shape functions for hexahedra
+# we use either 8-node mesh elements (bricks) or 27-node elements.
+# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
+# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
+# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
+NGNOD = 8
+
# models:
# available options are:
# default (model parameters described by mesh properties)
@@ -68,7 +75,6 @@ NTSTEP_BETWEEN_READ_ADJSRC = 0 # 0 = read the whole adjoint sources at
# use a force source (or several) 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, inclined force, impact etc.
-USE_FORCE_POINT_SOURCE = .false.
# If this flag is turned on, the FORCESOLUTION file must be edited by precising:
# - the corresponding time-shift parameter,
# - the half duration parameter of the source,
@@ -77,6 +83,7 @@ USE_FORCE_POINT_SOURCE = .false.
# - the components of a (non-unitary) direction vector for the force source in the E/N/Z_UP basis.
# The direction vector is made unitary internally in the code and thus only its direction matters here;
# its norm is ignored and the norm of the force used is the factor force source times the source time function.
+USE_FORCE_POINT_SOURCE = .false.
# set to use a Ricker source time function instead of a gaussian (when using CMTSOLUTION files)
# or a quasi-Heaviside (when using FORCESOLUTION files)
diff --git a/meshfem3D_examples/many_interfaces/Par_file b/meshfem3D_examples/many_interfaces/Par_file
index 2e04c4b..74e744f 100644
--- a/meshfem3D_examples/many_interfaces/Par_file
+++ b/meshfem3D_examples/many_interfaces/Par_file
@@ -18,6 +18,13 @@ NPROC = 81
NSTEP = 2000
DT = 0.03
+# number of nodes for 2D and 3D shape functions for hexahedra
+# we use either 8-node mesh elements (bricks) or 27-node elements.
+# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
+# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
+# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
+NGNOD = 8
+
# models:
# available options are:
# default (model parameters described by mesh properties)
@@ -68,7 +75,6 @@ NTSTEP_BETWEEN_READ_ADJSRC = 0 # 0 = read the whole adjoint sources at
# use a force source (or several) 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, inclined force, impact etc.
-USE_FORCE_POINT_SOURCE = .false.
# If this flag is turned on, the FORCESOLUTION file must be edited by precising:
# - the corresponding time-shift parameter,
# - the half duration parameter of the source,
@@ -77,6 +83,7 @@ USE_FORCE_POINT_SOURCE = .false.
# - the components of a (non-unitary) direction vector for the force source in the E/N/Z_UP basis.
# The direction vector is made unitary internally in the code and thus only its direction matters here;
# its norm is ignored and the norm of the force used is the factor force source times the source time function.
+USE_FORCE_POINT_SOURCE = .false.
# set to use a Ricker source time function instead of a gaussian (when using CMTSOLUTION files)
# or a quasi-Heaviside (when using FORCESOLUTION files)
diff --git a/meshfem3D_examples/simple_model/in_data_files/CMTSOLUTION b/meshfem3D_examples/simple_model/DATA/CMTSOLUTION
similarity index 100%
rename from meshfem3D_examples/simple_model/in_data_files/CMTSOLUTION
rename to meshfem3D_examples/simple_model/DATA/CMTSOLUTION
diff --git a/meshfem3D_examples/simple_model/in_data_files/FORCESOLUTION b/meshfem3D_examples/simple_model/DATA/FORCESOLUTION
similarity index 100%
rename from meshfem3D_examples/simple_model/in_data_files/FORCESOLUTION
rename to meshfem3D_examples/simple_model/DATA/FORCESOLUTION
diff --git a/meshfem3D_examples/simple_model/in_data_files/Par_file b/meshfem3D_examples/simple_model/DATA/Par_file
similarity index 87%
rename from meshfem3D_examples/simple_model/in_data_files/Par_file
rename to meshfem3D_examples/simple_model/DATA/Par_file
index 622785a..743a5da 100644
--- a/meshfem3D_examples/simple_model/in_data_files/Par_file
+++ b/meshfem3D_examples/simple_model/DATA/Par_file
@@ -18,6 +18,13 @@ NPROC = 4
NSTEP = 2000
DT = 0.03
+# number of nodes for 2D and 3D shape functions for hexahedra
+# we use either 8-node mesh elements (bricks) or 27-node elements.
+# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
+# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
+# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
+NGNOD = 8
+
# models:
# available options are:
# default (model parameters described by mesh properties)
@@ -68,7 +75,6 @@ NTSTEP_BETWEEN_READ_ADJSRC = 0 # 0 = read the whole adjoint sources at
# use a force source (or several) 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, inclined force, impact etc.
-USE_FORCE_POINT_SOURCE = .false.
# If this flag is turned on, the FORCESOLUTION file must be edited by precising:
# - the corresponding time-shift parameter,
# - the half duration parameter of the source,
@@ -77,6 +83,7 @@ USE_FORCE_POINT_SOURCE = .false.
# - the components of a (non-unitary) direction vector for the force source in the E/N/Z_UP basis.
# The direction vector is made unitary internally in the code and thus only its direction matters here;
# its norm is ignored and the norm of the force used is the factor force source times the source time function.
+USE_FORCE_POINT_SOURCE = .false.
# set to use a Ricker source time function instead of a gaussian (when using CMTSOLUTION files)
# or a quasi-Heaviside (when using FORCESOLUTION files)
diff --git a/meshfem3D_examples/simple_model/in_data_files/STATIONS b/meshfem3D_examples/simple_model/DATA/STATIONS
similarity index 100%
rename from meshfem3D_examples/simple_model/in_data_files/STATIONS
rename to meshfem3D_examples/simple_model/DATA/STATIONS
diff --git a/meshfem3D_examples/simple_model/in_data_files/meshfem3D_files/Mesh_Par_file b/meshfem3D_examples/simple_model/DATA/meshfem3D_files/Mesh_Par_file
similarity index 100%
rename from meshfem3D_examples/simple_model/in_data_files/meshfem3D_files/Mesh_Par_file
rename to meshfem3D_examples/simple_model/DATA/meshfem3D_files/Mesh_Par_file
diff --git a/meshfem3D_examples/simple_model/in_data_files/meshfem3D_files/interface1.dat b/meshfem3D_examples/simple_model/DATA/meshfem3D_files/interface1.dat
similarity index 100%
rename from meshfem3D_examples/simple_model/in_data_files/meshfem3D_files/interface1.dat
rename to meshfem3D_examples/simple_model/DATA/meshfem3D_files/interface1.dat
diff --git a/meshfem3D_examples/simple_model/in_data_files/meshfem3D_files/interface2.dat b/meshfem3D_examples/simple_model/DATA/meshfem3D_files/interface2.dat
similarity index 100%
rename from meshfem3D_examples/simple_model/in_data_files/meshfem3D_files/interface2.dat
rename to meshfem3D_examples/simple_model/DATA/meshfem3D_files/interface2.dat
diff --git a/meshfem3D_examples/simple_model/in_data_files/meshfem3D_files/interfaces.dat b/meshfem3D_examples/simple_model/DATA/meshfem3D_files/interfaces.dat
similarity index 100%
rename from meshfem3D_examples/simple_model/in_data_files/meshfem3D_files/interfaces.dat
rename to meshfem3D_examples/simple_model/DATA/meshfem3D_files/interfaces.dat
diff --git a/meshfem3D_examples/simple_model/process.sh b/meshfem3D_examples/simple_model/process.sh
index c6102fa..9f09831 100755
--- a/meshfem3D_examples/simple_model/process.sh
+++ b/meshfem3D_examples/simple_model/process.sh
@@ -24,15 +24,16 @@ echo " setting up example..."
echo
mkdir -p bin
-mkdir -p OUTPUT_FILES
-mkdir -p OUTPUT_FILES/DATABASES_MPI
+mkdir -p OUTPUT_FILES/
+mkdir -p DATABASES_MPI
-rm -rf OUTPUT_FILES/*
-rm -rf OUTPUT_FILES/DATABASES_MPI/*
+rm -f OUTPUT_FILES/*
+rm -rf DATABASES_MPI/*
# compiles executables in root directory
cd ../../../
+make clean
make > $currentdir/tmp.log
cd $currentdir
@@ -44,6 +45,7 @@ if [ ! -e bin/xspecfem3D ]; then echo "compilation failed, please check..."; exi
# stores setup
cp DATA/Par_file OUTPUT_FILES/
cp DATA/CMTSOLUTION OUTPUT_FILES/
+cp DATA/FORCESOLUTION OUTPUT_FILES/
cp DATA/STATIONS OUTPUT_FILES/
# creates and decomposes mesh
diff --git a/meshfem3D_examples/socal1D/in_data_files/Par_file b/meshfem3D_examples/socal1D/in_data_files/Par_file
index c9650f6..8363063 100644
--- a/meshfem3D_examples/socal1D/in_data_files/Par_file
+++ b/meshfem3D_examples/socal1D/in_data_files/Par_file
@@ -18,6 +18,13 @@ NPROC = 1
NSTEP = 840
DT = 0.012
+# number of nodes for 2D and 3D shape functions for hexahedra
+# we use either 8-node mesh elements (bricks) or 27-node elements.
+# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
+# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
+# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
+NGNOD = 8
+
# models:
# available options are:
# default (model parameters described by mesh properties)
@@ -68,7 +75,6 @@ NTSTEP_BETWEEN_READ_ADJSRC = 0 # 0 = read the whole adjoint sources at
# use a force source (or several) 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, inclined force, impact etc.
-USE_FORCE_POINT_SOURCE = .false.
# If this flag is turned on, the FORCESOLUTION file must be edited by precising:
# - the corresponding time-shift parameter,
# - the half duration parameter of the source,
@@ -77,6 +83,7 @@ USE_FORCE_POINT_SOURCE = .false.
# - the components of a (non-unitary) direction vector for the force source in the E/N/Z_UP basis.
# The direction vector is made unitary internally in the code and thus only its direction matters here;
# its norm is ignored and the norm of the force used is the factor force source times the source time function.
+USE_FORCE_POINT_SOURCE = .false.
# set to use a Ricker source time function instead of a gaussian (when using CMTSOLUTION files)
# or a quasi-Heaviside (when using FORCESOLUTION files)
diff --git a/tomographic_model/in_data_files/Par_file b/tomographic_model/in_data_files/Par_file
index 22cf46f..e7a10c8 100644
--- a/tomographic_model/in_data_files/Par_file
+++ b/tomographic_model/in_data_files/Par_file
@@ -18,6 +18,13 @@ NPROC = 4
NSTEP = 4000
DT = 0.01
+# number of nodes for 2D and 3D shape functions for hexahedra
+# we use either 8-node mesh elements (bricks) or 27-node elements.
+# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
+# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
+# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
+NGNOD = 8
+
# models:
# available options are:
# default (model parameters described by mesh properties)
@@ -68,7 +75,6 @@ NTSTEP_BETWEEN_READ_ADJSRC = 0 # 0 = read the whole adjoint sources at
# use a force source (or several) 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, inclined force, impact etc.
-USE_FORCE_POINT_SOURCE = .false.
# If this flag is turned on, the FORCESOLUTION file must be edited by precising:
# - the corresponding time-shift parameter,
# - the half duration parameter of the source,
@@ -77,6 +83,7 @@ USE_FORCE_POINT_SOURCE = .false.
# - the components of a (non-unitary) direction vector for the force source in the E/N/Z_UP basis.
# The direction vector is made unitary internally in the code and thus only its direction matters here;
# its norm is ignored and the norm of the force used is the factor force source times the source time function.
+USE_FORCE_POINT_SOURCE = .false.
# set to use a Ricker source time function instead of a gaussian (when using CMTSOLUTION files)
# or a quasi-Heaviside (when using FORCESOLUTION files)
diff --git a/waterlayered_halfspace/in_data_files/Par_file b/waterlayered_halfspace/in_data_files/Par_file
index d078829..103186e 100644
--- a/waterlayered_halfspace/in_data_files/Par_file
+++ b/waterlayered_halfspace/in_data_files/Par_file
@@ -18,6 +18,13 @@ NPROC = 4
NSTEP = 4500
DT = 0.005
+# number of nodes for 2D and 3D shape functions for hexahedra
+# we use either 8-node mesh elements (bricks) or 27-node elements.
+# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
+# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
+# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
+NGNOD = 8
+
# models:
# available options are:
# default (model parameters described by mesh properties)
@@ -68,7 +75,6 @@ NTSTEP_BETWEEN_READ_ADJSRC = 0 # 0 = read the whole adjoint sources at
# use a force source (or several) 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, inclined force, impact etc.
-USE_FORCE_POINT_SOURCE = .false.
# If this flag is turned on, the FORCESOLUTION file must be edited by precising:
# - the corresponding time-shift parameter,
# - the half duration parameter of the source,
@@ -77,6 +83,7 @@ USE_FORCE_POINT_SOURCE = .false.
# - the components of a (non-unitary) direction vector for the force source in the E/N/Z_UP basis.
# The direction vector is made unitary internally in the code and thus only its direction matters here;
# its norm is ignored and the norm of the force used is the factor force source times the source time function.
+USE_FORCE_POINT_SOURCE = .false.
# set to use a Ricker source time function instead of a gaussian (when using CMTSOLUTION files)
# or a quasi-Heaviside (when using FORCESOLUTION files)
More information about the CIG-COMMITS
mailing list