[cig-commits] [commit] devel, master: modify the process file for example of salt_dome (3ef80f5)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/specfem2d

On branches: devel,master
Link       : https://github.com/geodynamics/specfem2d/compare/fc67e6fd7ad890705b2b72b4b3c509accb22249e...e9ca46c40131588d89d7b0883250bc6584ce6b4c

>---------------------------------------------------------------

commit 3ef80f51b378a2d1f2d27649973e9b9c73af19b9
Author: Xie Zhinan <xiezhinan1984 at gmail.com>
Date:   Tue Nov 20 15:49:49 2012 +0000

    modify the process file for example of salt_dome


>---------------------------------------------------------------

3ef80f51b378a2d1f2d27649973e9b9c73af19b9
 salt_dome/Stacey_normal/Par_file   | 24 ++++++++++++------------
 salt_dome/Stacey_normal/process.sh |  6 ++++++
 2 files changed, 18 insertions(+), 12 deletions(-)

diff --git a/salt_dome/Stacey_normal/Par_file b/salt_dome/Stacey_normal/Par_file
index 6985739..7fc8ca1 100644
--- a/salt_dome/Stacey_normal/Par_file
+++ b/salt_dome/Stacey_normal/Par_file
@@ -23,7 +23,7 @@ freq0                           =  10            # frequency for viscous attenua
 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 @@ zfin                            = 1100.  # last receiver z in meters (ignored if
 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 @@ read_external_mesh              = .true.
 # 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 @@ PERIODIC_DETECT_TOL             = 3.3334d-6
 
 # 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
diff --git a/salt_dome/Stacey_normal/process.sh b/salt_dome/Stacey_normal/process.sh
index 10a977f..fe8af15 100755
--- a/salt_dome/Stacey_normal/process.sh
+++ b/salt_dome/Stacey_normal/process.sh
@@ -24,6 +24,12 @@ cd DATA/
 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