[cig-commits] r18053 - in seismo/2D/SPECFEM2D/trunk: DATA EXAMPLES/Abel_Brest EXAMPLES/M2_UPPA EXAMPLES/Tape2007 EXAMPLES/Tromp2005 EXAMPLES/canyon EXAMPLES/fluid_solid
carltape at geodynamics.org
carltape at geodynamics.org
Mon Mar 7 12:30:31 PST 2011
Author: carltape
Date: 2011-03-07 12:30:31 -0800 (Mon, 07 Mar 2011)
New Revision: 18053
Modified:
seismo/2D/SPECFEM2D/trunk/DATA/Par_file.in
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid
seismo/2D/SPECFEM2D/trunk/EXAMPLES/M2_UPPA/Par_file_M2_UPPA
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_132rec
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_onerec
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005_s100
seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/Par_file_canyon
seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/Par_file_fluid_solid
Log:
updated Par_file for all six examples to reflect recent changes; tested to make sure each one works
Modified: seismo/2D/SPECFEM2D/trunk/DATA/Par_file.in
===================================================================
--- seismo/2D/SPECFEM2D/trunk/DATA/Par_file.in 2011-03-07 20:28:33 UTC (rev 18052)
+++ seismo/2D/SPECFEM2D/trunk/DATA/Par_file.in 2011-03-07 20:30:31 UTC (rev 18053)
@@ -1,4 +1,3 @@
-
# title of job
title = Test for M2 UPPA
@@ -89,7 +88,7 @@
# PARAMETERS FOR EXTERNAL MESHING
# data concerning mesh, when generated using third-party app (more info in README)
-# (see also absorbing_conditions above)
+# (see also absorbing_conditions above)
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
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid 2011-03-07 20:28:33 UTC (rev 18052)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid 2011-03-07 20:30:31 UTC (rev 18053)
@@ -1,28 +1,14 @@
-# title of job, and file that contains interface data
+# title of job
title = Test for Abel Balanche UMR 6538 Brest
-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.
-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 partitioning
nproc = 1 # number of processes
-partitioning_method = 1 # SCOTCH = 3, ascending order (very bad idea) = 1
+partitioning_method = 3 # 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
-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
@@ -34,13 +20,6 @@
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
-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
@@ -90,8 +69,9 @@
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)
+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)
# velocity and density models
nbmodels = 2 # nb of different models
@@ -103,6 +83,42 @@
# 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
+
+# external mesh or not
+read_external_mesh = .false.
+
+# absorbing boundary active or not
+absorbing_conditions = .true.
+
+#-----------------------------------------------------------------------------
+# PARAMETERS FOR EXTERNAL MESHING
+
+# data concerning mesh, when generated using third-party app (more info in README)
+# (see also absorbing_conditions above)
+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 FOR INTERNAL MESHING
+
+# file containing interfaces for internal mesh
+interfacesfile = ../interfaces_Abel_Brest_bathymetry.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
+
+# absorbing boundary parameters (see absorbing_conditions above)
+absorbbottom = .true.
+absorbright = .true.
+absorbtop = .false.
+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
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/M2_UPPA/Par_file_M2_UPPA
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/M2_UPPA/Par_file_M2_UPPA 2011-03-07 20:28:33 UTC (rev 18052)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/M2_UPPA/Par_file_M2_UPPA 2011-03-07 20:30:31 UTC (rev 18053)
@@ -1,29 +1,14 @@
-
-# title of job, and file that contains interface data
+# title of job
title = Test for M2 UPPA
-interfacesfile = ../interfaces_M2_UPPA_curved.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 partitioning
nproc = 1 # number of processes
partitioning_method = 3 # 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
-xmax = 4000.d0 # abscissa of right side of the model
-nx = 80 # 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
@@ -35,13 +20,6 @@
freq0 = 10 # frequency for viscous attenuation
p_sv = .true. # set the type of calculation (P-SV or SH/membrane waves)
-# absorbing boundary parameters
-absorbing_conditions = .true. # absorbing boundary active or not
-absorbbottom = .true.
-absorbright = .true.
-absorbtop = .false.
-absorbleft = .true.
-
# time step parameters
nt = 1600 # total number of time steps
deltat = 1.d-3 # duration of a time step
@@ -61,7 +39,7 @@
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
+# 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
@@ -83,8 +61,9 @@
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)
+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)
# velocity and density models
nbmodels = 4 # nb of different models
@@ -98,6 +77,42 @@
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
+
+# external mesh or not
+read_external_mesh = .false.
+
+# absorbing boundary active or not
+absorbing_conditions = .true.
+
+#-----------------------------------------------------------------------------
+# PARAMETERS FOR EXTERNAL MESHING
+
+# data concerning mesh, when generated using third-party app (more info in README)
+# (see also absorbing_conditions above)
+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 FOR INTERNAL MESHING
+
+# file containing interfaces for internal mesh
+interfacesfile = ../interfaces_M2_UPPA_curved.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
+
+# absorbing boundary parameters (see absorbing_conditions above)
+absorbbottom = .true.
+absorbright = .true.
+absorbtop = .false.
+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
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_132rec
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_132rec 2011-03-07 20:28:33 UTC (rev 18052)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_132rec 2011-03-07 20:30:31 UTC (rev 18053)
@@ -1,29 +1,14 @@
-
-# title of job, and file that contains interface data
+# title of job
title = Tape-Liu-Tromp (GJI 2007)
-interfacesfile = ../interfaces_Tape2007.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/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
-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 # SCOTCH = 3, ascending order (very bad idea) = 1
+partitioning_method = 3 # 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
-xmax = 480000.d0 # abscissa of right side of the model
-nx = 40 # 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
@@ -35,13 +20,6 @@
freq0 = 10 # frequency for viscous attenuation
p_sv = .false. # set the type of calculation (P-SV or SH/membrane waves)
-# absorbing boundary parameters
-absorbing_conditions = .true. # absorbing boundary active or not
-absorbbottom = .true.
-absorbright = .true.
-absorbtop = .true.
-absorbleft = .true.
-
# time step parameters
nt = 4800 # total number of time steps
deltat = 6.0d-2 # duration of a time step
@@ -869,8 +847,9 @@
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)
+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)
# velocity and density models
nbmodels = 1 # nb of different models
@@ -881,6 +860,42 @@
# 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 3500.0d0 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
+
+# external mesh or not
+read_external_mesh = .false.
+
+# absorbing boundary active or not
+absorbing_conditions = .true.
+
+#-----------------------------------------------------------------------------
+# PARAMETERS FOR EXTERNAL MESHING
+
+# data concerning mesh, when generated using third-party app (more info in README)
+# (see also absorbing_conditions above)
+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
+tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
+
+#-----------------------------------------------------------------------------
+# PARAMETERS FOR INTERNAL MESHING
+
+# file containing interfaces for internal mesh
+interfacesfile = ../interfaces_Tape2007.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 = 480000.d0 # abscissa of right side of the model
+nx = 40 # number of elements along X
+
+# absorbing boundary parameters (see absorbing_conditions above)
+absorbbottom = .true.
+absorbright = .true.
+absorbtop = .true.
+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
1 40 1 40 1
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_onerec
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_onerec 2011-03-07 20:28:33 UTC (rev 18052)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_onerec 2011-03-07 20:30:31 UTC (rev 18053)
@@ -1,29 +1,14 @@
-
-# title of job, and file that contains interface data
+# title of job
title = Tape-Liu-Tromp (GJI 2007)
-interfacesfile = ../interfaces_Tape2007.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/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
-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 # SCOTCH = 3, ascending order (very bad idea) = 1
+partitioning_method = 3 # 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
-xmax = 480000.d0 # abscissa of right side of the model
-nx = 40 # 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
@@ -35,13 +20,6 @@
freq0 = 10 # frequency for viscous attenuation
p_sv = .false. # set the type of calculation (P-SV or SH/membrane waves)
-# absorbing boundary parameters
-absorbing_conditions = .true. # absorbing boundary active or not
-absorbbottom = .true.
-absorbright = .true.
-absorbtop = .true.
-absorbleft = .true.
-
# time step parameters
nt = 3000 # total number of time steps
deltat = 6.0d-2 # duration of a time step
@@ -84,8 +62,9 @@
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)
+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)
# velocity and density models
nbmodels = 1 # nb of different models
@@ -96,6 +75,42 @@
# 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 3500.0d0 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
+
+# external mesh or not
+read_external_mesh = .false.
+
+# absorbing boundary active or not
+absorbing_conditions = .true.
+
+#-----------------------------------------------------------------------------
+# PARAMETERS FOR EXTERNAL MESHING
+
+# data concerning mesh, when generated using third-party app (more info in README)
+# (see also absorbing_conditions above)
+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
+tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
+
+#-----------------------------------------------------------------------------
+# PARAMETERS FOR INTERNAL MESHING
+
+# file containing interfaces for internal mesh
+interfacesfile = ../interfaces_Tape2007.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 = 480000.d0 # abscissa of right side of the model
+nx = 40 # number of elements along X
+
+# absorbing boundary parameters (see absorbing_conditions above)
+absorbbottom = .true.
+absorbright = .true.
+absorbtop = .true.
+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
1 40 1 40 1
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005 2011-03-07 20:28:33 UTC (rev 18052)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005 2011-03-07 20:30:31 UTC (rev 18053)
@@ -1,29 +1,14 @@
-
-# title of job, and file that contains interface data
+# title of job
title = Tromp-Tape-Liu (GJI 2005)
-interfacesfile = ../interfaces_Tromp2005.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/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
-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 # SCOTCH = 3, ascending order (very bad idea) = 1
+partitioning_method = 3 # 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
-xmax = 200000.d0 # abscissa of right side of the model
-nx = 80 # 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
@@ -35,13 +20,6 @@
freq0 = 10 # frequency for viscous attenuation
p_sv = .true. # set the type of calculation (P-SV or SH/membrane waves)
-# absorbing boundary 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 = 2.0d-2 # duration of a time step
@@ -83,8 +61,9 @@
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)
+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)
# velocity and density models
nbmodels = 1 # nb of different models
@@ -95,6 +74,42 @@
# 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
+
+# external mesh or not
+read_external_mesh = .false.
+
+# absorbing boundary active or not
+absorbing_conditions = .true.
+
+#-----------------------------------------------------------------------------
+# PARAMETERS FOR EXTERNAL MESHING
+
+# data concerning mesh, when generated using third-party app (more info in README)
+# (see also absorbing_conditions above)
+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
+tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
+
+#-----------------------------------------------------------------------------
+# PARAMETERS FOR INTERNAL MESHING
+
+# file containing interfaces for internal mesh
+interfacesfile = ../interfaces_Tromp2005.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 = 200000.d0 # abscissa of right side of the model
+nx = 80 # number of elements along X
+
+# absorbing boundary parameters (see absorbing_conditions above)
+absorbbottom = .true.
+absorbright = .true.
+absorbtop = .false.
+absorbleft = .true.
+
# 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
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005_s100
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005_s100 2011-03-07 20:28:33 UTC (rev 18052)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005_s100 2011-03-07 20:30:31 UTC (rev 18053)
@@ -1,29 +1,14 @@
+# title of job
+title = Tromp-Tape-Liu (GJI 2005)
-# title of job, and file that contains interface data
-title = Tromp-Tape-Liu (GJI 2005) (scaled x100)
-interfacesfile = ../interfaces_Tromp2005_s100.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/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
-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 # SCOTCH = 3, ascending order (very bad idea) = 1
+partitioning_method = 3 # 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
-xmax = 2000.d0 # abscissa of right side of the model
-nx = 80 # 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
@@ -35,13 +20,6 @@
freq0 = 10 # frequency for viscous attenuation
p_sv = .true. # set the type of calculation (P-SV or SH/membrane waves)
-# absorbing boundary 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 = 2.0d-4 # duration of a time step
@@ -83,8 +61,9 @@
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)
+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)
# velocity and density models
nbmodels = 1 # nb of different models
@@ -95,6 +74,42 @@
# 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
+
+# external mesh or not
+read_external_mesh = .false.
+
+# absorbing boundary active or not
+absorbing_conditions = .true.
+
+#-----------------------------------------------------------------------------
+# PARAMETERS FOR EXTERNAL MESHING
+
+# data concerning mesh, when generated using third-party app (more info in README)
+# (see also absorbing_conditions above)
+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
+tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
+
+#-----------------------------------------------------------------------------
+# PARAMETERS FOR INTERNAL MESHING
+
+# file containing interfaces for internal mesh
+interfacesfile = ../interfaces_Tromp2005.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 = 2000.d0 # abscissa of right side of the model
+nx = 80 # number of elements along X
+
+# absorbing boundary parameters (see absorbing_conditions above)
+absorbbottom = .true.
+absorbright = .true.
+absorbtop = .false.
+absorbleft = .true.
+
# 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
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/Par_file_canyon
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/Par_file_canyon 2011-03-07 20:28:33 UTC (rev 18052)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/Par_file_canyon 2011-03-07 20:30:31 UTC (rev 18053)
@@ -1,29 +1,14 @@
-
-# title of job, and file that contains interface data
+# title of job
title = Calcul Mexique Alejandro
-interfacesfile = ../interfaces_no_canyon.dat
# forward or adjoint simulation
SIMULATION_TYPE = 1 # 1 = forward, 2 = adjoint + kernels
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 = .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
nproc = 1 # number of processes
-partitioning_method = 1 # SCOTCH = 3, ascending order (very bad idea) = 1
+partitioning_method = 3 # 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
-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 = .true. # use a plane wave as source or not
add_Bielak_conditions = .true. # add Bielak conditions or not if initial plane wave
@@ -35,13 +20,6 @@
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
-absorbbottom = .true.
-absorbright = .true.
-absorbtop = .false.
-absorbleft = .true.
-
# time step parameters
nt = 5000 # total number of time steps
deltat = 2.e-3 # duration of a time step
@@ -83,8 +61,9 @@
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)
+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)
# velocity and density models
nbmodels = 1 # nb of different models
@@ -95,6 +74,42 @@
# 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 10.d0 10.d0 0 0 0 0 0 0
-# define the different regions of the model in the (nx,nz) spectral element mesh
+
+# external mesh or not
+read_external_mesh = .true.
+
+# absorbing boundary active or not
+absorbing_conditions = .true.
+
+#-----------------------------------------------------------------------------
+# PARAMETERS FOR EXTERNAL MESHING
+
+# data concerning mesh, when generated using third-party app (more info in README)
+# (see also absorbing_conditions above)
+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 FOR INTERNAL MESHING
+
+# file containing interfaces for internal mesh
+interfacesfile = ../interfaces_no_canyon.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 = 19.d0 # abscissa of right side of the model
+nx = 95 # number of elements along X
+
+# absorbing boundary parameters (see absorbing_conditions above)
+absorbbottom = .true.
+absorbright = .true.
+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
1 95 1 45 1
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/Par_file_fluid_solid
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/Par_file_fluid_solid 2011-03-07 20:28:33 UTC (rev 18052)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/Par_file_fluid_solid 2011-03-07 20:30:31 UTC (rev 18053)
@@ -1,29 +1,14 @@
-
-# title of job, and file that contains interface data
+# title of job
title = External mesh with fluid and solid
-interfacesfile = interfaces_M2_UPPA_curved.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 = .true.
-mesh_file = ./mesh/mesh # file containing the mesh
-nodes_coords_file = ./mesh/nodes_coords # file containing the nodes coordinates
-materials_file = ./mesh/mat # file containing the material number for each element
-free_surface_file = ./mesh/surface_free # file containing the free surface
-absorbing_surface_file = ./mesh/surface_abs # 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 # SCOTCH = 3, ascending order (very bad idea) = 1
+partitioning_method = 3 # 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
-xmax = 4000.d0 # abscissa of right side of the model
-nx = 80 # 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
@@ -35,13 +20,6 @@
freq0 = 10 # frequency for viscous attenuation
p_sv = .true. # set the type of calculation (P-SV or SH/membrane waves)
-# absorbing boundary parameters
-absorbing_conditions = .true. # absorbing boundary active or not
-absorbbottom = .true.
-absorbright = .true.
-absorbtop = .false.
-absorbleft = .true.
-
# time step parameters
nt = 18000 # total number of time steps
deltat = 0.25d-3 # duration of a time step
@@ -83,8 +61,9 @@
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)
+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)
# velocity and density models
nbmodels = 5 # nb of different models
@@ -100,6 +79,42 @@
4 3 2200.d0 786.3d0 0.4 2.0 1d-11 0.0 1d-11 5.341d9 2d9 3d9 0.0d-4 3.204d9 10.d0
5 2 2500.d0 169.d9 122.d9 0.d0 169.d9 0.d0 75.3d9 0 0 0 0 0 0
#4 1 2200.d0 2200.d0 1343.375d0 0 0 10.d0 10.d0 0 0 0 0 0 0
+
+# external mesh or not
+read_external_mesh = .true.
+
+# absorbing boundary active or not
+absorbing_conditions = .true.
+
+#-----------------------------------------------------------------------------
+# PARAMETERS FOR EXTERNAL MESHING
+
+# data concerning mesh, when generated using third-party app (more info in README)
+# (see also absorbing_conditions above)
+mesh_file = ./mesh/mesh # file containing the mesh
+nodes_coords_file = ./mesh/nodes_coords # file containing the nodes coordinates
+materials_file = ./mesh/mat # file containing the material number for each element
+free_surface_file = ./mesh/surface_free # file containing the free surface
+absorbing_surface_file = ./mesh/surface_abs # file containing the absorbing surface
+tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
+
+#-----------------------------------------------------------------------------
+# PARAMETERS FOR INTERNAL MESHING
+
+# file containing interfaces for internal mesh
+interfacesfile = interfaces_M2_UPPA_curved.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
+
+# absorbing boundary parameters (see absorbing_conditions above)
+absorbbottom = .true.
+absorbright = .true.
+absorbtop = .false.
+absorbleft = .true.
+
# define the different regions of the model in the (nx,nz) spectral element mesh
nbregions = 6 # nb of regions and model number for each
1 80 1 20 1
More information about the CIG-COMMITS
mailing list