####################################
#
# 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=8000 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:
#	Window3d < ../RTM/sig_shot_prep_mad.H n4=1 min4=50000 squeeze=n > xxx.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=8000 nh=1 oh=0 dh=0 wbott=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:
#	< ../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 > $@

01sig_prep_cut.H:
	< ../OneWay/nfs-splitspread-shots.H Window3d min1=2 max1=6 max3=30000 | /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_cut.H:
	Window3d n4=1 < 01sig_prep_cut.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_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 > $@ 

###########################################################
#
# Testing - GF
#
999sig_mult_source.H:
	<../RTM/sig_illu_sou.H Pad beg1=17599 end1=17600 | /net/koko/bob/SEP_BOB/bin/Transf shotmig=y is2d=y wei=y f_min=3 f_min1=5 f_max1=60 f_max=80 maxsize=1000 | Transp plane=53 > $@ 

999sig_mult_illum.H:
	< 999input_sig_mult_source01.H ./bin/OWmig2d_par_DD.x sou=999sig_mult_source.H vel=velocidade1_pad_int25m.H nrefvel=8 xmin=-1200 xmax=82000 ymin=0 ymax=0 zmin=0 zmax=8000 nh=1 oh=0 dh=0 wbott=0 vsal=14700 G1=G2.H traces_in_taper=30> $@

#999sig_mult_illum.H:
#	< 999input_sig_mult_source00.H ./bin/OWmig2d_par_DD.x sou=999sig_mult_source.H vel=velocidade1_pad_int25m.H nrefvel=4 xmin=30025 #xmax=50025 ymin=0 ymax=0 zmin=0 zmax=8000 nh=1 oh=0 dh=0 wbott=0 vsal=14700 traces_in_taper=30 G1=G1.H > $@

999sig_mig01.H:
	<999input.H ./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> $@
############################################################

#############################################################
#
# 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
#
01sig_prep75.H:
	< nfs-splitspread-shots_75.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 > $@

5003areal_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=5013areal_synth_d$*_a.H shot_in=01sig_prep.H nphi=10 dphi=2. ophi=-70 encode=1 > $@

5003areal_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=5013areal_synth_d$*_b.H shot_in=01sig_prep.H nphi=10 dphi=2. ophi=-50 encode=1 > $@

5003areal_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=5013areal_synth_d$*_c.H shot_in=01sig_prep.H nphi=10 dphi=2. ophi=-30 encode=1 > $@

5003areal_synth_d%_d.H:
	< 01sig_prep75.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=5013areal_synth_d$*_d.H shot_in=01sig_prep.H nphi=10 dphi=2. ophi=-10 encode=1 > $@

5003areal_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=5013areal_synth_d$*_e.H shot_in=01sig_prep.H nphi=10 dphi=2. ophi=10 encode=1 > $@

5003areal_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=5013areal_synth_d$*_f.H shot_in=01sig_prep.H nphi=10 dphi=2. ophi=30 encode=1 > $@

5003areal_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=5013areal_synth_d$*_g.H shot_in=01sig_prep.H nphi=10 dphi=2. ophi=50 encode=1 > $@

5003areal_mig_d16000_%.H:
	< 5003areal_synth_d16000_$*.H Window3d min3=10000 max3=60000 squeeze=n | Transp plane=13 > 5003areal_synth_d16000_$*_transp.H
	< 5013areal_synth_d16000_$*.H Window3d min3=10000 max3=60000 squeeze=n | Transp plane=13 > 5013areal_synth_d16000_$*_transp.H
	< 5013areal_synth_d16000_$*_transp.H ./bin/OWmig2d_areal.x sou=5003areal_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> $@

5003tst_mig_d16000_d.H:
	< 5003areal_synth_d16000_d.H Window3d min3=10000 max3=60000 n4=1 squeeze=n | Transp plane=13 > 5003areal_tst_d16000_d_transp.H
	< 5013areal_synth_d16000_d.H Window3d min3=10000 max3=60000 n4=1 squeeze=n | Transp plane=13 > 5013areal_tst_d16000_d_transp.H
	< 5013areal_tst_d16000_d_transp.H ./bin/OWmig2d_areal.x sou=5003areal_tst_d16000_d_transp.H \
	vel=1000velmig_my.H nrefvel=4 xmin=10000 xmax=60000 ymin=0 ymax=0 zmin=0 zmax=20000 nh=1 33 oh=0 -1200 dh=75 wbott=0 vsal=14700 \
	traces_in_taper=0 > $@

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" > $@

5110areal_mig_d16000_50-50.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
	Math file1=5110areal_mig_d16000_b.H file2=5110areal_mig_d16000_c.H \
	file3=5110areal_mig_d16000_d.H file4=5110areal_mig_d16000_e.H file5=5110areal_mig_d16000_f.H \
	exp="file1+file2+file3+file4+file5" > $@ 

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> $@

5106areal_mig_d16000_%.H: # Verifying if imaging after summing frequencies solves crosstalk
#	< 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 ../Tese/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=1 33 oh=0 -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 > $@

#
#
#

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

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

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

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

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

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

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

5710areal_mig_d06000_%.H:
	< 5704areal_synth_d06000_$*.H Window3d min3=10000 max3=60000 squeeze=n | Transp plane=13 > 5700areal_synth_d06000_$*_transp.H
	< 5714areal_synth_d06000_$*.H Window3d min3=10000 max3=60000 squeeze=n | Transp plane=13 > 5701areal_synth_d06000_$*_transp.H
	< 5701areal_synth_d06000_$*_transp.H ./bin/OWmig2d_areal.x sou=5700areal_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=13000 \
	traces_in_taper=0 velout= velout.H wr= 6121Wr_d1.5_$*.H> $@

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

5720areal_mig_d16000_%.H:
	< 5704areal_synth_d16000_$*.H Window3d min3=10000 max3=60000 squeeze=n | Transp plane=13 > 5700areal_synth_d16000_$*_transp.H
	< 5714areal_synth_d16000_$*.H Window3d min3=10000 max3=60000 squeeze=n | Transp plane=13 > 5701areal_synth_d16000_$*_transp.H
	< 5701areal_synth_d16000_$*_transp.H ./bin/OWmig2d_areal.x sou=5700areal_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=13000 \
	traces_in_taper=0 velout= velout.H wr= 6121Wr_d1.5_$*.H> $@

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

5723areal_ang_velow.H:
	< 5721areal_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 > $@

5713areal_ang_velow.H:
	< 5711areal_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 > $@

5715semblance.H:
#	< 5501areal_ang_vel.H Window3d min3=15000 max3=55000 j3=10 > temp.H 
#	< 5504areal_ang_velow.H Window3d min3=15000 max3=55000 j3=25 > temp.H
	< 05sig_ang_brad_low.H Window3d min3=15000 max3=55000 j3=25 | Agc  detect=1 dwind=50 thresh=10 > temp.H
	<temp.H ./bin/OW_velan.x > $@

5780plane_wavemig_%.H:
	< 01sig_temp.H ./bin/OWmig2d_plane.x sou=02sig_sou1.H \
	vel=1000velmiglow_my.H nrefvel=4 xmin=10925 xmax=60000 ymin=0 ymax=0 zmin=0 zmax=20000 nh=33 oh=-1200 dh=75 wbott=0 vsal=14700 \
	traces_in_taper=50 padtraces=0 ray_par_n=10 prop_ang_o=$* prop_ang_d=2. planeshot= 5790planeshot.H  > $@

5790plane_wavemig_%.H:
#	< 01sig_prep.H Window3d max4=60000 min1=-26025 max1=26025 squeeze=n > 01sig_temp.H
	< 01sig_temp.H ./bin/OWmig2d_plane.x sou=02sig_sou1.H \
	vel=1000velmig_my.H nrefvel=4 xmin=10925 xmax=60000 ymin=0 ymax=0 zmin=0 zmax=20000 nh=33 oh=-1200 dh=75 wbott=0 vsal=14700 \
	traces_in_taper=50 padtraces=0 ray_par_n=10 prop_ang_o=$* prop_ang_d=2. planeshot= 5790planeshot.H  > $@

5791plane_wavemig_full.H:
	Window3d < 5790plane_wavemig_-70.H n4=1 f4=9 > 5791plane_wavemig_a.H
	Window3d < 5790plane_wavemig_-50.H n4=1 f4=9 > 5791plane_wavemig_b.H
	Window3d < 5790plane_wavemig_-30.H n4=1 f4=9 > 5791plane_wavemig_c.H
	Window3d < 5790plane_wavemig_-10.H n4=1 f4=9 > 5791plane_wavemig_d.H
	Window3d < 5790plane_wavemig_10.H n4=1 f4=9 > 5791plane_wavemig_e.H
	Window3d < 5790plane_wavemig_30.H n4=1 f4=9 > 5791plane_wavemig_f.H
	Window3d < 5790plane_wavemig_50.H n4=1 f4=9 > 5791plane_wavemig_g.H
	Math file1=5791plane_wavemig_a.H file2=5791plane_wavemig_b.H file3=5791plane_wavemig_c.H \
	file4=5791plane_wavemig_d.H file5=5791plane_wavemig_e.H file6=5791plane_wavemig_f.H \
	file7=5791plane_wavemig_g.H exp="file1+file2+file3+file4+file5+file6+file7" > $@

