###################################
#
# 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:
#	< ../RTM/sigsb_shot.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 | Pad beg1=200 end1=200 > $@

	< ../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=200 end1=200 > $@

#	< ../RTM/temp11.H Window3d j2=2 max1=8 | Bandpass flo=5 fhi=50 | Window3d j1=10 | Mute vmute=4720 tmute=1.4 | \
#	/net/koko/bob/SEP_BOB/bin/Transf shotmig=y is2d=y wei=y f_min=5 f_min1=8 f_max1=35 f_max=42 maxsize=1000 > j2.H
#	< j2.H Transp reshape=1,3,5 plane=23 | Pad beg1=200 end1=200 > $@

#	< ../RTM/temp11.H Window3d j2=2 max1=8 | Bandpass flo=5 fhi=50 |Window3d j1=10| 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 > $@

#	< sigsb2a.rec.time.H Bandpass flo=3 fhi=50 | Window3d j3=3 n4=1 f4=240 j1=2 max3=8| Transp plane=12| Transp plane=34 | \
#	/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 > $@

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&

02sig_source.H:
	<../../Vel-ill/nfs-splitspread-shots.H Window3d n2=1 min2=0 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 | \
	Window3d | Cabs |  Transp | Stack | Smooth rect1=5 > 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 > $@ 

#	Spike n1=711 n2=1 n3=1 n4=1 n5=1 o1=-3.0 o2=0 o3=0 o4=0 o5=0 d1=0.00845 d2=0 d3=0 d4=0 d5=0 nsp=1 k1=356 mag=1 | Bandpass flo=5 fhi=50 | \
#	/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 > $@ 

#############################################################
#
# Generating synthesis operator for Areal shots
#
#2000areal_synth.H:
#	< 01sig_prep.H Window3d squeeze=n n4=1 min4=40000 min1=-7500 max1=7500 | ./bin/OW_arealsht.x sou=02sig_sou1.H vel=1000velmig_my.H nrefvel=4 #xmin=32525 xmax=47525 ymin=0 ymax=0 zmin=0 zmax=20000 nh=1 oh=0 dh=0 wbott=0 vsal=14700 traces_in_taper=100 > $@

2000areal_synth.H:
	< 01sig_prep.H Window3d squeeze=n n4=401 | ./bin/OW_arealsht.x sou=02sig_sou1.H vel=1000velmig_my.H nrefvel=4 xmin=-30100 xmax=111950 ymin=0 ymax=0 zmin=0 zmax=20000 wbott=0 vsal=14700 traces_in_taper=0 synsh=2001Synth_Areal.H> $@

synth:
#	< 01sig_prep.H Window3d squeeze=n n4=10 f4=301 | ./bin/OW_arealsht.x sou=02sig_sou1.H vel=1000velmig_my.H nrefvel=4 xmin=15050 xmax=98600 #ymin=0 ymax=0 zmin=0 zmax=20000 nh=1 oh=0 dh=0 wbott=0 vsal=14700 traces_in_taper=600 synth_op=2000areal_synth.H synsh=2001Synth_Areal.H >/dev/null
#	< 01sig_prep.H Window3d squeeze=n n4=10 f4=301 | ./bin/OW_arealsht.x sou=02sig_sou1.H vel=1000velmig_my.H nrefvel=4 xmin=15050 xmax=98600 #ymin=0 ymax=0 zmin=0 zmax=20000 nh=1 oh=0 dh=0 wbott=0 vsal=14700 traces_in_taper=600 synth_op=2000areal_synth.H synsh=2001Synth_Areal.H >/dev/null
	< 01sig_prep.H Window3d squeeze=n n4=401 | ./bin/OW_arealsht.x sou=02sig_sou1.H vel=1000velmig_my.H nrefvel=4 xmin=-30100 xmax=111950 ymin=0 ymax=0 zmin=0 zmax=20000 nh=1 oh=0 dh=0 wbott=0 vsal=14700 traces_in_taper=600 synth_op=2000areal_synth.H synsh=2001Synth_Areal.H > /dev/null

#2000areal_synth.H:
#	<2000areal_synth ./bin/OWmig2d_par.x sou=02sig_sou1.H vel=1000velmig_my.H nrefvel=4 xmin=-2000 xmax=82000 ymin=0 ymax=0 zmin=0 zmax=20000 nh=1 #oh=0 dh=0 wbott=0 vsal=14700 wr= WS.H traces_in_taper=30> $@

2002Areal_mig.H:
	< 2001Synth_Areal_transp.H ./bin/OWmig2d_areal.x sou=2000areal_synth_transp.H vel=1000velmig_my_areal.H nrefvel=4 xmin=-30100 xmax=111950 ymin=0 ymax=0 zmin=0 zmax=20000 nh=33 oh=-1200 dh=75 wbott=0 vsal=14700 traces_in_taper=0 > $@
#####################################
#
# Sigsbee
#
5004areal_synth_d%_a.H:
	< 01sig_prep.H ./bin/OW_arealsht.x sou=02sig_sou1.H vel=1000velmig_my.H nrefvel=4 xmin=-30100 xmax=111950 ymin=0 ymax=0 zmin=0 zmax=20000 \
	wbott=0  depth_wav=$* vsal=14700 traces_in_taper=0 synsh=5014areal_synth_d$*_a.H shot_in=01sig_prep.H nphi=10 dphi=2. ophi=-70 encode=0 > $@

5004areal_synth_d%_b.H:
	< 01sig_prep.H ./bin/OW_arealsht.x sou=02sig_sou1.H vel=1000velmig_my.H nrefvel=4 xmin=-30100 xmax=111950 ymin=0 ymax=0 zmin=0 zmax=20000 \
	wbott=0  depth_wav=$* vsal=14700 traces_in_taper=0 synsh=5014areal_synth_d$*_b.H shot_in=01sig_prep.H nphi=10 dphi=2. ophi=-50 encode=0 > $@

5004areal_synth_d%_c.H:
	< 01sig_prep.H ./bin/OW_arealsht.x sou=02sig_sou1.H vel=1000velmig_my.H nrefvel=4 xmin=-30100 xmax=111950 ymin=0 ymax=0 zmin=0 zmax=20000 \
	wbott=0  depth_wav=$* vsal=14700 traces_in_taper=0 synsh=5014areal_synth_d$*_c.H shot_in=01sig_prep.H nphi=10 dphi=2. ophi=-30 encode=0 > $@

5004areal_synth_d%_d.H:
	< 01sig_prep.H ./bin/OW_arealsht.x sou=02sig_sou1.H vel=1000velmig_my.H nrefvel=4 xmin=-30100 xmax=111950 ymin=0 ymax=0 zmin=0 zmax=20000 \
	wbott=0  depth_wav=$* vsal=14700 traces_in_taper=0 synsh=5014areal_synth_d$*_d.H shot_in=01sig_prep.H nphi=10 dphi=2. ophi=-10 encode=0 > $@

5004areal_synth_d%_e.H:
	< 01sig_prep.H ./bin/OW_arealsht.x sou=02sig_sou1.H vel=1000velmig_my.H nrefvel=4 xmin=-30100 xmax=111950 ymin=0 ymax=0 zmin=0 zmax=20000 \
	wbott=0  depth_wav=$* vsal=14700 traces_in_taper=0 synsh=5014areal_synth_d$*_e.H shot_in=01sig_prep.H nphi=10 dphi=2. ophi=10 encode=0 > $@

5004areal_synth_d%_f.H:
	< 01sig_prep.H ./bin/OW_arealsht.x sou=02sig_sou1.H vel=1000velmig_my.H nrefvel=4 xmin=-30100 xmax=111950 ymin=0 ymax=0 zmin=0 zmax=20000 \
	wbott=0  depth_wav=$* vsal=14700 traces_in_taper=0 synsh=5014areal_synth_d$*_f.H shot_in=01sig_prep.H nphi=10 dphi=2. ophi=30 encode=0 > $@

5004areal_synth_d%_g.H:
	< 01sig_prep.H ./bin/OW_arealsht.x sou=02sig_sou1.H vel=1000velmig_my.H nrefvel=4 xmin=-30100 xmax=111950 ymin=0 ymax=0 zmin=0 zmax=20000 \
	wbott=0  depth_wav=$* vsal=14700 traces_in_taper=0 synsh=5014areal_synth_d$*_g.H shot_in=01sig_prep.H nphi=10 dphi=2. ophi=50 encode=0 > $@


5100areal_mig_d16000_%.H:
	< 5004areal_synth_d16000_$*.H Window3d min3=10000 max3=60000 squeeze=n | Transp plane=13 > 5100areal_synth_d16000_$*_transp.H
	< 5014areal_synth_d16000_$*.H Window3d min3=10000 max3=60000 squeeze=n | Transp plane=13 > 5101areal_synth_d16000_$*_transp.H
	< 5101areal_synth_d16000_$*_transp.H ./bin/OWmig2d_areal.x sou=5100areal_synth_d16000_$*_transp.H \
	vel=1000velmig_my.H nrefvel=4 xmin=10000 xmax=60000 ymin=0 ymax=0 zmin=0 zmax=20000 nh=33 oh=-1200 dh=75 wbott=0 vsal=14700 \
	traces_in_taper=0 velout= velout.H wr= 6121Wr_d1.5_$*.H> $@

5110areal_mig_d16000_full.H:
	Window3d < 5100areal_mig_d16000_a.H n4=1 f4=9 > 5110areal_mig_d16000_a.H
	Window3d < 5100areal_mig_d16000_b.H n4=1 f4=9 > 5110areal_mig_d16000_b.H
	Window3d < 5100areal_mig_d16000_c.H n4=1 f4=9 > 5110areal_mig_d16000_c.H
	Window3d < 5100areal_mig_d16000_d.H n4=1 f4=9 > 5110areal_mig_d16000_d.H
	Window3d < 5100areal_mig_d16000_e.H n4=1 f4=9 > 5110areal_mig_d16000_e.H
	Window3d < 5100areal_mig_d16000_f.H n4=1 f4=9 > 5110areal_mig_d16000_f.H
	Window3d < 5100areal_mig_d16000_g.H n4=1 f4=9 > 5110areal_mig_d16000_g.H
	Math file1=5110areal_mig_d16000_a.H file2=5110areal_mig_d16000_b.H file3=5110areal_mig_d16000_c.H \
	file4=5110areal_mig_d16000_d.H file5=5110areal_mig_d16000_e.H file6=5110areal_mig_d16000_f.H \
	file7=5110areal_mig_d16000_g.H exp="file1+file2+file3+file4+file5+file6+file7" > $@

5105areal_mig_d16000_%.H:
	< 5004areal_synth_d16000_$*.H Window3d min3=10000 max3=60000 squeeze=n | Transp plane=13 > 5100areal_synth_d16000_$*_transp.H
	< 5014areal_synth_d16000_$*.H Window3d min3=10000 max3=60000 squeeze=n | Transp plane=13 > 5101areal_synth_d16000_$*_transp.H
	< 5101areal_synth_d16000_$*_transp.H ./bin/OWmig2d_areal.x sou=5100areal_synth_d16000_$*_transp.H \
	vel=1000velmig_my.H nrefvel=8 xmin=10000 xmax=60000 ymin=0 ymax=0 zmin=0 zmax=20000 nh=33 oh=-1200 dh=75 wbott=0 vsal=14700 \
	traces_in_taper=0 velout= velout.H wr= 6121Wr_d1.5_$*.H> $@

5115areal_mig_d16000_full.H:
	Window3d < 5105areal_mig_d16000_a.H n4=1 f4=9 > 5115areal_mig_d16000_a.H
	Window3d < 5105areal_mig_d16000_b.H n4=1 f4=9 > 5115areal_mig_d16000_b.H
	Window3d < 5105areal_mig_d16000_c.H n4=1 f4=9 > 5115areal_mig_d16000_c.H
	Window3d < 5105areal_mig_d16000_d.H n4=1 f4=9 > 5115areal_mig_d16000_d.H
	Window3d < 5105areal_mig_d16000_e.H n4=1 f4=9 > 5115areal_mig_d16000_e.H
	Window3d < 5105areal_mig_d16000_f.H n4=1 f4=9 > 5115areal_mig_d16000_f.H
	Window3d < 5105areal_mig_d16000_g.H n4=1 f4=9 > 5115areal_mig_d16000_g.H
	Math file1=5115areal_mig_d16000_a.H file2=5115areal_mig_d16000_b.H file3=5115areal_mig_d16000_c.H \
	file4=5115areal_mig_d16000_d.H file5=5115areal_mig_d16000_e.H file6=5115areal_mig_d16000_f.H \
	file7=5115areal_mig_d16000_g.H exp="file1+file2+file3+file4+file5+file6+file7" > $@

5150areal_mig_d16000_%.H:
	< 5101areal_synth_d16000_$*_transp.H ./bin/OWmig2d_areal.x sou=5100areal_synth_d16000_$*_transp.H \
	vel=1000velmiglow_my.H nrefvel=4 xmin=10000 xmax=60000 ymin=0 ymax=0 zmin=0 zmax=20000 nh=33 oh=-1200 dh=75 wbott=0 vsal=14700 \
	traces_in_taper=0 > $@

5190areal_mig_d16000_full.H:
	Window3d < 5150areal_mig_d16000_a.H n4=1 f4=9 > 5190areal_mig_d16000_a.H
	Window3d < 5150areal_mig_d16000_b.H n4=1 f4=9 > 5190areal_mig_d16000_b.H
	Window3d < 5150areal_mig_d16000_c.H n4=1 f4=9 > 5190areal_mig_d16000_c.H
	Window3d < 5150areal_mig_d16000_d.H n4=1 f4=9 > 5190areal_mig_d16000_d.H
	Window3d < 5150areal_mig_d16000_e.H n4=1 f4=9 > 5190areal_mig_d16000_e.H
	Window3d < 5150areal_mig_d16000_f.H n4=1 f4=9 > 5190areal_mig_d16000_f.H
	Window3d < 5150areal_mig_d16000_g.H n4=1 f4=9 > 5190areal_mig_d16000_g.H
	Math file1=5190areal_mig_d16000_a.H file2=5190areal_mig_d16000_b.H file3=5190areal_mig_d16000_c.H \
	file4=5190areal_mig_d16000_d.H file5=5190areal_mig_d16000_e.H file6=5190areal_mig_d16000_f.H \
	file7=5190areal_mig_d16000_g.H exp="file1+file2+file3+file4+file5+file6+file7" > $@

