[cig-commits] [commit] devel, master: Added an acoustic-poroelastic example in EXAMPLES/acoustic_poroelastic (302d05c)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Wed Jun 18 15:19:57 PDT 2014
Repository : https://github.com/geodynamics/specfem2d
On branches: devel,master
Link : https://github.com/geodynamics/specfem2d/compare/fc67e6fd7ad890705b2b72b4b3c509accb22249e...e9ca46c40131588d89d7b0883250bc6584ce6b4c
>---------------------------------------------------------------
commit 302d05c039159f558b84f60db024ed95509d08c2
Author: Christina Morency <morency1 at llnl.gov>
Date: Fri Oct 14 16:58:30 2011 +0000
Added an acoustic-poroelastic example in EXAMPLES/acoustic_poroelastic
>---------------------------------------------------------------
302d05c039159f558b84f60db024ed95509d08c2
.../Par_file_acoustic_poroelastic | 44 +++++++++++-----------
{M2_UPPA => acoustic_poroelastic}/README | 14 +++----
.../SOURCE_acoustic_poroelastic | 6 +--
.../interfaces_acoustic_poroelastic.dat | 18 ++++++---
{M2_UPPA => acoustic_poroelastic}/process.sh | 6 +--
5 files changed, 48 insertions(+), 40 deletions(-)
diff --git a/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid b/acoustic_poroelastic/Par_file_acoustic_poroelastic
similarity index 81%
copy from Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid
copy to acoustic_poroelastic/Par_file_acoustic_poroelastic
index fce5c60..ab38c74 100644
--- a/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid
+++ b/acoustic_poroelastic/Par_file_acoustic_poroelastic
@@ -1,5 +1,5 @@
# title of job
-title = Test for Abel Balanche UMR 6538 Brest
+title = Test for 2 layers: acoustic/poroelastic
# forward or adjoint simulation
SIMULATION_TYPE = 1 # 1 = forward, 2 = adjoint + kernels
@@ -22,8 +22,8 @@ freq0 = 10 # frequency for viscous attenua
p_sv = .true. # set the type of calculation (P-SV or SH/membrane waves)
# time step parameters
-nt = 3000 # total number of time steps
-deltat = 1.d-3 # duration of a time step
+nt = 5000 # total number of time steps
+deltat = 3.d-4 # duration of a time step
# source parameters
NSOURCES = 1 # number of sources [source info read in CMTSOLUTION file]
@@ -34,32 +34,32 @@ N_SLS = 2 # number of standard li
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# receiver line parameters for seismograms
-seismotype = 4 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 2 # record 1=displ 2=veloc 3=accel 4=pressure
generate_STATIONS = .true. # creates a STATION file in ./DATA
nreceiverlines = 2 # number of receiver lines
anglerec = 0.d0 # angle to rotate components at receivers
rec_normal_to_surface = .false. # base anglerec normal to surface (external mesh and curve file needed)
-# first receiver line (in the ocean)
-nrec = 11 # number of receivers
+# first receiver line (in the acoustic layer)
+nrec = 1 # number of receivers
xdeb = 2000. # first receiver x in meters
-zdeb = -2000. # first receiver z in meters
-xfin = 8000. # last receiver x in meters (ignored if onlyone receiver)
-zfin = -2000. # last receiver z in meters (ignored if onlyone receiver)
+zdeb = 2933.33 # first receiver z in meters
+xfin = 3700. # last receiver x in meters (ignored if onlyone receiver)
+zfin = 2200. # last receiver z in meters (ignored if onlyone receiver)
enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface
-# second receiver line (in the solid)
-nrec = 11 # number of receivers
+# second receiver line (in the poroelastic layer)
+nrec = 1 # number of receivers
xdeb = 2000. # first receiver x in meters
-zdeb = -4500. # first receiver z in meters
-xfin = 8000. # last receiver x in meters (ignored if onlyone receiver)
-zfin = -4500. # last receiver z in meters (ignored if onlyone receiver)
+zdeb = 1866.67 # first receiver z in meters
+xfin = 3777 # last receiver x in meters (ignored if onlyone receiver)
+zfin = 1866.67 # last receiver z in meters (ignored if onlyone receiver)
enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface
# display parameters
NTSTEP_BETWEEN_OUTPUT_INFO = 300 # display frequency in time steps
output_postscript_snapshot = .true. # output Postscript snapshot of the results
-output_color_image = .true. # output color image of the results
+output_color_image = .false. # 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
@@ -82,8 +82,8 @@ nbmodels = 2 # nb of different models
# III: (model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu).
# 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 1020 1500 0 0 0 10.d0 10.d0 0 0 0 0 0 0
-2 1 2500 5500 3175.5196 0 0 10.d0 10.d0 0 0 0 0 0 0
+1 3 2500.d0 1020.d0 0.4d0 2.0 1d-11 0.d0 1d-11 1.60554d10 2.295d9 1.0d10 0.0d-4 9.63342d9 10.d0
+2 1 1020.d0 1500.d0 0.d0 0 0 10.d0 10.d0 0 0 0 0 0 0
# external mesh or not
read_external_mesh = .false.
@@ -107,12 +107,12 @@ tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the
# PARAMETERS FOR INTERNAL MESHING
# file containing interfaces for internal mesh
-interfacesfile = ../interfaces_Abel_Brest_bathymetry.dat
+interfacesfile = ../interfaces_acoustic_poroelastic.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 = 10000.d0 # abscissa of right side of the model
-nx = 134 # number of elements along X
+xmax = 4800.d0 # abscissa of right side of the model
+nx = 260 # number of elements along X
# absorbing boundary parameters (see absorbing_conditions above)
absorbbottom = .true.
@@ -122,5 +122,5 @@ absorbleft = .true.
# define the different regions of the model in the (nx,nz) spectral element mesh
nbregions = 2 # nb of regions and model number for each
-1 134 1 44 2
-1 134 45 84 1
+1 260 1 110 1
+1 260 111 220 2
diff --git a/M2_UPPA/README b/acoustic_poroelastic/README
similarity index 51%
copy from M2_UPPA/README
copy to acoustic_poroelastic/README
index f6c5383..e99df2d 100644
--- a/M2_UPPA/README
+++ b/acoustic_poroelastic/README
@@ -2,27 +2,27 @@
README
----------------------------------------------------------------------
-This default example creates an internal mesh with smooth topographic interfaces.
+This default example creates an internal mesh with a flat interface between
+an acoustic layer on top of a poroelastic layer.
TO RUN:
0. Read the user manual in SPECFEM2D/doc/manual_SPECFEM2D.pdf
-1. in SPECFEM2D root directory, configure, e.g.,
+1. in SPECFEM2D root directory, configure, e.g. for a serial calculation,
./configure FC=gfortran
2. compile:
make all
-3. cd EXAMPLES/M2_UPPA/
+3. cd EXAMPLES/acoustic_poroelastic
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
+5. check out the output files in the local directory OUTPUT_FILES,
+the seismograms should compare to Fig.8 in Morency and Tromp, Geophys. J. Int. (2008) 175, 301-345
+with a reversed polarity due to a reversed expression of the Ricker source time function in this paper.
-6. optional: set output_wavefield_snapshot = .true. in Par_file, re-run,
- then try plotting the wavefield using the script
- SPECFEM2D/UTILS/visualization/plot_wavefield.pl
----------------------------------------------------------------------
diff --git a/DATA_to_sort_older_examples/SOURCE_Abel_Balanche_bathy_source_solid b/acoustic_poroelastic/SOURCE_acoustic_poroelastic
similarity index 78%
copy from DATA_to_sort_older_examples/SOURCE_Abel_Balanche_bathy_source_solid
copy to acoustic_poroelastic/SOURCE_acoustic_poroelastic
index 9b93658..4fbe718 100644
--- a/DATA_to_sort_older_examples/SOURCE_Abel_Balanche_bathy_source_solid
+++ b/acoustic_poroelastic/SOURCE_acoustic_poroelastic
@@ -1,10 +1,10 @@
# source parameters
source_surf = .false. # source inside the medium or at the surface
-xs = 4000. # source location x in meters
-zs = -7000. # source location z in meters
+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
time_function_type = 1 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
-f0 = 8.0 # dominant source frequency (Hz) if not Dirac or Heaviside
+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)
angleforce = 0. # angle of the source (for a force only)
Mxx = 1. # Mxx component (for a moment tensor source only)
diff --git a/DATA_to_sort_older_examples/interfaces_no_canyon.dat b/acoustic_poroelastic/interfaces_acoustic_poroelastic.dat
similarity index 74%
copy from DATA_to_sort_older_examples/interfaces_no_canyon.dat
copy to acoustic_poroelastic/interfaces_acoustic_poroelastic.dat
index 0951555..c42ba85 100644
--- a/DATA_to_sort_older_examples/interfaces_no_canyon.dat
+++ b/acoustic_poroelastic/interfaces_acoustic_poroelastic.dat
@@ -1,7 +1,7 @@
#
# number of interfaces
#
- 2
+ 3
#
# for each interface below, we give the number of points and then x,y for each point
#
@@ -10,17 +10,25 @@
#
2
0 0
- 19 0
+ 4800 0
#
# interface number 2 (topography, top of the mesh)
#
2
- 0 9
- 19 9
+ 0 2400
+ 4800 2400
+#
+# interface number 3 (topography, top of the mesh)
+#
+ 2
+ 0 4800
+ 4800 4800
#
# for each layer, we give the number of spectral elements in the vertical direction
#
#
# layer number 1 (bottom layer)
#
- 45
+ 110
+#
+ 110
diff --git a/M2_UPPA/process.sh b/acoustic_poroelastic/process.sh
similarity index 89%
copy from M2_UPPA/process.sh
copy to acoustic_poroelastic/process.sh
index 4f94f18..5437bf4 100755
--- a/M2_UPPA/process.sh
+++ b/acoustic_poroelastic/process.sh
@@ -8,7 +8,7 @@ echo "running example: `date`"
currentdir=`pwd`
echo
-echo "(will take about 2 minutes)"
+echo "(will take about 30 minutes)"
echo
# sets up directory structure in current example directoy
@@ -21,8 +21,8 @@ mkdir -p DATA
# sets up local DATA/ directory
cd DATA/
-cp ../Par_file_M2_UPPA Par_file
-cp ../SOURCE_M2_UPPA SOURCE
+ln -s ../Par_file_acoustic_poroelastic Par_file
+ln -s ../SOURCE_acoustic_poroelastic SOURCE
cd ../
# cleans output files
More information about the CIG-COMMITS
mailing list