[cig-commits] [commit] devel, master: Added Abel_Brest and canyon as formal examples. (8ab5d5f)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Wed Jun 18 15:18:27 PDT 2014
Repository : https://github.com/geodynamics/specfem2d
On branches: devel,master
Link : https://github.com/geodynamics/specfem2d/compare/fc67e6fd7ad890705b2b72b4b3c509accb22249e...e9ca46c40131588d89d7b0883250bc6584ce6b4c
>---------------------------------------------------------------
commit 8ab5d5f84ed460029beff37861acdc4a687e2ec8
Author: Carl Tape <carltape at gi.alaska.edu>
Date: Fri Feb 25 19:48:08 2011 +0000
Added Abel_Brest and canyon as formal examples.
>---------------------------------------------------------------
8ab5d5f84ed460029beff37861acdc4a687e2ec8
.../Par_file_Abel_Balanche_bathy_source_solid | 46 +++++++++++----------
.../SOURCE_Abel_Balanche_bathy_source_solid | 0
.../interfaces_Abel_Brest_bathymetry.dat | 0
{M2_UPPA => Abel_Brest}/process.sh | 4 +-
.../Par_file_Tromp2005 => canyon/Par_file_canyon | 48 +++++++++++-----------
.../SOURCE_canyon | 2 +-
.../mesh}/canyon_absorbing_surface_file | 0
.../mesh}/canyon_free_surface_file | 0
.../mesh}/canyon_materials_file | 0
.../Mesh_canyon => canyon/mesh}/canyon_mesh_file | 0
.../mesh}/canyon_nodes_coords_file | 0
{fluid_solid => canyon}/process.sh | 4 +-
12 files changed, 53 insertions(+), 51 deletions(-)
diff --git a/DATA_to_sort_older_examples/Par_file_Abel_Balanche_bathy_source_solid b/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid
similarity index 73%
copy from DATA_to_sort_older_examples/Par_file_Abel_Balanche_bathy_source_solid
copy to Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid
index 62b70e3..9b26a88 100644
--- a/DATA_to_sort_older_examples/Par_file_Abel_Balanche_bathy_source_solid
+++ b/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid
@@ -1,7 +1,10 @@
-
# title of job, and file that contains interface data
title = Test for Abel Balanche UMR 6538 Brest
-interfacesfile = interfaces_Abel_Brest_bathymetry.dat
+interfacesfile = ../interfaces_Abel_Brest_bathymetry.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.
@@ -10,10 +13,11 @@ nodes_coords_file = ./DATA/Mesh_canyon/canyon_nodes_coords_file
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 partitioning
nproc = 1 # number of processes
-partitioning_method = 1 # ascending order = 1, Metis = 2, Scotch = 3
+partitioning_method = 1 # SCOTCH = 3, ascending order (very bad idea) = 1
# geometry of the model (origin lower-left corner = 0,0) and mesh description
xmin = 0.d0 # abscissa of left side of the model
@@ -23,8 +27,12 @@ ngnod = 9 # number of control nodes per e
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
+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_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 = .true. # absorbing boundary active or not
@@ -38,22 +46,11 @@ 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
+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
-Qp_attenuation = 136.4376068115 # quality factor P for attenuation
-Qs_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
@@ -61,6 +58,7 @@ seismotype = 4 # record 1=displ 2=veloc 3=acce
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
@@ -93,14 +91,18 @@ subsamp = 1 # subsampling of color snapshot
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 = 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 models as
+# I: (model_number 1 rho Vp Vs 0 0 Qp Qs 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 Qs).
+# 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
# 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
diff --git a/DATA_to_sort_older_examples/SOURCE_Abel_Balanche_bathy_source_solid b/Abel_Brest/SOURCE_Abel_Balanche_bathy_source_solid
similarity index 100%
copy from DATA_to_sort_older_examples/SOURCE_Abel_Balanche_bathy_source_solid
copy to Abel_Brest/SOURCE_Abel_Balanche_bathy_source_solid
diff --git a/DATA_to_sort_older_examples/interfaces_Abel_Brest_bathymetry.dat b/Abel_Brest/interfaces_Abel_Brest_bathymetry.dat
similarity index 100%
copy from DATA_to_sort_older_examples/interfaces_Abel_Brest_bathymetry.dat
copy to Abel_Brest/interfaces_Abel_Brest_bathymetry.dat
diff --git a/M2_UPPA/process.sh b/Abel_Brest/process.sh
similarity index 88%
copy from M2_UPPA/process.sh
copy to Abel_Brest/process.sh
index 7cd1651..84ad87f 100755
--- a/M2_UPPA/process.sh
+++ b/Abel_Brest/process.sh
@@ -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_Abel_Balanche_bathy_source_solid Par_file
+ln -s ../SOURCE_Abel_Balanche_bathy_source_solid SOURCE
cd ../
# cleans output files
diff --git a/Tromp2005/Par_file_Tromp2005 b/canyon/Par_file_canyon
similarity index 74%
copy from Tromp2005/Par_file_Tromp2005
copy to canyon/Par_file_canyon
index c45ac7c..7f0c90c 100644
--- a/Tromp2005/Par_file_Tromp2005
+++ b/canyon/Par_file_canyon
@@ -1,19 +1,19 @@
# title of job, and file that contains interface data
-title = Tromp-Tape-Liu (GJI 2005)
-interfacesfile = ../interfaces_Tromp2005.dat
+title = Calcul Mexique Alejandro
+interfacesfile = ../interfaces_no_canyon.dat
# forward or adjoint simulation
SIMULATION_TYPE = 1 # 1 = forward, 2 = adjoint + kernels
-SAVE_FORWARD = .false. # save the last frame, needed for adjoint simulation
+SAVE_FORWARD = .true. # 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/ice_water_rock_1D/ice_water_rock_1D_elements # file containing elements
-nodes_coords_file = ./DATA/ice_water_rock_1D/ice_water_rock_1D_nodes # file containing node coordinates
-materials_file = ./DATA/ice_water_rock_1D/ice_water_rock_1D_material # file containing material index for each element
-free_surface_file = ./DATA/ice_water_rock_1D/ice_water_rock_1D_surface_free # file containing free surface
-absorbing_surface_file = ./DATA/ice_water_rock_1D/ice_water_rock_1D_surface_absorb # file containing absorbing surface
+read_external_mesh = .true.
+mesh_file = ./mesh/canyon_mesh_file # file containing the mesh
+nodes_coords_file = ./mesh/canyon_nodes_coords_file # file containing the nodes coordinates
+materials_file = ./mesh/canyon_materials_file # file containing the material number for each element
+free_surface_file = ./mesh/canyon_free_surface_file # file containing the free surface
+absorbing_surface_file = ./mesh/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 partitioning
@@ -22,11 +22,11 @@ partitioning_method = 1 # SCOTCH = 3, ascending order (
# geometry of the model (origin lower-left corner = 0,0) and mesh description
xmin = 0.d0 # abscissa of left side of the model
-xmax = 200000.d0 # abscissa of right side of the model
-nx = 80 # number of elements along X
+xmax = 19.d0 # abscissa of right side of the model
+nx = 95 # 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
+initialfield = .true. # use a plane wave as source or not
+add_Bielak_conditions = .true. # 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
@@ -35,7 +35,7 @@ Q0 = 1 # quality factor for viscous at
freq0 = 10 # frequency for viscous attenuation
p_sv = .true. # set the type of calculation (P-SV or SH/membrane waves)
-# absorbing boundary parameters
+# absorbing boundaries parameters
absorbing_conditions = .true. # absorbing boundary active or not
absorbbottom = .true.
absorbright = .true.
@@ -43,8 +43,8 @@ absorbtop = .false.
absorbleft = .true.
# time step parameters
-nt = 3000 # total number of time steps
-deltat = 2.0d-2 # duration of a time step
+nt = 5000 # total number of time steps
+deltat = 2.e-3 # duration of a time step
# source parameters
NSOURCES = 1 # number of sources [source info read in CMTSOLUTION file]
@@ -62,15 +62,15 @@ anglerec = 0.d0 # angle to rotate components at
rec_normal_to_surface = .false. # base anglerec normal to surface (external mesh and curve file needed)
# first receiver line
-nrec = 1 # number of receivers
-xdeb = 150000. # first receiver x in meters
-zdeb = 40000. # first receiver z in meters
-xfin = 70000. # last receiver x in meters (ignored if onlyone receiver)
-zfin = 0. # last receiver z in meters (ignored if onlyone receiver)
+nrec = 20 # number of receivers
+xdeb = 1. # first receiver x in meters
+zdeb = 9. # first receiver z in meters
+xfin = 18. # last receiver x in meters (ignored if onlyone receiver)
+zfin = 9. # 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 = 400 # display frequency in time steps
+NTSTEP_BETWEEN_OUTPUT_INFO = 250 # 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
@@ -94,7 +94,7 @@ nbmodels = 1 # nb of different models
# III: (model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qs).
# 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 2600.d0 5800.d0 3198.6d0 0 0 10.d0 10.d0 0 0 0 0 0 0
+1 1 1.d0 2.d0 1.d0 0 0 10.d0 10.d0 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 80 1 32 1
+1 95 1 45 1
diff --git a/DATA_to_sort_older_examples/SOURCE_canyon b/canyon/SOURCE_canyon
similarity index 94%
copy from DATA_to_sort_older_examples/SOURCE_canyon
copy to canyon/SOURCE_canyon
index 30eb2d6..f3256d2 100644
--- a/DATA_to_sort_older_examples/SOURCE_canyon
+++ b/canyon/SOURCE_canyon
@@ -5,7 +5,7 @@ zs = 3. # 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 = 1.0 # dominant source frequency (Hz) if not Dirac or Heaviside
-t0 = 0.0 # time shift when multi sources (if one source, must be zero)
+tshift = 0.0 # time shift when multi sources (if one source, must be zero)
angleforce = 30. # 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/DATA_to_sort_older_examples/Mesh_canyon/canyon_absorbing_surface_file b/canyon/mesh/canyon_absorbing_surface_file
similarity index 100%
rename from DATA_to_sort_older_examples/Mesh_canyon/canyon_absorbing_surface_file
rename to canyon/mesh/canyon_absorbing_surface_file
diff --git a/DATA_to_sort_older_examples/Mesh_canyon/canyon_free_surface_file b/canyon/mesh/canyon_free_surface_file
similarity index 100%
rename from DATA_to_sort_older_examples/Mesh_canyon/canyon_free_surface_file
rename to canyon/mesh/canyon_free_surface_file
diff --git a/DATA_to_sort_older_examples/Mesh_canyon/canyon_materials_file b/canyon/mesh/canyon_materials_file
similarity index 100%
rename from DATA_to_sort_older_examples/Mesh_canyon/canyon_materials_file
rename to canyon/mesh/canyon_materials_file
diff --git a/DATA_to_sort_older_examples/Mesh_canyon/canyon_mesh_file b/canyon/mesh/canyon_mesh_file
similarity index 100%
rename from DATA_to_sort_older_examples/Mesh_canyon/canyon_mesh_file
rename to canyon/mesh/canyon_mesh_file
diff --git a/DATA_to_sort_older_examples/Mesh_canyon/canyon_nodes_coords_file b/canyon/mesh/canyon_nodes_coords_file
similarity index 100%
rename from DATA_to_sort_older_examples/Mesh_canyon/canyon_nodes_coords_file
rename to canyon/mesh/canyon_nodes_coords_file
diff --git a/fluid_solid/process.sh b/canyon/process.sh
similarity index 92%
copy from fluid_solid/process.sh
copy to canyon/process.sh
index 04e46ca..1d4f827 100755
--- a/fluid_solid/process.sh
+++ b/canyon/process.sh
@@ -21,8 +21,8 @@ mkdir -p DATA
# sets up local DATA/ directory
cd DATA/
-ln -s ../Par_file_fluid_solid Par_file
-ln -s ../SOURCE_fluid_solid SOURCE
+ln -s ../Par_file_canyon Par_file
+ln -s ../SOURCE_canyon SOURCE
cd ../
# cleans output files
More information about the CIG-COMMITS
mailing list