####################################
#
# Tests using Sigsbee data
#
xxx.H:
	Window3d squeeze=n n4=1 f4=11 <../RTM/sig_shot_prep_mad.H > $@

yyy.H:
	Window3d squeeze=n n3=10 <../RTM/sig_sou_mig_mad.H > $@

mig_tst.H:
	./bin/OWmig2d1.x < ../RTM/sig_shot_prep_mad.H sou=../RTM/sig_sou_mig_mad.H vel=../OneWay/sig_velmig_mad_transp.H nrefvel=4 xmin=-56250 xmax=105705 \
	ymin=0 ymax=0 zmin=0 zmax=16000 nh=1 oh=0 dh=0 >$@

#mig_tst_par.H:
#	./bin/OWmig2d_par.x < ../RTM/sig_shot_prep_mad.H sou=../RTM/sig_sou_mig_mad.H vel=../OneWay/sig_velmig_mad_transp.H nrefvel=4 xmin=-56250 #xmax=105705 \
#	ymin=0 ymax=0 zmin=0 zmax=16000 nh=1 oh=0 dh=0 >$@

mig_tst_par.H:
	./bin/OWmig2d_par.x < xxx.H sou=../RTM/sig_sou_mig_mad.H vel=../OneWay/sig_velmig_mad_transp.H nrefvel=4 xmin=-56250 xmax=105705 \
	ymin=0 ymax=0 zmin=0 zmax=16000 nh=1 oh=0 dh=0 >$@

ill_tst_par.H:
	./bin/OWmig2d_par_DD.x < ../RTM/sig_shot_prep_mad.H sou=../RTM/sig_sou_mig_mad.H vel=../OneWay/sig_velmig_mad_transp.H nrefvel=4 xmin=-56250 xmax=105705 \
	ymin=0 ymax=0 zmin=0 zmax=16000 nh=1 oh=0 dh=0 >$@

#mig_tst.H:
#	./bin/OWmig2d1.x < xxx.H sou=../RTM/sig_sou_mig_mad.H vel=../OneWay/sig_velmig_mad_transp.H nrefvel=4 xmin=-56250 xmax=105705 \
#	ymin=0 ymax=0 zmin=0 zmax=16000 nh=1 oh=0 dh=0 >$@

sig_prep.H:
	< ./data/nfs-splitspread-shots.H Window3d j2=2 j3=2 max3=70000 max1=8 | /net/koko/bob/SEP_BOB/bin/Transf shotmig=y is2d=y wei=y f_min=3 f_min1=5 \
	f_max1=32 f_max=40 maxsize=1000 > j2.H
	< j2.H Transp reshape=1,3,5 plane=23 | Pad beg1=200 end1=200 > ./data/$@

sig_prep1.H:
	< ./data/nfs-splitspread-shots.H Window3d n3=1 j2=2 f3=196 max1=8 squeeze=n | /net/koko/bob/SEP_BOB/bin/Transf shotmig=y is2d=y wei=y f_min=3 f_min1=5 \
	f_max1=32 f_max=40 maxsize=1000 > j2.H
	< j2.H Transp reshape=1,3,5 plane=23 | Pad beg1=200 end1=200 > ./data/$@

sig_source.H:
	Spike n1=1001 n2=1 n3=1 n4=1 n5=1 o1=0 o2=0 o3=0 o4=0 o5=0 d1=0.008 d2=0 d3=0 d4=0 d5=0 nsp=1 k1=501 mag=1 | Bandpass flo=5 fhi=32 | \
	/net/koko/bob/SEP_BOB/bin/Transf shotmig=y is2d=y wei=y f_min=3 f_min1=5 f_max1=32 f_max=40 maxsize=1000 > j2.H
	< j2.H Transp reshape=1,3,5 plane=23 > ./data/$@; echo o3=2.875 d3=0.125 >> ./data/$@

#sig_tst.H:
#	< ./data/sig_prep.H ./bin/OWmig2d_par.x sou=./data/sig_source.H vel=../OneWay/sig_velmig_mad_transp.H nrefvel=4 xmin=-48100 \
#	xmax=126500 ymin=0 ymax=0 zmin=0 zmax=17000 nh=1 oh=0 dh=0 > ./data/$@

sig_tst.H:
	< ./data/sig_prep1.H ./bin/OWmig2d_par.x sou=./data/sig_source.H vel=../OneWay/sig_velmig_mad_transp.H \
	nrefvel=4 xmin=-16000 xmax=96500 ymin=0 ymax=0 zmin=0 zmax=17000 nh=1 oh=0 dh=0 > ./data/$@
