############################################
#
# Migrating Gauss anomaly data
#
00gauss_noenc: # 
	< ../Tese/1004gauss_prep.H ./bin/OWmig2d_par.x sou=../Tese/1004gauss_source.H vel=../Tese/1003gaussvel_.02pad.H \
	image=../Tese/1005gauss_noenc.H encrec=encrec1.H encsou=encsou1.H \
	vsal=3 nrefvel=4 xmin=-5.1 xmax=9.1 noencode=1 zmin=0 zmax=1.15 nh=1 oh=0 dh=0 wbott=0 ngap=200 > /dev/null

00gauss_encall%: # 
	< ../Tese/1004gauss_prep.H ./bin/OWmig2d_par.x sou=../Tese/1004gauss_source.H vel=../Tese/1003gaussvel_.02pad.H \
	image=../Tese/1005gauss_encall$*.H noencode=0 \
	vsal=3 nrefvel=4 xmin=-5.1 xmax=9.1 random=1 randall=1 nreal=$* zmin=0 zmax=1.15 nh=1 oh=0 dh=0 wbott=0 ngap=1 > /dev/null

00v2000_noenc: # 
	< ../Tese/1004v2000_prep.H ./bin/OWmig2d_par.x sou=../Tese/1004v2000_source.H vel=../Tese/1000_v2000mig.H \
	image=../Tese/1005v2000_noenc.H \
	vsal=3 nrefvel=4 xmin=-5.1 xmax=9.1 noencode=1 zmin=0 zmax=1.15 nh=1 oh=0 dh=0 wbott=0 ngap=200 > /dev/null

00v2000_encall%:
	< ../Tese/1004v2000_prep.H ./bin/OWmig2d_par.x sou=../Tese/1004v2000_source.H vel=../Tese/1000_v2000mig.H \
	image=../Tese/1005v2000_encall$*.H noencode=0 encrec=encrec.H encsou=encsou.H \
	vsal=3 nrefvel=1 xmin=-5.1 xmax=9.1 random=1 randall=1 nreal=$* zmin=0 zmax=1.15 nh=1 oh=0 dh=0 wbott=0 ngap=1 > /dev/null

00v2000_goldall%: # 
	< ../Tese/1004v2000_prep.H ./bin/OWmig2d_par.x sou=../Tese/1004v2000_source.H vel=../Tese/1000_v2000mig.H \
	image=../Tese/1005v2000_encall$*.H noencode=0 encrec=encrec.H encsou=encsou.H \
	vsal=3 nrefvel=1 xmin=-5.1 xmax=9.1 gold=1 randall=1 nreal=$* shift1=$* zmin=0 zmax=1.15 nh=1 oh=0 dh=0 wbott=0 ngap=1 > /dev/null

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

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

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

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

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

########################################
#
# Testing random phases
#
001sig_rand_%: 
#	Window3d squeeze=n min4=12000 n4=400 < 01sig_prepw256.H > 01sig_prep_select.H
	< 01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1w256.H vel=00sig_velmig_fix.H image=001sig_mig_rand_gap$*.H \
	nrefvel=4 xmin=-38000 xmax=121000 random=1 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=$* > /dev/null

0015sig_rand_%: # new random function
#	Window3d squeeze=n min4=12000 n4=400 < 01sig_prepw256.H > 01sig_prep_select.H
	< 01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1w256.H vel=00sig_velmig_fix.H image=0015sig_mig_rand_gap$*.H \
	nrefvel=4 xmin=-38000 xmax=121000 random=1 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=$* > /dev/null

002sig_mseq_%: 
	< 01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1w256.H vel=00sig_velmig_fix.H image=004sig_mig_mseq_gap$*.H \
	nrefvel=4 xmin=-38000 xmax=121000 mseq=1 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=$* > /dev/null

003sig_golay_%: # one random code per shot
#	Window3d squeeze=n min4=12000 n4=400 < 01sig_prepw256.H > 01sig_prep_select.H
	< 01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1w256.H vel=00sig_velmig_fix.H image=003sig_mig_golay_gap$* \
	nrefvel=4 xmin=-38000 xmax=70000 golay=1 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=$* > /dev/null

004sig_gold_%: # one random code per shot
#	Window3d squeeze=n min4=12000 n4=400 < 01sig_prepw256.H > 01sig_prep_select.H
	< 01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1w256.H vel=00sig_velmig_fix.H image=005sig_mig_gold_gap$*.H \
	nrefvel=4 xmin=-38000 xmax=70000 gold=1 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=$* > /dev/null