5792plane_wavemig_%.H:# xmed is half-way
	< 01sig_temp.H ./bin/OWmig2d_plane.x sou=02sig_sou1.H \
	vel=1000velmig_my.H nrefvel=4 xmin=10925 xmax=60000 ymin=0 ymax=0 zmin=0 zmax=20000 nh=33 oh=-1200 dh=75 wbott=0 vsal=14700 \
	traces_in_taper=50 padtraces=0 ray_par_n=10 prop_ang_o=$* prop_ang_d=2. planeshot= 5790planeshot.H  > $@

5793plane_wavemig_full.H:
	Window3d < 5792plane_wavemig_-70.H n4=1 f4=9 > 5792plane_wavemig_a.H
	Window3d < 5792plane_wavemig_-50.H n4=1 f4=9 > 5792plane_wavemig_b.H
	Window3d < 5792plane_wavemig_-30.H n4=1 f4=9 > 5792plane_wavemig_c.H
	Window3d < 5792plane_wavemig_-10.H n4=1 f4=9 > 5792plane_wavemig_d.H
	Window3d < 5792plane_wavemig_10.H n4=1 f4=9 > 5792plane_wavemig_e.H
	Window3d < 5792plane_wavemig_30.H n4=1 f4=9 > 5792plane_wavemig_f.H
	Window3d < 5792plane_wavemig_50.H n4=1 f4=9 > 5792plane_wavemig_g.H
	Math file1=5792plane_wavemig_a.H file2=5792plane_wavemig_b.H file3=5792plane_wavemig_c.H \
	file4=5792plane_wavemig_d.H file5=5792plane_wavemig_e.H file6=5792plane_wavemig_f.H \
	file7=5792plane_wavemig_g.H exp="file1+file2+file3+file4+file5+file6+file7" > $@

5794plane_wavemig_%.H:# xmed is half-way and alias control and constant angles
	< 01sig_temp.H ./bin/OWmig2d_plane.x sou=02sig_sou1.H \
	vel=1000velmig_my.H nrefvel=4 xmin=10925 xmax=60000 ymin=0 ymax=0 zmin=0 zmax=20000 nh=33 oh=-1200 dh=75 wbott=0 vsal=14700 \
	traces_in_taper=50 padtraces=0 ray_par_n=10 prop_ang_o=$* prop_ang_d=2. planeshot= 5790planeshot.H  > $@

5795plane_wavemig_full.H:
	Window3d < 5794plane_wavemig_-70.H n4=1 f4=9 > 5794plane_wavemig_a.H
	Window3d < 5794plane_wavemig_-50.H n4=1 f4=9 > 5794plane_wavemig_b.H
	Window3d < 5794plane_wavemig_-30.H n4=1 f4=9 > 5794plane_wavemig_c.H
	Window3d < 5794plane_wavemig_-10.H n4=1 f4=9 > 5794plane_wavemig_d.H
	Window3d < 5794plane_wavemig_10.H n4=1 f4=9 > 5794plane_wavemig_e.H
	Window3d < 5794plane_wavemig_30.H n4=1 f4=9 > 5794plane_wavemig_f.H
	Window3d < 5794plane_wavemig_50.H n4=1 f4=9 > 5794plane_wavemig_g.H
	Math file1=5794plane_wavemig_a.H file2=5794plane_wavemig_b.H file3=5794plane_wavemig_c.H \
	file4=5794plane_wavemig_d.H file5=5794plane_wavemig_e.H file6=5794plane_wavemig_f.H \
	file7=5794plane_wavemig_g.H exp="file1+file2+file3+file4+file5+file6+file7" > $@

5796plane_wavemig_%.H:# xmed is half-way and alias control and constant Dp
	< 01sig_temp.H ./bin/OWmig2d_plane.x sou=02sig_sou1.H \
	vel=1000velmig_my.H nrefvel=4 xmin=10925 xmax=60000 ymin=0 ymax=0 zmin=0 zmax=20000 nh=1 33 oh=0 -1200 dh=75 wbott=0 vsal=14700 \
	traces_in_taper=50 padtraces=0 ray_par_n=20 prop_ang_o=$* prop_ang_d=0.000002 planeshot= 5790planeshot.H  > $@

5797plane_wavemig_full.H:
	Window3d < 5796plane_wavemig_-.00014.H n4=1 f4=16 19 > 5796plane_wavemig_a.H
	Window3d < 5796plane_wavemig_-.00010.H n4=1 f4=16 19 > 5796plane_wavemig_b.H
	Window3d < 5796plane_wavemig_-.00006.H n4=1 f4=16 19 > 5796plane_wavemig_c.H
	Window3d < 5796plane_wavemig_-.00002.H n4=1 f4=16 19 > 5796plane_wavemig_d.H
	Window3d < 5796plane_wavemig_.00002.H n4=1 f4=16 19 > 5796plane_wavemig_e.H
	Window3d < 5796plane_wavemig_.00006.H n4=1 f4=16 19 > 5796plane_wavemig_f.H
	Window3d < 5796plane_wavemig_.00010.H n4=1 f4=16 19 > 5796plane_wavemig_g.H
	Math file1=5796plane_wavemig_a.H file2=5796plane_wavemig_b.H file3=5796plane_wavemig_c.H \
	file4=5796plane_wavemig_d.H file5=5796plane_wavemig_e.H file6=5796plane_wavemig_f.H \
	file7=5796plane_wavemig_g.H exp="file1+file2+file3+file4+file5+file6+file7" > $@

5798plane_wavemig_%.H:# xmed is half-way and alias control and constant angle
	< 01sig_temp.H ./bin/OWmig2d_plane.x sou=02sig_sou1.H \
	vel=1000velmig_my.H nrefvel=4 xmin=10925 xmax=60000 ymin=0 ymax=0 zmin=0 zmax=20000 nh=1 33 oh=0 -1200 dh=75 wbott=0 vsal=14700 \
	traces_in_taper=50 padtraces=0 ray_par_n=20 prop_ang_o=$* prop_ang_d=0.5 > $@

57985plane_wavemig_full.H: # Compare with 5115areal_mig_d16000_full.H:
	Window3d < 5798plane_wavemig_-35.H n4=1 f4=19 > 57985plane_wavemig_a.H
	Window3d < 5798plane_wavemig_-25.H n4=1 f4=19 > 57985plane_wavemig_b.H
	Window3d < 5798plane_wavemig_-15.H n4=1 f4=19 > 57985plane_wavemig_c.H
	Window3d < 5798plane_wavemig_-05.H n4=1 f4=19 > 57985plane_wavemig_d.H
	Window3d < 5798plane_wavemig_05.H n4=1 f4=19 > 57985plane_wavemig_e.H
	Window3d < 5798plane_wavemig_15.H n4=1 f4=19 > 57985plane_wavemig_f.H
	Window3d < 5798plane_wavemig_25.H n4=1 f4=19 > 57985plane_wavemig_g.H
	Math file1=57985plane_wavemig_a.H file2=57985plane_wavemig_b.H file3=57985plane_wavemig_c.H \
	file4=57985plane_wavemig_d.H file5=57985plane_wavemig_e.H file6=57985plane_wavemig_f.H \
	file7=57985plane_wavemig_g.H exp="file1+file2+file3+file4+file5+file6+file7" > $@

5999plane_teste_%.H:
	< ../Tese/4003data_recip_prep.H ./bin/OWmig2d_plane.x sou=../Tese/4004data_recip_source.H \
	vel=../Tese/4000vel.H nrefvel=1 xmin=-4 xmax=8 ymin=0 ymax=0 zmin=0 zmax=1.4 nh=41 oh=-.4 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=50 padtraces=0 ray_par_n=10 prop_ang_o=$* prop_ang_d=2. planeshot= 5999planeshot.H ws= 5999ws_$*.H wr= 5999wr_$*.H > $@