#	< ./data/sig_prep1.H Window3d squeeze=n n4=1 f4=196 | ./bin/OWmig2d_par.x sou=./data/sig_source.H vel=../OneWay/sig_velmig_mad_transp.H \
#	nrefvel=4 xmin=48100 xmax=126500 ymin=0 ymax=0 zmin=0 zmax=17000 nh=1 oh=0 dh=0 > ./data/$@

###################################
#
# After testing, migrate a good 
# quantity of data
#

01sig_temp22.H:
	< temp22.H Window3d max1=6 | Mute vmute=4720 tmute=1.4 | /net/koko/bob/SEP_BOB/bin/Transf shotmig=y is2d=y wei=y f_min=3 f_min1=5 f_max1=40 f_max=50 maxsize=1000 > j2.H
	< j2.H Transp reshape=1,3,5 plane=23 | Pad beg1=200 end1=200 > $@

02sig_sou_temp.H:
	Window3d n3=1 n2=1 min2=5000 max1=2 < temp22.H | Pad n1out=710 > j1.H; echo o1=-1.016 >> j1.H
	<j1.H /net/koko/bob/SEP_BOB/bin/Transf shotmig=y is2d=y wei=y f_min=3 f_min1=5 f_max1=40 f_max=50 maxsize=1000 > j2.H
	< j2.H Transp reshape=1,3,5 plane=23 > $@

00sig_velmig.H:
	Pad < ~tang/Research.ln/Sigsb2a/sigsb2a.wei.vel.H beg1=1500 end1=1500 extend=1 > $@

00sig_velmig_res.H:
	Window3d j1=2 j3=2 < 00sig_velmig.H | Transp > $@

z_temp.H:
	< ../RTM/z.H Window3d max1=6 | /net/koko/bob/SEP_BOB/bin/Transf shotmig=y is2d=y wei=y f_min=3 f_min1=5 f_max1=40 f_max=50 maxsize=1000 > j2.H
	< j2.H Transp reshape=1,3,5 plane=23 | Pad beg1=200 end1=200 > $@
 
03sig_mig_temp.H:
	< z_temp.H ./bin/OWmig2d_par.x sou=02sig_sou_temp.H vel=00sig_velmig_res.H \
	nrefvel=4 xmin=-31000 xmax=142000 ymin=0 ymax=0 zmin=0 zmax=17000 nh=1 oh=0 dh=0 wbott=0 dat= datum.H> $@

#03sig_mig_temp.H:
#	< 01sig_temp22.H Window3d n4=1 f4=74 squeeze=n | ./bin/OWmig2d_par.x sou=02sig_sou_temp.H vel=00sig_velmig_res.H \
#	nrefvel=4 xmin=-24000 xmax=90000 ymin=0 ymax=0 zmin=0 zmax=20000 nh=1 oh=0 dh=0 wbott=6200 dat=datum.H> $@

#03sig_mig_temp.H:
#	< 01sig_temp22.H ./bin/OWmig2d_par.x sou=02sig_sou_temp.H vel=00sig_velmig_res.H \
#	nrefvel=4 xmin=-46000 xmax=90000 ymin=0 ymax=0 zmin=0 zmax=20000 nh=1 oh=0 dh=0 > $@

01sig_prep.H:
	< ../OneWay/nfs-splitspread-shots.H Window3d max1=8 | /net/koko/bob/SEP_BOB/bin/Transf shotmig=y is2d=y wei=y \
	f_min=3 f_min1=5 f_max1=36 f_max=42 maxsize=1000 > j2.H
	< j2.H Transp reshape=1,3,5 plane=23 | Pad beg1=300 end1=300 > $@

02sig_sou1.H:
	Window3d n4=1 < 01sig_prep.H | Transp | Cabs | Stack | Smooth rect1=9 > specsmoo.H
	Math file1=specsmoo.H exp="file1*0" > phsesmoo.H; Cmplx specsmoo.H phsesmoo.H | Transp reshape=1,3,5 plane=23| \
	Transp reshape=1,3,5 plane=12 > $@ 
	< $@ Window3d | Transp| Transf |Graph|Tube&

00sig_velmig.H:
	<../2Wpar/1000velmig_my_areal.H Pad beg2=1000 end2=1000 > $@

