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=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:
#	Window3d < ../RTM/sig_shot_prep_mad.H n4=1 min4=50000 squeeze=n > xxx.H
#	./bin/OWmig2d1_par.x < xxx.H sou=../RTM/sig_sou_mig_mad.H vel=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:
	Window3d < ../RTM/sig_shot_prep_mad.H n4=5 squeeze=n > xxx.H
	./bin/OWmig2d2.x < xxx.H sou=../RTM/sig_sou_mig_mad.H vel=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:
#	Window3d < ../RTM/sig_shot_prep_mad.H n4=1 min4=50000 squeeze=n > xxx.H
#	./bin/OWmig2d2_par.x < xxx.H sou=../RTM/sig_sou_mig_mad.H vel=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=yyy.H vel=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=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 >$@

##############################################################################################
#
# Test impulse response
#
999spk_inpt.H:
	Spike n1=501 d1=.004 o1=0 n2=1001 d2=10 o2=-5000 n3=1 d3=10 o3=0 nsp=1 mag=10 k1=101 k3=1 k2=501 | \
	Bandpass flo=4 fhi=30 > 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=30 f_max=40 maxsize=1000 > j2.H
	 < j2.H Transp reshape=1,3,5,6,7 plane=23 | Pad beg1=300 end1=300 > $@ 

999source.H:
	Spike n1=501 d1=.004 o1=0 n2=1 d2=10 o2=0 n3=1 d3=10 o3=0 nsp=1 mag=10 k1=251 k3=1 k2=1 | \
	Bandpass flo=4 fhi=30 > j1.H;echo o1=-1. >> 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=30 f_max=40 maxsize=1000 > j2.H
	 < j2.H Transp reshape=1,3,5,6,7 plane=23 > $@ 

999vcte.H:
	Vel n1=3001 d1=4 o1=0 n2=6001 d2=10 o2=-30000 n3=1 d3=1 o3=0 vc=2000 > $@

999vgrad.H:
	Vel n1=3001 d1=4 o1=0 n2=901 d2=20 o2=-9000 n3=1 d3=1 o3=0 vc=2000 vr1=2000 smooth=.5 > $@

999mig_cte.H:
	< 999spk_inpt.H ./bin/OWmig2d1.x sou=999source.H vel=999vcte.H nrefvel=1 xmin=-12000 xmax=12000 \
	ymin=0 ymax=0 zmin=0 zmax=600 nh=1 oh=0 dh=0 vsal=14700 >$@

999imp_vcte.H:

##############################################################################################
000data.H:
	Window3d squeeze=n < 02sig_sou.H n3=1 f3=120 >z.H 
	Window3d squeeze=n < 01sig_prep.H n4=25 n3=1 f3=120 | ./bin/greenxs4.x sou=z.H vel=00sig_velmig_my.H nrefvel=4 \
	xmin=-110000 xmax=128000 ymin=0 ymax=0 zmin=0 zmax=12000 >$@

#	Window3d squeeze=n < 01sig_prep.H n4=25 n3=1 f3=120 | ./bin/greenxs4.x sou=02sig_sou.H #vel=00sig_velmig_my.H nrefvel=4 \



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=32 f_max=40 maxsize=1000 > j2.H
	< j2.H Transp reshape=1,3,5 plane=23 | Pad beg1=200 end1=200 > $@

02sig_sou.H:
	Window3d < 01sig_prep.H n4=1 | Transp | Cabs | Stack | Smooth rect1=7 >  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 | Cabs |Graph|Tube&

00sig_velmig.H:
	Pad < sigsbee2b_strati_2pad.H beg2=5500 end2=5500 extend=1 > $@; echo n3=1 o3=0 d3=1

00sig_velmig_res.H:
	Window3d j1=2 j2=6 squeeze=n < 00sig_velmig.H | Pad extend=1 beg2=2000 end2=6000 | Transp reshape=1,3 plane=12 > $@

00sig_velmig_my.H:
	Window3d < 00sig_velmig_res.H | Transp > $@