59991plane_teste_full.H:
	Window3d < 5999plane_teste_-70.H n4=1 f4=9 > 59991plane_teste_a.H
	Window3d < 5999plane_teste_-50.H n4=1 f4=9 > 59991plane_teste_b.H
	Window3d < 5999plane_teste_-30.H n4=1 f4=9 > 59991plane_teste_c.H
	Window3d < 5999plane_teste_-10.H n4=1 f4=9 > 59991plane_teste_d.H
	Window3d < 5999plane_teste_10.H n4=1 f4=9 > 59991plane_teste_e.H
	Window3d < 5999plane_teste_30.H n4=1 f4=9 > 59991plane_teste_f.H
	Window3d < 5999plane_teste_50.H n4=1 f4=9 > 59991plane_teste_g.H
	Math file1=59991plane_teste_a.H file2=59991plane_teste_b.H file3=59991plane_teste_c.H \
	file4=59991plane_teste_d.H file5=59991plane_teste_e.H file6=59991plane_teste_f.H \
	file7=59991plane_teste_g.H exp="file1+file2+file3+file4+file5+file6+file7" > $@

59991plane_teste_ang.H:
	< 59991plane_teste_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=1.4 > $@

####################################
#
# 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=41 oh=-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

#################################
#
# Marmousi
#
#

8100marm_shot_recipr.H:
	Create3d < 8001marm_shot_fd.H keyname1=offset keyname2=sx verb=1 > 8001marm_shot_fd3d.H
	Headermath key1=cmp eqn1="sx+0.5*offset" < 8001marm_shot_fd3d.H > junk1.H
	Sort3d < junk1.H nkeys=2 key2=cmp ng2=860 og2=0 dg2=12 key1=offset ng1=68 og1=0 dg1=48 max_size=1000 > junk2.H
	<junk2.H Infill3d >junk3.H
	Interp d2out=24 o2out=0 n2out=134 maxsize=1000 <junk3.H >junk31.H
	Reverse < junk31.H opt=no which=2 | Window3d max2=-2 > 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=360 og2=0 dg2=24 key1=offset ng1=267 og1=-3192 dg1=24 max_size=1000 > junk8.H
	<junk8.H Infill3d > junk9.H
	<junk9.H Stack3d > $@
	 
8101marm_prep_areal_recip.H:
	< 8100marm_shot_recipr.H Pad end1=250 | Transf shotmig=y \
	is2d=y wei=y f_min=6 f_min1=10 f_max1=70 f_max=80 maxsize=1000 > j2.H
	< j2.H Transp reshape=1,3,5 plane=23 | Pad beg1=300 end1=300 > $@

8102marm_source_recip.H:
	 < 8101marm_prep_areal_recip.H Window3d n4=1 | Transp | Cabs | Stack | Smooth rect1=5 > 8102marm_spec.H;
	Math file1=8102marm_spec.H exp="file1*0" > 8102marm_phse.H; Cmplx 8102marm_spec.H 8102marm_phse.H |\
	Transp reshape=1,3,5 plane=23| Transp reshape=1,3,5 plane=12 > $@

#--------------------------
#
# Figures for thesis

8100marm_areal_synthesis.H:
	Wavelet n1=951 d1=0.004 wavelet=ricker2 phase=0 fund=20 order=0.2 tdelay=1.9 domain=time > j2.H;echo o1=-1.9 >> j2.H
	<j2.H Transf f_min=0 f_min1=4 f_max1=42 f_max=48 wei=y is2d=y > j1.H; <j1.H Window3d | Cabs > j2.H
	Math file1=j2.H exp="file1*0">j1.H;Cmplx j2.H j1.H|Transp plane=13 |Transp plane=24 > 11111.H
	echo n4=1 d4=1 o4=0 >> 11111.H
	< ../Tese/3000marmod_velsmoo.H Window3d squeeze=n | Pad beg3=8 end3=0 > 2222.H;echo o3=0 >> 2222.H
	< 2222.H ./bin/OW_arealsht_s.x sou=11111.H vel=../Tese/3000marm_vel_smooth1.tr.H nrefvel=4 xmin=-10392 xmax=19008 ymin=0 \
	ymax=0 zmin=0 zmax=3000 nh=1 oh=0 dh=0 vsal=14700 traces_in_taper=0 depth_wav=2300 shot_in=2222.H \
	synth_op= 8100marm_areal_synth.H synsh=8101marm_Synth_Areal_recip_d$*a.H velout= velout.H nphi=1 dphi=1. ophi=-20 > $@

8101marm_areal_mig:
	< 8101marm_Synth_Areal_recip_da.H Transp plane=13 > receiver.H
	< 8100marm_areal_synthesis.H Window3d squeeze=n min3=-6600 n3=926 | Transp plane=13 > source.H
	../Tese/bin/OWmig2d_areal.x sou=source.H rec=receiver.H  \
	vel=../Tese/3000marm_vel_smooth1.H bkgimg=8101marm_areal_mig.H \
	zmin=0 zmax=3000 nh=1 61 oh=0 -300 jh= 1 migration=1 areal=1 \
	padtraces=50 traces_in_taper=25 node=14 vsal=7000 nrefvel=4 > /dev/null	

#------------------------------------------------
# -30 to 30 every 1
#
8100marm_areal_synth_recip_d%a.H:
	< 8101marm_prep_areal_recip.H ./bin/OW_arealsht.x sou=8102marm_source_recip.H vel=8000marmvel_24m.H nrefvel=4 xmin=-10392 xmax=19008 ymin=0 \
	ymax=0 zmin=0 zmax=3000 nh=1 oh=0 dh=0 vsal=14700 traces_in_taper=0 depth_wav=$* shot_in=8101marm_prep_areal_recip.H \
	synth_op= 8100marm_areal_synth.H synsh=8101marm_Synth_Areal_recip_d$*a.H velout= velout.H nphi=10 dphi=1. ophi=-30 > $@

8100marm_areal_synth_recip_d%b.H:
	< 8101marm_prep_areal_recip.H ./bin/OW_arealsht.x sou=8102marm_source_recip.H vel=8000marmvel_24m.H nrefvel=4 xmin=-10392 xmax=19008 ymin=0 \
	ymax=0 zmin=0 zmax=3000 nh=1 oh=0 dh=0 vsal=14700 traces_in_taper=0 depth_wav=$* shot_in=8101marm_prep_areal_recip.H \
	synth_op= 8100marm_areal_synth.H synsh=8101marm_Synth_Areal_recip_d$*b.H velout= velout.H nphi=10 dphi=1. ophi=-20 > $@

8100marm_areal_synth_recip_d%c.H:
	< 8101marm_prep_areal_recip.H ./bin/OW_arealsht.x sou=8102marm_source_recip.H vel=8000marmvel_24m.H nrefvel=4 xmin=-10392 xmax=19008 ymin=0 \
	ymax=0 zmin=0 zmax=3000 nh=1 oh=0 dh=0 vsal=14700 traces_in_taper=0 depth_wav=$* shot_in=8101marm_prep_areal_recip.H \
	synth_op= 8100marm_areal_synth.H synsh=8101marm_Synth_Areal_recip_d$*c.H velout= velout.H nphi=10 dphi=1. ophi=-10 > $@

8100marm_areal_synth_recip_d%d.H:
	< 8101marm_prep_areal_recip.H ./bin/OW_arealsht.x sou=8102marm_source_recip.H vel=8000marmvel_24m.H nrefvel=4 xmin=-10392 xmax=19008 ymin=0 \
	ymax=0 zmin=0 zmax=3000 nh=1 oh=0 dh=0 vsal=14700 traces_in_taper=0 depth_wav=$* shot_in=8101marm_prep_areal_recip.H \
	synth_op= 8100marm_areal_synth.H synsh=8101marm_Synth_Areal_recip_d$*d.H velout= velout.H nphi=10 dphi=1. ophi=0> $@

8100marm_areal_synth_recip_d%e.H:
	< 8101marm_prep_areal_recip.H ./bin/OW_arealsht.x sou=8102marm_source_recip.H vel=8000marmvel_24m.H nrefvel=4 xmin=-10392 xmax=19008 ymin=0 \
	ymax=0 zmin=0 zmax=3000 nh=1 oh=0 dh=0 vsal=14700 traces_in_taper=0 depth_wav=$* shot_in=8101marm_prep_areal_recip.H \
	synth_op= 8100marm_areal_synth.H synsh=8101marm_Synth_Areal_recip_d$*e.H velout= velout.H nphi=10 dphi=1. ophi=10> $@

8100marm_areal_synth_recip_d%f.H:
	< 8101marm_prep_areal_recip.H ./bin/OW_arealsht.x sou=8102marm_source_recip.H vel=8000marmvel_24m.H nrefvel=4 xmin=-10392 xmax=19008 ymin=0 \
	ymax=0 zmin=0 zmax=3000 nh=1 oh=0 dh=0 vsal=14700 traces_in_taper=0 depth_wav=$* shot_in=8101marm_prep_areal_recip.H \
	synth_op= 8100marm_areal_synth.H synsh=8101marm_Synth_Areal_recip_d$*f.H velout= velout.H nphi=10 dphi=1. ophi=20 > $@

