[cig-commits] r17867 - in seismo/2D/SPECFEM2D/trunk: . EXAMPLES/M2_UPPA EXAMPLES/Tape2007 EXAMPLES/Tromp2005 doc
carltape at geodynamics.org
carltape at geodynamics.org
Mon Feb 14 23:25:22 PST 2011
Author: carltape
Date: 2011-02-14 23:25:21 -0800 (Mon, 14 Feb 2011)
New Revision: 17867
Added:
seismo/2D/SPECFEM2D/trunk/doc/manual_SPECFEM2D.pdf
Modified:
seismo/2D/SPECFEM2D/trunk/EXAMPLES/M2_UPPA/README
seismo/2D/SPECFEM2D/trunk/EXAMPLES/M2_UPPA/SOURCE_M2_UPPA
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/README
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/README
seismo/2D/SPECFEM2D/trunk/checkgrid.F90
Log:
updated README files for examples; added print line in checkgrid.F90 for USER_T0 ; added symbolic link for manual
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/M2_UPPA/README
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/M2_UPPA/README 2011-02-15 01:44:43 UTC (rev 17866)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/M2_UPPA/README 2011-02-15 07:25:21 UTC (rev 17867)
@@ -4,9 +4,19 @@
This default example creates an internal mesh with smooth topographic interfaces.
+TO RUN:
-To run the examples, type:
+0. Read the user manual in SPECFEM2D/doc/manual_SPECFEM2D.pdf
- ./process.sh
+1. in SPECFEM2D root directory, configure, e.g.,
+ ./configure FC=gfortran
-in this current directory.
+2. compile:
+ make all
+
+3. cd EXAMPLES/M2_UPPA/
+
+4. execute script to run mesher and solver for the PSV case:
+ ./process.sh
+
+5. check out the output files in the local directory OUTPUT_FILES
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/M2_UPPA/SOURCE_M2_UPPA
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/M2_UPPA/SOURCE_M2_UPPA 2011-02-15 01:44:43 UTC (rev 17866)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/M2_UPPA/SOURCE_M2_UPPA 2011-02-15 07:25:21 UTC (rev 17867)
@@ -4,7 +4,7 @@
zs = 2500. # source location z in meters
source_type = 1 # elastic force or acoustic pressure = 1 or moment tensor = 2
time_function_type = 1 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
-f0 = 5.0 # dominant source frequency (Hz) if not Dirac or Heaviside
+f0 = 10.0 # dominant source frequency (Hz) if not Dirac or Heaviside
t0 = 0.0 # time shift when multi sources (if one source, must be zero)
angleforce = 0.0 # angle of the source (for a force only)
Mxx = 1. # Mxx component (for a moment tensor source only)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/README
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/README 2011-02-15 01:44:43 UTC (rev 17866)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/README 2011-02-15 07:25:21 UTC (rev 17867)
@@ -2,10 +2,8 @@
README
----------------------------------------------------------------------
-Carl Tape, 05-Jan-2010
+Example for Tape-Liu-Tromp (GJI 2007)
-Examples for Tape-Liu-Tromp (GJI 2007)
-
Par_file_Tape2007_onerec
1 source, 1 receiver
Source index: 1
@@ -16,13 +14,27 @@
Source index: 5
Receiver index: 1-132
-In Tromp2005 there is a time shift of 48 s before the source starts.
-This is not yet implemented.
-
---------------------------
-To run the examples, type:
+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
- ./process.sh
+TO RUN:
-in this current directory.
+0. Read the user manual in SPECFEM2D/doc/manual_SPECFEM2D.pdf
+
+1. in SPECFEM2D root directory, configure, e.g.,
+ ./configure FC=gfortran
+
+2. modify constants.h file, setting USER_TO = 48 s
+
+3. compile:
+ make all
+
+4. cd EXAMPLES/Tape2007/
+
+5. 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
+
+---------------------------
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/README
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/README 2011-02-15 01:44:43 UTC (rev 17866)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/README 2011-02-15 07:25:21 UTC (rev 17867)
@@ -2,19 +2,48 @@
README
----------------------------------------------------------------------
-Examples for Tromp-Tape-Liu (GJI 2005)
+Example for Tromp-Tape-Liu (GJI 2005)
-The subscript s100 is the same set of simulations but scaled down in size
-by a factor of 100. This requires increasing the source frequency content
-and decreasing the time step.
+The example currently only generates the (synthetic) seismograms associated with the PSV example in Figure 7b and the SH example in Figure 2b.
-In Tromp2005 there is a time shift of 8 s before the source starts.
-This is not yet implemented.
+Demonstration of the kernel calculation is not yet added here.
+(The subscript s100 for the input files corresponds to the same set of simulations but scaled down in size by a factor of 100. This requires increasing the source frequency content and decreasing the time step; these changes are shown in the alternative files.)
+
---------------------------
-To run the examples, type:
+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
- ./process.sh
+TO RUN:
-in this current directory.
+0. Read the user manual in SPECFEM2D/doc/manual_SPECFEM2D.pdf
+
+1. in SPECFEM2D root directory, configure, e.g.,
+ ./configure FC=gfortran
+
+2. modify constants.h file, setting USER_TO = 8.0 s
+
+3. compile:
+ make all
+
+4. cd EXAMPLES/Tromp2005/
+
+5. 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
+
+ The seismograms (S0001.AA.BHX.semd,S0001.AA.BHZ.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.
+
+7. mv OUTPUT_FILES OUTPUT_FILES_PSV
+
+8. set p_sv = .false. in Par_file_Tromp2005
+
+9. 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
+
+ The seismogram (S0001.AA.BHY.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.
+
+---------------------------
Modified: seismo/2D/SPECFEM2D/trunk/checkgrid.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/checkgrid.F90 2011-02-15 01:44:43 UTC (rev 17866)
+++ seismo/2D/SPECFEM2D/trunk/checkgrid.F90 2011-02-15 07:25:21 UTC (rev 17867)
@@ -1609,6 +1609,7 @@
if(.not. initialfield) then
f0max = -HUGEVAL
f0min = HUGEVAL
+ write(IOUT,*) ' USER_T0 = ',USER_T0
do i = 1,NSOURCE
if(time_function_type(i) /= 4 .and. time_function_type(i) /= 5) then
Added: seismo/2D/SPECFEM2D/trunk/doc/manual_SPECFEM2D.pdf
===================================================================
--- seismo/2D/SPECFEM2D/trunk/doc/manual_SPECFEM2D.pdf (rev 0)
+++ seismo/2D/SPECFEM2D/trunk/doc/manual_SPECFEM2D.pdf 2011-02-15 07:25:21 UTC (rev 17867)
@@ -0,0 +1 @@
+link USER_MANUAL/manual_SPECFEM2D.pdf
\ No newline at end of file
Property changes on: seismo/2D/SPECFEM2D/trunk/doc/manual_SPECFEM2D.pdf
___________________________________________________________________
Name: svn:special
+ *
More information about the CIG-COMMITS
mailing list