# TECMakefile - tests for Chapter 2

#  . . Velocity - use BPVel.HH

BPVEL = /data/bpgom/
CBP.H:
	< $(BPVEL)/bpgom_velo.H Transp plane=12 | Window3d j2=2 | Transp plane=12|\
		 Pad beg3=200 end3=200 extend=1 > $@

# . . make Receiver Wavefield
TECrwf.H:
	Wavelet nt=2048 n1=2048 nx=1 n2=1 n3=1 np=1 domain=time wavelet=ricker2 \
			phase=0. fund=15 tdelay=2.0 > ta.H	
	Wavelet nt=2048 n1=2048 nx=1 n2=1 n3=1 np=1 domain=time wavelet=ricker2 \
			phase=0. fund=15 tdelay=2.5 > tb.H
	Wavelet nt=2048 n1=2048 nx=1 n2=1 n3=1 np=1 domain=time wavelet=ricker2 \
			phase=0. fund=15 tdelay=3.0 > tc.H
	Wavelet nt=2048 n1=2048 nx=1 n2=1 n3=1 np=1 domain=time wavelet=ricker2 \
			phase=0. fund=15 tdelay=3.5 > td.H
	Wavelet nt=2048 n1=2048 nx=1 n2=1 n3=1 np=1 domain=time wavelet=ricker2 \
			phase=0. fund=15 tdelay=4.0 > te.H
	Add ta.H tb.H tc.H td.H te.H > t0.H
	Pad end2=299 extend=1 < t0.H > t1.H
	Pad beg2=50 end2=50 < t1.H > t$@
	< t$@ Transf f_min=3 f_min1=6 f_max1=20 f_max=25 is2d=y shotmig=y > tt$@
	echo "n4=1 o4=0. d4=0.0005 o1=50000 d1=100 d2=100 n5=1 o5=75000 d5=250" >> tt$@
	< tt$@ Pad beg2=149 end2=150 > $@
	rm -f ta.H tb.H tc.H td.H te.H t0.H t1.H t$@ tt$@

RMIGPAR3D = norm=1 verb=1 verbose=1 newimage=1 eps=0.0001 nxtap=50 nytap=40 source_norm=0 timing=1\
                n1beg=0 n1end=0 overlap=5 nxint=1 planewave=1 n2beg=0 n2end=0 verb=1 verbose=1 nxint=16 
Cart.H:   TECrwf.H $(XBIN)/ConicalSRC3D.x CBP.H
	< TECrwf.H Window3d n3=160 > xxx.H
	echo "d2=200 o2=-30000" >> xxx.H
	time $(XBIN)/ConicalSRC3D.x rwf=xxx.H vel=CBP.H  $(RMIGPAR3D) image=$@ method=1 \
		ntilt=1 otilt=0 dtilt=0 wix=2 wiy=2 ne3=400 de3=100 oe3=0 nsx=2 nsz=2 nsy=2

CartRays.H:TECrwf.H $(XBIN)/ConicalSRC3Drays.H CBP.H
	<  TECrwf.H Window3d n3=160 > xxx.H
	echo "d2=200 o2=-30000" >> xxx.H
	time $(XBIN)/ConicalSRC3Drays.x rwf=xxx.H vel=CBP.H $(RMIGPAR3D) image=$@ method=1 \
		ntilt=1 otilt=0 dtilt=0 wix=2 wiy=2 ne3=400 de3=100 oe3=0 nsx=2 nsz=2 nsy=2

TEC.H:    TECrwf.H $(XBIN)/ConicalSRC3D.x CBP.H
	< TECrwf.H Window3d n3=160 > xxx.H
	time $(XBIN)/ConicalSRC3D.x rwf=xxx.H vel=CBP.H  $(RMIGPAR3D) image=$@ method=2 \
		ntilt=1 otilt=0 dtilt=0 wix=2 wiy=2 ne3=600 de3=0.0035 oe3=0 nsx=2 nsz=2 nsy=2

TECRays.H:TECrwf.H $(XBIN)/ConicalSRC3Drays.H CBP.H
	< TECrwf.H Window3d n3=160 > xxx.H
	time $(XBIN)/ConicalSRC3Drays.x rwf=xxx.H vel=CBP.H $(RMIGPAR3D) image=$@ method=2 \
		ntilt=1 otilt=0 dtilt=0 wix=2 wiy=2 ne3=600 de3=0.0035 oe3=0 nsx=2 nsz=2 nsy=2

PFRAME = nframe=1 frame1=180 frame2=300 frame3=250 point1=0.67 point2=0.73 popup=n \
		label1="Depth (m)" label2="Inline (m)" label3="Crossline (m)" title=" "\
		 $(dn) pclip=99 gainpanel=a labelsz=9 titlesz=11
vvv1.H:   Cart.H CBP.H
	< Cart.H Scale scale_to=1.  > t1.H
	Math file1="t1.H" file2="CBP.H" exp="40*file1+(file2-10000)/5000" |\
			Byte gainpanel=a pclip=98> $@
vvv2.H:   TEC.H CBP.H
	< TEC.H Scale scale_to=1.  > t2.H
	Math file1="t2.H" file2="CBP.H" exp="40*file1+(file2-10000)/5000" |\
			Byte gainpanel=a pclip=98> $@


NewPoint.v: vvv1.H vvv2.H
	< vvv1.H Window3d min2=50000 max2=90000 | Cubeplot $(PFRAME) out=v1.v 
	< vvv2.H Window3d min2=50000 max2=90000 | Cubeplot $(PFRAME) out=v2.v 
	vp_OverUnderAniso v2.v v1.v > $@