8111marm_areal_synth_recip_d%a.H:
	< 8101marm_prep_areal_recip.H ./bin/OW_arealsht.x sou=8102marm_source_recip.H vel=8000marmvel_24m.H nrefvel=4 xmin=-10392 xmax=19008 ymin=0 \
	ymax=0 zmin=0 zmax=3000 nh=1 oh=0 dh=0 vsal=14700 traces_in_taper=0 depth_wav=$* shot_in=8101marm_prep_areal_recip.H \
	synth_op= 8100marm_areal_synth.H synsh=8111marm_Synth_Areal_recip_d$*a.H velout= velout.H nphi=10 dphi=2. ophi=-50 encode=0 > $@

8111marm_areal_synth_recip_d%b.H:
	< 8101marm_prep_areal_recip.H ./bin/OW_arealsht.x sou=8102marm_source_recip.H vel=8000marmvel_24m.H nrefvel=4 xmin=-10392 xmax=19008 ymin=0 \
	ymax=0 zmin=0 zmax=3000 nh=1 oh=0 dh=0 vsal=14700 traces_in_taper=0 depth_wav=$* shot_in=8101marm_prep_areal_recip.H \
	synth_op= 8100marm_areal_synth.H synsh=8111marm_Synth_Areal_recip_d$*b.H velout= velout.H nphi=10 dphi=2. ophi=-30 encode=0 > $@

8111marm_areal_synth_recip_d%c.H:
	< 8101marm_prep_areal_recip.H ./bin/OW_arealsht.x sou=8102marm_source_recip.H vel=8000marmvel_24m.H nrefvel=4 xmin=-10392 xmax=19008 ymin=0 \
	ymax=0 zmin=0 zmax=3000 nh=1 oh=0 dh=0 vsal=14700 traces_in_taper=0 depth_wav=$* shot_in=8101marm_prep_areal_recip.H \
	synth_op= 8100marm_areal_synth.H synsh=8111marm_Synth_Areal_recip_d$*c.H velout= velout.H nphi=10 dphi=2. ophi=-10 encode=0 > $@

8111marm_areal_synth_recip_d%d.H:
	< 8101marm_prep_areal_recip.H ./bin/OW_arealsht.x sou=8102marm_source_recip.H vel=8000marmvel_24m.H nrefvel=4 xmin=-10392 xmax=19008 ymin=0 \
	ymax=0 zmin=0 zmax=3000 nh=1 oh=0 dh=0 vsal=14700 traces_in_taper=0 depth_wav=$* shot_in=8101marm_prep_areal_recip.H \
	synth_op= 8100marm_areal_synth.H synsh=8111marm_Synth_Areal_recip_d$*d.H velout= velout.H nphi=10 dphi=2. ophi=10 encode=0 > $@

8111marm_areal_synth_recip_d%e.H:
	< 8101marm_prep_areal_recip.H ./bin/OW_arealsht.x sou=8102marm_source_recip.H vel=8000marmvel_24m.H nrefvel=4 xmin=-10392 xmax=19008 ymin=0 \
	ymax=0 zmin=0 zmax=3000 nh=1 oh=0 dh=0 vsal=14700 traces_in_taper=0 depth_wav=$* shot_in=8101marm_prep_areal_recip.H \
	synth_op= 8100marm_areal_synth.H synsh=8111marm_Synth_Areal_recip_d$*e.H velout= velout.H nphi=10 dphi=2. ophi=30 encode=0 > $@

#-------------------------------------------------
#------------------------------------------------
# -50 to 50 every 2
#
8500marm_areal_synth_recip_d%a.H:
	< 8101marm_prep_areal_recip.H ./bin/OW_arealsht.x sou=8102marm_source_recip.H vel=8000marmvel_24m.H nrefvel=8 xmin=-10392 xmax=19008 ymin=0 \
	ymax=0 zmin=0 zmax=3000 nh=1 oh=0 dh=0 vsal=14700 traces_in_taper=0 depth_wav=$* shot_in=8101marm_prep_areal_recip.H \
	synth_op= 8100marm_areal_synth.H synsh=8501marm_Synth_Areal_recip_d$*a.H velout= velout.H nphi=10 dphi=2. ophi=-50 > $@

8500marm_areal_synth_recip_d%b.H:
	< 8101marm_prep_areal_recip.H ./bin/OW_arealsht.x sou=8102marm_source_recip.H vel=8000marmvel_24m.H nrefvel=8 xmin=-10392 xmax=19008 ymin=0 \
	ymax=0 zmin=0 zmax=3000 nh=1 oh=0 dh=0 vsal=14700 traces_in_taper=0 depth_wav=$* shot_in=8101marm_prep_areal_recip.H \
	synth_op= 8100marm_areal_synth.H synsh=8501marm_Synth_Areal_recip_d$*b.H velout= velout.H nphi=10 dphi=2. ophi=-30 > $@

8500marm_areal_synth_recip_d%c.H:
	< 8101marm_prep_areal_recip.H ./bin/OW_arealsht.x sou=8102marm_source_recip.H vel=8000marmvel_24m.H nrefvel=8 xmin=-10392 xmax=19008 ymin=0 \
	ymax=0 zmin=0 zmax=3000 nh=1 oh=0 dh=0 vsal=14700 traces_in_taper=0 depth_wav=$* shot_in=8101marm_prep_areal_recip.H \
	synth_op= 8100marm_areal_synth.H synsh=8501marm_Synth_Areal_recip_d$*c.H velout= velout.H nphi=10 dphi=2. ophi=-10 > $@

8500marm_areal_synth_recip_d%d.H:
	< 8101marm_prep_areal_recip.H ./bin/OW_arealsht.x sou=8102marm_source_recip.H vel=8000marmvel_24m.H nrefvel=8 xmin=-10392 xmax=19008 ymin=0 \
	ymax=0 zmin=0 zmax=3000 nh=1 oh=0 dh=0 vsal=14700 traces_in_taper=0 depth_wav=$* shot_in=8101marm_prep_areal_recip.H \
	synth_op= 8100marm_areal_synth.H synsh=8501marm_Synth_Areal_recip_d$*d.H velout= velout.H nphi=10 dphi=2. ophi=10> $@

8500marm_areal_synth_recip_d%e.H:
	< 8101marm_prep_areal_recip.H ./bin/OW_arealsht.x sou=8102marm_source_recip.H vel=8000marmvel_24m.H nrefvel=8 xmin=-10392 xmax=19008 ymin=0 \
	ymax=0 zmin=0 zmax=3000 nh=1 oh=0 dh=0 vsal=14700 traces_in_taper=0 depth_wav=$* shot_in=8101marm_prep_areal_recip.H \
	synth_op= 8100marm_areal_synth.H synsh=8501marm_Synth_Areal_recip_d$*e.H velout= velout.H nphi=10 dphi=2. ophi=30> $@

#--Migrating.....
85111marm_areal_mig_d2200_v6_dips%.H:
	< 8111marm_Synth_Areal_recip_d2200$*.H Window3d min3=-7200 n3=1093 squeeze=n | \
	Transp plane=13 > 8121marm_Synth_Areal_transp_dips$*.H
	Window3d squeeze=n < 8111marm_areal_synth_recip_d2200$*.H min3=-7200 n3=1093 | \
	Transp plane=13 > 8121marm_areal_synth_ltd_dips$*.H
	< 8121marm_Synth_Areal_transp_dips$*.H ./bin/OWmig2d_areal.x sou=8121marm_areal_synth_ltd_dips$*.H \
	vel=8000marmvel_24m.H nrefvel=6 xmin=-7200 xmax=19008 ymin=0 ymax=0 zmin=0 zmax=3000 nh=21 oh=-240 dh=24 wbott=0 vsal=14700 \
	traces_in_taper=0 velout= velout.H wr=851111Wr_dips$*.H > $@


85051marm_areal_mig_d2200_v8_dips%.H:
	< 8501marm_Synth_Areal_recip_d2200$*.H Window3d min3=-7200 n3=1093 squeeze=n | \
	Transp plane=13 > 8501marm_Synth_Areal_transp_dips$*.H
	Window3d squeeze=n < 8500marm_areal_synth_recip_d2200$*.H min3=-7200 n3=1093 | \
	Transp plane=13 > 8500marm_areal_synth_ltd_dips$*.H
	< 8501marm_Synth_Areal_transp_dips$*.H ./bin/OWmig2d_areal.x sou=8500marm_areal_synth_ltd_dips$*.H \
	vel=8000marmvel_24m.H nrefvel=8 xmin=-7200 xmax=19008 ymin=0 ymax=0 zmin=0 zmax=3000 nh=1 oh=0 dh=24 wbott=0 vsal=14700 \
	traces_in_taper=0 velout= velout.H > $@