03sig_mig.H:
	Window3d squeeze=n max4=60000 min4=30000 j4=20 n3= 2 < 01sig_prep.H > 01sig_prep_select.H
#	Window3d squeeze=n n3=2 < 02sig_sou1.H > 02sig_sou1_select.H
	< 01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1.H vel=00sig_velmig.H \
	nrefvel=4 xmin=-20000 xmax=110000 ymin=0 ymax=0 zmin=0 zmax=20000 nh=1 oh=0 dh=0 wbott=0 vsal=14000 > $@

tstsig_mig.H:
	Window3d squeeze=n n4=1 min4=30000 n3= 2 < 01sig_prep.H > 01sig_prep_select.H
#	Window3d squeeze=n n3=2 < 02sig_sou1.H > 02sig_sou1_select.H
	< 01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1.H _select.H vel=00sig_velmig.H \
	nrefvel=4 xmin=-20000 xmax=110000 ymin=0 ymax=0 zmin=0 zmax=20000 nh=1 oh=0 dh=0 wbott=0 vsal=14700 > $@

tstsig_mig1.H:
	Window3d squeeze=n min4=30000 n4=2 < 01sig_prep.H > xxx.H
#	Window3d squeeze=n n3=2 < 02sig_sou1.H > 02sig_sou1_select.H
	< xxx.H ./bin/OWmig2d_par.x sou=02sig_sou1.H vel=00sig_velmig.H \
	nrefvel=4 xmin=-20000 xmax=110000 ymin=0 ymax=0 zmin=0 zmax=20000 nh=1 oh=0 dh=0 wbott=0 > $@

tstsig_mig2:
	Window3d squeeze=n min4=12000 n4=2 < 01sig_prep.H > xxx.H
#	Window3d squeeze=n n3=2 < 02sig_sou1.H > 02sig_sou1_select.H
	< xxx.H ./bin/OWmig2d_par.x sou=02sig_sou1.H vel=00sig_velmig.H image=tstsig_mig2.H \
	nrefvel=4 xmin=-38000 xmax=70000 ymin=0 ymax=0 zmin=0 zmax=24000 nh=1 33 oh=0 -1200 dh=75 wbott=0 > /dev/null

# 	/net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX86_64/ShotMig \

tstsig_mig3.H:
	@date
	/net/koko/homes/sep/claudio/research.ln/Others/ShotMig/ShotMig \
	D=xxx.H D_d=02sig_sou1.H V=yyy.H R=$@ \
	pverb=2  nx_shots=2 ox_shots=29975 \
	par=go_sis.P taper=y taper_pcent=2 marine=y minx=-20000 maxx=110000 & #restart=y

04sig_mig.H:
#	Window3d squeeze=n max4=70000 min4=12000 < 01sig_prep.H > 01sig_prep_select.H
	< 01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1.H vel=00sig_velmig.H \
	nrefvel=4 xmin=-38000 xmax=120000 ymin=0 ymax=0 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 > $@

05sig_mig_grid2:
#	Window3d squeeze=n max4=70000 min4=68000 j4=2 j1=2 < 01sig_prep.H > zzz.H
#	Window3d squeeze=n max4=70000 min4=12000 j4=2 j1=2 < 01sig_prep.H > 01sig_prep_grid2.H
#	Window3d < 00sig_velmig_fix.H squeeze=n j2=2 > 00sig_velmig_grid2.H
	< 01sig_prep_grid2.H ./bin/OWmig2d_par.x sou=02sig_sou1.H vel=00sig_velmig_grid2.H image=05sig_mig_grid2.H  \
	velo=velo_check.H vsal=14500 nrefvel=4 xmin=-38000 xmax=120000 ymin=0 ymax=0 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 > /dev/null

05sig_mig_grid2_brad.H:
#	Window3d squeeze=n max4=70000 min4=12000 j4=2 j1=2 < 01sig_prep.H > 01sig_prep_grid2.H
	Window3d < 00sig_velmig_fix.H squeeze=n j2=2 |Transp plane=12 reshape=1,3> 00sig_velmig_grid2_brad.H
	@date
	/net/koko/homes/sep/claudio/research.ln/Others/ShotMig/ShotMig \
	D=01sig_prep_grid2.H D_d=02sig_sou1.H V=00sig_velmig_grid2_brad.H R=$@ \
	pverb=2 nx_shots=194 ox_shots=11975 \
	par=go_sis.P taper=y taper_pcent=2 marine=y minx=-38000 maxx=120000 & #restart=y