WS15 Joeri

Makefile

@@ -3,8 +3,10 @@
 all: mkdirs 
 	cd FFTlib		; $(MAKE)
 	cd fdelmodc		; $(MAKE) install
+	cd fdelmodc3D		; $(MAKE) install
 	cd utils		; $(MAKE) install
 	cd marchenko		; $(MAKE) install
+	cd marchenko3D		; $(MAKE) install
 	cd corrvir		; $(MAKE) install
 	cd raytime		; $(MAKE) install
 	cd MDD			; $(MAKE) install
@@ -17,8 +19,10 @@ mkdirs:
 clean:
 	cd FFTlib 		; $(MAKE) $@
 	cd fdelmodc		; $(MAKE) $@
+	cd fdelmodc3D		; $(MAKE) $@
 	cd utils		; $(MAKE) $@
 	cd marchenko		; $(MAKE) $@
+	cd marchenko3D		; $(MAKE) $@
 	cd corrvir		; $(MAKE) $@
 	cd raytime		; $(MAKE) $@
 	cd MDD			; $(MAKE) $@
@@ -26,8 +30,10 @@ clean:
 realclean:
 	cd FFTlib 		; $(MAKE) $@
 	cd fdelmodc		; $(MAKE) $@
+	cd fdelmodc3D		; $(MAKE) $@
 	cd utils		; $(MAKE) $@
 	cd marchenko		; $(MAKE) $@
+	cd marchenko3D		; $(MAKE) $@
 	cd corrvir		; $(MAKE) $@
 	cd raytime		; $(MAKE) $@
 	cd MDD			; $(MAKE) $@

marchenko/demo/twoD/marchenko_ray.scr

+#!/bin/bash -x
+
+export PATH=$HOME/src/OpenSource/bin:$PATH:
+
+makewave w=fw fmin=0 flef=5 frig=80 fmax=100 nt=1024 dt=0.004 t0=0.0 scale=0 scfft=1 file_out=wavesrc.su
+makewave w=g2 fp=25 nt=1024 dt=0.004 t0=0.0 scale=1 scfft=0 file_out=wavefp.su
+
+export OMP_NUM_THREADS=40
+
+#apply the Marchenko algorithm
+marchenko3D file_shot=shots/refl_rp.su verbose=1 \
+    tap=0 niter=6 hw=8 shift=7 smooth=3 geomspread=1 ampest=1 \
+    file_ray=raysrc_time.su file_amp=raysrc_amp.su file_wav=wavesrc.su \
+    file_green=greensrc.su
+
+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

marchenko/demo/twoD/rayvsp.scr

+#!/bin/bash
+
+export PATH=$HOME/src/OpenSource/bin:$PATH:
+
+Raytime \
+    verbose=1 useT2=1 geomspread=1 file_rcv=raysrc.su \
+    zsrc=1000 xrcv1=-3000 xrcv2=3000 dxrcv=10 file_cp=syncl_cp.su zrcv1=0 zrcv2=0 \
+    xsrc=0 nxshot=1 nzshot=1 dxshot=10 dzshot=10 nraystep=25 method=fd
+
+Raytime \
+    verbose=1 useT2=1 geomspread=1 file_rcv=rayz.su \
+    zsrc=100 xrcv1=-3000 xrcv2=3000 dxrcv=10 file_cp=syncl_cp.su zrcv1=0 zrcv2=0 \
+    xsrc=0 nxshot=1 nzshot=181 dxshot=10 dzshot=10 nraystep=25 method=fd