85052marm_areal_mig_d2200_v8_dips%.H:
	< 8501marm_Synth_Areal_recip_d2200$*.H Window3d min3=-7200 n3=1093 squeeze=n | \
	Transp plane=13 > 8502marm_Synth_Areal_transp_dips$*.H
	Window3d squeeze=n < 8500marm_areal_synth_recip_d2200$*.H min3=-7200 n3=1093 | \
	Transp plane=13 > 8500marm_areal_synth_ltd_dips$*.H
	< 8502marm_Synth_Areal_transp_dips$*.H ./bin/OWmig2d_areal.x sou=8500marm_areal_synth_ltd_dips$*.H \
	vel=8000marmvel_24m.H nrefvel=8 xmin=-7200 xmax=19008 ymin=0 ymax=0 zmin=0 zmax=3000 nh=21 oh=-240 dh=24 wbott=0 vsal=14700 \
	traces_in_taper=0 velout= velout.H > $@

#-------------------------------------------------


81051marm_areal_mig_d2200_dips%.H:
	< 8101marm_Synth_Areal_recip_d2200$*.H Window3d min3=-7200 n3=1093 squeeze=n | Transp plane=13 > 8101marm_Synth_Areal_transp_dips$*.H
	Window3d squeeze=n < 8100marm_areal_synth_recip_d2200$*.H min3=-7200 n3=1093 | Transp plane=13 > 8100marm_areal_synth_ltd_dips$*.H
	< 8101marm_Synth_Areal_transp_dips$*.H ./bin/OWmig2d_areal.x sou=8100marm_areal_synth_ltd_dips$*.H vel=8000marmvel_24m.H nrefvel=4 xmin=-7200 xmax=19008 ymin=0 ymax=0 zmin=0 zmax=3000 nh=1 oh=0 dh=24 wbott=0 vsal=14700 traces_in_taper=0 velout= velout.H > $@

81051marm_areal_mig_d2200_v12_dips%.H:
	< 8101marm_Synth_Areal_recip_d2200$*.H Window3d min3=-7200 n3=1093 squeeze=n | Transp plane=13 > 8101marm_Synth_Areal_transp_dips$*.H
	Window3d squeeze=n < 8100marm_areal_synth_recip_d2200$*.H min3=-7200 n3=1093 | Transp plane=13 > 8100marm_areal_synth_ltd_dips$*.H
	< 8101marm_Synth_Areal_transp_dips$*.H ./bin/OWmig2d_areal.x sou=8100marm_areal_synth_ltd_dips$*.H vel=8000marmvel_24m.H nrefvel=12 xmin=-7200 xmax=19008 ymin=0 ymax=0 zmin=0 zmax=3000 nh=1 oh=0 dh=24 wbott=0 vsal=14700 traces_in_taper=0 velout= velout.H > $@

8101marm_prep_areal.H:
	< 8001marm_shot_fd.H Pad end1=251| /net/koko/bob/SEP_BOB/bin/Transf shotmig=y \
	is2d=y wei=y f_min=6 f_min1=10 f_max1=70 f_max=80 maxsize=1000 > j2.H
	< j2.H Transp reshape=1,3,5 plane=23 | Pad beg1=300 end1=300 > $@

8102marm_source.H:
	 < 8101marm_prep_areal.H Window3d n4=1 | Transp | Cabs | Stack | Smooth rect1=5 > 8102marm_spec.H;
	Math file1=8102marm_spec.H exp="file1*0" > 8102marm_phse.H; Cmplx 8102marm_spec.H 8102marm_phse.H |\
	Transp reshape=1,3,5 plane=23| Transp reshape=1,3,5 plane=12 > $@

8102marm_spkw.H:
	Spike n1=501 n2=1 n3=1 n4=1 n5=1 o1=-2. o2=0 o3=0 o4=0 o5=0 d1=0.008 d2=0 d3=0 d4=0 d5=0 nsp=1 k1=251 mag=1 | \
	Bandpass flo=8 fhi=30 nplo=5 nphi=5 | /net/koko/bob/SEP_BOB/bin/Transf shotmig=y is2d=y wei=y f_min=2 f_min1=8 \
	f_max1=28 f_max=32 maxsize=1000 > j2.H; < j2.H Transp reshape=1,3,5 plane=23 > $@

8100marm_areal_synth_d%.H:
	< 8101marm_prep_areal.H ./bin/OW_arealsht.x sou=8102marm_source.H vel=8000marmvel_24m.H nrefvel=12 xmin=-7200 xmax=19032 ymin=0 \
	ymax=0 zmin=0 zmax=3000 nh=1 oh=0 dh=0 vsal=14700 traces_in_taper=0 depth_wav=$* shot_in=8101marm_prep_areal.H \
	synth_op= 8100marm_areal_synth.H synsh=8101marm_Synth_Areal_d$*.H velout=velout.H> $@

8100_marm_areal_synth_d%a.H:
	< 8101marm_prep_areal.H ./bin/OW_arealsht.x sou=8102marm_source.H vel=8000marmvel_24m.H nrefvel=12 xmin=-7200 xmax=19032 ymin=0 \
	ymax=0 zmin=0 zmax=3000 nh=1 oh=0 dh=0 vsal=14700 traces_in_taper=0 depth_wav=$* shot_in=8101marm_prep_areal.H \
	synth_op= 8100marm_areal_synth.H synsh=8101marm_Synth_Areal_d$*a.H velout=velout.H nphi=20 ophi=-40 dphi=1 > $@

8100_marm_areal_synth_d%b.H:
	< 8101marm_prep_areal.H ./bin/OW_arealsht.x sou=8102marm_source.H vel=8000marmvel_24m.H nrefvel=12 xmin=-7200 xmax=19032 ymin=0 \
	ymax=0 zmin=0 zmax=3000 nh=1 oh=0 dh=0 vsal=14700 traces_in_taper=0 depth_wav=$* shot_in=8101marm_prep_areal.H \
	synth_op= 8100marm_areal_synth.H synsh=8101marm_Synth_Areal_d$*b.H velout=velout.H nphi=20 ophi=21 dphi=1 > $@

8105marm_areal_mig_d%.H:
	< 8101marm_Synth_Areal_d$*.H Transp plane=13 > 8101marm_Synth_Areal_transp_d$*.H
	Window3d squeeze=n < 8100marm_areal_synth_d$*.H min3=-7200 n3=1093 | Transp plane=13 > 8100marm_areal_synth_ltd_d$*.H
	< 8101marm_Synth_Areal_transp_d$*.H ./bin/OWmig2d_areal.x sou=8100marm_areal_synth_ltd_d$*.H vel=8000marmvel_24m.H nrefvel=12 xmin=-7200 xmax=19008 ymin=0 ymax=0 zmin=0 zmax=3000 nh=1 oh=0 dh=24 wbott=0 vsal=14700 traces_in_taper=0 velout= velout.H > $@

81051marm_areal_mig_dips%.H:
	< 8101marm_Synth_Areal_d2300.H Window3d n4=1 f4=$* squeeze=n | Transp plane=13 > 8101marm_Synth_Areal_transp_dips$*.H
	Window3d squeeze=n < 8100marm_areal_synth_d2300.H n4=1 f4=$* min3=-7200 n3=1093 | Transp plane=13 > 8100marm_areal_synth_ltd_dips$*.H
	< 8101marm_Synth_Areal_transp_dips$*.H ./bin/OWmig2d_areal.x sou=8100marm_areal_synth_ltd_dips$*.H vel=8000marmvel_24m.H nrefvel=12 xmin=-7200 xmax=19008 ymin=0 ymax=0 zmin=0 zmax=3000 nh=1 oh=0 dh=24 wbott=0 vsal=14700 traces_in_taper=0 velout= velout.H > $@

81051marm_areal_mig_d2300_dips%.H:
	< 8101marm_Synth_Areal_d2300.H Window3d n4=1 f4=$* squeeze=n | Transp plane=13 > 8101marm_Synth_Areal_transp_dips$*.H
	Window3d squeeze=n < 8100marm_areal_synth_d2300.H n4=1 f4=$* min3=-7200 n3=1093 | Transp plane=13 > 8100marm_areal_synth_ltd_dips$*.H
	< 8101marm_Synth_Areal_transp_dips$*.H ./bin/OWmig2d_areal.x sou=8100marm_areal_synth_ltd_dips$*.H vel=8000marmvel_24m.H nrefvel=12 xmin=-7200 xmax=19008 ymin=0 ymax=0 zmin=0 zmax=3000 nh=1 oh=0 dh=24 wbott=0 vsal=14700 traces_in_taper=0 velout= velout.H > $@