5200areal_mig_d06000_%.H:
	< 5004areal_synth_d06000_$*.H Window3d min3=10000 max3=60000 squeeze=n | Transp plane=13 > 5200areal_synth_d06000_$*_transp.H
	< 5014areal_synth_d06000_$*.H Window3d min3=10000 max3=60000 squeeze=n | Transp plane=13 > 5201areal_synth_d06000_$*_transp.H
	< 5201areal_synth_d06000_$*_transp.H ./bin/OWmig2d_areal.x sou=5200areal_synth_d06000_$*_transp.H \
	vel=1000velmig_my.H nrefvel=4 xmin=10000 xmax=60000 ymin=0 ymax=0 zmin=0 zmax=20000 nh=33 oh=-1200 dh=75 wbott=0 vsal=14700 \
	traces_in_taper=0 velout= velout.H wr= 6121Wr_d1.5_$*.H> $@

5210areal_mig_d06000_full.H:
	Window3d < 5200areal_mig_d06000_a.H n4=1 f4=9 > 5210areal_mig_d06000_a.H
	Window3d < 5200areal_mig_d06000_b.H n4=1 f4=9 > 5210areal_mig_d06000_b.H
	Window3d < 5200areal_mig_d06000_c.H n4=1 f4=9 > 5210areal_mig_d06000_c.H
	Window3d < 5200areal_mig_d06000_d.H n4=1 f4=9 > 5210areal_mig_d06000_d.H
	Window3d < 5200areal_mig_d06000_e.H n4=1 f4=9 > 5210areal_mig_d06000_e.H
	Window3d < 5200areal_mig_d06000_f.H n4=1 f4=9 > 5210areal_mig_d06000_f.H
	Window3d < 5200areal_mig_d06000_g.H n4=1 f4=9 > 5210areal_mig_d06000_g.H
	Math file1=5210areal_mig_d06000_a.H file2=5210areal_mig_d06000_b.H file3=5210areal_mig_d06000_c.H \
	file4=5210areal_mig_d06000_d.H file5=5210areal_mig_d06000_e.H file6=5210areal_mig_d06000_f.H \
	file7=5210areal_mig_d06000_g.H exp="file1+file2+file3+file4+file5+file6+file7" > $@

5300areal_mig_d06000_%.H:
#	Math file1=1000velmig_my.H exp="file1*0.9" > 1000velmiglow_my.H
	< 5201areal_synth_d06000_$*_transp.H ./bin/OWmig2d_areal.x sou=5200areal_synth_d06000_$*_transp.H \
	vel=1000velmiglow_my.H nrefvel=4 xmin=10000 xmax=60000 ymin=0 ymax=0 zmin=0 zmax=20000 nh=33 oh=-1200 dh=75 wbott=0 vsal=14700 \
	traces_in_taper=0 velout= velout.H wr= 6121Wr_d1.5_$*.H> $@

5310areal_mig_d06000_full.H:
	Window3d < 5300areal_mig_d06000_a.H n4=1 f4=9 > 5310areal_mig_d06000_a.H
	Window3d < 5300areal_mig_d06000_b.H n4=1 f4=9 > 5310areal_mig_d06000_b.H
	Window3d < 5300areal_mig_d06000_c.H n4=1 f4=9 > 5310areal_mig_d06000_c.H
	Window3d < 5300areal_mig_d06000_d.H n4=1 f4=9 > 5310areal_mig_d06000_d.H
	Window3d < 5300areal_mig_d06000_e.H n4=1 f4=9 > 5310areal_mig_d06000_e.H
	Window3d < 5300areal_mig_d06000_f.H n4=1 f4=9 > 5310areal_mig_d06000_f.H
	Window3d < 5300areal_mig_d06000_g.H n4=1 f4=9 > 5310areal_mig_d06000_g.H
	Math file1=5310areal_mig_d06000_a.H file2=5310areal_mig_d06000_b.H file3=5310areal_mig_d06000_c.H \
	file4=5310areal_mig_d06000_d.H file5=5310areal_mig_d06000_e.H file6=5310areal_mig_d06000_f.H \
	file7=5310areal_mig_d06000_g.H exp="file1+file2+file3+file4+file5+file6+file7" > $@
###!!!
#
# Tem migracao 5300... rolando em glad - verificar
#
5501areal_ang_vel.H:
	< 5110areal_mig_d16000_full.H Transp plane=34 | Pad beg2=44 end2=43 end1=100 | Scale scale_to=1 |Reverse which=2 opt=ignore > xx.H
	< xx.H /net/koko/bob/SEP_BOB/bin/OFF2ANG amax=40 eps=10 maxmem=1000 verb=y > yy.H ; < yy.H Window3d max1=20000 > $@

5502areal_ang_vel.H:
	< 5210areal_mig_d06000_full.H Transp plane=34 | Pad beg2=44 end2=43 end1=100 | Scale scale_to=1 |Reverse which=2 opt=ignore > xx.H
	< xx.H /net/koko/bob/SEP_BOB/bin/OFF2ANG amax=40 eps=10 maxmem=1000 verb=y > yy.H ; < yy.H Window3d max1=20000 > $@

5503areal_ang_velow.H:
	< 5190areal_mig_d16000_full.H Transp plane=34 | Pad beg2=44 end2=43 end1=100 | Scale scale_to=1 |Reverse which=2 opt=ignore > xx.H
	< xx.H /net/koko/bob/SEP_BOB/bin/OFF2ANG amax=40 eps=10 maxmem=1000 verb=y > yy.H ; < yy.H Window3d max1=20000 > $@

5504areal_ang_velow.H:
	< 5310areal_mig_d06000_full.H Transp plane=34 | Pad beg2=44 end2=43 end1=100 | Scale scale_to=1 |Reverse which=2 opt=ignore > xx.H
	< xx.H /net/koko/bob/SEP_BOB/bin/OFF2ANG amax=40 eps=10 maxmem=1000 verb=y > yy.H ; < yy.H Window3d max1=20000 > $@

####################################
#
# Plane reflector
#
6001data_hor.H:
	susynlv ref="1:0,1.5;10.,1.5" nt=251 dt=0.008 dxs=0.02 nxs=201 fxs=0 dxo=0.02 nxo=201 fxo=0 v00=3 > junk1.su;< junk1.su Su2sep > junk1.H
	susynlv ref="1:0,-0.;10.,8." nt=251 dt=0.008 dxs=0.02 nxs=201 fxs=0 dxo=0.02 nxo=201 fxo=0 v00=3 > junk1.su;< junk1.su Su2sep > junk2.H
	Math file1=junk1.H file2=junk2.H exp="file1+file2" | Infill3d > $@; echo n2=201 o2=0 d2=0.02 n3=201 o3=0 d3=0.02 >> $@

6002data_recip.H:
	Create3d < 6001data_hor.H keyname1=offset keyname2=sx verb=1 > junk.H
	Headermath key1=cmp eqn1="sx+0.5*offset" < junk.H > junk1.H
	Sort3d < junk1.H nkeys=2 key2=cmp ng2=611 og2=0 dg2=0.01 key1=offset ng1=101 og1=0 dg1=0.04 max_size=1000 > junk2.H
	<junk2.H Infill3d >junk3.H;
	 <junk3.H Stack >junk30.H
	Interp d2out=0.02 o2out=0 n2out=201 maxsize=1000 <junk30.H >junk31.H
	Reverse < junk31.H opt=no which=2 | Window3d max2=-0.002 > junk4.H
	Merge axis=2 space=n junk4.H junk31.H > junk5.H
	Create3d < junk5.H keyname1=offset keyname2=cmp verb=1 > junk6.H
	Headermath key1=sx eqn1="cmp-0.5*offset" < junk6.H > junk7.H
	Sort3d < junk7.H nkeys=2 key2=sx ng2=201 og2=0 dg2=0.02 key1=offset ng1=401 og1=-4 dg1=0.02 max_size=1000 > junk8.H
	<junk8.H Infill3d > junk9.H
	<junk9.H Stack3d > $@

6003data_recip_prep.H:
	< 6002data_recip.H /net/koko/bob/SEP_BOB/bin/Transf shotmig=y \
	is2d=y wei=y f_min=6 f_min1=10 f_max1=50 f_max=60 maxsize=1000 > j2.H
	< j2.H Transp reshape=1,3,5 plane=23 | Pad beg1=300 end1=300 > $@

6004data_recip_source.H:
	 < 6003data_recip_prep.H Window3d n4=1 | Transp | Cabs | Stack | Smooth rect1=5 > 6004_spec.H;
	Math file1=6004_spec.H exp="file1*0" > 6004_phse.H; Cmplx 6004_spec.H 6004_phse.H |\
	Transp reshape=1,3,5 plane=23| Transp reshape=1,3,5 plane=12 > $@

6100areal_synth_d%_a.H:
	< 6003data_recip_prep.H ./bin/OW_arealsht.x sou=6004data_recip_source.H vel=7000vel_pad.H nrefvel=1 xmin=-10 xmax=14 ymin=0 \
	ymax=0 zmin=0 zmax=3 nh=1 oh=0 dh=0 vsal=6 traces_in_taper=0 depth_wav=$* shot_in=6003data_recip_prep.H \
	synth_op= 8100marm_areal_synth.H synsh=6101areal_synth_d$*_a.H nphi=10 dphi=2. ophi=-70 encode=0 > $@
 
6100areal_synth_d%_b.H:
	< 6003data_recip_prep.H ./bin/OW_arealsht.x sou=6004data_recip_source.H vel=7000vel_pad.H nrefvel=1 xmin=-10 xmax=14 ymin=0 \
	ymax=0 zmin=0 zmax=3. nh=1 oh=0 dh=0 vsal=6 traces_in_taper=0 depth_wav=$* shot_in=6003data_recip_prep.H \
	synth_op= 8100marm_areal_synth.H synsh=6101areal_synth_d$*_b.H nphi=10 dphi=2. ophi=-50 encode=0 > $@
 
6100areal_synth_d%_c.H:
	./bin/OW_arealsht.x < 6003data_recip_prep.H sou=6004data_recip_source.H vel=7000vel_pad.H nrefvel=1 xmin=-10 xmax=14 ymin=0 \
	ymax=0 zmin=0 zmax=3. nh=1 oh=0 dh=0 vsal=6 traces_in_taper=0 depth_wav=$* shot_in=6003data_recip_prep.H \
	synth_op= 8100marm_areal_synth.H synsh=6101areal_synth_d$*_c.H nphi=10 dphi=2. ophi=-30 encode=0 > $@
 
6100areal_synth_d%_d.H:
	< 6003data_recip_prep.H ./bin/OW_arealsht.x sou=6004data_recip_source.H vel=7000vel_pad.H nrefvel=1 xmin=-10 xmax=14 ymin=0 \
	ymax=0 zmin=0 zmax=3. nh=1 oh=0 dh=0 vsal=6 traces_in_taper=0 depth_wav=$* shot_in=6003data_recip_prep.H \
	synth_op= 8100marm_areal_synth.H synsh=6101areal_synth_d$*_d.H nphi=10 dphi=2. ophi=-10 encode=0 > $@
 
6100areal_synth_d%_e.H:
	< 6003data_recip_prep.H ./bin/OW_arealsht.x sou=6004data_recip_source.H vel=7000vel_pad.H nrefvel=1 xmin=-10 xmax=14 ymin=0 \
	ymax=0 zmin=0 zmax=3. nh=1 oh=0 dh=0 vsal=6 traces_in_taper=0 depth_wav=$* shot_in=6003data_recip_prep.H \
	synth_op= 8100marm_areal_synth.H synsh=6101areal_synth_d$*_e.H nphi=10 dphi=2. ophi=10 encode=0 > $@
 
6100areal_synth_d%_f.H:
	< 6003data_recip_prep.H ./bin/OW_arealsht.x sou=6004data_recip_source.H vel=7000vel_pad.H nrefvel=1 xmin=-10 xmax=14 ymin=0 \
	ymax=0 zmin=0 zmax=3. nh=1 oh=0 dh=0 vsal=6 traces_in_taper=0 depth_wav=$* shot_in=6003data_recip_prep.H \
	synth_op= 8100marm_areal_synth.H synsh=6101areal_synth_d$*_f.H nphi=10 dphi=2. ophi=30 encode=0 > $@
 
6100areal_synth_d%_g.H:
	< 6003data_recip_prep.H ./bin/OW_arealsht.x sou=6004data_recip_source.H vel=7000vel_pad.H nrefvel=1 xmin=-10 xmax=14 ymin=0 \
	ymax=0 zmin=0 zmax=3. nh=1 oh=0 dh=0 vsal=6 traces_in_taper=0 depth_wav=$* shot_in=6003data_recip_prep.H \
	synth_op= 8100marm_areal_synth.H synsh=6101areal_synth_d$*_g.H nphi=10 dphi=2. ophi=50 encode=0 > $@

6199areal_synth_d%_d.H:
	< 6003data_recip_prep.H ./bin/OW_arealsht.x sou=6004data_recip_source.H vel=7000vel.H nrefvel=1 xmin=-10 xmax=14 ymin=0 \
	ymax=0 zmin=0 zmax=2. nh=1 oh=0 dh=0 vsal=6 traces_in_taper=0 depth_wav=$* shot_in=6003data_recip_prep.H \
	synth_op= 8100marm_areal_synth.H synsh=6101areal_synth_d$*_d.H nphi=10 dphi=2. ophi=-10 encode=0 > $@

6120areal_mig_d1.5_%.H:
	< 6100areal_synth_d1.5_$*.H Window3d min3=-6 max3=8 squeeze=n n4=1 | Transp plane=13 > 6100areal_synth_d1.5_$*_transp.H
	< 6101areal_synth_d1.5_$*.H Window3d min3=-6 max3=8 squeeze=n n4=1 | Transp plane=13 > 6101areal_synth_d1.5_$*_transp.H
	< 6101areal_synth_d1.5_$*_transp.H ./bin/OWmig2d_areal.x sou=6100areal_synth_d1.5_$*_transp.H \
	vel=7000vel.H nrefvel=1 xmin=-6 xmax=8 ymin=0 ymax=0 zmin=0 zmax=2 nh=1 41 oh=0 -0.4 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 velout= velout.H wr= 6121Wr_d1.5_$*.H> $@

6130areal_mig_d0.500_ltd_full.H:
	Window3d < 6120areal_mig_d0.500_a.H n4=1 min3=0.5 max3=4 > 6130areal_mig_d0.500_ltd_a.H
	Window3d < 6120areal_mig_d0.500_b.H n4=1 min3=0.5 max3=4 > 6130areal_mig_d0.500_ltd_b.H
	Window3d < 6120areal_mig_d0.500_c.H n4=1 min3=0.5 max3=4 > 6130areal_mig_d0.500_ltd_c.H
	Window3d < 6120areal_mig_d0.500_d.H n4=1 min3=0.5 max3=4 > 6130areal_mig_d0.500_ltd_d.H
	Window3d < 6120areal_mig_d0.500_e.H n4=1 min3=0.5 max3=4 > 6130areal_mig_d0.500_ltd_e.H
	Math file1=6130areal_mig_d0.500_ltd_a.H file2=6130areal_mig_d0.500_ltd_b.H file3=6130areal_mig_d0.500_ltd_c.H \
	file4=6130areal_mig_d0.500_ltd_d.H file5=6130areal_mig_d0.500_ltd_e.H exp="file1+file2+file3+file4+file5" > $@

6140areal_mig_d0500_ltd_angle.H:
	Transp plane=34 < 6130areal_mig_d0.500_ltd_full.H | Scale scale_to=1 | Pad beg2=20 end2=20 end1=100 > xx.H
	< xx.H /net/koko/bob/SEP_BOB/bin/OFF2ANG amax=40 eps=10 maxmem=1000 verb=y > $@

7000velow.H:
	Math file1=7000vel.H exp="file1*0+2.7" > $@

7000velhi.H:
	Math file1=7000vel.H exp="file1*0+3.3" > $@

7000vel_pad.H:
	Pad

6220areal_mig_d1.5_%.H:
	< 6100areal_synth_d1.5_$*.H Window3d min3=-6 max3=8 n4=1 squeeze=n | Transp plane=13 > 6100areal_synth_d1.5_$*_transp.H
	< 6101areal_synth_d1.5_$*.H Window3d min3=-6 max3=8 n4=1 squeeze=n | Transp plane=13 > 6101areal_synth_d1.5_$*_transp.H
	< 6101areal_synth_d1.5_$*_transp.H ./bin/OWmig2d_areal.x sou=6100areal_synth_d1.5_$*_transp.H \
	vel=7000velow.H nrefvel=1 xmin=-6 xmax=8 ymin=0 ymax=0 zmin=0 zmax=2 nh=1 oh=0 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 velout= velout.H wr= 6121Wr_d1.5_$*.H> $@

6221areal_mig_d1.5_%.H:
	< 6100areal_synth_d1.5_$*.H Window3d min3=-6 max3=8 n4=1 squeeze=n | Transp plane=13 > 6100areal_synth_d1.5_$*_transp.H
	< 6101areal_synth_d1.5_$*.H Window3d min3=-6 max3=8 n4=1 squeeze=n | Transp plane=13 > 6101areal_synth_d1.5_$*_transp.H
	< 6101areal_synth_d1.5_$*_transp.H ./bin/OWmig2d_areal.x sou=6100areal_synth_d1.5_$*_transp.H \
	vel=7000vel.H nrefvel=1 xmin=-6 xmax=8 ymin=0 ymax=0 zmin=0 zmax=2 nh=1 oh=0 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 velout= velout.H wr= 6121Wr_d1.5_$*.H> $@

6230areal_mig_d0.500_ltd_full.H:
	Window3d < 6220areal_mig_d0.500_a.H n4=1 min3=0.5 max3=4 > 6230areal_mig_d0.500_ltd_a.H
	Window3d < 6220areal_mig_d0.500_b.H n4=1 min3=0.5 max3=4 > 6230areal_mig_d0.500_ltd_b.H
	Window3d < 6220areal_mig_d0.500_c.H n4=1 min3=0.5 max3=4 > 6230areal_mig_d0.500_ltd_c.H
	Window3d < 6220areal_mig_d0.500_d.H n4=1 min3=0.5 max3=4 > 6230areal_mig_d0.500_ltd_d.H
	Window3d < 6220areal_mig_d0.500_e.H n4=1 min3=0.5 max3=4 > 6230areal_mig_d0.500_ltd_e.H
	Math file1=6230areal_mig_d0.500_ltd_a.H file2=6230areal_mig_d0.500_ltd_b.H file3=6230areal_mig_d0.500_ltd_c.H \
	file4=6230areal_mig_d0.500_ltd_d.H file5=6230areal_mig_d0.500_ltd_e.H exp="file1+file2+file3+file4+file5" > $@

6240areal_mig_d0500_angle.H:
	Transp plane=34 < 6230areal_mig_d0.500_ltd_full.H | Scale scale_to=1 | Pad beg2=20 end2=20 end1=100 > xx.H
	< xx.H /net/koko/bob/SEP_BOB/bin/OFF2ANG amax=40 eps=10 maxmem=1000 verb=y > $@

6320areal_mig_d0.500_%.H:
	< 6101areal_synth_d0.500_$*_transp.H ./bin/OWmig2d_areal.x sou=6100areal_synth_d0.500_$*_transp.H \
	vel=7000velhi.H nrefvel=1 xmin=-6 xmax=10 ymin=0 ymax=0 zmin=0 zmax=1 nh=41 oh=-.4 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 velout= velout.H wr= 6121Wr_d0.500_$*.H> $@

6330areal_mig_d0.500_ltd_full.H:
	Window3d < 6320areal_mig_d0.500_a.H n4=1 min3=0.5 max3=4 > 6330areal_mig_d0.500_ltd_a.H
	Window3d < 6320areal_mig_d0.500_b.H n4=1 min3=0.5 max3=4 > 6330areal_mig_d0.500_ltd_b.H
	Window3d < 6320areal_mig_d0.500_c.H n4=1 min3=0.5 max3=4 > 6330areal_mig_d0.500_ltd_c.H
	Window3d < 6320areal_mig_d0.500_d.H n4=1 min3=0.5 max3=4 > 6330areal_mig_d0.500_ltd_d.H
	Window3d < 6320areal_mig_d0.500_e.H n4=1 min3=0.5 max3=4 > 6330areal_mig_d0.500_ltd_e.H
	Math file1=6330areal_mig_d0.500_ltd_a.H file2=6330areal_mig_d0.500_ltd_b.H file3=6330areal_mig_d0.500_ltd_c.H \
	file4=6330areal_mig_d0.500_ltd_d.H file5=6330areal_mig_d0.500_ltd_e.H exp="file1+file2+file3+file4+file5" > $@

6340areal_mig_d0500_ltd_angle.H:
	Transp plane=34 < 6330areal_mig_d0.500_ltd_full.H | Scale scale_to=1 | Pad beg2=20 end2=20 end1=100 > xx.H
	< xx.H /net/koko/bob/SEP_BOB/bin/OFF2ANG amax=40 eps=10 maxmem=1000 verb=y > $@

6420_mig_brad.H:
	@date
	/net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX86_64/ShotMig \
	D=6003data_recip_prep.H D_d=6004data_recip_source.H V=7000vel_brad.H R=$@ \
	pverb=2  nx_shots=201 ox_shots=0 \
	par=difr_go_sis.P taper=y taper_pcent=2 marine=y minx=-10.5 maxx=10.5 & #restart=y

6425_mig_brad.H:
	@date
	/net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX86_64/ShotMig \
	D=6003data_recip_prep.H D_d=6004data_recip_source.H V=7000velow_brad.H R=$@ \
	pverb=2  nx_shots=201 ox_shots=0 \
	par=difr_go_sis.P taper=y taper_pcent=2 marine=y minx=-10.5 maxx=10.5 & #restart=y

6501data_hor.H:
	susynlv ref="1:0,1.5;10.,1.5" nt=501 dt=0.008 dxs=0.02 nxs=201 fxs=0 dxo=0.02 nxo=401 fxo=0 v00=3 > junk1.su;< junk1.su Su2sep > junk1.H
	susynlv ref="1:0,-0.;10.,8." nt=501 dt=0.008 dxs=0.02 nxs=201 fxs=0 dxo=0.02 nxo=401 fxo=0 v00=3 > junk1.su;< junk1.su Su2sep > junk2.H
	Math file1=junk1.H file2=junk2.H exp="file1+file2" | Infill3d > $@; echo n2=401 o2=0 d2=0.02 n3=201 o3=0 d3=0.02 >> $@

6502data_recip.H:
	Create3d < 6501data_hor.H keyname1=offset keyname2=sx verb=1 > junk.H
	Headermath key1=cmp eqn1="sx+0.5*offset" < junk.H > junk1.H
	Sort3d < junk1.H nkeys=2 key2=cmp ng2=811 og2=0 dg2=0.01 key1=offset ng1=201 og1=0 dg1=0.04 max_size=1000 > junk2.H
	<junk2.H Infill3d >junk3.H;
	 <junk3.H Stack >junk30.H
	Interp d2out=0.02 o2out=0 n2out=401 maxsize=1000 <junk30.H >junk31.H
	Reverse < junk31.H opt=no which=2 | Window3d max2=-0.002 > junk4.H
	Merge axis=2 space=n junk4.H junk31.H > junk5.H
	Create3d < junk5.H keyname1=offset keyname2=cmp verb=1 > junk6.H
	Headermath key1=sx eqn1="cmp-0.5*offset" < junk6.H > junk7.H
	Sort3d < junk7.H nkeys=2 key2=sx ng2=201 og2=0 dg2=0.02 key1=offset ng1=801 og1=-8 dg1=0.02 max_size=1000 > junk8.H
	<junk8.H Infill3d > junk9.H
	<junk9.H Stack3d > $@

6503data_recip_prep.H:
	< 6002data_recip.H /net/koko/bob/SEP_BOB/bin/Transf shotmig=y \
	is2d=y wei=y f_min=6 f_min1=10 f_max1=50 f_max=60 maxsize=1000 > j2.H
	< j2.H Transp reshape=1,3,5 plane=23 | Pad beg1=300 end1=300 > $@

####################################
#
# Investigating cross-talk
#
7000vel.H:
	Spike n1=251 n2=2000 n3=1 o1=0 o2=-4 o3=0 d1=0.008 d2=0.02 d3=0.001 nsp=1 > temp.H ; Math file1=temp.H exp="file1*0+3" > $@

7001data.H:
#	susynlv ref="1:3,0.5;3.02,0.5,2:2,0.5;2.02,0.5" nt=501 dt=0.004 dxs=0.02 nxs=201 fxs=0 dxo=0.02 nxo=301 fxo=0 v00=3 > junk.su
#	susynlv ref="1:2.0,0.5;2.02,0.5,1:2.12,0.5;2.14,0.5,1:2.24,0.5;2.26.02,0.5,1:2.5,0.5;2.52,0.5,1:3,0.5;3.02,0.5" nt=501 dt=0.004 dxs=0.02 nxs=201 #fxs=0 dxo=0.02 nxo=301 fxo=0 v00=3 > junk.su
	susynlv ref="1:2.0,0.5;2.02,0.5" nt=501 dt=0.004 dxs=0.02 nxs=201 fxs=0 dxo=0.02 nxo=301 fxo=0 v00=3 > junk1.su;< junk1.su Su2sep > junk1.H
	susynlv ref="1:2.12,0.5;2.14,0.5" nt=501 dt=0.004 dxs=0.02 nxs=201 fxs=0 dxo=0.02 nxo=301 fxo=0 v00=3 > junk2.su;< junk2.su Su2sep > junk2.H
	susynlv ref="1:2.24,0.5;2.26,0.5" nt=501 dt=0.004 dxs=0.02 nxs=201 fxs=0 dxo=0.02 nxo=301 fxo=0 v00=3 > junk3.su;< junk3.su Su2sep > junk3.H
	susynlv ref="1:2.5,0.5;2.52,0.5" nt=501 dt=0.004 dxs=0.02 nxs=201 fxs=0 dxo=0.02 nxo=301 fxo=0 v00=3 > junk4.su;< junk4.su Su2sep > junk4.H
	susynlv ref="1:3.0,0.5;3.02,0.5" nt=501 dt=0.004 dxs=0.02 nxs=201 fxs=0 dxo=0.02 nxo=301 fxo=0 v00=3 > junk5.su;< junk5.su Su2sep > junk5.H
	Math file1=junk1.H file2=junk2.H file3=junk3.H file4=junk4.H file5=junk5.H exp="file1+file2+file3+file4+file5" | Infill3d > $@
	echo n2=301 o2=0 d2=0.02 n3=201 o3=0 d3=0.02 >> $@

7002data_recip.H:
	Create3d < 7001data.H keyname1=offset keyname2=sx verb=1 > junk.H
	Headermath key1=cmp eqn1="sx+0.5*offset" < junk.H > junk1.H
	Sort3d < junk1.H nkeys=2 key2=cmp ng2=711 og2=0 dg2=0.01 key1=offset ng1=151 og1=0 dg1=0.04 max_size=1000 > junk2.H
	<junk2.H Infill3d >junk3.H;
	 <junk3.H Stack >junk30.H
	Interp d2out=0.02 o2out=0 n2out=301 maxsize=1000 <junk30.H >junk31.H
	Reverse < junk31.H opt=no which=2 | Window3d max2=-0.002 > junk4.H
	Merge axis=2 space=n junk4.H junk31.H > junk5.H
	Create3d < junk5.H keyname1=offset keyname2=cmp verb=1 > junk6.H
	Headermath key1=sx eqn1="cmp-0.5*offset" < junk6.H > junk7.H
	Sort3d < junk7.H nkeys=2 key2=sx ng2=201 og2=0 dg2=0.02 key1=offset ng1=601 og1=-6 dg1=0.02 max_size=1000 > junk8.H
	<junk8.H Infill3d > junk9.H
	<junk9.H Stack3d > $@

7003data_recip_prep.H:
	< 7002data_recip.H /net/koko/bob/SEP_BOB/bin/Transf shotmig=y \
	is2d=y wei=y f_min=6 f_min1=10 f_max1=60 f_max=80 maxsize=1000 > j2.H
	< j2.H Transp reshape=1,3,5 plane=23 | Pad beg1=100 end1=100 > $@

7004data_recip_source.H:
	 < 7003data_recip_prep.H Window3d n4=1 | Transp | Cabs | Stack | Smooth rect1=5 > 7004_spec.H;
	Math file1=7004_spec.H exp="file1*0" > 7004_phse.H; Cmplx 7004_spec.H 7004_phse.H |\
	Transp reshape=1,3,5 plane=23| Transp reshape=1,3,5 plane=12 > $@

7100areal_synth_d%_a.H:
	< 7003data_recip_prep.H ./bin/OW_arealsht.x sou=7004data_recip_source.H vel=7000vel.H nrefvel=1 xmin=-8.2 xmax=12.2 ymin=0 \
	ymax=0 zmin=0 zmax=1. nh=1 oh=0 dh=0 vsal=6 traces_in_taper=0 depth_wav=$* shot_in=7003data_recip_prep.H \
	synth_op= 8100marm_areal_synth.H synsh=7101areal_synth_d$*_a.H nphi=10 dphi=2. ophi=-50 encode=0 > $@
 
7100areal_synth_d%_b.H:
	< 7003data_recip_prep.H ./bin/OW_arealsht.x sou=7004data_recip_source.H vel=7000vel.H nrefvel=1 xmin=-8.2 xmax=12.2 ymin=0 \
	ymax=0 zmin=0 zmax=1. nh=1 oh=0 dh=0 vsal=6 traces_in_taper=0 depth_wav=$* shot_in=7003data_recip_prep.H \
	synth_op= 8100marm_areal_synth.H synsh=7101areal_synth_d$*_b.H nphi=10 dphi=2. ophi=-30 encode=0 > $@
 
7100areal_synth_d%_c.H:
	< 7003data_recip_prep.H ./bin/OW_arealsht.x sou=7004data_recip_source.H vel=7000vel.H nrefvel=1 xmin=-8.2 xmax=12.2 ymin=0 \
	ymax=0 zmin=0 zmax=1. nh=1 oh=0 dh=0 vsal=6 traces_in_taper=0 depth_wav=$* shot_in=7003data_recip_prep.H \
	synth_op= 8100marm_areal_synth.H synsh=7101areal_synth_d$*_c.H nphi=10 dphi=2. ophi=-10 encode=0 > $@
 
7100areal_synth_d%_d.H:
	< 7003data_recip_prep.H ./bin/OW_arealsht.x sou=7004data_recip_source.H vel=7000vel.H nrefvel=1 xmin=-8.2 xmax=12.2 ymin=0 \
	ymax=0 zmin=0 zmax=1. nh=1 oh=0 dh=0 vsal=6 traces_in_taper=0 depth_wav=$* shot_in=7003data_recip_prep.H \
	synth_op= 8100marm_areal_synth.H synsh=7101areal_synth_d$*_d.H nphi=10 dphi=2. ophi=10 encode=0 > $@
 
7100areal_synth_d%_e.H:
	< 7003data_recip_prep.H ./bin/OW_arealsht.x sou=7004data_recip_source.H vel=7000vel.H nrefvel=1 xmin=-8.2 xmax=12.2 ymin=0 \
	ymax=0 zmin=0 zmax=1. nh=1 oh=0 dh=0 vsal=6 traces_in_taper=0 depth_wav=$* shot_in=7003data_recip_prep.H \
	synth_op= 8100marm_areal_synth.H synsh=7101areal_synth_d$*_e.H nphi=10 dphi=2. ophi=30 encode=0 > $@
 
7501marm_areal_mig_d0.500_%.H:
	< 7100areal_synth_d0.500_$*.H Window3d min3=-8 n3=1001 squeeze=n | \
	Transp plane=13 > 7100areal_synth_d0.500_$*_transp.H
	Window3d squeeze=n < 7101areal_synth_d0.500_$*.H min3=-8 n3=1001 | \
	Transp plane=13 > 7101areal_synth_d0.500_$*_transp.H
	< 7101areal_synth_d0.500_$*_transp.H ./bin/OWmig2d_areal.x sou=7100areal_synth_d0.500_$*_transp.H \
	vel=7000vel.H nrefvel=1 xmin=-8 xmax=12 ymin=0 ymax=0 zmin=0 zmax=1 nh=1 oh=0 dh=24 wbott=0 vsal=5 \
	traces_in_taper=0 velout= velout.H > $@

7200areal_synth_d%_a.H:
	< 7003data_recip_prep.H ./bin/OW_arealsht.x sou=7004data_recip_source.H vel=7000vel.H nrefvel=1 xmin=-8.2 xmax=12.2 ymin=0 \
	ymax=0 zmin=0 zmax=1. nh=1 oh=0 dh=0 vsal=6 traces_in_taper=0 depth_wav=$* shot_in=7003data_recip_prep.H \
	synth_op= 8100marm_areal_synth.H synsh=7201areal_synth_d$*_a.H nphi=10 dphi=2. ophi=-50 encode=1 > $@
 
7200areal_synth_d%_b.H:
	< 7003data_recip_prep.H ./bin/OW_arealsht.x sou=7004data_recip_source.H vel=7000vel.H nrefvel=1 xmin=-8.2 xmax=12.2 ymin=0 \
	ymax=0 zmin=0 zmax=1. nh=1 oh=0 dh=0 vsal=6 traces_in_taper=0 depth_wav=$* shot_in=7003data_recip_prep.H \
	synth_op= 8100marm_areal_synth.H synsh=7201areal_synth_d$*_b.H nphi=10 dphi=2. ophi=-30 encode=1 > $@
 
7200areal_synth_d%_c.H:
	< 7003data_recip_prep.H ./bin/OW_arealsht.x sou=7004data_recip_source.H vel=7000vel.H nrefvel=1 xmin=-8.2 xmax=12.2 ymin=0 \
	ymax=0 zmin=0 zmax=1. nh=1 oh=0 dh=0 vsal=6 traces_in_taper=0 depth_wav=$* shot_in=7003data_recip_prep.H \
	synth_op= 8100marm_areal_synth.H synsh=7201areal_synth_d$*_c.H nphi=10 dphi=2. ophi=-10 encode=1 > $@
 