03sig_mig_shot01.H:
	Window3d squeeze=n < 01sig_prep.H n4=1 min4=47000 | ./bin/OWmig2d1.x sou=02sig_sou.H vel=00sig_velmig_my.H nrefvel=4 xmin=-10000 xmax=110000 \
	ymin=0 ymax=0 zmin=0 zmax=20000 nh=1 oh=0 dh=0 vsal=14700 >$@

03sig_mig01.H:
	Window3d squeeze=n < 01sig_prep.H n4=25| ./bin/OWmig2d1.x sou=02sig_sou.H vel=00sig_velmig_res.H nrefvel=4 xmin=-56250 xmax=110000 \
	ymin=0 ymax=0 zmin=0 zmax=20000 nh=17 oh=-1200 dh=150 >$@

03sig_mig02.H:
	Window3d squeeze=n < 01sig_prep.H n4=25 f4=25| ./bin/OWmig2d1.x sou=02sig_sou.H vel=00sig_velmig_res.H nrefvel=4 xmin=-56250 xmax=110000 \
	ymin=0 ymax=0 zmin=0 zmax=20000 nh=17 oh=-1200 dh=150 >$@

03sig_mig03.H:
	Window3d squeeze=n < 01sig_prep.H n4=25 f4=50| ./bin/OWmig2d1.x sou=02sig_sou.H vel=00sig_velmig_res.H nrefvel=4 xmin=-56250 xmax=110000 \
	ymin=0 ymax=0 zmin=0 zmax=20000 nh=17 oh=-1200 dh=150 >$@

03sig_mig04.H:
	Window3d squeeze=n < 01sig_prep.H n4=25 f4=75| ./bin/OWmig2d1.x sou=02sig_sou.H vel=00sig_velmig_res.H nrefvel=4 xmin=-56250 xmax=110000 \
	ymin=0 ymax=0 zmin=0 zmax=20000 nh=17 oh=-1200 dh=150 >$@

03sig_mig05.H:
	Window3d squeeze=n < 01sig_prep.H n4=25 f4=100| ./bin/OWmig2d1.x sou=02sig_sou.H vel=00sig_velmig_res.H nrefvel=4 xmin=-56250 xmax=110000 \
	ymin=0 ymax=0 zmin=0 zmax=20000 nh=17 oh=-1200 dh=150 >$@

03sig_mig06.H:
	Window3d squeeze=n < 01sig_prep.H n4=25 f4=125| ./bin/OWmig2d1.x sou=02sig_sou.H vel=00sig_velmig_res.H nrefvel=4 xmin=-56250 xmax=110000 \
	ymin=0 ymax=0 zmin=0 zmax=20000 nh=17 oh=-1200 dh=150 >$@

03sig_mig01a.H:
	Window3d squeeze=n < 01sig_prep.H n4=25| ./bin/OWmig2d1.x sou=02sig_sou.H vel=00sig_velmig_res.H nrefvel=8 xmin=-56250 xmax=110000 \
	ymin=0 ymax=0 zmin=0 zmax=20000 nh=17 oh=-1200 dh=150 >$@

03sig_mig02a.H:
	Window3d squeeze=n < 01sig_prep.H n4=25 f4=25| ./bin/OWmig2d1.x sou=02sig_sou.H vel=00sig_velmig_res.H nrefvel=8 xmin=-56250 xmax=110000 \
	ymin=0 ymax=0 zmin=0 zmax=20000 nh=17 oh=-1200 dh=150 >$@

03sig_mig03a.H:
	Window3d squeeze=n < 01sig_prep.H n4=25 f4=50| ./bin/OWmig2d1.x sou=02sig_sou.H vel=00sig_velmig_res.H nrefvel=8 xmin=-56250 xmax=110000 \
	ymin=0 ymax=0 zmin=0 zmax=20000 nh=17 oh=-1200 dh=150 >$@

03sig_mig04a.H:
	Window3d squeeze=n < 01sig_prep.H n4=25 f4=75| ./bin/OWmig2d1.x sou=02sig_sou.H vel=00sig_velmig_res.H nrefvel=8 xmin=-56250 xmax=110000 \
	ymin=0 ymax=0 zmin=0 zmax=20000 nh=17 oh=-1200 dh=150 >$@