8106marm_areal_adcig_d%.H:
	< 8105marm_areal_mig_d$*.H Transp reshape=2,4,5 plane=13 |Pad end1=100 beg2=43 end2=44 > j1.H
	< j1.H /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y amax=60 eps=1 maxmem=1000 |Window3d max1=3000 > $@

8115marm_mig_areal_brad.H:
	< 8100marm_areal_synth_d2300.H Window3d squeeze=n n3=1093 min3=-7200 n4=4 | Transp plane=13 > 8110marm_areal_prep.H
	echo n4=4 d4=0 o4=5904 >> 8110marm_areal_prep.H
	< 8101marm_Synth_Areal_d2300.H Window3d squeeze=n n4=4 | Transp plane=13 > 8111marm_synthareal.H
	echo n4=4 d4=0 o4=5904 >> 8111marm_synthareal.H
	@date
	nohup /net/koko/claudio/Research/Off2ang/ShotMigNew \
	D=8111marm_synthareal.H D_d=8110marm_areal_prep.H V=8000marmvel_24m_brad.H R=$@ \
	nx_shots=4 ox_shots=5904 par=marm_areal_go_sis.P taper=n taper_pcent= 2 marine=y minx=-7200 maxx=19008 #restart=y

_d%.H:
	./bin/OW_shift.x dt=-1 < 8100marm_areal_synth_d$*.H > $@

8103marm_areal_synth_d%.H:
	< 8101marm_Synth_Areal_d$*.H Window3d | Transp reshape=1,3,4 plane=12 > 8101temp.H
	./bin/OW_shift.x dt=1.986 < 8101temp.H > $@

#
# End Areal
#
###################################################################

##################################################################
#
# Preparing Marmousi
#
# Pad < marmvel.H beg2=400 end2=400 extend=1 > marmvel_pad.H; echo n3=1 o3=0 d3=1 >> marmvel_pad.H; 
# 
# Math file1=marmvel_pad.H exp="file1*2" > marmvel_pad2x.H

8000marmvel_1m.H:
	Interp d1out=1 o2out=-10000 d2out=1 n2out=13001 n1out=3001 maxsize=1000 < marmvel_pad.H > $@

8000marmvel_2m.H:
	Interp d1out=2 o2out=-20000 d2out=2 n2out=30001 n1out=1501 maxsize=1000 < marmvel_pad.H > $@

8000marmvel_4m.H:
	Interp d1out=4 o2out=-20000 d2out=4 n2out=15001 n1out=751 maxsize=1000 < marmvel_pad.H > $@
###############################
#
# Modeling Marmousi
#
8000marmden_4m.H:
	Gpow < 8000marmvel_4m.H gpow=0.23 | Scale rscale=0.31 > $@
	< ../Tese/3000marm_vel_smooth1.H Window3d |Transp|Gpow gpow=-1|Pad end2=500 extend=1| \
	Interp d1out=4 d2out=4 type=1 maxsize=1000 > 3000marm_vel_smooth1.H

8001marm_shot_fd%:
	Fdmod intag=8000marmvel_4m.H  dfile= 8000marmden_4m.H no_stdout=1 oxs=$* dxs=24 nxs=1 25 tmax=4 ex=1650 bx=1650 jt=10 hsfile=temp$*.H > 111.H # /dev/null

8001marm_shots_fd.H:
	Cat temp0000.H temp0600.H temp1200.H temp1800.H temp2400.H temp3000.H temp3600.H temp4200.H temp4800.H temp5400.H \
	temp6000.H temp6600.H temp7200.H temp7800.H temp8400.H axis=3 | \
	Bandpass flo=4 | Window3d j2=6 | Interp d1out=0.004 |Mute vmute=1500 tmute=0.05  > $@

8000marmvel_4m.smooth%.H:
	Gpow gpow=-1 < 8000marmvel_24m.H | Smooth rect1=$* rect2=$* | Gpow gpow=-1 | Interp d1out=4 d2out=4 maxsize=100 | Pad end2=3000 extend=1 > $@

8000marmden_4m_log.H:
	Math file1=8000marmden_4m.H exp="@LOG(file1)" > $@

80013marmsmoo_shot_fd%:
	Fdmod intag=8000marmvel_4m.smooth3.H dfile=8000marmden_4m_log.H no_stdout=1 oxs=$* dxs=24 nxs=25 tmax=4 ex=1650 bx=1650 jt=10 hsfile=temp$*.H > /dev/null

80019marmsmoo_shot_fd%:
	Fdmod intag=8000marmvel_4m.smooth9.H dfile=8000marmden_4m_log.H no_stdout=1 oxs=$* dxs=24 nxs=25 tmax=4 ex=1650 bx=1650 jt=10 hsfile=temp$*.H > /dev/null

80019marm_shots_fd.H:
	Cat temp0000.H temp0600.H temp1200.H temp1800.H temp2400.H temp3000.H temp3600.H temp4200.H temp4800.H temp5400.H \
	temp6000.H temp6600.H temp7200.H temp7800.H temp8400.H axis=3 | \
	Bandpass flo=4 | Window3d j2=6 | Interp d1out=0.004 |Mute vmute=1500 tmute=0.05  > $@

###############################

8002marm_shot_fd.H:
	< 8001marm_shot_fd.H Mute vmute=1500 tmute=0.1 emute=1 eramp=0 tramp=0 > xxx.H
	< 8001marm_shot_fd.H Mute vmute=0 tmute=1 tramp=0 > yyy.H
	Math file1=xxx.H file2=yyy.H exp=file1+file2 > $@

8002marm_prep_fd.H:
	< 8002marm_shot_fd.H /net/koko/bob/SEP_BOB/bin/Transf \
	shotmig=y is2d=y wei=y f_min=6 f_min1=10 f_max1=70 f_max=80 maxsize=1000 > j2.H
	< j2.H Transp reshape=1,3,5 plane=23 | Pad beg1=300 end1=300 > $@

8001marm_prep_fd.H:
	< 8001marm_shot_fd.H /net/koko/bob/SEP_BOB/bin/Transf \
	shotmig=y is2d=y wei=y f_min=6 f_min1=10 f_max1=70 f_max=80 maxsize=1000 > j2.H
	< j2.H Transp reshape=1,3,5 plane=23 | Pad beg1=300 end1=300 > $@

8102marm_source_brad.H:
	 < 8001marm_prep_fd.H Window3d n4=1 | Transp | Cabs | Stack | Smooth rect1=5 > 8102marm_spec.H;
	Math file1=8102marm_spec.H exp="file1*0" > 8102marm_phse.H; Cmplx 8102marm_spec.H 8102marm_phse.H |\
	Transp reshape=1,3,5 plane=23| Transp reshape=1,3,5 plane=12 > $@

8004marm_mig_brad_fdmod.H:
	@date
	/net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX86_64/ShotMig \
	D=8001marm_prep_fd.H D_d=8102marm_source_brad.H V=8000marmvel_24m_brad.H R=$@ \
	pverb=2  nx_shots=360 ox_shots=0 \
	par=marm_go_sis.P taper=y taper_pcent=2 marine=y minx=-8000 maxx=20000 & #restart=y

80042marm_mig_brad_fdmod.H:
	@date
	nohup /net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX86_64/ShotMig \
	D=8002marm_prep_fd.H D_d=8102marm_source_brad.H V=8000marmvel_24m_brad.H R=$@ \
	pverb=2 nx_shots=360 ox_shots=0 \
	par=marm_go_sis.P taper=y taper_pcent=2 marine=y minx=-8000 maxx=20000 & #restart=y

8005marm_adcig_brad.H:
	< 8004marm_mig_brad_fdmod.H Transp reshape=2,4,5 plane=13 |Pad end1=100 beg2=43 end2=44 > j1.H
	< j1.H /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y amax=60 eps=1 maxmem=1000 |Window3d max1=3000 > $@

8002marm_source_fd001a.H:
	 < j2.H Window3d n4=1 | Transp | Cabs | Stack | Smooth rect1=5 > 8002marm_spec_fd001a.H;
	Math file1=8002marm_spec_fd001a.H exp="file1*0" > 8002marm_phse_fd001a.H; Cmplx 8002marm_spec_fd001a.H \
	8002marm_phse_fd001a.H | Transp reshape=1,3,5 plane=23| Transp reshape=1,3,5 plane=12 > $@

8004marm_mig_brad_fd.H:
	@date
	/net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX86_64/ShotMig \
	D=8001marm_prep_fd001a.H D_d=8102marm_source.H V=8000marmvel_24m_brad.H R=$@ \
	pverb=2  nx_shots=360 ox_shots=0 \
	par=marm_go_sis.P taper=y taper_pcent=2 marine=y minx=-8000 maxx=20000 & #restart=y

