[cig-commits] [commit] devel, master: Updated Par_file of fault-examples to new format (25d4c46)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Wed Jun 18 16:50:46 PDT 2014
Repository : https://github.com/geodynamics/specfem3d
On branches: devel,master
Link : https://github.com/geodynamics/specfem3d/compare/6026e367984905ab133865f62fa6293b343759b9...47f703851338234f96397e7da9fbff63d8178b8a
>---------------------------------------------------------------
commit 25d4c4647951456afdb8d68a11caaa1f696d1727
Author: Surendra Somala <surendra at caltech.edu>
Date: Sat Dec 8 23:07:30 2012 +0000
Updated Par_file of fault-examples to new format
>---------------------------------------------------------------
25d4c4647951456afdb8d68a11caaa1f696d1727
splay_faults/DATA/Par_file | 79 ++++++++++++++++++++++++++++++++++++-------
tpv102/DATA/Par_file | 83 ++++++++++++++++++++++++++++++++++++++--------
tpv103/DATA/Par_file | 81 ++++++++++++++++++++++++++++++++++++--------
tpv15/DATA/Par_file | 81 ++++++++++++++++++++++++++++++++++++--------
tpv16/DATA/Par_file | 83 ++++++++++++++++++++++++++++++++++++++--------
tpv5/DATA/Par_file | 81 ++++++++++++++++++++++++++++++++++++--------
6 files changed, 409 insertions(+), 79 deletions(-)
diff --git a/splay_faults/DATA/Par_file b/splay_faults/DATA/Par_file
index 8c6ef2a..4c4c946 100644
--- a/splay_faults/DATA/Par_file
+++ b/splay_faults/DATA/Par_file
@@ -1,6 +1,10 @@
-
+# simulation input parameters
+#
# forward or adjoint simulation
-SIMULATION_TYPE = 1 # 1 = forward, 2 = adjoint, 3 = both simultaneously
+# 1 = forward, 2 = adjoint, 3 = both simultaneously
+SIMULATION_TYPE = 1
+# 0 = earthquake simulation, 1/2/3 = three steps in noise simulation
+NOISE_TOMOGRAPHY = 0
SAVE_FORWARD = .false.
# UTM projection parameters
@@ -14,29 +18,59 @@ NPROC = 128
NSTEP = 10000
DT = 0.00325d0
+# 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)
+# 1D models available are:
+# 1d_prem,1d_socal,1d_cascadia
+# 3D models available are:
+# aniso,external,gll,salton_trough,tomo
+MODEL = default
+
# parameters describing the model
OCEANS = .false.
+TOPOGRAPHY = .false.
ATTENUATION = .false.
-USE_OLSEN_ATTENUATION = .false.
ANISOTROPY = .false.
+GRAVITY = .false.
+
+# path for external tomographic models files
+TOMOGRAPHY_PATH = ../DATA/tomo_files/
+
+# Olsen's constant for Q_mu = constant * v_s attenuation rule
+USE_OLSEN_ATTENUATION = .false.
+OLSEN_ATTENUATION_RATIO = 0.05
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .true.
-ABSORB_FREE_SURFACE = .false.
+
+# absorbing top surface (defined in mesh as 'free_surface_file')
+ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
+
# save AVS or OpenDX movies
-MOVIE_SURFACE = .true.
-MOVIE_VOLUME = .false.
-NTSTEP_BETWEEN_FRAMES = 16500
+# MOVIE_TYPE = 1 to show the top surface
+# MOVIE_TYPE = 2 to show all the external faces of the mesh
CREATE_SHAKEMAP = .false.
-SAVE_DISPLACEMENT = .true.
+MOVIE_SURFACE = .false.
+MOVIE_TYPE = 1
+MOVIE_VOLUME = .false.
+SAVE_DISPLACEMENT = .false.
USE_HIGHRES_FOR_MOVIES = .false.
+NTSTEP_BETWEEN_FRAMES = 200
HDUR_MOVIE = 0.0
# save AVS or OpenDX mesh files to check the mesh
SAVE_MESH_FILES = .false.
# path to store the local database file on each node
-LOCAL_PATH = /cluster/work/scr1/gpercy
+LOCAL_PATH = ../OUTPUT_FILES/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 500
@@ -44,10 +78,31 @@ NTSTEP_BETWEEN_OUTPUT_INFO = 500
# interval in time steps for writing of seismograms
NTSTEP_BETWEEN_OUTPUT_SEISMOS = 10000
+# interval in time steps for reading adjoint traces
+# 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.
+# 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:
+# - the corresponding time-shift parameter,
+# - the half duration parameter of the source,
+# - the coordinates of the source,
+# - the magnitude of the force source,
+# - 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 true to use a Ricker source time function instead of the source time functions set by default
+# to represent a (tilted) FORCESOLUTION force point source or a CMTSOLUTION moment-tensor source.
+USE_RICKER_TIME_FUNCTION = .false.
+
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .false.
-## Here defined tags that will be taken for faults.
-## which friction laws .
-
+# set to true to use GPUs
+GPU_MODE = .false.
diff --git a/tpv102/DATA/Par_file b/tpv102/DATA/Par_file
index 1c28149..64ca619 100644
--- a/tpv102/DATA/Par_file
+++ b/tpv102/DATA/Par_file
@@ -1,6 +1,10 @@
-
+# simulation input parameters
+#
# forward or adjoint simulation
-SIMULATION_TYPE = 1 # 1 = forward, 2 = adjoint, 3 = both simultaneously
+# 1 = forward, 2 = adjoint, 3 = both simultaneously
+SIMULATION_TYPE = 1
+# 0 = earthquake simulation, 1/2/3 = three steps in noise simulation
+NOISE_TOMOGRAPHY = 0
SAVE_FORWARD = .false.
# UTM projection parameters
@@ -8,46 +12,97 @@ UTM_PROJECTION_ZONE = 11
SUPPRESS_UTM_PROJECTION = .true.
# number of MPI processors
-NPROC = 48
+NPROC = 192
# time step parameters
NSTEP = 4800
DT = 0.0025d0
+# 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)
+# 1D models available are:
+# 1d_prem,1d_socal,1d_cascadia
+# 3D models available are:
+# aniso,external,gll,salton_trough,tomo
+MODEL = default
+
# parameters describing the model
OCEANS = .false.
+TOPOGRAPHY = .false.
ATTENUATION = .false.
-USE_OLSEN_ATTENUATION = .false.
ANISOTROPY = .false.
+GRAVITY = .false.
+
+# path for external tomographic models files
+TOMOGRAPHY_PATH = ../DATA/tomo_files/
+
+# Olsen's constant for Q_mu = constant * v_s attenuation rule
+USE_OLSEN_ATTENUATION = .false.
+OLSEN_ATTENUATION_RATIO = 0.05
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .true.
-ABSORB_FREE_SURFACE = .false.
+
+# absorbing top surface (defined in mesh as 'free_surface_file')
+ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
+
# 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
+CREATE_SHAKEMAP = .false.
MOVIE_SURFACE = .false.
+MOVIE_TYPE = 1
MOVIE_VOLUME = .false.
-NTSTEP_BETWEEN_FRAMES = 100
-CREATE_SHAKEMAP = .false.
-SAVE_DISPLACEMENT = .true.
+SAVE_DISPLACEMENT = .false.
USE_HIGHRES_FOR_MOVIES = .false.
+NTSTEP_BETWEEN_FRAMES = 200
HDUR_MOVIE = 0.0
# save AVS or OpenDX mesh files to check the mesh
SAVE_MESH_FILES = .false.
# path to store the local database file on each node
-LOCAL_PATH = /home/surendra/SCEC_benchmarks/TPV102_DATABASES_MPI_48proc
+LOCAL_PATH = ../OUTPUT_FILES/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
-NTSTEP_BETWEEN_OUTPUT_INFO = 200
+NTSTEP_BETWEEN_OUTPUT_INFO = 500
# interval in time steps for writing of seismograms
-NTSTEP_BETWEEN_OUTPUT_SEISMOS = 20000
+NTSTEP_BETWEEN_OUTPUT_SEISMOS = 10000
+
+# interval in time steps for reading adjoint traces
+# 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.
+# 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:
+# - the corresponding time-shift parameter,
+# - the half duration parameter of the source,
+# - the coordinates of the source,
+# - the magnitude of the force source,
+# - 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 true to use a Ricker source time function instead of the source time functions set by default
+# to represent a (tilted) FORCESOLUTION force point source or a CMTSOLUTION moment-tensor source.
+USE_RICKER_TIME_FUNCTION = .false.
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .false.
-## Here defined tags that will be taken for faults.
-## which friction laws .
-
+# set to true to use GPUs
+GPU_MODE = .false.
diff --git a/tpv103/DATA/Par_file b/tpv103/DATA/Par_file
index b8d14b5..64ca619 100644
--- a/tpv103/DATA/Par_file
+++ b/tpv103/DATA/Par_file
@@ -1,6 +1,10 @@
-
+# simulation input parameters
+#
# forward or adjoint simulation
-SIMULATION_TYPE = 1 # 1 = forward, 2 = adjoint, 3 = both simultaneously
+# 1 = forward, 2 = adjoint, 3 = both simultaneously
+SIMULATION_TYPE = 1
+# 0 = earthquake simulation, 1/2/3 = three steps in noise simulation
+NOISE_TOMOGRAPHY = 0
SAVE_FORWARD = .false.
# UTM projection parameters
@@ -14,40 +18,91 @@ NPROC = 192
NSTEP = 4800
DT = 0.0025d0
+# 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)
+# 1D models available are:
+# 1d_prem,1d_socal,1d_cascadia
+# 3D models available are:
+# aniso,external,gll,salton_trough,tomo
+MODEL = default
+
# parameters describing the model
OCEANS = .false.
+TOPOGRAPHY = .false.
ATTENUATION = .false.
-USE_OLSEN_ATTENUATION = .false.
ANISOTROPY = .false.
+GRAVITY = .false.
+
+# path for external tomographic models files
+TOMOGRAPHY_PATH = ../DATA/tomo_files/
+
+# Olsen's constant for Q_mu = constant * v_s attenuation rule
+USE_OLSEN_ATTENUATION = .false.
+OLSEN_ATTENUATION_RATIO = 0.05
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .true.
-ABSORB_FREE_SURFACE = .true.
+
+# absorbing top surface (defined in mesh as 'free_surface_file')
+ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
+
# 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
+CREATE_SHAKEMAP = .false.
MOVIE_SURFACE = .false.
+MOVIE_TYPE = 1
MOVIE_VOLUME = .false.
-NTSTEP_BETWEEN_FRAMES = 100
-CREATE_SHAKEMAP = .false.
-SAVE_DISPLACEMENT = .true.
+SAVE_DISPLACEMENT = .false.
USE_HIGHRES_FOR_MOVIES = .false.
+NTSTEP_BETWEEN_FRAMES = 200
HDUR_MOVIE = 0.0
# save AVS or OpenDX mesh files to check the mesh
SAVE_MESH_FILES = .false.
# path to store the local database file on each node
-LOCAL_PATH = /home/surendra/SCEC_benchmarks/TPV101_DATABASES_MPI_serial_Fault_EXT2ABS_larger
+LOCAL_PATH = ../OUTPUT_FILES/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
-NTSTEP_BETWEEN_OUTPUT_INFO = 200
+NTSTEP_BETWEEN_OUTPUT_INFO = 500
# interval in time steps for writing of seismograms
-NTSTEP_BETWEEN_OUTPUT_SEISMOS = 20000
+NTSTEP_BETWEEN_OUTPUT_SEISMOS = 10000
+
+# interval in time steps for reading adjoint traces
+# 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.
+# 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:
+# - the corresponding time-shift parameter,
+# - the half duration parameter of the source,
+# - the coordinates of the source,
+# - the magnitude of the force source,
+# - 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 true to use a Ricker source time function instead of the source time functions set by default
+# to represent a (tilted) FORCESOLUTION force point source or a CMTSOLUTION moment-tensor source.
+USE_RICKER_TIME_FUNCTION = .false.
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .false.
-## Here defined tags that will be taken for faults.
-## which friction laws .
-
+# set to true to use GPUs
+GPU_MODE = .false.
diff --git a/tpv15/DATA/Par_file b/tpv15/DATA/Par_file
index 644b597..0182a3f 100644
--- a/tpv15/DATA/Par_file
+++ b/tpv15/DATA/Par_file
@@ -1,6 +1,10 @@
-
+# simulation input parameters
+#
# forward or adjoint simulation
-SIMULATION_TYPE = 1 # 1 = forward, 2 = adjoint, 3 = both simultaneously
+# 1 = forward, 2 = adjoint, 3 = both simultaneously
+SIMULATION_TYPE = 1
+# 0 = earthquake simulation, 1/2/3 = three steps in noise simulation
+NOISE_TOMOGRAPHY = 0
SAVE_FORWARD = .false.
# UTM projection parameters
@@ -14,29 +18,59 @@ NPROC = 100
NSTEP = 12000
DT = 0.001d0
+# 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)
+# 1D models available are:
+# 1d_prem,1d_socal,1d_cascadia
+# 3D models available are:
+# aniso,external,gll,salton_trough,tomo
+MODEL = default
+
# parameters describing the model
OCEANS = .false.
+TOPOGRAPHY = .false.
ATTENUATION = .false.
-USE_OLSEN_ATTENUATION = .false.
ANISOTROPY = .false.
+GRAVITY = .false.
+
+# path for external tomographic models files
+TOMOGRAPHY_PATH = ../DATA/tomo_files/
+
+# Olsen's constant for Q_mu = constant * v_s attenuation rule
+USE_OLSEN_ATTENUATION = .false.
+OLSEN_ATTENUATION_RATIO = 0.05
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .true.
-ABSORB_FREE_SURFACE = .false.
+
+# absorbing top surface (defined in mesh as 'free_surface_file')
+ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
+
# save AVS or OpenDX movies
-MOVIE_SURFACE = .true.
-MOVIE_VOLUME = .false.
-NTSTEP_BETWEEN_FRAMES = 6000
+# MOVIE_TYPE = 1 to show the top surface
+# MOVIE_TYPE = 2 to show all the external faces of the mesh
CREATE_SHAKEMAP = .false.
-SAVE_DISPLACEMENT = .true.
+MOVIE_SURFACE = .false.
+MOVIE_TYPE = 1
+MOVIE_VOLUME = .false.
+SAVE_DISPLACEMENT = .false.
USE_HIGHRES_FOR_MOVIES = .false.
+NTSTEP_BETWEEN_FRAMES = 200
HDUR_MOVIE = 0.0
# save AVS or OpenDX mesh files to check the mesh
-SAVE_MESH_FILES = .true.
+SAVE_MESH_FILES = .false.
# path to store the local database file on each node
-LOCAL_PATH = /cluster/work/erdw/gpercy
+LOCAL_PATH = ../OUTPUT_FILES/DATABASES_MPI
# interval at which we output time step info and max of norm of displacement
NTSTEP_BETWEEN_OUTPUT_INFO = 500
@@ -44,10 +78,31 @@ NTSTEP_BETWEEN_OUTPUT_INFO = 500
# interval in time steps for writing of seismograms
NTSTEP_BETWEEN_OUTPUT_SEISMOS = 10000
+# interval in time steps for reading adjoint traces
+# 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.
+# 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:
+# - the corresponding time-shift parameter,
+# - the half duration parameter of the source,
+# - the coordinates of the source,
+# - the magnitude of the force source,
+# - 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 true to use a Ricker source time function instead of the source time functions set by default
+# to represent a (tilted) FORCESOLUTION force point source or a CMTSOLUTION moment-tensor source.
+USE_RICKER_TIME_FUNCTION = .false.
+
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .false.
-## Here defined tags that will be taken for faults.
-## which friction laws .
-
+# set to true to use GPUs
+GPU_MODE = .false.
diff --git a/tpv16/DATA/Par_file b/tpv16/DATA/Par_file
index 43f3fd3..d576134 100644
--- a/tpv16/DATA/Par_file
+++ b/tpv16/DATA/Par_file
@@ -1,6 +1,10 @@
-
+# simulation input parameters
+#
# forward or adjoint simulation
-SIMULATION_TYPE = 1 # 1 = forward, 2 = adjoint, 3 = both simultaneously
+# 1 = forward, 2 = adjoint, 3 = both simultaneously
+SIMULATION_TYPE = 1
+# 0 = earthquake simulation, 1/2/3 = three steps in noise simulation
+NOISE_TOMOGRAPHY = 0
SAVE_FORWARD = .false.
# UTM projection parameters
@@ -8,46 +12,97 @@ UTM_PROJECTION_ZONE = 11
SUPPRESS_UTM_PROJECTION = .true.
# number of MPI processors
-NPROC = 400
+NPROC = 200
# time step parameters
NSTEP = 6000
DT = 0.0025d0
+# 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)
+# 1D models available are:
+# 1d_prem,1d_socal,1d_cascadia
+# 3D models available are:
+# aniso,external,gll,salton_trough,tomo
+MODEL = default
+
# parameters describing the model
OCEANS = .false.
+TOPOGRAPHY = .false.
ATTENUATION = .false.
-USE_OLSEN_ATTENUATION = .false.
ANISOTROPY = .false.
+GRAVITY = .false.
+
+# path for external tomographic models files
+TOMOGRAPHY_PATH = ../DATA/tomo_files/
+
+# Olsen's constant for Q_mu = constant * v_s attenuation rule
+USE_OLSEN_ATTENUATION = .false.
+OLSEN_ATTENUATION_RATIO = 0.05
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .true.
-ABSORB_FREE_SURFACE = .false.
+
+# absorbing top surface (defined in mesh as 'free_surface_file')
+ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
+
# 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
+CREATE_SHAKEMAP = .false.
MOVIE_SURFACE = .false.
+MOVIE_TYPE = 1
MOVIE_VOLUME = .false.
-NTSTEP_BETWEEN_FRAMES = 100
-CREATE_SHAKEMAP = .false.
-SAVE_DISPLACEMENT = .true.
+SAVE_DISPLACEMENT = .false.
USE_HIGHRES_FOR_MOVIES = .false.
+NTSTEP_BETWEEN_FRAMES = 200
HDUR_MOVIE = 0.0
# save AVS or OpenDX mesh files to check the mesh
SAVE_MESH_FILES = .false.
# path to store the local database file on each node
-LOCAL_PATH = 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 = 200
+NTSTEP_BETWEEN_OUTPUT_INFO = 500
# interval in time steps for writing of seismograms
-NTSTEP_BETWEEN_OUTPUT_SEISMOS = 20000
+NTSTEP_BETWEEN_OUTPUT_SEISMOS = 10000
+
+# interval in time steps for reading adjoint traces
+# 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.
+# 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:
+# - the corresponding time-shift parameter,
+# - the half duration parameter of the source,
+# - the coordinates of the source,
+# - the magnitude of the force source,
+# - 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 true to use a Ricker source time function instead of the source time functions set by default
+# to represent a (tilted) FORCESOLUTION force point source or a CMTSOLUTION moment-tensor source.
+USE_RICKER_TIME_FUNCTION = .false.
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .false.
-## Here defined tags that will be taken for faults.
-## which friction laws .
-
+# set to true to use GPUs
+GPU_MODE = .false.
diff --git a/tpv5/DATA/Par_file b/tpv5/DATA/Par_file
index ad0ba74..04cfdb3 100644
--- a/tpv5/DATA/Par_file
+++ b/tpv5/DATA/Par_file
@@ -1,6 +1,10 @@
-
+# simulation input parameters
+#
# forward or adjoint simulation
-SIMULATION_TYPE = 1 # 1 = forward, 2 = adjoint, 3 = both simultaneously
+# 1 = forward, 2 = adjoint, 3 = both simultaneously
+SIMULATION_TYPE = 1
+# 0 = earthquake simulation, 1/2/3 = three steps in noise simulation
+NOISE_TOMOGRAPHY = 0
SAVE_FORWARD = .false.
# UTM projection parameters
@@ -8,35 +12,65 @@ UTM_PROJECTION_ZONE = 11
SUPPRESS_UTM_PROJECTION = .true.
# number of MPI processors
-NPROC = 100
+NPROC = 120
# time step parameters
-NSTEP = 12000
-DT = 0.001d0
+NSTEP = 4000
+DT = 0.002d0
+
+# 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)
+# 1D models available are:
+# 1d_prem,1d_socal,1d_cascadia
+# 3D models available are:
+# aniso,external,gll,salton_trough,tomo
+MODEL = default
# parameters describing the model
OCEANS = .false.
+TOPOGRAPHY = .false.
ATTENUATION = .false.
-USE_OLSEN_ATTENUATION = .false.
ANISOTROPY = .false.
+GRAVITY = .false.
+
+# path for external tomographic models files
+TOMOGRAPHY_PATH = ../DATA/tomo_files/
+
+# Olsen's constant for Q_mu = constant * v_s attenuation rule
+USE_OLSEN_ATTENUATION = .false.
+OLSEN_ATTENUATION_RATIO = 0.05
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .true.
-ABSORB_FREE_SURFACE = .false.
+
+# absorbing top surface (defined in mesh as 'free_surface_file')
+ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
+
# 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
+CREATE_SHAKEMAP = .false.
MOVIE_SURFACE = .false.
+MOVIE_TYPE = 1
MOVIE_VOLUME = .false.
-NTSTEP_BETWEEN_FRAMES = 200
-CREATE_SHAKEMAP = .false.
SAVE_DISPLACEMENT = .false.
USE_HIGHRES_FOR_MOVIES = .false.
+NTSTEP_BETWEEN_FRAMES = 200
HDUR_MOVIE = 0.0
# save AVS or OpenDX mesh files to check the mesh
SAVE_MESH_FILES = .false.
# path to store the local database file on each node
-LOCAL_PATH = /home/galvez/FAULT_SOURCE/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 = 500
@@ -44,10 +78,31 @@ NTSTEP_BETWEEN_OUTPUT_INFO = 500
# interval in time steps for writing of seismograms
NTSTEP_BETWEEN_OUTPUT_SEISMOS = 10000
+# interval in time steps for reading adjoint traces
+# 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.
+# 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:
+# - the corresponding time-shift parameter,
+# - the half duration parameter of the source,
+# - the coordinates of the source,
+# - the magnitude of the force source,
+# - 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 true to use a Ricker source time function instead of the source time functions set by default
+# to represent a (tilted) FORCESOLUTION force point source or a CMTSOLUTION moment-tensor source.
+USE_RICKER_TIME_FUNCTION = .false.
+
# print source time function
PRINT_SOURCE_TIME_FUNCTION = .false.
-## Here defined tags that will be taken for faults.
-## which friction laws .
-
+# set to true to use GPUs
+GPU_MODE = .false.
More information about the CIG-COMMITS
mailing list