03sig_mig05a.H:
	Window3d squeeze=n < 01sig_prep.H n4=25 f4=100| ./bin/OWmig2d1.x sou=02sig_sou.H vel=00sig_velmig_res.H nrefvel=8 xmin=-56250 xmax=110000 \
	ymin=0 ymax=0 zmin=0 zmax=20000 nh=17 oh=-1200 dh=150 >$@

03sig_mig06a.H:
	Window3d squeeze=n < 01sig_prep.H n4=25 f4=125| ./bin/OWmig2d1.x sou=02sig_sou.H vel=00sig_velmig_res.H nrefvel=8 xmin=-56250 xmax=110000 \
	ymin=0 ymax=0 zmin=0 zmax=20000 nh=17 oh=-1200 dh=150 >$@

03sig_mig_tota.H:
	Math file1=03sig_mig02a.H file2=03sig_mig03a.H file3=03sig_mig04a.H file4=03sig_mig05a.H file5=03sig_mig06a.H \
	exp="file1+file2+file3+file4+file5" > $@
######################################################################
#
03sig_mig01b.H:
	Window3d squeeze=n < 01sig_prep.H n4=25 | ./bin/OWmig2d1.x sou=02sig_sou.H vel=00sig_velmig_my.H nrefvel=4 xmin=-56250 xmax=110000 \
	ymin=0 ymax=0 zmin=0 zmax=20000 nh=17 oh=-1200 dh=150 vsal=14700>$@

03sig_mig02b.H:
	Window3d squeeze=n < 01sig_prep.H n4=25 f4=25 | ./bin/OWmig2d1.x sou=02sig_sou.H vel=00sig_velmig_my.H nrefvel=4 xmin=-56250 xmax=110000 \
	ymin=0 ymax=0 zmin=0 zmax=20000 nh=17 oh=-1200 dh=150 vsal=14700>$@

03sig_mig03b.H:
	Window3d squeeze=n < 01sig_prep.H n4=25 f4=50| ./bin/OWmig2d1.x sou=02sig_sou.H vel=00sig_velmig_my.H nrefvel=4 xmin=-56250 xmax=110000 \
	ymin=0 ymax=0 zmin=0 zmax=20000 nh=17 oh=-1200 dh=150 vsal=14700>$@

03sig_mig04b.H:
	Window3d squeeze=n < 01sig_prep.H n4=25 f4=75| ./bin/OWmig2d1.x sou=02sig_sou.H vel=00sig_velmig_my.H nrefvel=4 xmin=-56250 xmax=110000 \
	ymin=0 ymax=0 zmin=0 zmax=20000 nh=17 oh=-1200 dh=150 vsal=14700>$@

03sig_mig05b.H:
	Window3d squeeze=n < 01sig_prep.H n4=25 f4=100| ./bin/OWmig2d1.x sou=02sig_sou.H vel=00sig_velmig_my.H nrefvel=4 xmin=-56250 xmax=110000 \
	ymin=0 ymax=0 zmin=0 zmax=20000 nh=17 oh=-1200 dh=150 vsal=14700>$@

03sig_mig06b.H:
	Window3d squeeze=n < 01sig_prep.H n4=25 f4=125| ./bin/OWmig2d1.x sou=02sig_sou.H vel=00sig_velmig_my.H nrefvel=4 xmin=-56250 xmax=110000 \
	ymin=0 ymax=0 zmin=0 zmax=20000 nh=17 oh=-1200 dh=150 vsal=14700>$@

03sig_mig_totb.H:
	Math file1=03sig_mig01b.H file2=03sig_mig02b.H file3=03sig_mig03b.H file4=03sig_mig04b.H file5=03sig_mig05b.H file6=03sig_mig06b.H exp="file1+file2+file3+file4+file5+file6" > $@

03sig_mig_par.H:
	nohup Parallel command="./bin/OWmig2d1.x" files="01sig_prep.H,02sig_sou.H,00sig_velmig_my.H,$@" \
	tags="D,D_d,V,R" usage="INPUT,INPUT,INPUT,OUTPUT" \
	nblock=126 axis=4,4,4,6 pverb=2 file_type="DISTRIBUTE,COPY,COPY,COPY" mach_mfile=machinefile_200 \
	nrefvel=4 xmin=-56250 xmax=110000 ymin=0 ymax=0 zmin=0 zmax=20000 nh=17 oh=-1200 dh=150 vsal=14700 

