#############################################
#
# . . BP Velocity Benchmark
#
nz=801
oz=0
dz=37.5
nx=8000
ox=-80000
dx=37.5 
Cartpar=a1n=$(nz) a1d=$(dz) a1o=$(oz) a2n=$(nx) a2o=$(ox) a2d=$(dx) 
VPAR = n1=1201 d1=10 o1=0 n2=2401 d2=10 o2=-12000 vc=2000 n3=1 d3=20 o3=0
bppar = flo=1.0 fhi=32.5
BPweipar= min5=5 n5=283 j5=1 squeeze=n
sidecut=50
nfreq=250
B = ./Bin/
dn = /dev/null
RESDIR=./Fig/
#################################################
#
# . . Velocity
Vtest.H:
	Vel n1=600 n2=1000 n3=1 o1=0. d3=1. o3=0. o2=-5000. o3=0. d1=10. d2=10. vc=2000. > $@

#################################################
#
# . . Source and Receiver Wavefields
#
add1=511
add2=512

bppar = flo=2.0 fhi=30 
Test_WPAR = min5=5 max5=35 j5=2
XBIN = ./Bin/
################
#
# . . Data
add3=399
add4=624
RWFtest2.H:
	Wavelet n1=4096 n2=1 n3=1 d1=0.004 d2=10 d3=10 d4=10 d5=10 \
		o1=-0.08 tdelay=2.08 o2=0. o3=0. \
		domain=time fund=25 wavelet=ricker2> ta.H
	Wavelet n1=4096 n2=1 n3=1 d1=0.004 d2=10 d3=10 d4=10 d5=10 \
		o1=-0.08 tdelay=2.58 o2=0. o3=0. \
		domain=time fund=25 wavelet=ricker2> tb.H
	Wavelet n1=4096 n2=1 n3=1 d1=0.004 d2=10 d3=10 d4=10 d5=10 \
		o1=-0.08 tdelay=3.08 o2=0. o3=0. \
		domain=time fund=25 wavelet=ricker2> tc.H
	Wavelet n1=4096 n2=1 n3=1 d1=0.004 d2=10 d3=10 d4=10 d5=10 \
		o1=-0.08 tdelay=3.58 o2=0. o3=0. \
		domain=time fund=25 wavelet=ricker2> td.H
	Math file1=ta.H file2=tb.H file3=tc.H file4=td.H exp="file1+file2+file3+file4" > t0.H
	< t0.H Pad  beg2=$(add3) end2=$(add4)  |\
		 Smooth rect2=1 rect3=1 rect1=0 tridiag=1  > t1.H
	echo "o2=-5100 d2=10" >> t1.H
	< t1.H Bandpass $(bppar) | Transf wei=y verb=y  > t2.H
	< t2.H Window3d $(Test_WPAR) | Window3d > $@

SWFtest2.H:
	Wavelet n1=4096 n2=1 n3=1 d1=0.004 d2=10 d3=10 d4=10 d5=10 \
		o1=-0.08 tdelay=0.08 o2=0. o3=0. \
		domain=time fund=25 wavelet=ricker2 > t0.H
	< t0.H Pad beg2=$(add4) end2=$(add3)  |\
		 Smooth rect2=1 rect3=1 rect1=0 tridiag=1> t1.H
	echo "o2=-5100 d2=10" >> t1.H
	< t1.H Bandpass $(bppar) | Transf wei=y verb=y > t2.H
	< t2.H Window3d $(Test_WPAR) | Window3d > $@

########################
#
# . . Migrations
#
# . . Elliptic Test - Correct Velocity
MIGPARS =  eps=0.001 verbose=1 source_norm=0 newimage=1 Aeps=1. nxtap=50 norm=1 nsx=2 nsz=1 rscale=1. nh=64

SHOTtest3.H: RWFtest2.H SWFtest2.H Vtest.H
	$(XBIN)/RWE2D.x rwf=RWFtest2.H swf=SWFtest2.H vel=Vtest.H image=t$@ rays=Rrays.H \
		$(MIGPARS) method=2 oe3=0. de3=0.002 ne3=512 minang=1 maxang=179
	< t$@ Window3d min2=-3000 max2=3000 max1=4000  > $@
	rm -f t$@

# . . Cartesian Test - Correct Velocity
SHOTtest4.H: RWFtest2.H SWFtest2.H Vtest.H
	$(XBIN)/RWE2D.x rwf=RWFtest2.H swf=SWFtest2.H vel=Vtest.H image=t$@ rays=Rrays.H\
		$(MIGPARS) method=1 oe3=0. de3=10 ne3=512 
	< t$@ Window3d min2=-3000 max2=3000  max1=4000 > $@
	rm -f t$@

# . . Elliptic Test - Wrong Velocity
ELLIP-95test.H: RWFtest2.H SWFtest2.H Vtest.H
	Math file1="Vtest.H" exp="0.95*file1" > VVV.H 
	$(XBIN)/RWE2D.x rwf=RWFtest2.H swf=SWFtest2.H vel=VVV.H image=t$@ rays=Rrays.H \
		$(MIGPARS) method=2 oe3=0. de3=0.002 ne3=512 minang=1 maxang=179
	< t$@ Window3d min2=-3000 max2=3000 max1=4000  > $@
	rm -f t$@ VVV.H

# . . Cartesian Test - Wrong Velocity
CART-95test.H: RWFtest2.H SWFtest2.H Vtest.H
	Math file1="Vtest.H" exp="0.95*file1" > VVV.H 
	$(XBIN)/RWE2D.x rwf=RWFtest2.H swf=SWFtest2.H vel=VVV.H image=t$@ rays=Rrays.H $(MIGPARS)\
		method=1 oe3=0. de3=10 ne3=512 
	< t$@ Window3d min2=-3000 max2=3000  max1=4000 > $@
	rm -f t$@ VVV.H

# . . Figures for paper
Figure5.ps: SHOTtest3.H SHOTtest4.H
	matlab < ELtest.m
Figure6.ps: ELLIP-95test.H CART-95test.H
	matlab < WELtest.m