Skip to content

Commit

Permalink
corrected Par_files and MPI
Browse files Browse the repository at this point in the history
  • Loading branch information
@cmorency1
cmorency1 committed May 23, 2024
1 parent 09d698a commit 38dd3fd
Show file tree
Hide file tree
Showing 84 changed files with 2,582 additions and 2,472 deletions.
62 changes: 32 additions & 30 deletions DATA/Par_file
View file
Original file line number Diff line number Diff line change
@@ -1,8 +1,8 @@
#-----------------------------------------------------------------------------
#-----------------------------------------------------------
#
# simulation input parameters
# Simulation input parameters
#
#-----------------------------------------------------------------------------
#-----------------------------------------------------------

# title of job
title = Test of SPECFEM2D with curved interfaces
Expand All @@ -17,7 +17,7 @@ NOISE_TOMOGRAPHY = 0
SAVE_FORWARD = .false.

# parameters concerning partitioning
NPROC = 1
NPROC = 1 # number of processes

# time step parameters
# total number of time steps
Expand All @@ -30,7 +30,7 @@ DT = 1.1d-3
# 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta
time_stepping_scheme = 1

# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM))
# set the type of calculation (P-SV or SH/membrane waves)
P_SV = .true.

# axisymmetric (2.5D) or Cartesian planar (2D) simulation
Expand Down Expand Up @@ -72,20 +72,20 @@ SAVE_MODEL = default
# (used for tomography materials with negative material ids and/or MODEL==tomo settings)
TOMOGRAPHY_FILE = dummy

#-----------------------------------------------------------------------------
#-----------------------------------------------------------
#
# Attenuation
#
#-----------------------------------------------------------------------------
#-----------------------------------------------------------

# attenuation parameters
ATTENUATION_VISCOELASTIC = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
ATTENUATION_VISCOACOUSTIC = .false. # turn attenuation (viscoacousticity) on or off for non-poroelastic fluid parts of the model

# for viscoelastic attenuation
# for viscoelastic or viscoacoustic attenuation
N_SLS = 3 # number of standard linear solids for attenuation (3 is usually the minimum)
ATTENUATION_f0_REFERENCE = 5.196 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file
READ_VELOCITIES_AT_f0 = .false. # shift velocities to account for physical dispersion (see user manual for more information)
ATTENUATION_f0_REFERENCE = 5.196 # in case of attenuation, reference frequency in Hz at which the velocity values in the velocity model are given (unused otherwise); relevant only if source is a Dirac or a Heaviside, otherwise it is automatically set to f0 the dominant frequency of the source in the DATA/SOURCE file
READ_VELOCITIES_AT_f0 = .false. # read seismic velocities at ATTENUATION_f0_REFERENCE instead of at infinite frequency (see user manual for more information)
USE_SOLVOPT = .false. # use more precise but much more expensive way of determining the Q factor relaxation times, as in https://doi.org/10.1093/gji/ggw024

# for poroelastic attenuation
Expand All @@ -98,8 +98,8 @@ ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener mode
ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material
f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file

# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD
# use the flag below. It performs undoing of attenuation in an exact way for sensitivity kernel calculations
# to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD
# use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations
# but requires disk space for temporary storage, and uses a significant amount of memory used as buffers for temporary storage.
# When that option is on the second parameter indicates how often the code dumps restart files to disk (if in doubt, use something between 100 and 1000).
UNDO_ATTENUATION_AND_OR_PML = .false.
Expand All @@ -109,11 +109,11 @@ NT_DUMP_ATTENUATION = 500
# Outperforms conventional mode using a value of NTSTEP_BETWEEN_COMPUTE_KERNELS high enough.
NO_BACKWARD_RECONSTRUCTION = .false.

#-----------------------------------------------------------------------------
#-----------------------------------------------------------
#
# Sources
#
#-----------------------------------------------------------------------------
#-----------------------------------------------------------

# source parameters
NSOURCES = 1 # number of sources (source information is then read from the DATA/SOURCE file)
Expand Down Expand Up @@ -142,14 +142,14 @@ noise_source_time_function_type = 4
# This is only useful for GPU version (for now)
write_moving_sources_database = .false.

