[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