04sig_mig_brad.H:
	@date
	nohup Parallel command=/net/koko/claudio/Research/Off2ang/ShotMigNew \
	files="01sig_prep.H,02sig_sou.H,00sig_velmig_brad.H,$@" tags="D,D_d,V,R" usage="INPUT,INPUT,INPUT,OUTPUT" \
	nblock=200 axis=4,4,4,6 pverb=2 file_type="DISTRIBUTE,COPY,COPY,COPY" mach_mfile=machinefile_200 \
	nx_shots=200 ox_shots=10000 par=go_sis.P taper=y taper_pcent=2 marine=y minx=-56250 maxx=110000 #restart=y

04sig_mig_brad1.H:
	@date
	Window3d squeeze=n < 01sig_prep.H n4=126 > zprep.H
	nohup Parallel command=/net/koko/claudio/Research/Off2ang/ShotMigNew \
	files="zprep.H,02sig_sou.H,00sig_velmig_brad.H,$@" tags="D,D_d,V,R" usage="INPUT,INPUT,INPUT,OUTPUT" \
	nblock=126 axis=4,4,4,6 pverb=2 file_type="DISTRIBUTE,COPY,COPY,COPY" mach_mfile=machinefile_200 \
	nx_shots=126 ox_shots=10000 par=go_sis.P taper=y taper_pcent=2 marine=y minx=-56250 maxx=110000 #restart=y

04sig_mig_brad2.H:
	@date
	Window3d squeeze=n < 01sig_prep.H n4=126 > zprep.H
	nohup Parallel command=/net/koko/claudio/Research/Off2ang/ShotMigNew \
	files="zprep.H,02sig_sou.H,00sig_velmig_brad.H,$@" tags="D,D_d,V,R" usage="INPUT,INPUT,INPUT,OUTPUT" \
	nblock=126 axis=4,4,4,6 pverb=2 file_type="DISTRIBUTE,COPY,COPY,COPY" mach_mfile=machinefile_200 \
	nx_shots=126 ox_shots=10000 par=go_sis.P taper=y taper_pcent=2 marine=y minx=-56250 maxx=110000 #restart=y

05sig_ang_totb.H:
	< 03sig_mig_totb.H Transp reshape=1,2,4,5 plane=34 | Pad end1=100 beg2=52 end2=52 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 amax=40 \
	> $@

05sig_ang_brad.H:
	< 04sig_mig_brad1.H Transp reshape=2,4,5 plane=13 | Pad end1=100 beg2=52 end2=52 | /net/koko/bob/SEP_BOB/bin/OFF2ANG verb=y eps=30 amax=40 \
	> $@


###############################################################
#
# Original sigsbee - 1000's
#
#   Window3d min2=10000 n2=534 max1=30000 < 00sig_velmig_my.H 
1999prep.H:
	< nfs-splitspread-shots.H Window3d max1=5 n3=4 f3=30 | /net/koko/bob/SEP_BOB/bin/Transf shotmig=y is2d=y wei=y f_min=3 f_min1=5 f_max1=38 f_max=42 maxsize=1000 > j2.H;
	< j2.H Transp reshape=1,3,5 plane=23 | Pad beg1=400 end1=400 > $@

2000source.H:
	Window3d n4=1 < 1999prep.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&

2001mig.H:
	< 1999prep.H ./bin/OWmig2d2p.x sou=2000source.H vel=1000velmig_my.H nrefvel=4 xmin=-56250 xmax=110000 \
	ymin=0 ymax=0 zmin=0 zmax=20000 nh=33 oh=-1200 dh=75 vsal=14700 wbott=6000 >$@

#
#
1000velmig_my.H:
	< sigsbee2b_migration_velocity.H Window3d j2=2 > veltemp1.H;echo o2=10025 d2=75 >> veltemp1.H; <veltemp1.H Pad beg2=5000 end2=5000 extend=1 > $@; echo n3=1 o3=0 d3=1 >> $@

1000velmig_transp.H:
	< 1000velmig_my.H Transp reshape=1,3 plane=12 > $@

