From b108cc3088c338b1e8ec9fc97d57650a733eaf57 Mon Sep 17 00:00:00 2001
From: Jbrackenhoff <joeri_scw@live.nl>
Date: Thu, 6 Jun 2019 12:32:26 -0500
Subject: [PATCH] WS15 Joeri

---
 marchenko/demo/WS15/README.5          |  2 +-
 marchenko/demo/twoD/homgpng.scr       |  2 ++
 marchenko/demo/twoD/homgview.scr      |  1 +
 marchenko/demo/twoD/marchenko_ray.scr | 13 ++++++++++++-
 4 files changed, 16 insertions(+), 2 deletions(-)

diff --git a/marchenko/demo/WS15/README.5 b/marchenko/demo/WS15/README.5
index cec5287..ff43109 100644
--- a/marchenko/demo/WS15/README.5
+++ b/marchenko/demo/WS15/README.5
@@ -4,7 +4,7 @@ rayvsp.scr
 
 After the raytimes have been determined you can create the source function of the homogeneous Green's function by running:
 
-marchenko_ray.scr
+marchenko_ray.pbs
 
 A reference Green's function can be modeled by running
 
diff --git a/marchenko/demo/twoD/homgpng.scr b/marchenko/demo/twoD/homgpng.scr
index dfd55d0..05ba0f5 100755
--- a/marchenko/demo/twoD/homgpng.scr
+++ b/marchenko/demo/twoD/homgpng.scr
@@ -2,5 +2,7 @@
 
 suwind < homg_rayz.su key=fldr min=513 | supsimage title='Retrieved' perc=99 > homg_retrieved.eps
 /usr/bin/convert -quality 90 -antialias homg_retrieved.eps homg_retrieved.png
+suwind < homg_class_rayz.su key=fldr min=513 | supsimage title='Retrieved classical' perc=99 > homg_class_retrieved.eps
+/usr/bin/convert -quality 90 -antialias homg_class_retrieved.eps homg_class_retrieved.png
 supsimage < homg_ref_sp.su title='Reference' perc=99 > homg_reference.eps
 /usr/bin/convert -quality 90 -antialias homg_reference.eps homg_reference.png
diff --git a/marchenko/demo/twoD/homgview.scr b/marchenko/demo/twoD/homgview.scr
index d5efb59..c5bc8d9 100755
--- a/marchenko/demo/twoD/homgview.scr
+++ b/marchenko/demo/twoD/homgview.scr
@@ -1,4 +1,5 @@
 #!/bin/bash -x
 
 suwind < homg_rayz.su key=fldr min=513 | suximage title='Retrieved' perc=99 &
+suwind < homg_class_rayz.su key=fldr min=513 | suximage title='Retrieved Classical' perc=99 &
 suximage < homg_ref_sp.su title='Reference' perc=99 &
diff --git a/marchenko/demo/twoD/marchenko_ray.scr b/marchenko/demo/twoD/marchenko_ray.scr
index ad12049..5df653c 100755
--- a/marchenko/demo/twoD/marchenko_ray.scr
+++ b/marchenko/demo/twoD/marchenko_ray.scr
@@ -13,10 +13,21 @@ marchenko3D file_shot=shots/refl_rp.su verbose=1 \
     file_ray=raysrc_time.su file_amp=raysrc_amp.su file_wav=wavesrc.su \
     file_green=greensrc.su
 
+#Create the homogeneous Green's function
 marchenko3D file_shot=shots/refl_rp.su verbose=1 file_inp=greensrc.su \
     tap=0 niter=6 hw=8 shift=7 smooth=3 geomspread=1 ampest=1 \
     file_ray=rayz_time.su file_amp=rayz_amp.su file_wav=wavefp.su \
-    file_homg=homg_rayz.su file_green=green_rayz.su
+    file_homg=homg_rayz.su scheme=0 file_green=green_rayz.su
 
+#Create the homogeneous Green's function
+marchenko3D file_shot=shots/refl_rp.su verbose=1 file_inp=greensrc.su \
+    tap=0 niter=6 hw=8 shift=7 smooth=3 geomspread=1 ampest=1 \
+    file_ray=rayz_time.su file_amp=rayz_amp.su file_wav=wavefp.su \
+    file_homg=homg_class_rayz.su scheme=1 file_green=green_class_rayz.su
+
+#Change info in the header for display in SU
 sushw < homg_rayz.su key=ntr a=0 > tmp.su
 mv tmp.su homg_rayz.su
+
+sushw < homg_class_rayz.su key=ntr a=0 > tmp.su
+mv tmp.su homg_class_rayz.su
-- 
GitLab