[cig-commits] [commit] devel, master: added EXAMPLES/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom (5d5aa2d)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Wed Jun 18 15:20:40 PDT 2014
Repository : https://github.com/geodynamics/specfem2d
On branches: devel,master
Link : https://github.com/geodynamics/specfem2d/compare/fc67e6fd7ad890705b2b72b4b3c509accb22249e...e9ca46c40131588d89d7b0883250bc6584ce6b4c
>---------------------------------------------------------------
commit 5d5aa2d52151edb80ea21c7b70400e1089f04fae
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date: Mon Jan 9 02:01:59 2012 +0000
added EXAMPLES/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom
>---------------------------------------------------------------
5d5aa2d52151edb80ea21c7b70400e1089f04fae
M2_UPPA/topo_cours_M2_UPPA_curved.dat | 10 ----
fluid_solid/fluid_solid_external_mesh/README | 8 ++-
.../Par_file_fluid_solid | 61 +++++++++++-----------
.../README | 12 +++++
.../SOURCE_fluid_solid | 6 +--
.../interfaces_fluid_flat.dat | 20 +++----
.../process.sh | 6 ++-
7 files changed, 59 insertions(+), 64 deletions(-)
diff --git a/M2_UPPA/topo_cours_M2_UPPA_curved.dat b/M2_UPPA/topo_cours_M2_UPPA_curved.dat
deleted file mode 100644
index 96f3b94..0000000
--- a/M2_UPPA/topo_cours_M2_UPPA_curved.dat
+++ /dev/null
@@ -1,10 +0,0 @@
-9
-0 3000
-500 3000
-1000 3100
-1500 3400
-2000 3300
-2500 3200
-3000 3100
-3500 3000
-5000 3000
diff --git a/fluid_solid/fluid_solid_external_mesh/README b/fluid_solid/fluid_solid_external_mesh/README
index 86607f9..056f59a 100644
--- a/fluid_solid/fluid_solid_external_mesh/README
+++ b/fluid_solid/fluid_solid_external_mesh/README
@@ -2,12 +2,10 @@
README
----------------------------------------------------------------------
-This example creates and external mesh, with fluid and solid domains,
-and topography.
+This example creates an external mesh with fluid and solid domains, and topography.
-
-To run the examples, type:
+To run the example, type:
./process.sh
-in this current directory.
+in this directory.
diff --git a/M2_UPPA/Par_file_M2_UPPA b/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/Par_file_fluid_solid
similarity index 80%
copy from M2_UPPA/Par_file_M2_UPPA
copy to fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/Par_file_fluid_solid
index 2a4295e..a9e141b 100644
--- a/M2_UPPA/Par_file_M2_UPPA
+++ b/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/Par_file_fluid_solid
@@ -1,5 +1,5 @@
-# title of job r19201
-title = Test for M2 UPPA
+# title of job
+title = Flat fluid/solid interface
# forward or adjoint simulation
SIMULATION_TYPE = 1 # 1 = forward, 2 = adjoint + kernels
@@ -11,11 +11,11 @@ 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
-ngnod = 9 # number of control nodes per element (4 or 9)
+ngnod = 4 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
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
+READ_EXTERNAL_SEP_FILE = .false. # Read external model from DATA/model_velocity.dat_input, or use routine
TURN_ATTENUATION_ON = .false. # turn attenuation on or off for solid medium
TURN_VISCATTENUATION_ON = .false. # turn viscous attenuation on or off
Q0 = 1 # quality factor for viscous attenuation
@@ -23,8 +23,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 = 1600 # total number of time steps
-deltat = 1.d-3 # duration of a time step
+nt = 4500 # total number of time steps
+deltat = 0.70d-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
@@ -36,7 +36,7 @@ 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 set parameters for seismograms
-seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 2 # record 1=displ 2=veloc 3=accel 4=pressure 6=potential
generate_STATIONS = .true. # creates a STATION file in ./DATA
nreceiversets = 1 # number of receiver sets
anglerec = 0.d0 # angle to rotate components at receivers
@@ -45,24 +45,24 @@ SU_FORMAT = .false. # output seismograms in Seismic
# first receiver set (repeat these 6 lines and adjust nreceiversets accordingly)
nrec = 11 # number of receivers
-xdeb = 300. # first receiver x in meters
-zdeb = 2200. # 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 = .true. # receivers inside the medium or at the surface
+xdeb = 2500.d0 # first receiver x in meters
+zdeb = 2933.33333d0 # first receiver z in meters
+xfin = 6000.d0 # last receiver x in meters (ignored if onlyone receiver)
+zfin = 2933.33333d0 # 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 = 100 # 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 = .true. # output color image of the results
-imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
+output_color_image = .true. # output JPEG color image of the results
+imagetype = 2 # 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
modelvect = .false. # display velocity model on vector plots
boundvect = .true. # display boundary conditions on plots
interpol = .true. # interpolation of the display or not
pointsdisp = 6 # points for interpolation of display (set to 1 for lower-left corner only)
-subsamp_postscript = 1 # subsampling of color snapshots
+subsamp_postscript = 1 # subsampling of background velocity model in PostScript snapshots
factor_subsample_image = 1 # factor to subsample color images output by the code (useful for very large models)
POWER_DISPLAY_COLOR = 0.30d0 # non linear display to enhance small amplitudes in color images
DRAW_WATER_CONSTANT_BLUE_IN_JPG = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water
@@ -72,20 +72,19 @@ USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the fi
gnuplot = .false. # generate a GNUPLOT file for the grid
output_grid = .false. # save the grid in a text file or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
+output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
-nbmodels = 4 # nb of different models
+nbmodels = 2 # 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
# 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 2700.d0 3000.d0 1732.051d0 0 0 9999 9999 0 0 0 0 0 0
-2 1 2500.d0 2700.d0 0 0 0 9999 9999 0 0 0 0 0 0
-3 1 2200.d0 2500.d0 1443.375d0 0 0 9999 9999 0 0 0 0 0 0
-4 1 2200.d0 2200.d0 1343.375d0 0 0 9999 9999 0 0 0 0 0 0
+1 1 2500.d0 3400.d0 1963.d0 0 0 9999 9999 0 0 0 0 0 0
+2 1 1020.d0 1500.d0 0.d0 0 0 9999 9999 0 0 0 0 0 0
+
# external mesh or not
read_external_mesh = .false.
@@ -118,22 +117,22 @@ tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the
# PARAMETERS FOR INTERNAL MESHING
# file containing interfaces for internal mesh
-interfacesfile = ../interfaces_M2_UPPA_curved.dat
+interfacesfile = ../interfaces_fluid_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 = 4000.d0 # abscissa of right side of the model
-nx = 80 # number of elements along X
+xmax = 6400.d0 # abscissa of right side of the model
+## DK DK the original 2000 Geophysics paper used nx = 120 but NGLLX = 6
+## DK DK here I rescale it to nx = 144 and NGLLX = 5 because nowadays we almost always use NGLLX = 5
+nx = 144 # number of elements along X
# absorbing boundary parameters (see absorbing_conditions above)
absorbbottom = .true.
absorbright = .true.
-absorbtop = .false.
+absorbtop = .true.
absorbleft = .true.
# define the different regions of the model in the (nx,nz) spectral element mesh
-nbregions = 4 # nb of regions and model number for each
-1 80 1 20 1
-1 80 21 40 2
-1 80 41 60 3
-60 70 21 40 4
+nbregions = 2 # nb of regions and model number for each
+1 144 1 54 1
+1 144 55 108 2
diff --git a/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/README b/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/README
new file mode 100644
index 0000000..b8211cc
--- /dev/null
+++ b/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/README
@@ -0,0 +1,12 @@
+----------------------------------------------------------------------
+README
+----------------------------------------------------------------------
+
+This example creates the fluid-solid example with flat ocean bottom from the 2000 Geophysics paper of Komatitsch, Barnes and Tromp.
+
+To run the example, type:
+
+ ./process.sh
+
+in this directory.
+
diff --git a/M2_UPPA/SOURCE_M2_UPPA b/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/SOURCE_fluid_solid
similarity index 82%
copy from M2_UPPA/SOURCE_M2_UPPA
copy to fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/SOURCE_fluid_solid
index 73cd2c0..b1073c3 100644
--- a/M2_UPPA/SOURCE_M2_UPPA
+++ b/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/SOURCE_fluid_solid
@@ -1,12 +1,12 @@
#source 1. 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.
source_surf = .false. # source inside the medium or at the surface
-xs = 2500. # source location x in meters
-zs = 2500. # source location z in meters
+xs = 1575. # 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 = 10.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.0 # angle of the source (for a force only)
+angleforce = 0.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)
Mxz = 0. # Mxz component (for a moment tensor source only)
diff --git a/acoustic_poroelastic/interfaces_acoustic_poroelastic.dat b/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/interfaces_fluid_flat.dat
similarity index 61%
copy from acoustic_poroelastic/interfaces_acoustic_poroelastic.dat
copy to fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/interfaces_fluid_flat.dat
index c42ba85..8a87500 100644
--- a/acoustic_poroelastic/interfaces_acoustic_poroelastic.dat
+++ b/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/interfaces_fluid_flat.dat
@@ -1,34 +1,28 @@
#
# number of interfaces
#
- 3
+ 2
#
-# for each interface below, we give the number of points and then x,y for each point
+# for each interface below, we give the number of points and then x,z for each point
#
#
# interface number 1 (bottom of the mesh)
#
2
0 0
- 4800 0
+ 6400 0
#
# interface number 2 (topography, top of the mesh)
#
2
- 0 2400
- 4800 2400
-#
-# interface number 3 (topography, top of the mesh)
-#
- 2
0 4800
- 4800 4800
+ 6400 4800
#
# for each layer, we give the number of spectral elements in the vertical direction
#
#
# layer number 1 (bottom layer)
#
- 110
-#
- 110
+## DK DK the original 2000 Geophysics paper used nz = 90 but NGLLZ = 6
+## DK DK here I rescale it to nz = 108 and NGLLZ = 5 because nowadays we almost always use NGLLZ = 5
+ 108
diff --git a/fluid_solid/fluid_solid_external_mesh/process.sh b/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/process.sh
similarity index 90%
copy from fluid_solid/fluid_solid_external_mesh/process.sh
copy to fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/process.sh
index 004adb5..b5cb380 100755
--- a/fluid_solid/fluid_solid_external_mesh/process.sh
+++ b/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/process.sh
@@ -46,13 +46,15 @@ cp DATA/SOURCE OUTPUT_FILES/
echo
echo " running mesher..."
echo
-./xmeshfem2D > OUTPUT_FILES/output_mesher.txt
+#./xmeshfem2D > OUTPUT_FILES/output_mesher.txt
+./xmeshfem2D
# runs simulation
echo
echo " running solver..."
echo
-./xspecfem2D > OUTPUT_FILES/output_solver.txt
+#./xspecfem2D > OUTPUT_FILES/output_solver.txt
+./xspecfem2D
# stores output
cp DATA/SOURCE_xz.dat OUTPUT_FILES/
More information about the CIG-COMMITS
mailing list