05sig_mig_grid2: # one random code per shot
#	Window3d squeeze=n min4=12000 n4=400 < 01sig_prepw256.H > 01sig_prep_select.H
	< ../1Wpar/01sig_prep_grid2.H ./bin/OWmig2d_par.x sou=../1Wpar/02sig_sou1.H 02sig_sou1w256.H vel=../1Wpar/00sig_velmig_grid2.H 00sig_velmig_fix.H image=05sig_mig_grid2.H\
	encsou= 001rnd_encode_sou.H vsal=14500 nrefvel=4 xmin=-38000 xmax=121000 random=1 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=194 400 > /dev/null

05sig_mig: # 
#	Window3d squeeze=n min4=12000 n4=400 < 01sig_prepw256.H > 01sig_prep_select.H
	< 01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1w256.H vel=00sig_velmig_fix.H image=05sig_mig_noencode.H\
	vsal=14500 nrefvel=4 xmin=-38000 xmax=121000 noencode=1 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=400 > /dev/null

004tst: # one random code per shot
#	Window3d squeeze=n min4=12000 n4=400 < 01sig_prepw256.H > 01sig_prep_select.H
	< 01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1w256.H vel=00sig_velmig_fix.H image=004tst.H \
	nrefvel=4 xmin=-38000 xmax=121000 mseq=1 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=40 > /dev/null

z: 
	< 01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1w256.H vel=00sig_velmig_fix.H image=z.H \
	nrefvel=4 xmin=-38000 xmax=121000 mseq=1 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=1 > /dev/null

#####################################
#
# Phase encoding all the shots 
#      in n realizations
000sig_rndall_tst_%:
#	Window3d squeeze=n min4=12000 n4=40 < 01sig_prepw256.H > z01sig_prep_select.H
	< 01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1w256.H vel=00sig_velmig_fix.H image=000sig_rndall_nre$*.H \
	nrefvel=4 xmin=-38000 xmax=121000 random=1 randall=1 nreal=$* zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=4 node= 6 > /dev/null

000sig_goldall_tst_%:
#	Window3d squeeze=n min4=12000 n4=40 < 01sig_prepw256.H > z01sig_prep_select.H
	< 01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1w256.H vel=00sig_velmig_fix.H image=000sig_goldall_nre$*.H \
	nrefvel=4 xmin=-38000 xmax=121000 gold=1 randall=1 nreal=$* shift=$* shift1=5 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=4 node= 6 > /dev/null

writephase_%:
#	Window3d squeeze=n min4=12000 n4=40 < 01sig_prepw256.H > z01sig_prep_select.H
	< 01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1w256.H vel=00sig_velmig_fix.H image=wrt$*.H \
	nrefvel=4 xmin=-38000 xmax=121000 random= 1 gold=1 randall=1 nreal=$* shift=$* shift1=5 zmin=0 zmax=24000 \
	nh=1 oh=0 dh=0 wbott=0 ngap=4 node= 6 code=gold_all.H random_all.H > /dev/null

##################################
#
# Xiao's frequency independent
#       phase encoding
001sig_xiao_tst:
	Window3d squeeze=n min4=12000 n4=40 < 01sig_prepw256.H > z01sig_prep_select.H
	< z01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1w256.H vel=00sig_velmig_fix.H image=001sig_mig_xiao_tst.H \
	nrefvel=4 xmin=-38000 xmax=121000 xiao=1 random=1 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=4 node= 6 > /dev/null

001sig_xiao_%:
	< 01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1w256.H vel=00sig_velmig_fix.H image=001sig_mig_xiao_gap$*.H \
	nrefvel=4 xmin=-38000 xmax=121000 xiao=1 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=$* node=6 > /dev/null


###################################
#
# Test of fftw + omp
#
tst_omp:
	Window3d squeeze=n n4=4 < ../Tese/1004gauss_prep.H > z1004gauss_prep.H
	< z1004gauss_prep.H ./bin/OWmig2d_par.x sou=../Tese/1004gauss_source.H vel=../Tese/1003gaussvel_.02pad.H \
	image=1.H \
	vsal=3 nrefvel=4 xmin=-5.1 xmax=9.1 noencode=1 zmin=0 zmax=1.15 nh=1 oh=0 dh=0 wbott=0 ngap=4 > /dev/null

tst_omp_sig:
	Window3d squeeze=n min4=12000 n4=2 < 01sig_prepw256.H > tst_omp_sig.H
	< tst_omp_sig.H ./bin/OWmig2d_par.x sou=02sig_sou1w256.H vel=00sig_velmig_fix.H image=1.H \
	nrefvel=4 xmin=-38000 xmax=121000 noencode=1 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=6000 ngap=2 > /dev/null