02sig_sou1w128.H:
	< ../OneWay/nfs-splitspread-shots.H Window3d min1=4.28 max1=8 n3=1 | /net/koko/bob/SEP_BOB/bin/Transf shotmig=y is2d=y wei=y \
	f_min=2 f_min1=6 f_max1=30 f_max=36 maxsize=1000 > j2.H
	Window3d n4=1 < j2.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&

####################################
#
# Xwell encoding
#
000no_encode: # changed code=0 in the program
	Pad <../1Wpar_phenc/02sig_sou1w256.H beg1=400 end1=400 n4out=3 > 02sig_sou1w256_pad.H; echo o1=0 d1=75 d3=0.05 o3=2 d4=9975. o4=19950. >> 02sig_sou1w256_pad.H
	Math file1=../1Wpar_phenc/00sig_velmig_fix.H exp="file1*0+5000" > 00velmig_5000.H
	< 02sig_sou1w256_pad.H ./bin/OW_encode.x vel=00velmig_5000.H 00sig_velmig_fix.H wvfld=000no_encode.H decoded=000no_decoded.H \
	nrefvel=1 zmax=4000 > /dev/null

001rnd_encode: # one random code per shot
#	Pad <../1Wpar_phenc/02sig_sou1w256.H beg1=400 end1=400 n4out=3 > 02sig_sou1w256_pad.H; echo o1=0 d1=75 d3=0.05 o3=2 d4=9975. o4=19950. >> 02sig_sou1w256_pad.H
#	Math file1=../1Wpar_phenc/00sig_velmig_fix.H exp="file1*0+5000" > 00velmig_5000.H
	< 02sig_sou1w256_pad.H ./bin/OW_encode1.x vel=00velmig_5000.H 00sig_velmig_fix.H wvfld=001rnd_encode.H decoded=001rnd_decoded.H \
	nrefvel=1 zmax=4000 random=1 nreal=1 > /dev/null

002mseq_encode: # one random code per shot
#	Pad <../1Wpar_phenc/02sig_sou1w256.H beg1=400 end1=400 n4out=3 > 02sig_sou1w256_pad.H; echo o1=0 d1=75 d3=0.05 o3=2 d4=9975. o4=19950. >> 02sig_sou1w256_pad.H
#	Math file1=../1Wpar_phenc/00sig_velmig_fix.H exp="file1*0+5000" > 00velmig_5000.H
	< 02sig_sou1w256_pad.H ./bin/OW_encode.x vel=00velmig_5000.H 00sig_velmig_fix.H wvfld=002mseq_encode.H decoded=002mseq_decoded.H \
	nrefvel=1 zmax=4000 mseq=1 nreal=1> /dev/null

002mseq_encode1: # one random code per shot
#	Pad <../1Wpar_phenc/02sig_sou1w256.H beg1=400 end1=400 n4out=3 > 02sig_sou1w256_pad.H; echo o1=0 d1=75 d3=0.05 o3=2 d4=9975. o4=19950. >> 02sig_sou1w256_pad.H
#	Math file1=../1Wpar_phenc/00sig_velmig_fix.H exp="file1*0+5000" > 00velmig_5000.H
	< 02sig_sou1w256_pad.H ./bin/OW_encode1.x vel=00velmig_5000.H 00sig_velmig_fix.H wvfld=002mseq_encode1.H decoded=002mseq_decoded1.H \
	nrefvel=1 zmax=4000 mseq=1 nreal=1> /dev/null

003gold_encode: # one random code per shot
	Pad <../1Wpar_phenc/02sig_sou1w256.H beg1=400 end1=400 n4out=3 > 02sig_sou1w256_pad.H; echo o1=0 d1=75 d3=0.05 o3=2 d4=9975. o4=19950. >> 02sig_sou1w256_pad.H
#	Math file1=../1Wpar_phenc/00sig_velmig_fix.H exp="file1*0+5000" > 00velmig_5000.H
	< 02sig_sou1w256_pad.H ./bin/OW_encode1.x vel=00velmig_5000.H 00sig_velmig_fix.H wvfld=003gold_encode.H decoded=003gold_decoded.H \
	nrefvel=1 zmax=4000 gold=1 nreal=1> /dev/null

003gold_encode128: # one random code per shot d3=0.05 o3=2
	Pad <02sig_sou1w128.H beg1=400 end1=400 n4out=3 > 02sig_sou1w128_pad.H; echo o1=0 d1=75 d4=9975. o4=19950. >> 02sig_sou1w128_pad.H
#	Math file1=../1Wpar_phenc/00sig_velmig_fix.H exp="file1*0+5000" > 00velmig_5000.H
	< 02sig_sou1w128_pad.H ./bin/OW_encode1.x vel=00velmig_5000.H 00sig_velmig_fix.H wvfld=003gold_encode128.H decoded=003gold_decoded128.H \
	nrefvel=1 zmax=4000 gold=1 nreal=1> /dev/null

#######################################
#
# Using outcode => Biondo's autocorrelation 
#
001rnd: # one random code per shot
	Pad <../1Wpar_phenc/02sig_sou1w256.H beg1=400 end1=400 n4out=200 > 02sig_sou1w256_pad200.H; echo o1=0 d1=75 d3=0.2 o3=2 d4=75. o4=19950. >> 02sig_sou1w256_pad200.H
#	Math file1=../1Wpar_phenc/00sig_velmig_fix.H exp="file1*0+5000" > 00velmig_5000.H
	< 02sig_sou1w256_pad200.H ./bin/OW_encode1.x vel=00velmig_5000.H wvfld=001rnd.H outcode=1 code=001rnd_code.H \
	nrefvel=1 zmax=4000 random=1 nreal=1 > /dev/null

002mseq: # one random code per shot
	< 02sig_sou1w256_pad200.H ./bin/OW_encode1.x vel=00velmig_5000.H wvfld=002mseq.H outcode=1 code=002mseq_code.H \
	nrefvel=1 zmax=4000 mseq=1 nreal=1> /dev/null

########################################
#
# Testing random phases
#
001sig_mig_rand: # one random code per shot in the set (4 codes in total)
	Window3d squeeze=n min4=12000 n4=40 < 01sig_prep.H > 01sig_prep_select.H
	< 01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1.H vel=00sig_velmig_fix.H image=001sig_mig_rand.H encrec=encrec.H \
	encsou=encsou.H GS= GS.H GR= GR.H nrefvel=4 xmin=-38000 xmax=70000 random=1 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=10 > /dev/null

001sig_rand_%: # one random code per shot
	< 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 mseq=1 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=$* > /dev/null

003sig_mig_golay: # one random code per shot
	Window3d squeeze=n min4=12000 n4=40 < 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.H encrec= encrec2.H \
	encsou= encsou2.H nrefvel=4 xmin=-38000 xmax=70000 golay=1 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=10 > /dev/null

004sig_mig_mseq: # 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=004sig_mig_mseq.H \
	nrefvel=4 xmin=-38000 xmax=121000 mseq=1 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=10 > /dev/null

004sig_mig_mseqint: # 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=004sig_mig_mseqint.H \
	nrefvel=4 xmin=-38000 xmax=121000 mseq=1 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=10 > /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

004sig_mseq_%: # one random code per shot
	< 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

005sig_mig_gold: # one random code per shot
	Window3d squeeze=n min4=12000 n4=40 < 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.H \
	nrefvel=4 xmin=-38000 xmax=70000 gold=1 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=10 > /dev/null
z: # one random code per shot
#	Window3d squeeze=n min4=12000 n4=40 < 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=z.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