#############################################
#
# . . BP Velocity Benchmark
#
bppar = flo=1.0 fhi=35.
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=35.
BPweipar= min5=5 max5=28 j5=1 squeeze=n
sidecut=50
nfreq=250
B = /net/koko/jeff/Bin/LINUX/
#################################################
#
# . . Velocity

BP.Vel.H:#BPAIT-vel.HH is released with data
	< BPAIT-vel.HH Window3d j2=1 j1=2 > t1.H
	< t1.H Pad beg2=3000 extend=1 end2=2200 end1=1000 > $@	
	\rm -f t1.H

# . . model test point 1
SL1.H:
	< /net/koko/jeff/BPAIT-vel.HH Window3d min2=20000 max2=30000 max1=10000 j1=4 j2=2 > t1.H
	Fdmod intag=t1.H tmax=1 xs=25000 zs=25 nt=1000 abs=0,1,1,1 outtag=$@

#################################################
#
# . . Source and Receiver Wavefields
#
BP.rwf.H:#shots.small.H is included with the data release
	< /data/bpait/shots.small.H Window3d | Pad end1=47 > t1.H
	< t1.H Interp n2out=2048 o2out=-16994.30 d2out=9.287 type=2 maxsize=400 > tt1.H  #-16003.13 16.618
	< tt1.H		Window3d min2=-15000 max2=0 | Pad beg2=215 end2=218 > t2.H 
	< t2.H Mute vmute=2000 tmute=.5 > t3.H
	< t3.H Transf wei=y verb=y > t4.H
	< t4.H Window3d $(BPweipar) | Transp plane=25  |\
		 Window3d | Window3d n3=1 f3=0 > $@
	\rm -f t1.H tt1.H t2.H t3.H t4.H
BP.Sdata.H:
	Spike n1=2048 n2=1 n3=1 d1=0.006 d2=9.287 d3=1 o1=-0.024 o2=0. o3=0. \
		nsp=1 k1=5 k3=1 | Bandpass fhi=45> t0.H
	< t0.H Pad beg2=1830 end2=217  > t1.H
	< t1.H Bandpass $(bppar) | Transf wei=y verb=y > t2.H
	< t2.H Window3d $(BPweipar) | Transp plane=25 | Window3d > $@
	\rm -f t0.H t1.H t2.H

#################################################
#
# . . Migration
#
BP.Mig.H: BP.Sdata.H BP.Vel.H $(B)/RWE2D_GOGI.x
	< BP.rwf.H Window3d n3=320 f3=400 j2=1 |\
		 Window3d f1=50 max1=-25 |\
		 Pad beg1=50 end1=146 > rrwf.H
	< BP.Sdata.H  Window3d j2=1  > sswf.H
	< BP.Vel.H Window3d max1=24400 >  BPV.H
	Parallel command=$(B)/RWE2D_GOGI.x \
		files="rrwf.H,sswf.H,BPV.H,$@" \
		tags="rwf,swf,vel,image" \
		usage="INPUT,INPUT,INPUT,OUTPUT"\
		axis="3,3,3,3" verb=0 nblock=68 \
		file_type="DISTRIBUTE,COPY,COPY,COPY" \
		mach_mfile=mfile verbose=1  \
		ntap=50 norm=1 nsx=2 nsz=2 \
		ot=0. dt=0.002 nt=875 minang=5 maxang=175
#################################################
#
# . . Image
#
xxxcut1=10000
xxxcut2=40000
zzzcut2=12000
aspect =  crowd1=0.85 crowd2=0.8 

Figure6.v: BP.Mig.H
	<  BP.Mig.H Window3d min2=$(xxxcut1) max2=$(xxxcut2) max1=$(zzzcut2) min1=0   |\
		Tpow tpow=1.| Grey gainpanel=e pclip=99 out=t1.v $(dn) min2=$(xxxcut1)\
		max2=$(xxxcut2) max1=$(zzzcut2)	min1=0 label2="Distance (m)" 
		label1="Depth (m)" title=" " $(aspect)
	vp_OverUnderAniso t1.v $(R)/CARTMIG.v > $@

clean: jclean