7200areal_synth_d%_d.H:
	< 7003data_recip_prep.H ./bin/OW_arealsht.x sou=7004data_recip_source.H vel=7000vel.H nrefvel=1 xmin=-8.2 xmax=12.2 ymin=0 \
	ymax=0 zmin=0 zmax=1. nh=1 oh=0 dh=0 vsal=6 traces_in_taper=0 depth_wav=$* shot_in=7003data_recip_prep.H \
	synth_op= 8100marm_areal_synth.H synsh=7201areal_synth_d$*_d.H nphi=10 dphi=2. ophi=10 encode=1 > $@
 
7200areal_synth_d%_e.H:
	< 7003data_recip_prep.H ./bin/OW_arealsht.x sou=7004data_recip_source.H vel=7000vel.H nrefvel=1 xmin=-8.2 xmax=12.2 ymin=0 \
	ymax=0 zmin=0 zmax=1. nh=1 oh=0 dh=0 vsal=6 traces_in_taper=0 depth_wav=$* shot_in=7003data_recip_prep.H \
	synth_op= 8100marm_areal_synth.H synsh=7201areal_synth_d$*_e.H nphi=10 dphi=2. ophi=30 encode=1 > $@
 
7300areal_synth_d%_a.H:
	< 7003data_recip_prep.H Window3d squeeze=n j4=10 > tempa.H
	< tempa.H ./bin/OW_arealsht.x sou=7004data_recip_source.H vel=7000vel.H nrefvel=1 xmin=-8.2 xmax=12.2 ymin=0 \
	ymax=0 zmin=0 zmax=1. nh=1 oh=0 dh=0 vsal=6 traces_in_taper=0 depth_wav=$* shot_in=tempa.H \
	synth_op=7100areal_synth_d$*_a.H synsh=7301areal_synth_d$*_a.H nphi=10 dphi=2. ophi=-50 encode=0 > $@

7502areal_mig_d0.500_%.H:
	< 7200areal_synth_d0.500_$*.H Window3d min3=-8 n3=1001 squeeze=n | \
	Transp plane=13 > 7200areal_synth_d0.500_$*_transp.H
	Window3d squeeze=n < 7201areal_synth_d0.500_$*.H min3=-8 n3=1001 | \
	Transp plane=13 > 7201areal_synth_d0.500_$*_transp.H
	< 7201areal_synth_d0.500_$*_transp.H ./bin/OWmig2d_areal.x sou=7200areal_synth_d0.500_$*_transp.H \
	vel=7000vel.H nrefvel=1 xmin=-8 xmax=12 ymin=0 ymax=0 zmin=0 zmax=1 nh=1 oh=0 dh=24 wbott=0 vsal=5 \
	traces_in_taper=0 velout= velout.H > $@

7503areal_mig_d0.500_%.H:
	< 7300areal_synth_d0.500_$*.H Window3d min3=-8 n3=1001 squeeze=n | \
	Transp plane=13 > 7300areal_synth_d0.500_$*_transp.H
	Window3d squeeze=n < 7301areal_synth_d0.500_$*.H min3=-8 n3=1001 | \
	Transp plane=13 > 7301areal_synth_d0.500_$*_transp.H
	< 7301areal_synth_d0.500_$*_transp.H ./bin/OWmig2d_areal.x sou=7300areal_synth_d0.500_$*_transp.H \
	vel=7000vel.H nrefvel=1 xmin=-8 xmax=12 ymin=0 ymax=0 zmin=0 zmax=1 nh=1 oh=0 dh=24 wbott=0 vsal=5 \
	traces_in_taper=0 velout= velout.H > $@

7505mig_brad.H:
	@date
	/net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX86_64/ShotMig \
	D=7003data_recip_prep.H D_d=7004data_recip_source.H V=7000vel_brad.H R=$@ \
	pverb=2  nx_shots=201 ox_shots=0 \
	par=difr_go_sis.P taper=y taper_pcent=2 marine=y minx=-8 maxx=12 & #restart=y

############################################################
04sig_mig_brad.H:
#	< 01sig_prep.H Window3d squeeze=n min4=12000 max4=72000 > z01sig_prep.H
#	Math file1=1000velmig_transp.H exp="file1*0.9" |Window3d squeeze=n max3=20000 > 1000velmiglow_transp.H 
	@date
	nohup /net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX86_64/ShotMig \
	D=z01sig_prep.H D_d=02sig_sou1.H V=1000velmig_transp.H R=$@ \
	pverb=2  nx_shots=400 ox_shots=12000 \
	par=go_sis.P taper=y taper_pcent=2 marine=y minx=-34000 maxx=120000 & #restart=y

z02sig_sou1.H:
	Window3d n4=1 < ../OneWay/sigsb2a.rec.freq.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&


04sig2a_mig_brad.H:
#	< ../OneWay/sigsb2a.rec.freq.H Pad beg1=200 end1=200| Window3d squeeze=n min4=12000 max4=72000 > z02sig_prep.H
#	Math file1=1000velmig_transp.H exp="file1*0.9" |Window3d squeeze=n max3=20000 > 1000velmiglow_transp.H 
	@date
	nohup /net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX86_64/ShotMig \
	D=z02sig_prep.H D_d=z02sig_sou1.H V=1000velmig_transp.H R=$@ \
	pverb=2  nx_shots=400 ox_shots=12000 \
	par=go_sis.P taper=y taper_pcent=2 marine=y minx=-34000 maxx=120000 & #restart=y

05sig_ang_brad.H:
	< 04sig_mig_brad.H Window3d min1=10025 n1=666 squeeze=n | Transp reshape=2,4,5 plane=13  | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG amax=40 eps=10 maxmem=1000 verb=y  > junk1.H; <junk1.H Window3d max1=20000 > $@

04sig_mig_brad_low.H:
	< 01sig_prep.H Window3d squeeze=n min4=12000 max4=72000 > z01sig_prep.H
	Math file1=1000velmig_transp.H exp="file1*0.9" |Window3d squeeze=n max3=20000 > 1000velmiglow_transp.H 
	@date
	nohup /net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX86_64/ShotMig \
	D=z01sig_prep.H D_d=02sig_sou1.H V=1000velmiglow_transp.H R=$@ \
	pverb=2  nx_shots=400 ox_shots=12000 \
	par=go_sis.P taper=y taper_pcent=2 marine=y minx=-34000 maxx=120000 & #restart=y

05sig_ang_brad_low.H:
	< 04sig_mig_brad_low.H Window3d min1=10025 n1=666 squeeze=n | Transp reshape=2,4,5 plane=13  | Pad beg2=43 end2=44 end1=100 | /net/koko/bob/SEP_BOB/bin/OFF2ANG amax=40 eps=10 maxmem=1000 verb=y  > junk1.H; <junk1.H Window3d max1=20000 > $@