[cig-commits] r19955 - in seismo/2D/SPECFEM2D/trunk: . EXAMPLES/Tape2007 EXAMPLES/Tromp2005 UTILS/visualization
carltape at geodynamics.org
carltape at geodynamics.org
Fri Apr 20 10:32:10 PDT 2012
Author: carltape
Date: 2012-04-20 10:32:10 -0700 (Fri, 20 Apr 2012)
New Revision: 19955
Modified:
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_132rec_checker
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_onerec
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/README
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005_s100
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/README
seismo/2D/SPECFEM2D/trunk/README
seismo/2D/SPECFEM2D/trunk/UTILS/visualization/plot_wavefield.pl
Log:
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
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_132rec_checker
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_132rec_checker 2012-04-19 11:40:25 UTC (rev 19954)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_132rec_checker 2012-04-20 17:32:10 UTC (rev 19955)
@@ -25,7 +25,7 @@
# 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
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_onerec
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_onerec 2012-04-19 11:40:25 UTC (rev 19954)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_onerec 2012-04-20 17:32:10 UTC (rev 19955)
@@ -25,7 +25,7 @@
# 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
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/README
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/README 2012-04-19 11:40:25 UTC (rev 19954)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/README 2012-04-20 17:32:10 UTC (rev 19955)
@@ -18,7 +18,7 @@
---------------------------
-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 @@
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
+
---------------------------
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005 2012-04-19 11:40:25 UTC (rev 19954)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005 2012-04-20 17:32:10 UTC (rev 19955)
@@ -20,12 +20,12 @@
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
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005_s100
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005_s100 2012-04-19 11:40:25 UTC (rev 19954)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005_s100 2012-04-20 17:32:10 UTC (rev 19955)
@@ -25,7 +25,7 @@
# 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
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/README
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/README 2012-04-19 11:40:25 UTC (rev 19954)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/README 2012-04-20 17:32:10 UTC (rev 19955)
@@ -12,7 +12,7 @@
---------------------------
-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 @@
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
---------------------------
Modified: seismo/2D/SPECFEM2D/trunk/README
===================================================================
--- seismo/2D/SPECFEM2D/trunk/README 2012-04-19 11:40:25 UTC (rev 19954)
+++ seismo/2D/SPECFEM2D/trunk/README 2012-04-20 17:32:10 UTC (rev 19955)
@@ -1,6 +1,11 @@
-
Instructions on how to install and use SPECFEM2D are
available in the PDF manual located in directory doc/USER_MANUAL.
Main developers: Dimitri Komatitsch and Jeroen Tromp
+For a quick test, run the default example with these commands:
+ ./configure FC=gfortran
+ make all
+ ./bin/xmeshfem2D
+ ./bin/xspecfem2D
+Check the output files in ./OUTPUT_FILES/
\ No newline at end of file
Modified: seismo/2D/SPECFEM2D/trunk/UTILS/visualization/plot_wavefield.pl
===================================================================
--- seismo/2D/SPECFEM2D/trunk/UTILS/visualization/plot_wavefield.pl 2012-04-19 11:40:25 UTC (rev 19954)
+++ seismo/2D/SPECFEM2D/trunk/UTILS/visualization/plot_wavefield.pl 2012-04-20 17:32:10 UTC (rev 19955)
@@ -7,7 +7,8 @@
#
# This perl script writes a shell script that calls GMT to plot kernels or snapshots of the wavefield.
# Most of the input commands are associated with plotting.
-# The script could use quite a bit of cleaning up.
+# Anyone is welcome/encouraged to improve this script.
+# It has been tested for the Tromp2005 and Tape2007 examples.
#
# NOTES:
# 1. must have output_wavefield_snapshot = .true. in Par_file to generate wavefield snapshots
@@ -16,7 +17,7 @@
# WAVEFIELD EXAMPLES:
# plot_wavefield.pl 100/1600/400 200/800/200 0/4/0/4 1/0.1/1/0.1 -3/-3/-3 6/6/6 0.0/0.001 1/0/1/1/1 1.7/1/1 1/0/1/120 M2_UPPA PSV
# 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
-# plot_wavefield.pl 400/2800/400 400/2000/400 0/200/0/80 50/10/40/10 -4/-4/-4 1/1/1 -8.0/0.02 0/1/0/0/1 3.0/1/0 1/0/1/200 Tromp2005 SH_homo
+# 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
# 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
# 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
#
@@ -46,8 +47,8 @@
#$numf = @frames;
$numf = ($pend - $pfirst)/$pint + 1;
-# directory with data files
-$bdir = "/data/svn/seismo/2D/SPECFEM2D_work";
+# directory with data files (USER CHANGE THIS)
+$bdir = "/data2/SVN/seismo/2D/SPECFEM2D_20120420";
#$idir1 = "$bdir/OUTPUT_FILES"; # if running from the default directory
$idir1 = "$bdir/EXAMPLES/$tlab/OUTPUT_FILES"; # if running from an examples directory
if (not -e $idir1) {die("check if idir1 $idir1 exist or not\n");}
More information about the CIG-COMMITS
mailing list