Commit 63a24ee1 authored by thomas.forbriger's avatar thomas.forbriger
Browse files

[WP][DOC] (issue10): add step-by-step examples

End-user documentation for lisousi shall be improved on this branch.
Add a step-by-step example for the sparse receiver array example, testcase
Bietigheim.
parent 8d7894a6
......@@ -4,11 +4,13 @@ Toy examples and test cases for lisousi
---------------------------------------
This subdirectory contains test cases used by Forbriger et al. (2014). The
processing of the test cases is controlled by make (see Makefile). Synthetic
waveform are calculated with the reflectivity method (Fuchs and Müller, 1971;
Müller, 1985) as implemented by Ungerer (1990) and as provided by the program
refmet in Seitosh. Several refmet control files are provided in this
subdirectory. Their purpose can be displayed by
processing of the test cases is controlled by make (see Makefile).
Step-by-step examples are given at the end of this file.
Synthetic waveforms are calculated with the reflectivity method (Fuchs and
Müller, 1971; Müller, 1985) as implemented by Ungerer (1990) and as provided
by the program refmet in Seitosh. Several refmet control files are provided in
this subdirectory. Their purpose can be displayed by
make reportcontrol
......@@ -99,4 +101,71 @@ Ungerer, J., 1990. Berechnung von Nahfeldseismogrammen mit der
Reflektivitätsmethode, Diplomarbeit, Institut für Geophysik der Universität
Stuttgart.
=============================================================================
Step-by-step example
====================
The sparse receiver array test cases can be run manually by the following
steps:
Test-case Bietigheim:
---------------------
1. Create subsurface model in file format required by refmet
mocon P biet000601.mod biet000601.dmod
2. Calculate synthetic seismograms for a vertical line force
refmetx -v -4 -c -o biet.zforce.sparse -ty ascii -l \
biet.zforce.sparse.main
3. Swap sign of signals of vertical component, rename radial component file
printf "fac -1.\nend\n" | tidofi -v -o -cs -type ascii -Type ascii \
biet.zforce.sparse.LZ.ascii biet.zforce.sparse.TZ
/bin/mv -fv biet.zforce.sparse.TR biet.zforce.sparse.LR.ascii
4. Calculate synthetic seismograms for a vertical point force
refmetx -v -4 -c -o biet.zforce.sparse -ty ascii biet.zforce.sparse.main
5. Swap sign of signals of vertical component, rename radial component file
printf "fac -1.\nend\n" | tidofi -v -o -cs -type ascii -Type ascii \
biet.zforce.sparse.PZ.ascii biet.zforce.sparse.TZ
/bin/mv -fv biet.zforce.sparse.TR biet.zforce.sparse.PR.ascii
6. Apply direct wave transformation to vertical component
lisousi -v -o -type ascii -Type ascii -limitlength \
-integ 0.5 -tapdel 0.00625 -pad 2 -tapslo 10. \
biet.zforce.sparse.PZ.LS.ascii
biet.zforce.sparse.PZ.ascii
7. Apply reflected wave transformation to vertical component
lisousi -v -o -type ascii -Type ascii -sqrttaper -limitlength \
-integ 0.5 -tapdel 0.00625 -pad 2 -tapslo 10. \
biet.zforce.sparse.PZ.LSP.ascii \
biet.zforce.sparse.PZ.ascii
8. Apply single velocity to vertical component
lisousi -v -o -type ascii -Type ascii -fredomain -velocity 0.18 \
biet.zforce.sparse.PZ.LSG.ascii \
biet.zforce.sparse.PZ.ascii
9. Plot vertical component
refractx -d biet.zforce.sparse.Z.ps/vcps \
-C -Eg -Sm 2 -Se 1.1 -Sa 1. -Sc 100. \
-Sn -St -0.0,0.7 -ty ascii -Lc -Sx -5.,105. -Ts 1.0 -TL 1 \
-Si -Ta "true amplitudes scaled by (offset/1m)\u1.1\d" -TR \
-Tt "1D testcase Bietigheim; vertical component"\
biet.zforce.sparse.LZ.ascii o:0.0 s:4,1,10 r:0.6,0.6,0.6 b:F \
n:"line source (reference)" o:0.0 \
biet.zforce.sparse.PZ.LS.ascii o:-2.0 s:2,1,4 r:0.,0.0,1.00 \
n:"direct wave" o:0.0 \
biet.zforce.sparse.PZ.LSG.ascii o:-4.0 s:5,1,2 r:0.,0.59,0.58 \
n:"single velocity" o:-1.0 \
biet.zforce.sparse.PZ.LSP.ascii o:-6.0 s:6,1,2 r:0.82,0.55,0.01 \
n:"reflected wave" o:-2.0 '
----- END OF README -----
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment