[cig-commits] r21112 - in seismo/3D/SPECFEM3D/trunk/examples: splay_faults/DATA tpv102/DATA tpv103/DATA tpv15/DATA tpv16/DATA tpv5/DATA
surendra at geodynamics.org
surendra at geodynamics.org
Sat Dec 8 15:07:30 PST 2012
Author: surendra
Date: 2012-12-08 15:07:30 -0800 (Sat, 08 Dec 2012)
New Revision: 21112
Modified:
seismo/3D/SPECFEM3D/trunk/examples/splay_faults/DATA/Par_file
seismo/3D/SPECFEM3D/trunk/examples/tpv102/DATA/Par_file
seismo/3D/SPECFEM3D/trunk/examples/tpv103/DATA/Par_file
seismo/3D/SPECFEM3D/trunk/examples/tpv15/DATA/Par_file
seismo/3D/SPECFEM3D/trunk/examples/tpv16/DATA/Par_file
seismo/3D/SPECFEM3D/trunk/examples/tpv5/DATA/Par_file
Log:
Updated Par_file of fault-examples to new format
Modified: seismo/3D/SPECFEM3D/trunk/examples/splay_faults/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/splay_faults/DATA/Par_file 2012-12-08 22:52:00 UTC (rev 21111)
+++ seismo/3D/SPECFEM3D/trunk/examples/splay_faults/DATA/Par_file 2012-12-08 23:07:30 UTC (rev 21112)
@@ -1,42 +1,76 @@
-
+# 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
UTM_PROJECTION_ZONE = 11
SUPPRESS_UTM_PROJECTION = .true.
-# number of MPI processors
-NPROC = 128
+# number of MPI processors
+NPROC = 128
# time step parameters
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_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 = 16500
-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 = /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 @@
# 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.
-
Modified: seismo/3D/SPECFEM3D/trunk/examples/tpv102/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/tpv102/DATA/Par_file 2012-12-08 22:52:00 UTC (rev 21111)
+++ seismo/3D/SPECFEM3D/trunk/examples/tpv102/DATA/Par_file 2012-12-08 23:07:30 UTC (rev 21112)
@@ -1,53 +1,108 @@
-
+# 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
UTM_PROJECTION_ZONE = 11
SUPPRESS_UTM_PROJECTION = .true.
-# number of MPI processors
-NPROC = 48
+# number of MPI processors
+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.
-
Modified: seismo/3D/SPECFEM3D/trunk/examples/tpv103/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/tpv103/DATA/Par_file 2012-12-08 22:52:00 UTC (rev 21111)
+++ seismo/3D/SPECFEM3D/trunk/examples/tpv103/DATA/Par_file 2012-12-08 23:07:30 UTC (rev 21112)
@@ -1,53 +1,108 @@
-
+# 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
UTM_PROJECTION_ZONE = 11
SUPPRESS_UTM_PROJECTION = .true.
-# number of MPI processors
+# number of MPI processors
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 = .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.
-
Modified: seismo/3D/SPECFEM3D/trunk/examples/tpv15/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/tpv15/DATA/Par_file 2012-12-08 22:52:00 UTC (rev 21111)
+++ seismo/3D/SPECFEM3D/trunk/examples/tpv15/DATA/Par_file 2012-12-08 23:07:30 UTC (rev 21112)
@@ -1,42 +1,76 @@
-
+# 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
UTM_PROJECTION_ZONE = 11
SUPPRESS_UTM_PROJECTION = .true.
-# number of MPI processors
-NPROC = 100
+# number of MPI processors
+NPROC = 100
# time step parameters
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_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 = 6000
-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 = .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 @@
# 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.
-
Modified: seismo/3D/SPECFEM3D/trunk/examples/tpv16/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/tpv16/DATA/Par_file 2012-12-08 22:52:00 UTC (rev 21111)
+++ seismo/3D/SPECFEM3D/trunk/examples/tpv16/DATA/Par_file 2012-12-08 23:07:30 UTC (rev 21112)
@@ -1,53 +1,108 @@
-
+# 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
UTM_PROJECTION_ZONE = 11
SUPPRESS_UTM_PROJECTION = .true.
-# number of MPI processors
-NPROC = 400
+# number of MPI processors
+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.
-
Modified: seismo/3D/SPECFEM3D/trunk/examples/tpv5/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/tpv5/DATA/Par_file 2012-12-08 22:52:00 UTC (rev 21111)
+++ seismo/3D/SPECFEM3D/trunk/examples/tpv5/DATA/Par_file 2012-12-08 23:07:30 UTC (rev 21112)
@@ -1,42 +1,76 @@
-
+# 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
UTM_PROJECTION_ZONE = 11
SUPPRESS_UTM_PROJECTION = .true.
-# number of MPI processors
-NPROC = 100
+# number of MPI processors
+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 @@
# 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