[cig-commits] [commit] devel,master: added an example with an incident Rayleigh wave propagating towards a crack; modified the source code in order to "manually" add a crack (i.e., some parameters are hardwired and the crack detection and handling code I added is not very flexible yet) (7deafb5)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Wed Jun 18 15:20:30 PDT 2014
Repository : https://github.com/geodynamics/specfem2d
On branches: devel,master
Link : https://github.com/geodynamics/specfem2d/compare/fc67e6fd7ad890705b2b72b4b3c509accb22249e...e9ca46c40131588d89d7b0883250bc6584ce6b4c
>---------------------------------------------------------------
commit 7deafb50efd0e2e90619c96c87aae429d625f5e8
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date: Thu Dec 29 00:40:50 2011 +0000
added an example with an incident Rayleigh wave propagating towards a crack; modified the source code in order to "manually" add a crack (i.e., some parameters are hardwired and the crack detection and handling code I added is not very flexible yet)
>---------------------------------------------------------------
7deafb50efd0e2e90619c96c87aae429d625f5e8
.../Par_file_Rayleigh_2D | 0
{Rayleigh_wave => Rayleigh_wave_no_crack}/README | 2 +-
.../SOURCE_Rayleigh_2D | 0
.../interfaces_Rayleigh_flat.dat | 0
{Rayleigh_wave => Rayleigh_wave_no_crack}/process.sh | 0
.../Par_file_Rayleigh_2D | 18 ++++++++++++------
{Rayleigh_wave => Rayleigh_wave_with_crack}/README | 4 ++--
.../SOURCE_Rayleigh_2D | 0
.../interfaces_Rayleigh_flat.dat | 0
{Rayleigh_wave => Rayleigh_wave_with_crack}/process.sh | 0
10 files changed, 15 insertions(+), 9 deletions(-)
diff --git a/Rayleigh_wave/Par_file_Rayleigh_2D b/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D
similarity index 100%
copy from Rayleigh_wave/Par_file_Rayleigh_2D
copy to Rayleigh_wave_no_crack/Par_file_Rayleigh_2D
diff --git a/Rayleigh_wave/README b/Rayleigh_wave_no_crack/README
similarity index 96%
copy from Rayleigh_wave/README
copy to Rayleigh_wave_no_crack/README
index ed34818..fa74ae8 100644
--- a/Rayleigh_wave/README
+++ b/Rayleigh_wave_no_crack/README
@@ -14,7 +14,7 @@ TO RUN:
2. compile:
make all
-3. cd EXAMPLES/Rayleigh_wave
+3. cd EXAMPLES/Rayleigh_wave_no_crack
4. execute script to run mesher and solver for the PSV case:
./process.sh
diff --git a/Rayleigh_wave/SOURCE_Rayleigh_2D b/Rayleigh_wave_no_crack/SOURCE_Rayleigh_2D
similarity index 100%
copy from Rayleigh_wave/SOURCE_Rayleigh_2D
copy to Rayleigh_wave_no_crack/SOURCE_Rayleigh_2D
diff --git a/Rayleigh_wave/interfaces_Rayleigh_flat.dat b/Rayleigh_wave_no_crack/interfaces_Rayleigh_flat.dat
similarity index 100%
copy from Rayleigh_wave/interfaces_Rayleigh_flat.dat
copy to Rayleigh_wave_no_crack/interfaces_Rayleigh_flat.dat
diff --git a/Rayleigh_wave/process.sh b/Rayleigh_wave_no_crack/process.sh
similarity index 100%
copy from Rayleigh_wave/process.sh
copy to Rayleigh_wave_no_crack/process.sh
diff --git a/Rayleigh_wave/Par_file_Rayleigh_2D b/Rayleigh_wave_with_crack/Par_file_Rayleigh_2D
similarity index 92%
rename from Rayleigh_wave/Par_file_Rayleigh_2D
rename to Rayleigh_wave_with_crack/Par_file_Rayleigh_2D
index a01c0b0..3d35385 100644
--- a/Rayleigh_wave/Par_file_Rayleigh_2D
+++ b/Rayleigh_wave_with_crack/Par_file_Rayleigh_2D
@@ -1,5 +1,5 @@
# title of job
-title = Rayleigh wave in a homogeneous 2D medium
+title = Rayleigh wave in a homogeneous medium with crack
# forward or adjoint simulation
SIMULATION_TYPE = 1 # 1 = forward, 2 = adjoint + kernels
@@ -9,11 +9,11 @@ SAVE_FORWARD = .false. # save the last frame, needed for adj
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
ngnod = 4 # number of control nodes per element (4 or 9)
initialfield = .true. # use a plane wave as source or not
-add_Bielak_conditions = .true. # add Bielak conditions or not if initial plane wave
+add_Bielak_conditions = .false. # add Bielak conditions or not if initial plane wave
assign_external_model = .false. # define external earth model or not
READ_EXTERNAL_SEP_FILE = .false. # Read external SEP file from DATA/model_velocity.dat_input, or use routine
TURN_ATTENUATION_ON = .false. # turn attenuation on or off for solid medium
@@ -75,7 +75,7 @@ output_energy = .false. # compute and output acoustic a
output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
-nbmodels = 1 # nb of different models
+nbmodels = 3 # nb of different models
# define models as
# I: (model_number 1 rho Vp Vs 0 0 QKappa Qmu 0 0 0 0 0 0) or
# II: (model_number 2 rho c11 c13 c15 c33 c35 c55 0 0 0 0 0 0) or
@@ -83,6 +83,10 @@ nbmodels = 1 # nb of different models
# For istropic elastic/acoustic material use I and set Vs to zero to make a given model acoustic, for anisotropic elastic use II,
# and for isotropic poroelastic material use III. The mesh can contain acoustic, elastic, and poroelastic models simultaneously.
1 1 1.d0 2.d0 1.d0 0 0 9999.d0 9999.d0 0 0 0 0 0 0
+## DK DK create three identical media (thus it is still a homogeneous medium)
+## DK DK just to have a way of detecting the crack in the mesh and to indicate with colors where the crack is located in the PostScript displays
+2 1 1.d0 2.d0 1.d0 0 0 9999.d0 9999.d0 0 0 0 0 0 0
+3 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 = .false.
@@ -115,7 +119,7 @@ tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the
# PARAMETERS FOR INTERNAL MESHING
# file containing interfaces for internal mesh
-interfacesfile = interfaces_Rayleigh_flat.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
@@ -129,5 +133,7 @@ absorbtop = .false.
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
+nbregions = 3 # nb of regions and model number for each
1 60 1 28 1
+30 30 26 28 2
+31 31 26 28 3
diff --git a/Rayleigh_wave/README b/Rayleigh_wave_with_crack/README
similarity index 83%
rename from Rayleigh_wave/README
rename to Rayleigh_wave_with_crack/README
index ed34818..245146a 100644
--- a/Rayleigh_wave/README
+++ b/Rayleigh_wave_with_crack/README
@@ -2,7 +2,7 @@
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.
+This example propagates an initial (incident) Rayleigh wave along the surface of a 2D homogeneous medium with a crack, similar to Figure 3 of Komatitsch and Tromp (1999) but with a crack instead of a canyon.
TO RUN:
@@ -14,7 +14,7 @@ TO RUN:
2. compile:
make all
-3. cd EXAMPLES/Rayleigh_wave
+3. cd EXAMPLES/Rayleigh_wave_with_crack
4. execute script to run mesher and solver for the PSV case:
./process.sh
diff --git a/Rayleigh_wave/SOURCE_Rayleigh_2D b/Rayleigh_wave_with_crack/SOURCE_Rayleigh_2D
similarity index 100%
rename from Rayleigh_wave/SOURCE_Rayleigh_2D
rename to Rayleigh_wave_with_crack/SOURCE_Rayleigh_2D
diff --git a/Rayleigh_wave/interfaces_Rayleigh_flat.dat b/Rayleigh_wave_with_crack/interfaces_Rayleigh_flat.dat
similarity index 100%
rename from Rayleigh_wave/interfaces_Rayleigh_flat.dat
rename to Rayleigh_wave_with_crack/interfaces_Rayleigh_flat.dat
diff --git a/Rayleigh_wave/process.sh b/Rayleigh_wave_with_crack/process.sh
similarity index 100%
rename from Rayleigh_wave/process.sh
rename to Rayleigh_wave_with_crack/process.sh
More information about the CIG-COMMITS
mailing list