working on Marchenko planewave demo

marchenko/demo/planewave/README

@@ -10,16 +10,19 @@ Description of files:
 
 3) ./direct.scr compute 1 shot that contains direct wave only
    
-4) ./remove_direct.scr removes the direct wave from the data created in 2 and places all shots in one file shots/refl_rp.su 
+4) ./remove_direct.scr removes the direct wave from the data created in 3 and places all shots in one file shots/refl_rp.su 
 
    Figures of the model and the middle shot are generated by
    ./epsModel.scr model
 
-5) IniFocus.scr compute plane wave responses at 600 m depth with angles -5 0 and 5 degrees.
+5) ./IniFocus.scr compute plane wave responses at 600 m depth with angles -5 0 and 5 degrees.
+
    To compute the plane wave respone with the sources all starting at t=0 on a slanted line in the model run
-   ./IniFocus.scr
+
    ./Ini2Focus.scr
-   To generate the 12 snapshot figures
+
+   To generate 12 snapshot figures based on the results of Ini(2)Focus.scr
+
    ./epsModel.scr snapshots
 
 *** Steps to generate the Marchenko plane wave results.

marchenko/demo/planewave/marchenkoPlaneIterMuteWindow.scr

-#!/bin/bash -x
-
-ROOT=/vardim/home/thorbcke/src/
-export OMP_NUM_THREADS=4
-exe=$ROOT/OpenSource/marchenko/marchenko
-
-src_angle=-5
-inifile=SourceNew${src_angle}a120_rp.su
-base=${inifile%_rp.su}
-base=MuteW
-
-#change source code marchenko.c to generate mute windows. /* for testing time-windows with dipping plane waves */
-fmute file_shot=$inifile file_out=fdplus${base}.su above=-1 shift=-6 verbose=1 check=0 hw=2
-#for testing mute window in manual  (B/W figures)
-$exe file_shot=shots/refl_rp.su file_tinv=fdplus${base}.su nshots=601 verbose=2 \
-        tap=3 niter=2 hw=2 shift=7 smooth=3 plane_wave=1 src_angle=$src_angle rotate=1 \
-        file_green=green${base}.su file_gplus=Gplus${base}.su file_gmin=Gmin${base}.su  \
-        file_f1plus=f1plus${base}.su file_f1min=f1min${base}.su file_f2=f2${base}.su file_iter=iter${base}.su 
-
-
-
-