[cig-commits] [commit] devel,master: modified USER_T0 in the Tromp2005 and Tape2007 examples; also updated the instructions for these examples. added a few 'quick test' lines to the main README (4a2914d)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Wed Jun 18 15:21:06 PDT 2014
Repository : https://github.com/geodynamics/specfem2d
On branches: devel,master
Link : https://github.com/geodynamics/specfem2d/compare/fc67e6fd7ad890705b2b72b4b3c509accb22249e...e9ca46c40131588d89d7b0883250bc6584ce6b4c
>---------------------------------------------------------------
commit 4a2914d4321e209d36ddfcdfd03090d6cb188260
Author: Carl Tape <carltape at gi.alaska.edu>
Date: Fri Apr 20 17:32:10 2012 +0000
modified USER_T0 in the Tromp2005 and Tape2007 examples; also updated the instructions for these examples. added a few 'quick test' lines to the main README
>---------------------------------------------------------------
4a2914d4321e209d36ddfcdfd03090d6cb188260
Tape2007/Par_file_Tape2007_132rec_checker | 2 +-
Tape2007/Par_file_Tape2007_onerec | 2 +-
Tape2007/README | 24 +++++++++++++---------
Tromp2005/Par_file_Tromp2005 | 4 ++--
Tromp2005/Par_file_Tromp2005_s100 | 2 +-
Tromp2005/README | 33 +++++++++++++++++--------------
6 files changed, 38 insertions(+), 29 deletions(-)
diff --git a/Tape2007/Par_file_Tape2007_132rec_checker b/Tape2007/Par_file_Tape2007_132rec_checker
index ea927a7..dfc556b 100644
--- a/Tape2007/Par_file_Tape2007_132rec_checker
+++ b/Tape2007/Par_file_Tape2007_132rec_checker
@@ -25,7 +25,7 @@ p_sv = .false. # set the type of calculation (
# time step parameters
nt = 4800 # total number of time steps
deltat = 6.0d-2 # duration of a time step
-USER_T0 = 0.0d0 # use this t0 as earliest starting time rather than the automatically calculated one
+USER_T0 = 48.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 4th-order 4-stage Runge-Kutta
# source parameters
diff --git a/Tape2007/Par_file_Tape2007_onerec b/Tape2007/Par_file_Tape2007_onerec
index b39ecdc..f0b208d 100644
--- a/Tape2007/Par_file_Tape2007_onerec
+++ b/Tape2007/Par_file_Tape2007_onerec
@@ -25,7 +25,7 @@ p_sv = .false. # set the type of calculation (
# time step parameters
nt = 3000 # total number of time steps
deltat = 6.0d-2 # duration of a time step
-USER_T0 = 0.0d0 # use this t0 as earliest starting time rather than the automatically calculated one
+USER_T0 = 48.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 4th-order 4-stage Runge-Kutta
# source parameters
diff --git a/Tape2007/README b/Tape2007/README
index ec96937..8890f7d 100644
--- a/Tape2007/README
+++ b/Tape2007/README
@@ -18,7 +18,7 @@ checkerboard model (external file model_velocity.dat_checker)
---------------------------
-The examples in Tape-Liu-Tromp (GJI 2007) assume a simulation start time of t=0, and a source start time of t = t0 = 48 s. SPECFEM2D assumes that the simulation starts at t = -t0 and the source start time is t = 0. In order to generate the same time series as in Tape2007, the user must first modify USER_T0 in constants.h
+The examples in Tape-Liu-Tromp (GJI 2007) assume a simulation start time of t=0, and a source start time of t = t0 = 48 s. SPECFEM2D assumes that the simulation starts at t = -t0 and the source start time is t = 0. Thus, USER_T0 is set to 48.0 in the Par_file
TO RUN:
@@ -27,24 +27,30 @@ TO RUN:
1. in SPECFEM2D root directory, configure, e.g.,
./configure FC=gfortran
-2. modify setup/constants.h file, setting USER_TO = 48.0 s
-
-3. compile:
+2. compile:
make all
-4. cd EXAMPLES/Tape2007/
+3. cd EXAMPLES/Tape2007/
-5. execute script to run mesher and solver for the SH/membrane case:
+4. execute script to run mesher and solver for the SH/membrane case:
./process.sh
-6. check out the output files in the local directory OUTPUT_FILES
+5. check out the output files in the local directory OUTPUT_FILES
- optional: try plotting the wavefield using the script
- SPECFEM2D/UTILS/visualization/plot_wavefield.pl
+ optional: try plotting the wavefield using the script plot_wavefield.pl
+ (note: you will need to change a path variable in plot_wavefield.pl)
+ cd OUTPUT_FILES
+ ln -s ../../../UTILS/visualization/plot_wavefield.pl .
+ plot_wavefield.pl 400/4800/400 400/2800/800 0/480/0/480 120/20/120/20 -3/-3/-3 6/6/6 -48.0/0.06 0/1/0/0/1 1.7/1/0 1/0/1/120 Tape2007 onerec_homo
---------------------------
After running the example above, try the second example:
./process_checker.sh
+ optional: try plotting the wavefield using the script plot_wavefield.pl
+ cd OUTPUT_FILES
+ ln -s ../../../UTILS/visualization/plot_wavefield.pl .
+ plot_wavefield.pl 400/4800/400 800/2000/400 0/480/0/480 120/20/120/20 -3/-3/-3 6/6/6 -48.0/0.06 0/1/0/0/1 1.7/1/0 1/0/1/120 Tape2007 132rec_checker
+
---------------------------
diff --git a/Tromp2005/Par_file_Tromp2005 b/Tromp2005/Par_file_Tromp2005
index b027f25..d78cf3e 100644
--- a/Tromp2005/Par_file_Tromp2005
+++ b/Tromp2005/Par_file_Tromp2005
@@ -20,12 +20,12 @@ ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasti
ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
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)
+p_sv = .false. # set the type of calculation (P-SV or SH/membrane waves)
# time step parameters
nt = 3000 # total number of time steps
deltat = 2.0d-2 # duration of a time step
-USER_T0 = 0.0d0 # use this t0 as earliest starting time rather than the automatically calculated one
+USER_T0 = 8.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 4th-order 4-stage Runge-Kutta
# source parameters
diff --git a/Tromp2005/Par_file_Tromp2005_s100 b/Tromp2005/Par_file_Tromp2005_s100
index c6ea35f..3a2e3aa 100644
--- a/Tromp2005/Par_file_Tromp2005_s100
+++ b/Tromp2005/Par_file_Tromp2005_s100
@@ -25,7 +25,7 @@ p_sv = .true. # set the type of calculation (
# time step parameters
nt = 3000 # total number of time steps
deltat = 2.0d-4 # duration of a time step
-USER_T0 = 0.0d0 # use this t0 as earliest starting time rather than the automatically calculated one
+USER_T0 = 8.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 4th-order 4-stage Runge-Kutta
# source parameters
diff --git a/Tromp2005/README b/Tromp2005/README
index 7d500d0..1b488e9 100644
--- a/Tromp2005/README
+++ b/Tromp2005/README
@@ -12,7 +12,7 @@ Demonstration of the kernel calculation is not yet added here.
---------------------------
-The examples in Tromp-Tape-Liu (GJI 2005) assume a simulation start time of t=0, and a source start time of t = t0 = 8.0 s. SPECFEM2D assumes that the simulation starts at t = -t0 and the source start time is t = 0. In order to generate the same time series as in Tromp2005, the user must first modify USER_T0 in constants.h
+The examples in Tromp-Tape-Liu (GJI 2005) assume a simulation start time of t=0, and a source start time of t = t0 = 8.0 s. SPECFEM2D assumes that the simulation starts at t = -t0 and the source start time is t = 0. Thus, USER_T0 is set to 8.0 in the Par_file
TO RUN:
@@ -21,35 +21,38 @@ TO RUN:
1. in SPECFEM2D root directory, configure, e.g.,
./configure FC=gfortran
-2. modify constants.h file, setting USER_TO = 8.0 s
-
-3. compile:
+2. compile:
make all
-4. cd EXAMPLES/Tromp2005/
+3. cd EXAMPLES/Tromp2005/
-5. execute script to run mesher and solver for the PSV case:
+4. execute script to run mesher and solver for the PSV case:
./process.sh
-6. check out the output files in the local directory OUTPUT_FILES
+5. check out the output files in the local directory OUTPUT_FILES
The seismograms (S0001.AA.BXX.semd,S0001.AA.BXZ.semd) will be similar to those in Tromp2005, Figure 7b, but with an 8-second time shift due to the difference in origin time convention.
- optional: try plotting the wavefield using the script
- SPECFEM2D/UTILS/visualization/plot_wavefield.pl
+ optional: try plotting the wavefield using the script plot_wavefield.pl
+ (note: you will need to change a path variable in plot_wavefield.pl)
+ cd OUTPUT_FILES
+ ln -s ../../../UTILS/visualization/plot_wavefield.pl .
+ plot_wavefield.pl 400/2800/400 400/2000/400 0/200/0/80 50/10/40/10 -3/-3/-3 4/4/4 -8.0/0.02 1/0/1/1/1 3.0/1/0 1/0/1/200 Tromp2005 PSV_homo
-7. mv OUTPUT_FILES OUTPUT_FILES_PSV
+6. mv OUTPUT_FILES OUTPUT_FILES_PSV
-8. set p_sv = .false. in Par_file_Tromp2005
+7. set p_sv = .false. in Par_file_Tromp2005
-9. execute script to run mesher and solver for the SH case:
+8. execute script to run mesher and solver for the SH case:
./process.sh
-10. check out the output files in the local directory OUTPUT_FILES
+9. check out the output files in the local directory OUTPUT_FILES
The seismogram (S0001.AA.BXY.semd) will be similar to those in Tromp2005, Figure 2b, but with an 8-second time shift due to the difference in origin time convention.
- optional: try plotting the wavefield using the script
- SPECFEM2D/UTILS/visualization/plot_wavefield.pl
+ optional: try plotting the wavefield using the script plot_wavefield.pl
+ cd OUTPUT_FILES
+ ln -s ../../../UTILS/visualization/plot_wavefield.pl .
+ plot_wavefield.pl 400/2800/400 400/2000/400 0/200/0/80 50/10/40/10 -2/-2/-2 1/1/1 -8.0/0.02 0/1/0/0/1 3.0/1/0 1/0/1/200 Tromp2005 SH_homo
---------------------------
More information about the CIG-COMMITS
mailing list