[cig-commits] [commit] devel: updated utils/DSM_FOR_SPECFEM3D/example/DATA/Par_file to the current format of Par_files (it was obsolete) (a11d1f6)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Tue Jul 8 12:05:49 PDT 2014
Repository : https://github.com/geodynamics/specfem3d
On branch : devel
Link : https://github.com/geodynamics/specfem3d/compare/7869d8b26275f52d5adcb3c1cf3f7054dafe3ec3...a11d1f6d6edb57032ef169be59b4b317d40e9c4b
>---------------------------------------------------------------
commit a11d1f6d6edb57032ef169be59b4b317d40e9c4b
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date: Tue Jul 8 20:54:59 2014 +0200
updated utils/DSM_FOR_SPECFEM3D/example/DATA/Par_file to the current format of Par_files (it was obsolete)
>---------------------------------------------------------------
a11d1f6d6edb57032ef169be59b4b317d40e9c4b
utils/DSM_FOR_SPECFEM3D/example/DATA/Par_file | 101 ++++++++++++++++++++++----
1 file changed, 88 insertions(+), 13 deletions(-)
diff --git a/utils/DSM_FOR_SPECFEM3D/example/DATA/Par_file b/utils/DSM_FOR_SPECFEM3D/example/DATA/Par_file
index af4ed5a..b2d2eb0 100644
--- a/utils/DSM_FOR_SPECFEM3D/example/DATA/Par_file
+++ b/utils/DSM_FOR_SPECFEM3D/example/DATA/Par_file
@@ -1,7 +1,10 @@
-
+# simulation input parameters
+#
# forward or adjoint simulation
-SIMULATION_TYPE = 1 # 1 = forward, 2 = adjoint, 3 = both simultaneously
-NOISE_TOMOGRAPHY = 0 # 0 = earthquake simulation, 1/2/3 = three steps in noise simulation
+# 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 = .true.
# UTM projection parameters
@@ -15,23 +18,66 @@ NPROC = 24
NSTEP = 12000
DT = 0.1
+# 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.
+APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
-USE_OLSEN_ATTENUATION = .false.
+FULL_ATTENUATION_SOLID = .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
+
+# C-PML boundary conditions for a regional simulation
+PML_CONDITIONS = .false.
+
+# C-PML top surface
+PML_INSTEAD_OF_FREE_SURFACE = .false.
+
+# C-PML dominant frequency
+f0_FOR_PML = 0.05555
+
+# parameters used to rotate C-PML boundary conditions by a given angle (not completed yet)
+# ROTATE_PML_ACTIVATE = .false.
+# ROTATE_PML_ANGLE = 0.
# absorbing boundary conditions for a regional simulation
-ABSORBING_CONDITIONS = .true.
+STACEY_ABSORBING_CONDITIONS = .true.
+
+# absorbing top surface (defined in mesh as 'free_surface_file')
+STACEY_INSTEAD_OF_FREE_SURFACE = .false.
# save AVS or OpenDX movies
-MOVIE_SURFACE = .false.
-MOVIE_VOLUME = .true.
-NTSTEP_BETWEEN_FRAMES = 25
+# 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.
SAVE_DISPLACEMENT = .false.
-USE_HIGHRES_FOR_MOVIES = .true.
+USE_HIGHRES_FOR_MOVIES = .false.
+NTSTEP_BETWEEN_FRAMES = 200
HDUR_MOVIE = 0.0
# save AVS or OpenDX mesh files to check the mesh
@@ -47,11 +93,40 @@ NTSTEP_BETWEEN_OUTPUT_INFO = 25
NTSTEP_BETWEEN_OUTPUT_SEISMOS = 100
# interval in time steps for reading adjoint traces
-NTSTEP_BETWEEN_READ_ADJSRC = 0 # 0 = read the whole adjoint sources at the same time
+# 0 = read the whole adjoint sources at the same time
+NTSTEP_BETWEEN_READ_ADJSRC = 0
+
+# use a (tilted) FORCESOLUTION force point source (or several) 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.
-# path to store the local database file on each node
-TRAC_PATH = ../OUTPUT_FILES/Tractions
+# to couple with DSM
+COUPLE_WITH_DSM = .true.
+TRACTION_PATH = ../OUTPUT_FILES/Tractions
+
+# set to true to use GPUs
+GPU_MODE = .false.
+
+# ADIOS Options for I/Os
+ADIOS_ENABLED = .false.
+ADIOS_FOR_DATABASES = .false.
+ADIOS_FOR_MESH = .false.
+ADIOS_FOR_FORWARD_ARRAYS = .false.
+ADIOS_FOR_KERNELS = .false.
More information about the CIG-COMMITS
mailing list