[cig-commits] [commit] devel, master: added EXAMPLES/Rayleigh_wave and fixed a small problem in the horizontal offset of Rayleigh sources (344d449)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Wed Jun 18 15:20:28 PDT 2014
Repository : https://github.com/geodynamics/specfem2d
On branches: devel,master
Link : https://github.com/geodynamics/specfem2d/compare/fc67e6fd7ad890705b2b72b4b3c509accb22249e...e9ca46c40131588d89d7b0883250bc6584ce6b4c
>---------------------------------------------------------------
commit 344d4493ccfff78496bae05dd84ba35fa53c5f75
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date: Wed Dec 28 14:36:37 2011 +0000
added EXAMPLES/Rayleigh_wave and fixed a small problem in the horizontal offset of Rayleigh sources
>---------------------------------------------------------------
344d4493ccfff78496bae05dd84ba35fa53c5f75
.../Par_file_Rayleigh_2D | 18 +++++++-------
Rayleigh_wave/README | 29 ++++++++++++++++++++++
.../SOURCE_Rayleigh_2D | 10 ++++----
.../interfaces_Rayleigh_flat.dat | 2 +-
{attenuation => Rayleigh_wave}/process.sh | 6 ++---
attenuation/README | 4 +--
6 files changed, 49 insertions(+), 20 deletions(-)
diff --git a/canyon/Par_file_canyon b/Rayleigh_wave/Par_file_Rayleigh_2D
similarity index 94%
copy from canyon/Par_file_canyon
copy to Rayleigh_wave/Par_file_Rayleigh_2D
index 5817ae7..a01c0b0 100644
--- a/canyon/Par_file_canyon
+++ b/Rayleigh_wave/Par_file_Rayleigh_2D
@@ -1,5 +1,5 @@
-# title of job r19201
-title = Canyon Mexico Paco P-SV
+# title of job
+title = Rayleigh wave in a homogeneous 2D medium
# forward or adjoint simulation
SIMULATION_TYPE = 1 # 1 = forward, 2 = adjoint + kernels
@@ -24,7 +24,7 @@ p_sv = .true. # set the type of calculation (
# time step parameters
nt = 5000 # total number of time steps
-deltat = 2.e-3 # duration of a time step
+deltat = 3.90625e-3 # duration of a time step
USER_T0 = 0.0d0 # use this t0 as earliest starting time rather than the automatically calculated one
# source parameters
@@ -52,9 +52,9 @@ zfin = 9. # last receiver z in meters (ig
enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 250 # display frequency in time steps
+NTSTEP_BETWEEN_OUTPUT_INFO = 200 # display frequency in time steps
output_postscript_snapshot = .true. # output Postscript snapshot of the results
-output_color_image = .false. # output color image of the results
+output_color_image = .true. # output color image of the results
imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
cutsnaps = 1. # minimum amplitude in % for snapshots
meshvect = .true. # display mesh on vector plots or not
@@ -85,7 +85,7 @@ nbmodels = 1 # nb of different models
1 1 1.d0 2.d0 1.d0 0 0 9999.d0 9999.d0 0 0 0 0 0 0
# external mesh or not
-read_external_mesh = .true.
+read_external_mesh = .false.
# absorbing boundary active or not
absorbing_conditions = .true.
@@ -115,12 +115,12 @@ tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the
# PARAMETERS FOR INTERNAL MESHING
# file containing interfaces for internal mesh
-interfacesfile = ../interfaces_no_canyon.dat
+interfacesfile = interfaces_Rayleigh_flat.dat
# geometry of the model (origin lower-left corner = 0,0) and mesh description
xmin = 0.d0 # abscissa of left side of the model
xmax = 19.d0 # abscissa of right side of the model
-nx = 95 # number of elements along X
+nx = 60 # number of elements along X
# absorbing boundary parameters (see absorbing_conditions above)
absorbbottom = .true.
@@ -130,4 +130,4 @@ absorbleft = .true.
# define the different regions of the model in the (nx,nz) spectral element mes
nbregions = 1 # nb of regions and model number for each
-1 95 1 45 1
+1 60 1 28 1
diff --git a/Rayleigh_wave/README b/Rayleigh_wave/README
new file mode 100644
index 0000000..ed34818
--- /dev/null
+++ b/Rayleigh_wave/README
@@ -0,0 +1,29 @@
+----------------------------------------------------------------------
+README
+----------------------------------------------------------------------
+
+This example propagates an initial (incident) Rayleigh wave along the surface of a 2D homogeneous medium, similar to Figure 3 of Komatitsch and Tromp (1999) but without the canyon.
+
+TO RUN:
+
+0. Read the user manual in SPECFEM2D/doc/manual_SPECFEM2D.pdf
+
+1. in the SPECFEM2D root directory, configure, e.g.,
+ ./configure FC=gfortran
+
+2. compile:
+ make all
+
+3. cd EXAMPLES/Rayleigh_wave
+
+4. execute script to run mesher and solver for the PSV case:
+ ./process.sh
+
+5. check out the output files in the local directory OUTPUT_FILES
+
+References:
+-----------
+
+Dimitri Komatitsch, Jean-Pierre Vilotte, Rossana Vai, José M. Castillo-Covarrubias and Francisco J. Sánchez-Sesma, The Spectral Element method for elastic wave equations: application to 2D and 3D seismic problems, International Journal for Numerical Methods in Engineering, vol. 45, p. 1139-1164 (1999).
+
+----------------------------------------------------------------------
diff --git a/acoustic_poroelastic/SOURCE_acoustic_poroelastic b/Rayleigh_wave/SOURCE_Rayleigh_2D
similarity index 68%
copy from acoustic_poroelastic/SOURCE_acoustic_poroelastic
copy to Rayleigh_wave/SOURCE_Rayleigh_2D
index 4fbe718..9bc3aba 100644
--- a/acoustic_poroelastic/SOURCE_acoustic_poroelastic
+++ b/Rayleigh_wave/SOURCE_Rayleigh_2D
@@ -1,11 +1,11 @@
# source parameters
source_surf = .false. # source inside the medium or at the surface
-xs = 1600. # source location x in meters
-zs = 2900. # source location z in meters
-source_type = 1 # elastic force or acoustic pressure = 1 or moment tensor = 2
+xs = 3.5 # source location x in meters
+zs = 9. # source location z in meters
+source_type = 3 # 1 for plane P waves, 2 for plane SV waves, 3 for Rayleigh wave
time_function_type = 1 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
-f0 = 15.0 # dominant source frequency (Hz) if not Dirac or Heaviside
-t0 = 0.0 # time shift when multi sources (if one source, must be zero)
+f0 = 1.0 # dominant source frequency (Hz) if not Dirac or Heaviside
+tshift = 0.0 # time shift when multi sources (if one source, must be zero)
angleforce = 0. # angle of the source (for a force only)
Mxx = 1. # Mxx component (for a moment tensor source only)
Mzz = 1. # Mzz component (for a moment tensor source only)
diff --git a/unused_older_examples_DATA_to_sort/interfaces_no_canyon.dat b/Rayleigh_wave/interfaces_Rayleigh_flat.dat
similarity index 98%
copy from unused_older_examples_DATA_to_sort/interfaces_no_canyon.dat
copy to Rayleigh_wave/interfaces_Rayleigh_flat.dat
index 0951555..1a9c7e6 100644
--- a/unused_older_examples_DATA_to_sort/interfaces_no_canyon.dat
+++ b/Rayleigh_wave/interfaces_Rayleigh_flat.dat
@@ -23,4 +23,4 @@
#
# layer number 1 (bottom layer)
#
- 45
+ 28
diff --git a/attenuation/process.sh b/Rayleigh_wave/process.sh
similarity index 90%
copy from attenuation/process.sh
copy to Rayleigh_wave/process.sh
index b62a643..b18512f 100755
--- a/attenuation/process.sh
+++ b/Rayleigh_wave/process.sh
@@ -21,9 +21,9 @@ mkdir -p DATA
# sets up local DATA/ directory
cd DATA/
-cp ../Par_file_attenuation_2D Par_file
-cp ../interfaces_attenuation_analytic.dat .
-cp ../SOURCE_attenuation_2D SOURCE
+cp ../Par_file_Rayleigh_2D Par_file
+cp ../interfaces_Rayleigh_flat.dat .
+cp ../SOURCE_Rayleigh_2D SOURCE
cd ../
# cleans output files
diff --git a/attenuation/README b/attenuation/README
index df4a869..677d97e 100644
--- a/attenuation/README
+++ b/attenuation/README
@@ -2,7 +2,7 @@
README
----------------------------------------------------------------------
-This default example creates the 2D attenuation example of Komatitsch and Tromp (1999, Figure 16).
+This example creates the 2D attenuation benchmark of Komatitsch and Tromp (1999, Figure 16).
TO RUN:
@@ -14,7 +14,7 @@ TO RUN:
2. compile:
make all
-3. cd EXAMPLES/attenuation/
+3. cd EXAMPLES/attenuation
4. execute script to run mesher and solver for the PSV case:
./process.sh
More information about the CIG-COMMITS
mailing list