[cig-commits] r21059 - seismo/2D/SPECFEM2D/trunk/EXAMPLES/salt_dome/Stacey_normal
xie.zhinan at geodynamics.org
xie.zhinan at geodynamics.org
Tue Nov 20 07:49:49 PST 2012
Author: xie.zhinan
Date: 2012-11-20 07:49:49 -0800 (Tue, 20 Nov 2012)
New Revision: 21059
Modified:
seismo/2D/SPECFEM2D/trunk/EXAMPLES/salt_dome/Stacey_normal/Par_file
seismo/2D/SPECFEM2D/trunk/EXAMPLES/salt_dome/Stacey_normal/process.sh
Log:
modify the process file for example of salt_dome
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/salt_dome/Stacey_normal/Par_file
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/salt_dome/Stacey_normal/Par_file 2012-11-20 15:11:35 UTC (rev 21058)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/salt_dome/Stacey_normal/Par_file 2012-11-20 15:49:49 UTC (rev 21059)
@@ -23,7 +23,7 @@
p_sv = .true. # set the type of calculation (P-SV or SH/membrane waves)
# time step parameters
-nt = 10000 # total number of time steps
+nt = 5000 # total number of time steps
deltat = 4.25e-4 # duration of a time step (see section "How to choose the time step" of the manual for how to do this)
USER_T0 = 0.0d0 # use this t0 as earliest starting time rather than the automatically calculated one
time_stepping_scheme = 1 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta
@@ -58,9 +58,9 @@
enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface
# display parameters
-NSTEP_BETWEEN_OUTPUT_INFO = 500 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_INFO = 100 # every how many time steps we display information about the simulation (costly, do not use a very small value)
####
-NSTEP_BETWEEN_OUTPUT_IMAGES = 500 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
cutsnaps = 1. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
#### for JPEG color images ####
output_color_image = .true. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
@@ -123,8 +123,8 @@
# absorbing boundary active or not
PML_BOUNDARY_CONDITIONS = .false.
NELEM_PML_THICKNESS = 3
-STACEY_ABSORBING_CONDITIONS = .true.
-ADD_SPRING_TO_STACEY = .true.
+STACEY_ABSORBING_CONDITIONS = .false.
+ADD_SPRING_TO_STACEY = .false.
# for horizontal periodic conditions: detect common points between left and right edges
ADD_PERIODIC_CONDITIONS = .false.
@@ -140,13 +140,13 @@
# data concerning mesh, when generated using third-party app (more info in README)
# (see also absorbing_conditions above)
-mesh_file = modelY1_mesh_file # file containing the mesh
-nodes_coords_file = modelY1_nodes_coords_file # file containing the nodes coordinates
-materials_file = modelY1_materials_file # file containing the material number for each element
-free_surface_file = modelY1_free_surface_file # file containing the free surface
-absorbing_surface_file = modelY1_absorbing_surface_file # file containing the absorbing surface
-CPML_element_file = Elements_CPML_list # file containing the CPML element numbers
-tangential_detection_curve_file = courbe_eros_nodes # file containing the curve delimiting the velocity model
+mesh_file = ./DATA/modelY1_mesh_file # file containing the mesh
+nodes_coords_file = ./DATA/modelY1_nodes_coords_file # file containing the nodes coordinates
+materials_file = ./DATA/modelY1_materials_file # file containing the material number for each element
+free_surface_file = ./DATA/modelY1_free_surface_file # file containing the free surface
+absorbing_surface_file = ./DATA/modelY1_absorbing_surface_file # file containing the absorbing surface
+CPML_element_file = ./DATA/Elements_CPML_list # file containing the CPML element numbers
+tangential_detection_curve_file = ./DATA/courbe_eros_nodes # file containing the curve delimiting the velocity model
#-----------------------------------------------------------------------------
# PARAMETERS FOR INTERNAL MESHING
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/salt_dome/Stacey_normal/process.sh
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/salt_dome/Stacey_normal/process.sh 2012-11-20 15:11:35 UTC (rev 21058)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/salt_dome/Stacey_normal/process.sh 2012-11-20 15:49:49 UTC (rev 21059)
@@ -24,6 +24,12 @@
cp ../Par_file Par_file
#cp ../interfaces_elastic_analytic.dat .
cp ../SOURCE SOURCE
+cp ../modelY1_mesh_file modelY1_mesh_file
+cp ../modelY1_nodes_coords_file modelY1_nodes_coords_file
+cp ../modelY1_materials_file modelY1_materials_file
+cp ../modelY1_free_surface_file modelY1_free_surface_file
+cp ../modelY1_absorbing_surface_file modelY1_absorbing_surface_file
+cp ../courbe_eros_nodes courbe_eros_nodes
cd ../
# cleans output files
More information about the CIG-COMMITS
mailing list