8001marm_prep_fd001.H:
	< temp10.H Bandpass flo=12 fhi=90 | Window3d j2=12 | /net/koko/bob/SEP_BOB/bin/Transf shotmig=y \
	is2d=y wei=y f_min=6 f_min1=12 f_max1=80 f_max=90 maxsize=1000 > j2.H
	< j2.H Transp reshape=1,3,5 plane=23 | Pad beg1=300 end1=300 > $@

8002marm_source_fd001.H:
	 < j2.H Window3d n4=1 | Transp | Cabs | Stack | Smooth rect1=5 > 8002marm_spec_fd001.H;
	Math file1=8002marm_spec_fd001.H exp="file1*0" > 8002marm_phse_fd001.H; Cmplx 8002marm_spec_fd001.H \
	8002marm_phse_fd001.H | Transp reshape=1,3,5 plane=23| Transp reshape=1,3,5 plane=12 > $@

8004marm_mig_brad_fd001.H:
	@date
	/net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX86_64/ShotMig \
	D=8001marm_prep_fd001.H D_d=8002marm_source_fd001.H V=8000marmvel_24m_brad.H R=$@ \
	pverb=2  nx_shots=1 ox_shots=0 \
	par=marm_go_sis.P taper=y taper_pcent=2 marine=y minx=-8000 maxx=8000 & #restart=y

8000marmvel_24m.H:
	Interp d1out=8 o2out=-20000 d2out=24 n2out=2001 maxsize=1000 < marmvel_pad.H > $@

8000marmvel_25m.H:
	Interp d1out=8 o2out=-10000 d2out=25 n2out=1801 < marmvel_pad2x.H > $@

8004marm_mig.H:
	<8003marm_fd_prep.H Window3d squeeze=n j4=20 | ./bin/OWmig2d_par.x sou=8002marm_source.H vel=8000marmvel_24m.H nrefvel=8 xmin=-8000 xmax=22000 ymin=0 ymax=0 zmin=0 zmax=3000 nh=1 oh=0 dh=0 wbott=0 vsal=000 wr= WS.H traces_in_taper=30> $@

9000marmvel_25m.H:
	Interp d1out=8 o2out=-10000 d2out=25 n2out=1201 < marmvel.H > $@

9001marm_prep.H:
	< marmdata.H Scale scale_to=1. | /net/koko/bob/SEP_BOB/bin/Transf shotmig=y is2d=y wei=y \
	f_min=1 f_min1=4 f_max1=50 f_max=60 maxsize=1000 > j2.H
	< j2.H Transp reshape=1,3,5 plane=23 | Pad beg1=300 end1=300 > $@; echo d1=25 o1=-10075 >> $@

9002marm_source.H:
	 < marmdata.H /net/koko/bob/SEP_BOB/bin/Transf shotmig=y is2d=y wei=y f_min=1 f_min1=3 f_max1=50 f_max=60 maxsize=1000 \
	| Window3d | Transp plane=12 | Cabs | Stack | Transp | Stack | Smooth rect1=15 > 9002marm_spec.H;
	Math file1=9002marm_spec.H exp="file1*0" > 9002marm_phse.H; Cmplx 9002marm_spec.H 9002marm_phse.H | Transp reshape=1,3,5 plane=23| \
	Transp reshape=1,3,5 plane=12 > $@

9003marm_mig.H:
	<8101marm_prep_areal_recip.H Window3d squeeze=n j4=40 | ./bin/OWmig2d_par.x sou=8102marm_source_recip.H vel=8000marmvel_24m.H nrefvel=1 xmin=-11000 xmax=19100 ymin=0 ymax=0 zmin=0 zmax=500 nh=1 oh=0 dh=0 wbott=0 vsal=6000 wr= WS.H traces_in_taper=30> $@

9004marm_mig_brad.H:
	@date
	nohup Parallel command=/net/koko/claudio/Research/Off2ang/ShotMigNew \
	files="9001marm_prep.H,9002marm_source.H,9000marmvel_25m_brad.H,$@" tags="D,D_d,V,R" \
	usage="INPUT,INPUT,INPUT,OUTPUT" \
	nblock=240 axis=4,4,4,6 pverb=2 file_type="DISTRIBUTE,COPY,COPY,COPY" mach_mfile=machinefile_200 \
	nx_shots=240 ox_shots=3000 par=marm_go_sis.P taper=y taper_pcent=2 marine=y minx=-8000 maxx=18000 #restart=y

#03sig_mig.H:
#	< ../OneWay/xxx.H ./bin/OWmig2d1.x sou=../RTM/sig_sou_mig_mad.H vel=00sig_velmig_res.H nrefvel=4 xmin=-2000 xmax=95000 ymin=0 ymax=0 zmin=0 #zmax=20000 nh=1 oh=0 dh=0 wbott=0 > $@

#03sig_mig.H:
#	< ../OneWay/xxx.H ./bin/OWmig2d2.x sou=../RTM/sig_sou_mig_mad.H vel=00sig_velmig_res.H nrefvel=4 xmin=4000 xmax=95000 ymin=0 ymax=0 #zmin=0 zmax=20000 nh=1 oh=0 dh=0 wbott=0 > $@

03sig_mig.H:
	< ../RTM/sig_shot_prep_mad.H Window3d squeeze=n n4=1 | ./bin/OWmig2d_par.x sou=../RTM/sig_sou_mig_mad.H vel=00sig_velmig_res.H \
	nrefvel=4 xmin=-56250 xmax=105705 wr=wr.H \
	ymin=0 ymax=0 zmin=0 zmax=8000 nh=1 oh=0 dh=0 wbott=0 > $@

#03sig_mig.H:
#	< 01sig_prep.H ./bin/OWmig2d_par.x sou=02sig_source.H vel=00sig_velmig_res.H \
#	nrefvel=4 xmin=10000 xmax=128000 ymin=0 ymax=0 zmin=0 zmax=20000 nh=1 oh=0 dh=0 wbott=0 > $@

zprep.H:
	< ../OneWay/nfs-splitspread-shots.H Window3d n3=1 min3=40000 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 > $@


031sig_mig.H:
	< zprep.H ./bin/OWmig2d_par.x sou=02sig_sou1.H vel=1000velmig_my.H nrefvel=4 xmin=-32000 xmax=128000 ymin=0 ymax=0 zmin=0 zmax=16000 nh=1 oh=0 dh=0 wbott=0 vsal=14700 G1= G1.H> $@

031sig_mig01.H:
	< 01sig_prep.H Window3d squeeze=n n4=50 f4=0 | ./bin/OWmig2d_par.x sou=02sig_sou1.H vel=1000velmig_my.H nrefvel=4 xmin=-32000 xmax=128000 ymin=0 ymax=0 zmin=0 zmax=20000 nh=1 oh=0 dh=0 wbott=0 vsal=14700 G1= G1.H> $@

031sig_mig02.H:
	< 01sig_prep.H Window3d squeeze=n n4=50 f4=50 | ./bin/OWmig2d_par.x sou=02sig_sou1.H vel=1000velmig_my.H nrefvel=4 xmin=-32000 xmax=128000 ymin=0 ymax=0 zmin=0 zmax=20000 nh=1 oh=0 dh=0 wbott=0 vsal=14700 G1= G1.H> $@

031sig_mig03.H:
	< 01sig_prep.H Window3d squeeze=n n4=50 f4=100 | ./bin/OWmig2d_par.x sou=02sig_sou1.H vel=1000velmig_my.H nrefvel=4 xmin=-32000 xmax=128000 ymin=0 ymax=0 zmin=0 zmax=20000 nh=1 oh=0 dh=0 wbott=0 vsal=14700 G1= G1.H> $@

031sig_mig04.H:
	< 01sig_prep.H Window3d squeeze=n n4=50 f4=150 | ./bin/OWmig2d_par.x sou=02sig_sou1.H vel=1000velmig_my.H nrefvel=4 xmin=-32000 xmax=128000 ymin=0 ymax=0 zmin=0 zmax=20000 nh=1 oh=0 dh=0 wbott=0 vsal=14700 G1= G1.H> $@

031sig_mig05.H:
	< 01sig_prep.H Window3d squeeze=n n4=50 f4=200 | ./bin/OWmig2d_par.x sou=02sig_sou1.H vel=1000velmig_my.H nrefvel=4 xmin=-32000 xmax=128000 ymin=0 ymax=0 zmin=0 zmax=20000 nh=1 oh=0 dh=0 wbott=0 vsal=14700 G1= G1.H> $@