#-----------------------------------------------------------------------------
#-----------------------------------------------------------
#
# Receivers
#
#-----------------------------------------------------------------------------
#-----------------------------------------------------------

# receiver set parameters for recording stations (i.e. recording points)
# record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential 7=electric field
# seismotype : record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential 7=electric field
seismotype = 1 # several values can be chosen. For example : 1,2,4

# interval in time steps for writing of seismograms
Expand Down Expand Up @@ -201,11 +201,11 @@ zfin = 0. # last receiver z in meters (ig
record_at_surface_same_vertical = .false. # receivers inside the medium or at the surface (z values are ignored if this is set to true, they are replaced with the topography height)


#-----------------------------------------------------------------------------
#-----------------------------------------------------------
#
# Adjoint kernel outputs
# adjoint kernel outputs
#
#-----------------------------------------------------------------------------
#-----------------------------------------------------------

# save sensitivity kernels in ASCII format (much bigger files, but compatible with current GMT scripts) or in binary format
save_ASCII_kernels = .true.
Expand All @@ -216,11 +216,11 @@ NTSTEP_BETWEEN_COMPUTE_KERNELS = 1
# outputs approximate Hessian for preconditioning
APPROXIMATE_HESS_KL = .false.

#-----------------------------------------------------------------------------
#-----------------------------------------------------------
#
# Boundary conditions
#
#-----------------------------------------------------------------------------
#-----------------------------------------------------------

# Perfectly Matched Layer (PML) boundaries
# absorbing boundary active or not
Expand All @@ -245,11 +245,11 @@ STACEY_ABSORBING_CONDITIONS = .false.
ADD_PERIODIC_CONDITIONS = .false.
PERIODIC_HORIZ_DIST = 4000.d0

#-----------------------------------------------------------------------------
#-----------------------------------------------------------
#
# MESHING - Velocity and density models
#
#-----------------------------------------------------------------------------
#-----------------------------------------------------------

# use an external mesh created by an external meshing tool or use the internal mesher
read_external_mesh = .false.
Expand Down Expand Up @@ -278,6 +278,7 @@ tangential_detection_curve_file = dummy # file
# PARAMETERS FOR INTERNAL MESHING
#
#-----------------------------------------------------------

# material properties
# number of model materials
nbmodels = 4
Expand All @@ -287,6 +288,7 @@ nbmodels = 4
# anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu
# anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu
# poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu
# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0
# tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0
#
# note: When viscoelasticity or viscoacousticity is turned on,
Expand Down Expand Up @@ -324,11 +326,11 @@ nbregions = 5 # then set below the different
1 80 41 60 3
60 70 21 40 4

#-----------------------------------------------------------------------------
#-----------------------------------------------------------
#
# Display parameters
#
#-----------------------------------------------------------------------------
#-----------------------------------------------------------

# interval at which we output time step info and max of norm of displacement
# (every how many time steps we display information about the simulation. costly, do not use a very small value)
Expand All @@ -349,11 +351,11 @@ NTSTEP_BETWEEN_OUTPUT_ENERGY = 10
# the script utils/visualisation/plotIntegratedEnergyFile.py to watch. It is refreshed at the same time than the seismograms
COMPUTE_INTEGRATED_ENERGY_FIELD = .false.

#-----------------------------------------------------------------------------
#-----------------------------------------------------------
#
# Movies/images/snaphots visualizations
#
#-----------------------------------------------------------------------------
#-----------------------------------------------------------

# every how many time steps we draw JPEG or PostScript pictures of the simulation
# and/or we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
Expand Down Expand Up @@ -406,9 +408,9 @@ use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision
# When that option is on, of course the number of processor cores used to start
# the code in the batch system must be a multiple of NUMBER_OF_SIMULTANEOUS_RUNS,
# all the individual runs must use the same number of processor cores,
# which as usual is proc in the Par_file,
# which as usual is NPROC in the Par_file,
# and thus the total number of processor cores to request from the batch system
# should be NUMBER_OF_SIMULTANEOUS_RUNS * proc.
# should be NUMBER_OF_SIMULTANEOUS_RUNS * NPROC.
# All the runs to perform must be placed in directories called run0001, run0002, run0003 and so on
# (with exactly four digits).
#
Expand Down
Loading

0 comments on commit 38dd3fd

Please sign in to comment.