Skip to content

Commit

Permalink
Merge pull request #1204 from mnagaso/devel
Browse files Browse the repository at this point in the history
Gmsh example updated for the latest Gmsh and Specfem2D
  • Loading branch information
danielpeter authored Feb 20, 2024
2 parents c7133e5 + c759143 commit b735263
Show file tree
Hide file tree
Showing 14 changed files with 115,425 additions and 0 deletions.
438 changes: 438 additions & 0 deletions EXAMPLES/Gmsh_updated_example/DATA/Par_file

Large diffs are not rendered by default.

57 changes: 57 additions & 0 deletions EXAMPLES/Gmsh_updated_example/DATA/SOURCE
Original file line number Diff line number Diff line change
@@ -0,0 +1,57 @@
## Source 1
source_surf = .false. # source inside the medium, or source automatically moved exactly at the surface by the solver
xs = 7000.0
zs = 10000.0
## Source type parameters:
# 1 = elastic force or acoustic pressure
# 2 = moment tensor
# or Initial field type (when initialfield set in Par_file):
# For a plane wave including converted and reflected waves at the free surface:
# 1 = P wave,
# 2 = S wave,
# 3 = Rayleigh wave
# For a plane wave without converted nor reflected waves at the free surface, i.e. with the incident wave only:
# 4 = P wave,
# 5 = S wave
# For initial mode displacement:
# 6 = mode (2,3) of a rectangular membrane
source_type = 1
# Source time function:
# In the case of a source located in an acoustic medium,
# to get pressure for a Ricker in the seismograms, here we need to select a Gaussian for the potential Chi
# used as a source, rather than a Ricker, because pressure = - Chi_dot_dot.
# This is true both when USE_TRICK_FOR_BETTER_PRESSURE is set to .true. or to .false.
# Options:
# 1 = second derivative of a Gaussian (a.k.a. Ricker),
# 2 = first derivative of a Gaussian,
# 3 = Gaussian,
# 4 = Dirac,
# 5 = Heaviside (4 and 5 will produce noisy recordings because of frequencies above the mesh resolution limit),
# 6 = ocean acoustics type I,
# 7 = ocean acoustics type II,
# 8 = external source time function = 8 (source read from file),
# 9 = burst,
# 10 = Sinus source time function,
# 11 = Marmousi Ormsby wavelet
time_function_type = 1
# If time_function_type == 8, enter below the custom source file to read (two columns file with time and amplitude) :
# (For the moment dt must be equal to the dt of the simulation. File name cannot exceed 150 characters)
# IMPORTANT: do NOT put quote signs around the file name, just put the file name itself otherwise the run will stop
name_of_source_file = DUMMY
burst_band_width = 0. # Only for option 9 : band width of the burst
f0 = 0.678
tshift = 0.0 # time shift when multi sources (if one source, must be zero)
## Force source
# angle of the source (for a force only); for a plane wave, this is the incidence angle; for moment tensor sources this is unused
anglesource = 90.
## Moment tensor
# The components of a moment tensor source must be given in N.m, not in dyne.cm as in the DATA/CMTSOLUTION source file of the 3D version of the code.
Mxx = 1. # Mxx component (for a moment tensor source only)
Mzz = 1. # Mzz component (for a moment tensor source only)
Mxz = 0. # Mxz component (for a moment tensor source only)
## Amplification (factor to amplify source time function)
factor = 1.d10 # amplification factor
## Moving source parameters
vx = 0.0 # Horizontal source velocity (m/s)
vz = 0.0 # Vertical source velocity (m/s)