031sig_mig06.H:
	< 01sig_prep.H Window3d squeeze=n n4=50 f4=250 | ./bin/OWmig2d_par.x sou=02sig_sou1.H vel=1000velmig_my.H nrefvel=4 xmin=-32000 xmax=128000 ymin=0 ymax=0 zmin=0 zmax=20000 nh=1 oh=0 dh=0 wbott=0 vsal=14700 G1= G1.H> $@

031sig_mig07.H:
	< 01sig_prep.H Window3d squeeze=n n4=50 f4=300 | ./bin/OWmig2d_par.x sou=02sig_sou1.H vel=1000velmig_my.H nrefvel=4 xmin=-32000 xmax=128000 ymin=0 ymax=0 zmin=0 zmax=20000 nh=1 oh=0 dh=0 wbott=0 vsal=14700 G1= G1.H> $@

031sig_mig08.H:
	< 01sig_prep.H Window3d squeeze=n n4=50 f4=350 | ./bin/OWmig2d_par.x sou=02sig_sou1.H vel=1000velmig_my.H nrefvel=4 xmin=-32000 xmax=128000 ymin=0 ymax=0 zmin=0 zmax=20000 nh=1 oh=0 dh=0 wbott=0 vsal=14700 G1= G1.H> $@

############################################################
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

04sig_mig_brad1.H:
#	< 01sig_prep.H Window3d squeeze=n min4=12000 n4=400 | Pad beg1=100 end1=100 > 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/claudio/research.ln/Others/ShotMig/ShotMig \
	D=../1Wpar_phenc/01sig_prep_select.H D_d=../1Wpar_phenc/02sig_sou1w256.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=-38000 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&

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

04sig2a_mig_brad.H:
	< ../OneWay/sigsb2a.rec.freq.H Pad beg1=300 end1=300| Window3d squeeze=n min4=12000 n4=400 > 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/claudio/research.ln/Others/ShotMig/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=20000 max1=55000 min5=6000 max5=20000 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=wem.sigsb2a.migvwrg.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 > $@

04sig_mig_brad_cut.H:
	Window3d max3=16000 squeeze=n <1000velmig_transp.H > 1000velmig_transp_z16000.H 
	@date
	nohup /net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX86_64/ShotMig \
	D=01sig_prep_cut.H D_d=02sig_sou1_cut.H V=1000velmig_transp_z16000.H R=$@ \
	pverb=2  nx_shots=127 ox_shots=10925 watercolumn=5500 \
	par=go_sis.P taper=y taper_pcent=2 marine=y minx=-30100 maxx=72000 & #restart=y

#####
# Sigsbee2a
#
sigsb2a.shots.time.H:
	 Window3d < sigsb2a.rec.time.H | Transp plane=12 | Bandpass flo=2 fhi=60 > x.H ; < x.H Interp d1out=0.008 > $@

9900sigsb2a.shots.recip.H:
#	Create3d < sigsb2a.shots.time.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=2400 og2=10025 dg2=37.5 key1=offset ng1=88 og1=0 dg1=300 max_size=1000 > junk2.H
#	<junk2.H Infill3d >junk3.H;
##	 <junk3.H Stack >junk30.H;
#	Interp d2out=75 o2out=0 n2out=348 maxsize=1000 <junk3.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=1001 og2=10025 dg2=75 key1=offset ng1=697 og1=-26025 dg1=75 max_size=1000 > junk8.H
#	<junk8.H Infill3d > junk9.H
#	<junk9.H Stack3d > $@
	Transp plane=23 < sigsb2a.shots.time.H | Interp d2out=75 >x.H ; <x.H Transp plane=23 >junk0.H
	Create3d < junk0.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=2500 og2=10025 dg2=37.5 key1=offset ng1=177 og1=0 dg1=150 max_size=1000 > junk2.H
	<junk2.H Infill3d >junk3.H;
#	 <junk3.H Stack >junk30.H;
	Interp d2out=75 o2out=0 n2out=348 maxsize=1000 <junk3.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=1000 og2=10025 dg2=75 key1=offset ng1=697 og1=-26025 dg1=75 max_size=1000 > junk8.H
	<junk8.H Infill3d > junk9.H
	<junk9.H Stack3d > junk10.H
	Window3d <junk10.H min3=10025 j3=2 > $@

#	Window3d <9900sigsb2a.shots.recip.H n2=695 min3=10925 > junk1_new.H
#	Window3d <9900sigsb2a.shots.recip.H n2=695 min3=10925 > junk1_new.H
#	Math file1=junk1_new.H file2=../OneWay/nfs-splitspread-shots.H exp="file1-file2" > junk4_new.H
#	time ~claudio/research.ln/ResMult/Bin/Ultimatch-new-all.x < junk1_new.H \
#		primaries_in=../OneWay/nfs-splitspread-shots.H multiples_in=junk4_new.H mu=0.001 \

9901mult_est.H:
	Window3d < 9900sigsb2a.shots.recip.H n2=695 n3=1 min3=35925 > junk1_new.H
	Window3d <../OneWay/nfs-splitspread-shots.H n3=1 min3=35925 > junk2_new.H
	Math file1=junk2_new.H file2=junk1_new.H exp="file1-file2" > junk4_new.H
	time ~claudio/research.ln/ResMult/Bin/Ultimatch-new-all.x < junk2_new.H \
		primaries_in=junk1_new.H multiples_in=junk4_new.H mu=0.00 \
		primaries_out= 9902prim_est.H \
		mode=2 nsmode=2 wmode=1 nfilt=21,21 npatch=31,31 \
		op=1 niter=30 2D=1 eps=15. > $@

9901sigsb2a_prep.H:
	< 9900sigsb2a.shots.recip.H Window3d max1=8 min3=10925 n2=695 max3=72000 | /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 > $@

9902sigsb2a_source.H:
	Window3d n4=1 < 9901sigsb2a_prep.H | Transp | Cabs | Stack | Smooth rect1=19 tridiag=1 > 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&

#  max4=72000
#	< 9901sigsb2a_prep.H Window3d squeeze=n min4=12000 > z02sig_prep.H

9903sigsb2a_mig.H:
	@date
	nohup /net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX86_64/ShotMig \
	D=z02sig_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=-46000 maxx=140000 & #restart=y

1000velmig_transp_slow5.H:
	Math file1=1000velmig_transp_ltdz.H exp="file1*0.95" > $@

9903sigsb2a_mig_slow.H:
	< 9901sigsb2a_prep.H Window3d squeeze=n min4=12000 | Pad beg1=200 end1=200 > z02sig_prep.H
	@date
	nohup /net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX86_64/ShotMig \
	D=z02sig_prep.H D_d=02sig_sou1.H V=../Tese/wem.sigsb2a.migvwrg.H R=$@ \
	pverb=2  nx_shots=400 ox_shots=12000 \
	par=go_sis.P taper=y taper_pcent=2 marine=y minx=-46000 maxx=140000 & #restart=y

#	< 9901sigsb2a_prep.H Window3d squeeze=n min4=12000 | Pad beg1=200 end1=200 > z02sig_prep.H

9903sigsb2a_mig_slow5.H:
	@date
	nohup /net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX86_64/ShotMig \
	D=z02sig_prep.H D_d=02sig_sou1.H V=1000velmig_transp_slow5.H R=$@ \
	pverb=2  nx_shots=200 ox_shots=11975 \
	par=go_sis.P taper=y taper_pcent=2 marine=y minx=-46000 maxx=100000 & #restart=y

9903sigsb2a_mig_slow_cluster.H:
	< 9901sigsb2a_prep.H Window3d squeeze=n min4=12000 | Pad beg1=200 end1=200 > z02sig_prep.H
	@date
	nohup Parallel command=/net/koko/claudio/Research/Off2ang/ShotMigNew \
	files="z02sig_prep.H,02sig_sou1.H,../Tese/wem.sigsb2a.migvwrg.H,$@" tags="D,D_d,V,R" \
	usage="INPUT,INPUT,INPUT,OUTPUT" \
	nblock=400 axis=4,4,4,6 pverb=2 file_type="DISTRIBUTE,COPY,COPY,COPY" mach_mfile=machinefile_200 \
	nx_shots=400 ox_shots=12000 \
	par=go_sis.P taper=y taper_pcent=2 marine=y minx=-46000 maxx=140000 & #restart=y

9903sigsb2a_ang_slow.H:
	< 9903sigsb2a_mig_slow.H Window3d min1=12000 max1=55000 squeeze=n | Transp reshape=2,4,5 plane=13 | Pad beg2=83 end2=84 end1=200 > x.H
	< x.H /net/koko/bob/SEP_BOB/bin/OFF2ANG amax=40 eps=20 maxmem=1000 verb=y  > junk1.H; <junk1.H Window3d max1=24000 > $@


