Commit 44eb521b authored by thomas.forbriger's avatar thomas.forbriger
Browse files

[TASK][DOC] (issue10): add homogeneous full space example

parent b3355389
......@@ -159,6 +159,7 @@ Transformation of vertical component data:
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 \
......@@ -212,4 +213,84 @@ Transformation of radial component data:
n:"reflected wave" o:-2.0 '
Test-case homogeneous full space:
---------------------------------
1. rename subsurface model
/bin/cp -vpd fullspace.d.mod fullspace.dmod
2. Calculate synthetic seismograms for a vertical line force
refmetx -v -4 -c -o fs.zforce0.sparse -ty ascii -l fs.zforce0.sparse.main
3. Rename vertical component file
/bin/mv -fv fs.zforce0.sparse.TZ fs.zforce0.sparse.LZ.inv.ascii
4. Calculate synthetic seismograms for a vertical point force
refmetx -v -4 -c -o fs.zforce0.sparse -ty ascii fs.zforce0.sparse.main
5. Swap sign of signals of vertical component
printf "fac -1.\nend\n" | tidofi -v -o -cs -type ascii -Type ascii \
fs.zforce0.sparse.PZ.ascii fs.zforce0.sparse.TZ
Transformation of vertical component data:
6. Apply direct wave transformation to vertical component
and swap sign of result
lisousi -v -o -type ascii -Type ascii -limitlength \
-integ 0.5 -tapdel 0.00625 -pad 2 -tapslo 10. \
fs.zforce0.sparse.PZ.LS.ascii \
fs.zforce0.sparse.PZ.ascii
printf "fac -1.\nend\n" | tidofi -v -o -cs -type ascii -Type ascii \
fs.zforce0.sparse.PZ.LS.inv.ascii \
fs.zforce0.sparse.PZ.LS.ascii
7. Apply reflected wave transformation to vertical component
and swap sign of result
lisousi -v -o -type ascii -Type ascii -sqrttaper -velocity 1. \
-limitlength -integ 0.5 -tapdel 0.00625 -pad 2 -tapslo 10. \
fs.zforce0.sparse.PZ.LSP.ascii \
fs.zforce0.sparse.PZ.ascii
printf "fac -1.\nend\n" | tidofi -v -o -cs -type ascii -Type ascii \
fs.zforce0.sparse.PZ.LSP.inv.ascii \
fs.zforce0.sparse.PZ.LSP.ascii
8. Apply single velocity transformation to vertical component
and swap sign of result
lisousi -v -o -type ascii -Type ascii -fredomain -velocity 1. \
fs.zforce0.sparse.PZ.LSG.ascii \
fs.zforce0.sparse.PZ.ascii
printf "fac -1.\nend\n" | tidofi -v -o -cs -type ascii -Type ascii \
fs.zforce0.sparse.PZ.LSG.inv.ascii \
fs.zforce0.sparse.PZ.LSG.ascii
9. Plot vertical component
refractx -d fs.zforce0.sparse.Z.ps/vcps -C -Eg -Sm 2 -Se 0.5 \
-Sa 8.0 -Sc 100. -TL 1 -Ts 1.2 \
-Fp lzline.pck -Ep T \
-Sn -Sr 1. \
-St -0.03,0.045 -ty ascii -Lc -Sx -5.,115. \
-Ta "true amplitudes scaled by (offset/1m)\u0.5\d v\dred\u=1. km/s" \
-Tx "reduced time / s" -TR \
-Tt "homogeneous full space testcase"\
fs.zforce0.sparse.LZ.inv.ascii o:0.0 s:4,1,10 r:0.6,0.6,0.6 b:F \
n:"line source (reference)" o:0.0 \
fs.zforce0.sparse.PZ.LS.inv.ascii o:-2.0 s:2,1,4 r:0.,0.0,1.00 \
n:"direct wave" o:0.0 \
fs.zforce0.sparse.PZ.LSG.inv.ascii o:-4.0 s:5,1,2 r:0.,0.59,0.58 \
n:"single velocity" o:-1.0 \
fs.zforce0.sparse.PZ.LSP.inv.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