1001prep.H:
	< nfs-splitspread-shots.H Window3d max1=8.5 | /net/koko/bob/SEP_BOB/bin/Transf shotmig=y is2d=y wei=y f_min=3 f_min1=5 f_max1=38 f_max=42 maxsize=1000 > j2.H;
	< j2.H Transp reshape=1,3,5 plane=23 | Pad beg1=400 end1=400 > $@

1002source.H:
	Window3d n4=1 < 1001prep.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&

#	Window3d squeeze=n < 1001prep.H n4=1 f4=25| ./bin/OWmig2d1.x sou=1002source.H vel=1000velmig_my.H nrefvel=4 xmin=-56250 xmax=10000 \
#	ymin=0 ymax=0 zmin=0 zmax=10000 nh=17 oh=-1200 dh=150 vsal=14700>$@

1003mig_velmig_my01.H:
	Window3d squeeze=n < 1001prep.H n4=50 j4=2 f4=0| ./bin/OWmig2d2.x sou=1002source.H vel=1000velmig_my.H nrefvel=4 xmin=-56250 xmax=140000 \
	ymin=0 ymax=0 zmin=0 zmax=20000 nh=33 oh=-1200 dh=75 vsal=14700 wbott= 6000 >$@

1003mig_velmig_my02.H:
	Window3d squeeze=n < 1001prep.H n4=50 j4=2 f4=100| ./bin/OWmig2d2.x sou=1002source.H vel=1000velmig_my.H nrefvel=4 xmin=-56250 xmax=140000 \
	ymin=0 ymax=0 zmin=0 zmax=20000 nh=33 oh=-1200 dh=75 vsal=14700 wbott= 6000 >$@

1003mig_velmig_my03.H:
	Window3d squeeze=n < 1001prep.H n4=50 j4=2 f4=200| ./bin/OWmig2d2.x sou=1002source.H vel=1000velmig_my.H nrefvel=4 xmin=-56250 xmax=140000 \
	ymin=0 ymax=0 zmin=0 zmax=20000 nh=33 oh=-1200 dh=75 vsal=14700 wbott= 6000 >$@

1003mig_velmig_my04.H:
	Window3d squeeze=n < 1001prep.H n4=50 j4=2 f4=300| ./bin/OWmig2d2.x sou=1002source.H vel=1000velmig_my.H nrefvel=4 xmin=-56250 xmax=140000 \
	ymin=0 ymax=0 zmin=0 zmax=20000 nh=33 oh=-1200 dh=75 vsal=14700 wbott= 6000 >$@

1003mig_velmig_my05.H:
	Window3d squeeze=n < 1001prep.H n4=27 j4=2 f4=400| ./bin/OWmig2d2.x sou=1002source.H vel=1000velmig_my.H nrefvel=4 xmin=-56250 xmax=140000 \
	ymin=0 ymax=0 zmin=0 zmax=20000 nh=33 oh=-1200 dh=75 vsal=14700 wbott= 6000 >$@

#004mig_brad.H:
#@date
#nohup Parallel command=/net/koko/claudio/Research/Off2ang/ShotMigNew \
#files="1001prep.H,1002source.H,1000velmig_transp.H,$@" tags="D,D_d,V,R" usage="INPUT,INPUT,INPUT,OUTPUT" \
#nblock=454 axis=4,4,4,6 pverb=2 file_type="DISTRIBUTE,COPY,COPY,COPY" mach_mfile=machinefile_200 \
#nx_shots=454 ox_shots=10925 par=go_sis.P taper=y taper_pcent=2 marine=y minx=-56250 maxx=110000 #restart=y

1004mig_brad.H:
	Window3d squeeze=n n4=2 j4=20 f4=300 <1001prep.H >z.H
	date
	nohup Parallel command=/net/koko/claudio/Research/Off2ang/ShotMigNew \
	files="z.H,1002source.H,1000velmig_transp.H,$@" tags="D,D_d,V,R" usage="INPUT,INPUT,INPUT,OUTPUT" \
	nblock=2 axis=4,4,4,6 pverb=2 file_type="DISTRIBUTE,COPY,COPY,COPY" mach_mfile=machinefile_200 \
	nx_shots=2 ox_shots=10925 par=go_sis.P taper=y taper_pcent=2 marine=y minx=-56250 maxx=110000 #restart=y