[cig-commits] [commit] devel, master: added a cleaned version of the 2D attenuation example from our 1999 GJI paper (Figure 16) (fd5fc2e)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Wed Jun 18 15:20:16 PDT 2014
Repository : https://github.com/geodynamics/specfem2d
On branches: devel,master
Link : https://github.com/geodynamics/specfem2d/compare/fc67e6fd7ad890705b2b72b4b3c509accb22249e...e9ca46c40131588d89d7b0883250bc6584ce6b4c
>---------------------------------------------------------------
commit fd5fc2e7cdcf85b125b5519b41a5a50009b269ef
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date: Tue Dec 27 21:24:26 2011 +0000
added a cleaned version of the 2D attenuation example from our 1999 GJI paper (Figure 16)
>---------------------------------------------------------------
fd5fc2e7cdcf85b125b5519b41a5a50009b269ef
.../Par_file_attenuation_2D | 58 +++--
{M2_UPPA => attenuation}/README | 13 +-
.../SOURCE_attenuation_2D | 0
.../interfaces_attenuation_analytic.dat | 0
{Tape2007_kernel => attenuation}/process.sh | 12 +-
.../Par_file_Abel_Balanche_bathy_source_solid | 107 ----------
.../Par_file_attenuation_2D | 100 ---------
.../SOURCE_Abel_Balanche_bathy_source_solid | 13 --
.../interfaces_Abel_Brest_bathymetry.dat | 234 ---------------------
.../interfaces_Abel_Brest_flat.dat | 36 ----
.../topo_rouen_aniso.dat | 3 -
11 files changed, 40 insertions(+), 536 deletions(-)
diff --git a/M2_UPPA/Par_file_M2_UPPA b/attenuation/Par_file_attenuation_2D
similarity index 78%
copy from M2_UPPA/Par_file_M2_UPPA
copy to attenuation/Par_file_attenuation_2D
index 2801c95..f8531b1 100644
--- a/M2_UPPA/Par_file_M2_UPPA
+++ b/attenuation/Par_file_attenuation_2D
@@ -1,5 +1,5 @@
-# title of job r19201
-title = Test for M2 UPPA
+# title of job
+title = Test for 2D attenuation 1999 GJI paper
# forward or adjoint simulation
SIMULATION_TYPE = 1 # 1 = forward, 2 = adjoint + kernels
@@ -11,20 +11,20 @@ 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
-TURN_ATTENUATION_ON = .false. # turn attenuation on or off for solid medium
+TURN_ATTENUATION_ON = .true. # 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
freq0 = 10 # frequency for viscous attenuation
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 = 1500 # total number of time steps
+deltat = 7.5e-4 # duration of a time step
USER_T0 = 0.0d0 # use this t0 as earliest starting time rather than the automatically calculated one
# source parameters
@@ -43,13 +43,13 @@ anglerec = 0.d0 # angle to rotate components at
rec_normal_to_surface = .false. # base anglerec normal to surface (external mesh and curve file needed)
SU_FORMAT = .false. # output seismograms in Seismic Unix format (adjoint traces will be read in the same format)
-# first receiver line (repeat these 6 lines and adjust nreceiverlines 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
+# first receiver line
+nrec = 2 # number of receivers
+xdeb = 1400.d0
+zdeb = 1500.d0
+xfin = 1500.d0
+zfin = 1500.d0
+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
@@ -70,28 +70,25 @@ sizemax_arrows = 1.d0 # maximum size of arrows on vec
US_LETTER = .false. # US letter paper or European A4
USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
gnuplot = .false. # generate a GNUPLOT file for the grid
-output_grid = .false. # save the grid in a text file or not
+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 = 1 # 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 2000.d0 3000.d0 2000.d0 0 0 27. 20. 0 0 0 0 0 0
# external mesh or not
read_external_mesh = .false.
# absorbing boundary active or not
-absorbing_conditions = .true.
+absorbing_conditions = .false.
# for horizontal periodic conditions: detect common points between left and right edges
ADD_PERIODIC_CONDITIONS = .false.
@@ -118,22 +115,19 @@ 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_attenuation_analytic.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 = 2000.d0 # abscissa of right side of the model
+nx = 44 # number of elements along X
# absorbing boundary parameters (see absorbing_conditions above)
-absorbbottom = .true.
-absorbright = .true.
+absorbbottom = .false.
+absorbright = .false.
absorbtop = .false.
-absorbleft = .true.
+absorbleft = .false.
# 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 = 1 # nb of regions and model number for each
+1 44 1 44 1
diff --git a/M2_UPPA/README b/attenuation/README
similarity index 55%
copy from M2_UPPA/README
copy to attenuation/README
index f6c5383..cb4da90 100644
--- a/M2_UPPA/README
+++ b/attenuation/README
@@ -2,27 +2,28 @@
README
----------------------------------------------------------------------
-This default example creates an internal mesh with smooth topographic interfaces.
+This default example creates the 2D attenuation example of Komatitsch and Tromp (1999, Figure 16).
TO RUN:
0. Read the user manual in SPECFEM2D/doc/manual_SPECFEM2D.pdf
-1. in SPECFEM2D root directory, configure, e.g.,
+1. in the SPECFEM2D root directory, configure, e.g.,
./configure FC=gfortran
2. compile:
make all
-3. cd EXAMPLES/M2_UPPA/
+3. cd EXAMPLES/attenuation/
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
-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
+Reference:
+----------
+
+Dimitri Komatitsch and Jeroen Tromp, Introduction to the spectral-element method for 3-D seismic wave propagation, Geophysical Journal International, vol. 139, p. 806-822 (1999)
----------------------------------------------------------------------
diff --git a/unused_older_examples_DATA_to_sort/SOURCE_attenuation_2D b/attenuation/SOURCE_attenuation_2D
similarity index 100%
rename from unused_older_examples_DATA_to_sort/SOURCE_attenuation_2D
rename to attenuation/SOURCE_attenuation_2D
diff --git a/unused_older_examples_DATA_to_sort/interfaces_attenuation_analytic.dat b/attenuation/interfaces_attenuation_analytic.dat
similarity index 100%
rename from unused_older_examples_DATA_to_sort/interfaces_attenuation_analytic.dat
rename to attenuation/interfaces_attenuation_analytic.dat
diff --git a/Tape2007_kernel/process.sh b/attenuation/process.sh
similarity index 81%
copy from Tape2007_kernel/process.sh
copy to attenuation/process.sh
index 852ed3a..b62a643 100755
--- a/Tape2007_kernel/process.sh
+++ b/attenuation/process.sh
@@ -8,7 +8,7 @@ echo "running example: `date`"
currentdir=`pwd`
echo
-echo "(will take about 1 minute)"
+echo "(will take a few minutes)"
echo
# sets up directory structure in current example directoy
@@ -18,12 +18,12 @@ echo
mkdir -p OUTPUT_FILES
mkdir -p DATA
-mkdir -p SEM
# sets up local DATA/ directory
cd DATA/
-ln -s ../Par_file_Tape2007_onerec Par_file
-ln -s ../SOURCE_001 SOURCE
+cp ../Par_file_attenuation_2D Par_file
+cp ../interfaces_attenuation_analytic.dat .
+cp ../SOURCE_attenuation_2D SOURCE
cd ../
# cleans output files
@@ -48,12 +48,14 @@ echo
echo " running mesher..."
echo
./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/
diff --git a/unused_older_examples_DATA_to_sort/Par_file_Abel_Balanche_bathy_source_solid b/unused_older_examples_DATA_to_sort/Par_file_Abel_Balanche_bathy_source_solid
deleted file mode 100644
index 1d43bb7..0000000
--- a/unused_older_examples_DATA_to_sort/Par_file_Abel_Balanche_bathy_source_solid
+++ /dev/null
@@ -1,107 +0,0 @@
-
-# title of job, and file that contains interface data
-title = Test for Abel Balanche UMR 6538 Brest
-interfacesfile = interfaces_Abel_Brest_bathymetry.dat
-
-# data concerning mesh, when generated using third-party app (more info in README)
-read_external_mesh = .false.
-mesh_file = ./DATA/Mesh_canyon/canyon_mesh_file # file containing the mesh
-nodes_coords_file = ./DATA/Mesh_canyon/canyon_nodes_coords_file # file containing the nodes coordinates
-materials_file = ./DATA/Mesh_canyon/canyon_materials_file # file containing the material number for each element
-free_surface_file = ./DATA/Mesh_canyon/canyon_free_surface_file # file containing the free surface
-absorbing_surface_file = ./DATA/Mesh_canyon/canyon_absorbing_surface_file # file containing the absorbing surface
-
-# parameters concerning partitioning
-nproc = 1 # number of processes
-partitioning_method = 1 # ascending order = 1, Metis = 2, Scotch = 3
-
-# 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
-ngnod = 9 # 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
-TURN_ANISOTROPY_ON = .false. # turn anisotropy on or off for solid medium
-TURN_ATTENUATION_ON = .false. # turn attenuation on or off for solid medium
-
-# absorbing boundaries parameters
-absorbing_conditions = .true. # absorbing boundary active or not
-absorbbottom = .true.
-absorbright = .true.
-absorbtop = .false.
-absorbleft = .true.
-
-# time step parameters
-nt = 3000 # total number of time steps
-deltat = 1.d-3 # duration of a time step
-
-# 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
-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
-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)
-Mxz = 0. # Mxz component (for a moment tensor source only)
-factor = 1.d10 # amplification factor
-
-# constants for attenuation
-N_SLS = 2 # number of standard linear solids for attenuation
-QKappa_attenuation = 136.4376068115 # quality factor P for attenuation
-Qmu_attenuation = 136.4376068115 # quality factor S for attenuation
-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
-generate_STATIONS = .true. # creates a STATION file in ./DATA
-nreceiverlines = 2 # number of receiver lines
-anglerec = 0.d0 # angle to rotate components at receivers
-
-# first receiver line (in the ocean)
-nrec = 11 # 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)
-enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface
-
-# second receiver line (in the solid)
-nrec = 11 # 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)
-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
-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
-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 = 1 # subsampling of color snapshots
-sizemax_arrows = 1.d0 # maximum size of arrows on vector plots in cm
-gnuplot = .false. # generate a GNUPLOT file for the grid
-outputgrid = .false. # save the grid in a text file or not
-
-# velocity and density models
-nbmodels = 2 # nb of different models
-# define models as (model_number,1,rho,vp,vs,0,0) or (model_number,2,rho,c11,c13,c33,c44)
-# set vs to zero to make a given model acoustic
-# the mesh can contain both acoustic and elastic models simultaneously
-1 1 1020 1500 0 0 0
-2 1 2500 5500 3175.5196 0 0
-# 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
diff --git a/unused_older_examples_DATA_to_sort/Par_file_attenuation_2D b/unused_older_examples_DATA_to_sort/Par_file_attenuation_2D
deleted file mode 100644
index 7add240..0000000
--- a/unused_older_examples_DATA_to_sort/Par_file_attenuation_2D
+++ /dev/null
@@ -1,100 +0,0 @@
-
-# title of job, and file that contains interface data
-title = Test for M2 UPPA
-interfacesfile = interfaces_attenuation_analytic.dat
-
-# forward or adjoint simulation
-SIMULATION_TYPE = 1 # 1 = forward, 2 = adjoint + kernels
-SAVE_FORWARD = .false. # save the last frame, needed for adjoint simulation
-
-# data concerning mesh, when generated using third-party app (more info in README)
-read_external_mesh = .false.
-mesh_file = ./DATA/Mesh_canyon/canyon_mesh_file # file containing the mesh
-nodes_coords_file = ./DATA/Mesh_canyon/canyon_nodes_coords_file # file containing the nodes coordinates
-materials_file = ./DATA/Mesh_canyon/canyon_materials_file # file containing the material number for each element
-free_surface_file = ./DATA/Mesh_canyon/canyon_free_surface_file # file containing the free surface
-absorbing_surface_file = ./DATA/Mesh_canyon/canyon_absorbing_surface_file # file containing the absorbing surface
-tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
-
-# parameters concerning partitionning
-nproc = 1 # number of processes
-partionning_method = 3 # ascending order = 1, Metis = 2, Scotch = 3
-
-# geometry of the model (origin lower-left corner = 0,0) and mesh description
-xmin = 0.d0 # abscissa of left side of the model
-xmax = 2000.d0 # abscissa of right side of the model
-nx = 44 # number of elements along X
-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
-TURN_ATTENUATION_ON = .true. # 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
-freq0 = 10 # frequency for viscous attenuation
-p_sv = .true. # set the type of calculation (P-SV or SH/membrane waves)
-
-# absorbing boundaries parameters
-absorbing_conditions = .false. # absorbing boundary active or not
-absorbbottom = .false.
-absorbright = .false.
-absorbtop = .false.
-absorbleft = .false.
-
-# time step parameters
-nt = 1500 # total number of time steps
-deltat = 7.5e-4 # duration of a time step
-
-# source parameters
-NSOURCES = 1 # number of sources [source info read in CMTSOLUTION file]
-force_normal_to_surface = .false. # angleforce normal to surface (external mesh and curve file needed)
-
-# constants for attenuation
-N_SLS = 2 # number of standard linear solids for attenuation
-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 = 1 # record 1=displ 2=veloc 3=accel 4=pressure
-generate_STATIONS = .true. # creates a STATION file in ./DATA
-nreceiverlines = 1 # 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
-nrec = 2 # number of receivers
-xdeb =1400.d0
-zdeb =1500.d0
-xfin =1500.d0
-zfin =1500.d0
-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
-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
-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 = 1 # subsampling of color snapshots
-sizemax_arrows = 1.d0 # maximum size of arrows on vector plots in cm
-gnuplot = .false. # generate a GNUPLOT file for the grid
-outputgrid = .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)
-
-# velocity and density models
-nbmodels = 1 # 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 2000.d0 3000.d0 2000.d0 0 0 27. 20. 0 0 0 0 0 0
-# define the different regions of the model in the (nx,nz) spectral element mesh
-nbregions = 1 # nb of regions and model number for each
-1 44 1 44 1
diff --git a/unused_older_examples_DATA_to_sort/SOURCE_Abel_Balanche_bathy_source_solid b/unused_older_examples_DATA_to_sort/SOURCE_Abel_Balanche_bathy_source_solid
deleted file mode 100644
index 9b93658..0000000
--- a/unused_older_examples_DATA_to_sort/SOURCE_Abel_Balanche_bathy_source_solid
+++ /dev/null
@@ -1,13 +0,0 @@
-# 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
-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
-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)
-Mzz = 1. # Mzz component (for a moment tensor source only)
-Mxz = 0. # Mxz component (for a moment tensor source only)
-factor = 1.d10 # amplification factor
diff --git a/unused_older_examples_DATA_to_sort/interfaces_Abel_Brest_bathymetry.dat b/unused_older_examples_DATA_to_sort/interfaces_Abel_Brest_bathymetry.dat
deleted file mode 100644
index 9b54369..0000000
--- a/unused_older_examples_DATA_to_sort/interfaces_Abel_Brest_bathymetry.dat
+++ /dev/null
@@ -1,234 +0,0 @@
-#
-# number of interfaces
-#
- 3
-#
-# for each interface below, we give the number of points and then x,y for each point
-#
-#
-# interface number 1 (bottom of the mesh)
-#
- 2
- 0 -10000
- 10000 -10000
-#
-# interface number 2 (shape of the ocean bottom, "middle" of the mesh)
-#
- 200
- 0.000000000000000E+000 -3000.00000000000
- 50.2512562814070 -3000.00000000000
- 100.502512562814 -3000.00000000000
- 150.753768844221 -3000.00000000000
- 201.005025125628 -3000.00000000000
- 251.256281407035 -3000.00000000000
- 301.507537688442 -3000.00000000000
- 351.758793969849 -3000.00000000000
- 402.010050251256 -3000.00000000000
- 452.261306532663 -3000.00000000000
- 502.512562814070 -3000.00000000000
- 552.763819095477 -3000.00000000000
- 603.015075376884 -3000.00000000000
- 653.266331658292 -3000.00000000000
- 703.517587939699 -3000.00000000000
- 753.768844221106 -3000.00000000000
- 804.020100502513 -3000.00000000000
- 854.271356783920 -3000.00000000000
- 904.522613065327 -3000.00000000000
- 954.773869346734 -3000.00000000000
- 1005.02512562814 -2999.99999999999
- 1055.27638190955 -2999.99999999999
- 1105.52763819095 -2999.99999999997
- 1155.77889447236 -2999.99999999994
- 1206.03015075377 -2999.99999999987
- 1256.28140703518 -2999.99999999973
- 1306.53266331658 -2999.99999999943
- 1356.78391959799 -2999.99999999882
- 1407.03517587940 -2999.99999999755
- 1457.28643216080 -2999.99999999498
- 1507.53768844221 -2999.99999998982
- 1557.78894472362 -2999.99999997957
- 1608.04020100503 -2999.99999995939
- 1658.29145728643 -2999.99999992011
- 1708.54271356784 -2999.99999984439
- 1758.79396984925 -2999.99999969997
- 1809.04522613065 -2999.99999942731
- 1859.29648241206 -2999.99999891786
- 1909.54773869347 -2999.99999797578
- 1959.79899497487 -2999.99999625159
- 2010.05025125628 -2999.99999312854
- 2060.30150753769 -2999.99998753006
- 2110.55276381910 -2999.99997759768
- 2160.80402010050 -2999.99996015856
- 2211.05527638191 -2999.99992985604
- 2261.30653266332 -2999.99987774717
- 2311.55778894472 -2999.99978906879
- 2361.80904522613 -2999.99963972340
- 2412.06030150754 -2999.99939082149
- 2462.31155778894 -2999.99898031417
- 2512.56281407035 -2999.99831033156
- 2562.81407035176 -2999.99722827654
- 2613.06532663317 -2999.99549897317
- 2663.31658291457 -2999.99276419917
- 2713.56783919598 -2999.98848470935
- 2763.81909547739 -2999.98185836178
- 2814.07035175879 -2999.97170619229
- 2864.32160804020 -2999.95631628237
- 2914.57286432161 -2999.93323311997
- 2964.82412060302 -2999.89897802111
- 3015.07537688442 -2999.84868431144
- 3065.32663316583 -2999.77562971313
- 3115.57788944724 -2999.67064821178
- 3165.82914572864 -2999.52140516632
- 3216.08040201005 -2999.31152324080
- 3266.33165829146 -2999.01955359034
- 3316.58291457286 -2998.61779730828
- 3366.83417085427 -2998.07099698650
- 3417.08542713568 -2997.33493763956
- 3467.33668341709 -2996.35502005601
- 3517.58793969849 -2995.06489714776
- 3567.83919597990 -2993.38529361747
- 3618.09045226131 -2991.22315897430
- 3668.34170854271 -2988.47133045741
- 3718.59296482412 -2985.00890186138
- 3768.84422110553 -2980.70250213323
- 3819.09547738693 -2975.40867928914
- 3869.34673366834 -2968.97755639549
- 3919.59798994975 -2961.25787379407
- 3969.84924623116 -2952.10345386256
- 4020.10050251256 -2941.38102223754
- 4070.35175879397 -2928.97919641565
- 4120.60301507538 -2914.81831610901
- 4170.85427135678 -2898.86065004476
- 4221.10552763819 -2881.12038422109
- 4271.35678391960 -2861.67269193906
- 4321.60804020101 -2840.66112160371
- 4371.85929648241 -2818.30252836067
- 4422.11055276382 -2794.88883085730
- 4472.36180904523 -2770.78500039549
- 4522.61306532663 -2746.42288541601
- 4572.86432160804 -2722.29073097341
- 4623.11557788945 -2698.91855430598
- 4673.36683417085 -2676.85986065549
- 4723.61809045226 -2656.67050006428
- 4773.86934673367 -2638.88574560987
- 4824.12060301508 -2623.99688702652
- 4874.37185929648 -2612.42875590047
- 4924.62311557789 -2604.51961216496
- 4974.87437185930 -2600.50471905680
- 5025.12562814070 -2600.50471905680
- 5075.37688442211 -2604.51961216496
- 5125.62814070352 -2612.42875590047
- 5175.87939698492 -2623.99688702652
- 5226.13065326633 -2638.88574560987
- 5276.38190954774 -2656.67050006428
- 5326.63316582915 -2676.85986065549
- 5376.88442211055 -2698.91855430598
- 5427.13567839196 -2722.29073097341
- 5477.38693467337 -2746.42288541601
- 5527.63819095477 -2770.78500039549
- 5577.88944723618 -2794.88883085730
- 5628.14070351759 -2818.30252836067
- 5678.39195979900 -2840.66112160371
- 5728.64321608040 -2861.67269193906
- 5778.89447236181 -2881.12038422109
- 5829.14572864322 -2898.86065004476
- 5879.39698492462 -2914.81831610901
- 5929.64824120603 -2928.97919641565
- 5979.89949748744 -2941.38102223754
- 6030.15075376884 -2952.10345386256
- 6080.40201005025 -2961.25787379407
- 6130.65326633166 -2968.97755639549
- 6180.90452261307 -2975.40867928914
- 6231.15577889447 -2980.70250213323
- 6281.40703517588 -2985.00890186138
- 6331.65829145729 -2988.47133045741
- 6381.90954773869 -2991.22315897430
- 6432.16080402010 -2993.38529361747
- 6482.41206030151 -2995.06489714776
- 6532.66331658292 -2996.35502005601
- 6582.91457286432 -2997.33493763956
- 6633.16582914573 -2998.07099698650
- 6683.41708542714 -2998.61779730828
- 6733.66834170854 -2999.01955359034
- 6783.91959798995 -2999.31152324080
- 6834.17085427136 -2999.52140516632
- 6884.42211055276 -2999.67064821178
- 6934.67336683417 -2999.77562971313
- 6984.92462311558 -2999.84868431144
- 7035.17587939699 -2999.89897802111
- 7085.42713567839 -2999.93323311997
- 7135.67839195980 -2999.95631628237
- 7185.92964824121 -2999.97170619229
- 7236.18090452261 -2999.98185836178
- 7286.43216080402 -2999.98848470935
- 7336.68341708543 -2999.99276419917
- 7386.93467336683 -2999.99549897317
- 7437.18592964824 -2999.99722827654
- 7487.43718592965 -2999.99831033156
- 7537.68844221106 -2999.99898031417
- 7587.93969849246 -2999.99939082149
- 7638.19095477387 -2999.99963972340
- 7688.44221105528 -2999.99978906879
- 7738.69346733668 -2999.99987774717
- 7788.94472361809 -2999.99992985604
- 7839.19597989950 -2999.99996015856
- 7889.44723618090 -2999.99997759768
- 7939.69849246231 -2999.99998753006
- 7989.94974874372 -2999.99999312854
- 8040.20100502513 -2999.99999625159
- 8090.45226130653 -2999.99999797578
- 8140.70351758794 -2999.99999891786
- 8190.95477386935 -2999.99999942731
- 8241.20603015075 -2999.99999969997
- 8291.45728643216 -2999.99999984439
- 8341.70854271357 -2999.99999992011
- 8391.95979899498 -2999.99999995939
- 8442.21105527638 -2999.99999997957
- 8492.46231155779 -2999.99999998982
- 8542.71356783920 -2999.99999999498
- 8592.96482412060 -2999.99999999755
- 8643.21608040201 -2999.99999999882
- 8693.46733668342 -2999.99999999943
- 8743.71859296482 -2999.99999999973
- 8793.96984924623 -2999.99999999987
- 8844.22110552764 -2999.99999999994
- 8894.47236180905 -2999.99999999997
- 8944.72361809045 -2999.99999999999
- 8994.97487437186 -2999.99999999999
- 9045.22613065327 -3000.00000000000
- 9095.47738693467 -3000.00000000000
- 9145.72864321608 -3000.00000000000
- 9195.97989949749 -3000.00000000000
- 9246.23115577889 -3000.00000000000
- 9296.48241206030 -3000.00000000000
- 9346.73366834171 -3000.00000000000
- 9396.98492462312 -3000.00000000000
- 9447.23618090452 -3000.00000000000
- 9497.48743718593 -3000.00000000000
- 9547.73869346734 -3000.00000000000
- 9597.98994974874 -3000.00000000000
- 9648.24120603015 -3000.00000000000
- 9698.49246231156 -3000.00000000000
- 9748.74371859296 -3000.00000000000
- 9798.99497487437 -3000.00000000000
- 9849.24623115578 -3000.00000000000
- 9899.49748743719 -3000.00000000000
- 9949.74874371859 -3000.00000000000
- 10000.0000000000 -3000.00000000000
-#
-# interface number 3 (ocean surface, top of the mesh)
-#
- 2
- 0 0
- 10000 0
-#
-# for each layer, we give the number of spectral elements in the vertical direction
-#
-#
-# layer number 1 (bottom layer)
-#
- 44
-#
-# layer number 2 (top layer)
-#
- 40
diff --git a/unused_older_examples_DATA_to_sort/interfaces_Abel_Brest_flat.dat b/unused_older_examples_DATA_to_sort/interfaces_Abel_Brest_flat.dat
deleted file mode 100644
index f400336..0000000
--- a/unused_older_examples_DATA_to_sort/interfaces_Abel_Brest_flat.dat
+++ /dev/null
@@ -1,36 +0,0 @@
-#
-# number of interfaces
-#
- 3
-#
-# for each interface below, we give the number of points and then x,y for each point
-#
-#
-# interface number 1 (bottom of the mesh)
-#
- 2
- 0 -10000
- 10000 -10000
-#
-# interface number 2 (shape of the ocean bottom, "middle" of the mesh)
-#
- 2
- 0 -3000
- 10000 -3000
-#
-# interface number 3 (ocean surface, top of the mesh)
-#
- 2
- 0 0
- 10000 0
-#
-# for each layer, we give the number of spectral elements in the vertical direction
-#
-#
-# layer number 1 (bottom layer)
-#
- 44
-#
-# layer number 2 (top layer)
-#
- 40
diff --git a/unused_older_examples_DATA_to_sort/topo_rouen_aniso.dat b/unused_older_examples_DATA_to_sort/topo_rouen_aniso.dat
deleted file mode 100644
index c60db61..0000000
--- a/unused_older_examples_DATA_to_sort/topo_rouen_aniso.dat
+++ /dev/null
@@ -1,3 +0,0 @@
-2
-0 0.06d0
-1 0.06d0
More information about the CIG-COMMITS
mailing list