[TASK] (lisousi): add radial component example

src/ts/lisousi/examples/README

@@ -136,6 +136,8 @@ Test-case Bietigheim:
        biet.zforce.sparse.PZ.ascii biet.zforce.sparse.TZ 
      /bin/mv -fv biet.zforce.sparse.TR biet.zforce.sparse.PR.ascii 
 
+Transformation of vertical component data:
+
 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. \
@@ -143,12 +145,12 @@ Test-case Bietigheim:
       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. \
+  lisousi -v -o -type ascii -Type ascii -sqrttaper -velocity 0.18 \ 
+      -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
+8. Apply single velocity transformation 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
@@ -168,4 +170,38 @@ Test-case Bietigheim:
       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 '
 
+Transformation of radial component data:
+
+10. Apply direct wave transformation to radial component
+  lisousi -v -o -type ascii -Type ascii -limitlength \
+     -integ 0.5 -tapdel 0.00625 -pad 2 -tapslo 10. \
+     biet.zforce.sparse.PR.LS.ascii \
+     biet.zforce.sparse.PR.ascii
+
+11. Apply reflected wave transformation to radial component
+  lisousi -v -o -type ascii -Type ascii -sqrttaper -velocity 0.18  \
+     -limitlength -integ 0.5 -tapdel 0.00625 -pad 2 -tapslo 10. \
+     biet.zforce.sparse.PR.LSP.ascii \
+     biet.zforce.sparse.PR.ascii
+
+12. Apply single velocity transformation to radial component
+  lisousi -v -o -type ascii -Type ascii -fredomain -velocity 0.18 \
+    biet.zforce.sparse.PR.LSG.ascii \
+    biet.zforce.sparse.PR.ascii
+
+13. Plot vertical component
+  refractx -d biet.zforce.sparse.R.ps/vcps \
+     -C -Eg -Sm 2 -Se 1.1 -Sa 1.0 -Sc 100. \
+     -St 0.,0.7 -ty ascii -Lc -Sx -5.,105. -Ts 1.0 -TL 1\
+     -Ta "true amplitudes scaled by (offset/1m)\u1.1\d" -TR \
+     -Tt "1D testcase Bietigheim; radial component"\
+       biet.zforce.sparse.LR.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.PR.LS.ascii o:-2.0 s:2,1,4 r:0.,0.0,1.00 \
+         n:"direct wave" o:0.0 \
+       biet.zforce.sparse.PR.LSG.ascii o:-4.0 s:5,1,2 r:0.,0.59,0.58 \
+         n:"single velocity"  o:-1.0 \
+       biet.zforce.sparse.PR.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 -----