245 changes: 245 additions & 0 deletions EXAMPLES/Gmsh_updated_example/DATA/STATIONS
Original file line number Diff line number Diff line change
@@ -0,0 +1,245 @@
x_0 z_0 2000.0 16000.0 0.0 0.0
x_0 z_1 2000.0 17000.0 0.0 0.0
x_0 z_2 2000.0 18000.0 0.0 0.0
x_0 z_3 2000.0 19000.0 0.0 0.0
x_0 z_4 2000.0 20000.0 0.0 0.0
x_1 z_0 4000.0 16000.0 0.0 0.0
x_1 z_1 4000.0 17000.0 0.0 0.0
x_1 z_2 4000.0 18000.0 0.0 0.0
x_1 z_3 4000.0 19000.0 0.0 0.0
x_1 z_4 4000.0 20000.0 0.0 0.0
x_2 z_0 6000.0 16000.0 0.0 0.0
x_2 z_1 6000.0 17000.0 0.0 0.0
x_2 z_2 6000.0 18000.0 0.0 0.0
x_2 z_3 6000.0 19000.0 0.0 0.0
x_2 z_4 6000.0 20000.0 0.0 0.0
x_3 z_0 8000.0 16000.0 0.0 0.0
x_3 z_1 8000.0 17000.0 0.0 0.0
x_3 z_2 8000.0 18000.0 0.0 0.0
x_3 z_3 8000.0 19000.0 0.0 0.0
x_3 z_4 8000.0 20000.0 0.0 0.0
x_4 z_0 10000.0 16000.0 0.0 0.0
x_4 z_1 10000.0 17000.0 0.0 0.0
x_4 z_2 10000.0 18000.0 0.0 0.0
x_4 z_3 10000.0 19000.0 0.0 0.0
x_4 z_4 10000.0 20000.0 0.0 0.0
x_5 z_0 12000.0 16000.0 0.0 0.0
x_5 z_1 12000.0 17000.0 0.0 0.0
x_5 z_2 12000.0 18000.0 0.0 0.0
x_5 z_3 12000.0 19000.0 0.0 0.0
x_5 z_4 12000.0 20000.0 0.0 0.0
x_6 z_0 14000.0 16000.0 0.0 0.0
x_6 z_1 14000.0 17000.0 0.0 0.0
x_6 z_2 14000.0 18000.0 0.0 0.0
x_6 z_3 14000.0 19000.0 0.0 0.0
x_6 z_4 14000.0 20000.0 0.0 0.0
x_7 z_0 16000.0 16000.0 0.0 0.0
x_7 z_1 16000.0 17000.0 0.0 0.0
x_7 z_2 16000.0 18000.0 0.0 0.0
x_7 z_3 16000.0 19000.0 0.0 0.0
x_7 z_4 16000.0 20000.0 0.0 0.0
x_8 z_0 18000.0 16000.0 0.0 0.0
x_8 z_1 18000.0 17000.0 0.0 0.0
x_8 z_2 18000.0 18000.0 0.0 0.0
x_8 z_3 18000.0 19000.0 0.0 0.0
x_8 z_4 18000.0 20000.0 0.0 0.0
x_9 z_0 20000.0 16000.0 0.0 0.0
x_9 z_1 20000.0 17000.0 0.0 0.0
x_9 z_2 20000.0 18000.0 0.0 0.0
x_9 z_3 20000.0 19000.0 0.0 0.0
x_9 z_4 20000.0 20000.0 0.0 0.0
x_10 z_0 22000.0 16000.0 0.0 0.0
x_10 z_1 22000.0 17000.0 0.0 0.0
x_10 z_2 22000.0 18000.0 0.0 0.0
x_10 z_3 22000.0 19000.0 0.0 0.0
x_10 z_4 22000.0 20000.0 0.0 0.0
x_11 z_0 24000.0 16000.0 0.0 0.0
x_11 z_1 24000.0 17000.0 0.0 0.0
x_11 z_2 24000.0 18000.0 0.0 0.0
x_11 z_3 24000.0 19000.0 0.0 0.0
x_11 z_4 24000.0 20000.0 0.0 0.0
x_12 z_0 26000.0 16000.0 0.0 0.0
x_12 z_1 26000.0 17000.0 0.0 0.0
x_12 z_2 26000.0 18000.0 0.0 0.0
x_12 z_3 26000.0 19000.0 0.0 0.0
x_12 z_4 26000.0 20000.0 0.0 0.0
x_13 z_0 28000.0 16000.0 0.0 0.0
x_13 z_1 28000.0 17000.0 0.0 0.0
x_13 z_2 28000.0 18000.0 0.0 0.0
x_13 z_3 28000.0 19000.0 0.0 0.0
x_13 z_4 28000.0 20000.0 0.0 0.0
x_14 z_0 30000.0 16000.0 0.0 0.0
x_14 z_1 30000.0 17000.0 0.0 0.0
x_14 z_2 30000.0 18000.0 0.0 0.0
x_14 z_3 30000.0 19000.0 0.0 0.0
x_14 z_4 30000.0 20000.0 0.0 0.0
x_15 z_0 32000.0 16000.0 0.0 0.0
x_15 z_1 32000.0 17000.0 0.0 0.0
x_15 z_2 32000.0 18000.0 0.0 0.0
x_15 z_3 32000.0 19000.0 0.0 0.0
x_15 z_4 32000.0 20000.0 0.0 0.0
x_16 z_0 34000.0 16000.0 0.0 0.0
x_16 z_1 34000.0 17000.0 0.0 0.0
x_16 z_2 34000.0 18000.0 0.0 0.0
x_16 z_3 34000.0 19000.0 0.0 0.0
x_16 z_4 34000.0 20000.0 0.0 0.0
x_17 z_0 36000.0 16000.0 0.0 0.0
x_17 z_1 36000.0 17000.0 0.0 0.0
x_17 z_2 36000.0 18000.0 0.0 0.0
x_17 z_3 36000.0 19000.0 0.0 0.0
x_17 z_4 36000.0 20000.0 0.0 0.0
x_18 z_0 38000.0 16000.0 0.0 0.0
x_18 z_1 38000.0 17000.0 0.0 0.0
x_18 z_2 38000.0 18000.0 0.0 0.0
x_18 z_3 38000.0 19000.0 0.0 0.0
x_18 z_4 38000.0 20000.0 0.0 0.0
x_19 z_0 40000.0 16000.0 0.0 0.0
x_19 z_1 40000.0 17000.0 0.0 0.0
x_19 z_2 40000.0 18000.0 0.0 0.0
x_19 z_3 40000.0 19000.0 0.0 0.0
x_19 z_4 40000.0 20000.0 0.0 0.0
x_20 z_0 42000.0 16000.0 0.0 0.0
x_20 z_1 42000.0 17000.0 0.0 0.0
x_20 z_2 42000.0 18000.0 0.0 0.0
x_20 z_3 42000.0 19000.0 0.0 0.0
x_20 z_4 42000.0 20000.0 0.0 0.0
x_21 z_0 44000.0 16000.0 0.0 0.0
x_21 z_1 44000.0 17000.0 0.0 0.0
x_21 z_2 44000.0 18000.0 0.0 0.0
x_21 z_3 44000.0 19000.0 0.0 0.0
x_21 z_4 44000.0 20000.0 0.0 0.0
x_22 z_0 46000.0 16000.0 0.0 0.0
x_22 z_1 46000.0 17000.0 0.0 0.0
x_22 z_2 46000.0 18000.0 0.0 0.0
x_22 z_3 46000.0 19000.0 0.0 0.0
x_22 z_4 46000.0 20000.0 0.0 0.0
x_23 z_0 48000.0 16000.0 0.0 0.0
x_23 z_1 48000.0 17000.0 0.0 0.0
x_23 z_2 48000.0 18000.0 0.0 0.0
x_23 z_3 48000.0 19000.0 0.0 0.0
x_23 z_4 48000.0 20000.0 0.0 0.0
x_24 z_0 50000.0 16000.0 0.0 0.0
x_24 z_1 50000.0 17000.0 0.0 0.0
x_24 z_2 50000.0 18000.0 0.0 0.0
x_24 z_3 50000.0 19000.0 0.0 0.0
x_24 z_4 50000.0 20000.0 0.0 0.0
x_25 z_0 52000.0 16000.0 0.0 0.0
x_25 z_1 52000.0 17000.0 0.0 0.0
x_25 z_2 52000.0 18000.0 0.0 0.0
x_25 z_3 52000.0 19000.0 0.0 0.0
x_25 z_4 52000.0 20000.0 0.0 0.0
x_26 z_0 54000.0 16000.0 0.0 0.0
x_26 z_1 54000.0 17000.0 0.0 0.0
x_26 z_2 54000.0 18000.0 0.0 0.0
x_26 z_3 54000.0 19000.0 0.0 0.0
x_26 z_4 54000.0 20000.0 0.0 0.0
x_27 z_0 56000.0 16000.0 0.0 0.0
x_27 z_1 56000.0 17000.0 0.0 0.0
x_27 z_2 56000.0 18000.0 0.0 0.0
x_27 z_3 56000.0 19000.0 0.0 0.0
x_27 z_4 56000.0 20000.0 0.0 0.0
x_28 z_0 58000.0 16000.0 0.0 0.0
x_28 z_1 58000.0 17000.0 0.0 0.0
x_28 z_2 58000.0 18000.0 0.0 0.0
x_28 z_3 58000.0 19000.0 0.0 0.0
x_28 z_4 58000.0 20000.0 0.0 0.0
x_29 z_0 60000.0 16000.0 0.0 0.0
x_29 z_1 60000.0 17000.0 0.0 0.0
x_29 z_2 60000.0 18000.0 0.0 0.0
x_29 z_3 60000.0 19000.0 0.0 0.0
x_29 z_4 60000.0 20000.0 0.0 0.0
x_30 z_0 62000.0 16000.0 0.0 0.0
x_30 z_1 62000.0 17000.0 0.0 0.0
x_30 z_2 62000.0 18000.0 0.0 0.0
x_30 z_3 62000.0 19000.0 0.0 0.0
x_30 z_4 62000.0 20000.0 0.0 0.0
x_31 z_0 64000.0 16000.0 0.0 0.0
x_31 z_1 64000.0 17000.0 0.0 0.0
x_31 z_2 64000.0 18000.0 0.0 0.0
x_31 z_3 64000.0 19000.0 0.0 0.0
x_31 z_4 64000.0 20000.0 0.0 0.0
x_32 z_0 66000.0 16000.0 0.0 0.0
x_32 z_1 66000.0 17000.0 0.0 0.0
x_32 z_2 66000.0 18000.0 0.0 0.0
x_32 z_3 66000.0 19000.0 0.0 0.0
x_32 z_4 66000.0 20000.0 0.0 0.0
x_33 z_0 68000.0 16000.0 0.0 0.0
x_33 z_1 68000.0 17000.0 0.0 0.0
x_33 z_2 68000.0 18000.0 0.0 0.0
x_33 z_3 68000.0 19000.0 0.0 0.0
x_33 z_4 68000.0 20000.0 0.0 0.0
x_34 z_0 70000.0 16000.0 0.0 0.0
x_34 z_1 70000.0 17000.0 0.0 0.0
x_34 z_2 70000.0 18000.0 0.0 0.0
x_34 z_3 70000.0 19000.0 0.0 0.0
x_34 z_4 70000.0 20000.0 0.0 0.0
x_35 z_0 72000.0 16000.0 0.0 0.0
x_35 z_1 72000.0 17000.0 0.0 0.0
x_35 z_2 72000.0 18000.0 0.0 0.0
x_35 z_3 72000.0 19000.0 0.0 0.0
x_35 z_4 72000.0 20000.0 0.0 0.0
x_36 z_0 74000.0 16000.0 0.0 0.0
x_36 z_1 74000.0 17000.0 0.0 0.0
x_36 z_2 74000.0 18000.0 0.0 0.0
x_36 z_3 74000.0 19000.0 0.0 0.0
x_36 z_4 74000.0 20000.0 0.0 0.0
x_37 z_0 76000.0 16000.0 0.0 0.0
x_37 z_1 76000.0 17000.0 0.0 0.0
x_37 z_2 76000.0 18000.0 0.0 0.0
x_37 z_3 76000.0 19000.0 0.0 0.0
x_37 z_4 76000.0 20000.0 0.0 0.0
x_38 z_0 78000.0 16000.0 0.0 0.0
x_38 z_1 78000.0 17000.0 0.0 0.0
x_38 z_2 78000.0 18000.0 0.0 0.0
x_38 z_3 78000.0 19000.0 0.0 0.0
x_38 z_4 78000.0 20000.0 0.0 0.0
x_39 z_0 80000.0 16000.0 0.0 0.0
x_39 z_1 80000.0 17000.0 0.0 0.0
x_39 z_2 80000.0 18000.0 0.0 0.0
x_39 z_3 80000.0 19000.0 0.0 0.0
x_39 z_4 80000.0 20000.0 0.0 0.0
x_40 z_0 82000.0 16000.0 0.0 0.0
x_40 z_1 82000.0 17000.0 0.0 0.0
x_40 z_2 82000.0 18000.0 0.0 0.0
x_40 z_3 82000.0 19000.0 0.0 0.0
x_40 z_4 82000.0 20000.0 0.0 0.0
x_41 z_0 84000.0 16000.0 0.0 0.0
x_41 z_1 84000.0 17000.0 0.0 0.0
x_41 z_2 84000.0 18000.0 0.0 0.0
x_41 z_3 84000.0 19000.0 0.0 0.0
x_41 z_4 84000.0 20000.0 0.0 0.0
x_42 z_0 86000.0 16000.0 0.0 0.0
x_42 z_1 86000.0 17000.0 0.0 0.0
x_42 z_2 86000.0 18000.0 0.0 0.0
x_42 z_3 86000.0 19000.0 0.0 0.0
x_42 z_4 86000.0 20000.0 0.0 0.0
x_43 z_0 88000.0 16000.0 0.0 0.0
x_43 z_1 88000.0 17000.0 0.0 0.0
x_43 z_2 88000.0 18000.0 0.0 0.0
x_43 z_3 88000.0 19000.0 0.0 0.0
x_43 z_4 88000.0 20000.0 0.0 0.0
x_44 z_0 90000.0 16000.0 0.0 0.0
x_44 z_1 90000.0 17000.0 0.0 0.0
x_44 z_2 90000.0 18000.0 0.0 0.0
x_44 z_3 90000.0 19000.0 0.0 0.0
x_44 z_4 90000.0 20000.0 0.0 0.0
x_45 z_0 92000.0 16000.0 0.0 0.0
x_45 z_1 92000.0 17000.0 0.0 0.0
x_45 z_2 92000.0 18000.0 0.0 0.0
x_45 z_3 92000.0 19000.0 0.0 0.0
x_45 z_4 92000.0 20000.0 0.0 0.0
x_46 z_0 94000.0 16000.0 0.0 0.0
x_46 z_1 94000.0 17000.0 0.0 0.0
x_46 z_2 94000.0 18000.0 0.0 0.0
x_46 z_3 94000.0 19000.0 0.0 0.0
x_46 z_4 94000.0 20000.0 0.0 0.0
x_47 z_0 96000.0 16000.0 0.0 0.0
x_47 z_1 96000.0 17000.0 0.0 0.0
x_47 z_2 96000.0 18000.0 0.0 0.0
x_47 z_3 96000.0 19000.0 0.0 0.0
x_47 z_4 96000.0 20000.0 0.0 0.0
x_48 z_0 98000.0 16000.0 0.0 0.0
x_48 z_1 98000.0 17000.0 0.0 0.0
x_48 z_2 98000.0 18000.0 0.0 0.0
x_48 z_3 98000.0 19000.0 0.0 0.0
x_48 z_4 98000.0 20000.0 0.0 0.0
Loading

0 comments on commit b735263

Please sign in to comment.