#############################################################
#
#sep136
###############################################################
9000data_hor.H:
	susynlv ref="1:0,.5.;5.,.5." nt=2001 dt=0.002 dxs=0.02 nxs=401 fpeak=20 fxs=0 dxo=0.02 nxo=401 fxo=-4 v00=2 > junk1.su
	< junk1.su Su2sep > junk1.H
	<junk1.H Infill3d |Window3d j1=4 > $@; echo n2=401 o2=-4 d2=0.02 n3=401 o3=0 d3=0.02 >> $@

9000data_prep1.H:
	< 9000data_hor.H /net/koko/bob/SEP_BOB/bin/Transf shotmig=y \
	is2d=y wei=y f_min=2 f_min1=6 f_max1=40 f_max=50 maxsize=1000 > j2.H
	< j2.H Transp reshape=1,3,5 plane=23 | Pad beg1=200 end1=200 > $@

9000data_source1.H:
	< 9000data_hor.H Window3d n3=1 n2=1 min2=0 max1=.6 | Pad n1out=501 | /net/koko/bob/SEP_BOB/bin/Transf shotmig=y \
	is2d=y wei=y f_min=2 f_min1=6 f_max1=40 f_max=50 maxsize=1000 > y.H
	< y.H Window3d | Cabs > spec.H
	Math file1=spec.H exp="file1*0" > phse.H; Cmplx spec.H phse.H |\
	Transp reshape=1,3,5 plane=23| Transp reshape=1,3,5 plane=12 > $@

9000_mig_vlow.H: # With 401 shots separated 0.01km correct velocity
	@date
	/net/koko/claudio/Research/Others/ShotMig/ShotMig \
	D=9000data_prep1.H D_d=9000data_source1.H V=InvVel/44413vlow1.9_transp.H R=$@ \
	pverb=2  nx_shots=401 ox_shots=0 \
	par=go_sis.P taper=y taper_pcent=2 marine=y minx=-8 maxx=17  #restart=y

44413vlow1.9.H:
	Spike n1=600 o1=0 d1=0.004 n2=4000 o2=-8 d2=0.02 n3=1 o3=0 d3=0.001 > 1.H
	Math file1=1.H exp="file1*0+1.9" > $@

44413vlow1.9_transp.H:
	Transp plane=12 reshape=1,3 > $@

44413enc0: # correct velocity areal shot random encode in w
	< 44413vlow1.9.H
	< 9000_mig_vlow.H Window3d min1=1 max1=3 n5=256 | Transp plane=13 > junk.H
	./bin/OW_arealsht.x rec_cig=junk.H sou_cig=junk.H sou=4004data_recip_source_pad.H \
	vel=44413vlow1.9.H rec_out=44413areal_rec_enc0_velow1.H sou_out=44413areal_sou_enc0_velow1.H \
	ngap=101 nrefvel=1 vsal=6 padtraces=41 encode=0 injection=1 nreal= 5 > /dev/null 

44413mig_enc0_velow1.H: # correct velocity areal shot random encode in w
	< 44413areal_rec_enc0_velow1.H Transp plane=13 > 44413areal_rec_tempw1.H
	< 44413areal_sou_enc0_velow1.H Transp plane=13 > 44413areal_sou_tempw1.H
	< 44413areal_rec_tempw1.H InvVel/bin/OWmig2d_areal.x sou=44413areal_sou_tempw1.H \
	vel=44413vlow1.9.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=21 oh=-.20 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 migration=1 > $@ 
	Rm 44413*tempw1.H

44413mig_enc0_vel.H: # correct velocity areal shot random encode in w
	< 44413areal_rec_enc0_velow1.H Transp plane=13 > 44413areal_rec_tempw1.H
	< 44413areal_sou_enc0_velow1.H Transp plane=13 > 44413areal_sou_tempw1.H
	< 44413areal_rec_tempw1.H InvVel/bin/OWmig2d_areal.x sou=44413areal_sou_tempw1.H \
	vel=4000vel.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=21 oh=-.20 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 migration=1 > $@ 
	Rm 44413*tempw1.H

44413enca: # correct velocity areal shot random encode in w
	< InvVel/44413vlow1.9_transp.H Window3d | Transp > 44413vlow1.9.H
	< 9000_mig_vlow.H Window3d min1=1 max1=3 n5=256 | Transp plane=13 > junk.H
	./bin/OW_arealsht.x rec_cig=junk.H sou_cig=junk.H sou=4004data_recip_source_pad.H \
	vel=44413vlow1.9.H rec_out=44413areal_rec_enca_velow1.H sou_out=44413areal_sou_enca_velow1.H \
	ngap=10 nrefvel=1 vsal=6 padtraces=41 encode=0 random=1 encw=1 eps=1 injection=1 nreal=1 > /dev/null 

44413mig_enca_velow1.H: # correct velocity areal shot random encode in w
	< 44413areal_rec_enca_velow1.H Transp plane=13 > 44413areal_rec_tempw1.H
	< 44413areal_sou_enca_velow1.H Transp plane=13 > 44413areal_sou_tempw1.H
	< 44413areal_rec_tempw1.H InvVel/bin/OWmig2d_areal.x sou=44413areal_sou_tempw1.H \
	vel=44413vlow1.9.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=21 oh=-.20 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 migration=1 > $@ 
	Rm 44413*tempw1.H

44413genca: # correct velocity areal shot random encode in w
	< InvVel/44413vlow1.9_transp.H Window3d | Transp > 44413vlow1.9.H
	< 9000_mig_vlow.H Window3d min1=1 max1=3 n5=256 | Transp plane=13 > junk.H
	./bin/OW_arealsht.x rec_cig=junk.H sou_cig=junk.H sou=4004data_recip_source_pad.H \
	vel=44413vlow1.9.H rec_out=44413areal_rec_genca_velow1.H sou_out=44413areal_sou_genca_velow1.H \
	ngap=1 nrefvel=1 vsal=6 padtraces=41 encode=1 gold=1 perfect=1 encw=1 eps=1 injection=1 nreal=5 > /dev/null 

44413mig_genca_velow1.H: # correct velocity areal shot random encode in w
	< 44413areal_rec_genca_velow1.H Transp plane=13 > 44413areal_rec_tempw1.H
	< 44413areal_sou_genca_velow1.H Transp plane=13 > 44413areal_sou_tempw1.H
	< 44413areal_rec_tempw1.H InvVel/bin/OWmig2d_areal.x sou=44413areal_sou_tempw1.H \
	vel=44413vlow1.9.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=21 oh=-.20 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 migration=1 > $@ 
	Rm 44413*tempw1.H

44413optia: # correct velocity areal shot random encode in w
	< 9000_mig_vlow.H Window3d min1=1 max1=3 n5=256 | Transp plane=13 > junk.H
	< ~biondo/Code/Anneal/GenAnneal/Test/Claudio/Anneal_Corr_indfreq_wide.H Window3d n1=100 f1=50 n6=5| Transp > \
	xx.H ; Math file1=xx.H exp="file1*2" > Anneal_Corr_indfreq_wide.H
	./bin/OW_arealsht.x rec_cig=junk.H sou_cig=junk.H sou=4004data_recip_source_pad.H InCode=Anneal_Corr_indfreq_wide.H \
	vel=44413vlow1.9.H rec_out=44413areal_rec_optia_velow1.H sou_out=44413areal_sou_optia_velow1.H ngap=1 rndcte=1231 nrefvel=1 \
	vsal=6 padtraces=41 encode=1 encw=1 eps=1 injection=1 nreal=5 node=12 8 > /dev/null 

44413mig_optia_velow1.H: # correct velocity areal shot random encode in w
	< 44413areal_rec_optia_velow1.H Transp plane=13 > 44403areal_rec_tempw.H
	< 44413areal_sou_optia_velow1.H Transp plane=13 > 44403areal_sou_tempw.H
	< 44403areal_rec_tempw.H InvVel/bin/OWmig2d_areal.x sou=44403areal_sou_tempw.H \
	vel=44413vlow1.9.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=21 oh=-.20 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 ws= xxx_GS.H wr= xxx_GR.H migration=1> $@ 
	Rm 44404*tempz.H

44413encw11: # correct velocity areal shot random encode in w
	< InvVel/44413vlow1.9_transp.H Window3d | Transp > 44413vlow1.9.H
	< 9000_mig_vlow.H Window3d min1=1 max1=3 n5=256 | Transp plane=13 > junk.H
	./bin/OW_arealsht.x rec_cig=junk.H sou_cig=junk.H sou=4004data_recip_source_pad.H \
	vel=44413vlow1.9.H rec_out=44413areal_rec_encw_velow11.H sou_out=44413areal_sou_encw_velow11.H \
	OutCode=Random_conventional.H ngap=10 rndcte=1231 nrefvel=1 \
	vsal=6 padtraces=41 encode=1 random=1 encw=1 eps=1 injection=1 nreal=1 correlations= 44413_rnd_correl.H > /dev/null 

44413mig_encw_velow11.H: # correct velocity areal shot random encode in w
	< 44413areal_rec_encw_velow11.H Transp plane=13 > 44413areal_rec_tempw1.H
	< 44413areal_sou_encw_velow11.H Transp plane=13 > 44413areal_sou_tempw1.H
	< 44413areal_rec_tempw1.H InvVel/bin/OWmig2d_areal.x sou=44413areal_sou_tempw1.H \
	vel=44413vlow1.9.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=21 oh=-.20 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 migration=1 > $@ 
	Rm 44413*tempw1.H

44413encw1: # correct velocity areal shot random encode in w
	< InvVel/44413vlow1.9_transp.H Window3d | Transp > 44413vlow1.9.H
	< 9000_mig_vlow.H Window3d min1=1 max1=3 n5=256 | Transp plane=13 > junk.H
	./bin/OW_arealsht.x rec_cig=junk.H sou_cig=junk.H sou=4004data_recip_source_pad.H \
	vel=44413vlow1.9.H rec_out= 44413areal_rec_encw_velow1.H sou_out= 44413areal_sou_encw_velow1.H \
	OutCode=Random_conventional.H ngap=10 rndcte=1231 nrefvel=1 \
	vsal=6 padtraces=41 encode=1 random=1 encw=1 eps=1 injection=1 nreal=5 correlations=44413_rnd_correl.H > /dev/null 

44413mig_encw_velow1.H: # correct velocity areal shot random encode in w
	< 44413areal_rec_encw_velow1.H Transp plane=13 > 44413areal_rec_tempw1.H
	< 44413areal_sou_encw_velow1.H Transp plane=13 > 44413areal_sou_tempw1.H
	< 44413areal_rec_tempw1.H InvVel/bin/OWmig2d_areal.x sou=44413areal_sou_tempw1.H \
	vel=44413vlow1.9.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=21 oh=-.20 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 migration=1 > $@ 
	Rm 44413*tempw1.H

44413mig_encw_vel.H: # correct velocity areal shot random encode in w
	< 44413areal_rec_encw_velow1.H Transp plane=13 > 44413areal_rec_tempw1.H
	< 44413areal_sou_encw_velow1.H Transp plane=13 > 44413areal_sou_tempw1.H
	< 44413areal_rec_tempw1.H InvVel/bin/OWmig2d_areal.x sou=44413areal_sou_tempw1.H \
	vel=4000vel.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=21 oh=-.20 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 migration=1 > $@ 
	Rm 44413*tempw1.H

44413gencw1: # correct velocity areal shot random encode in w
	< InvVel/44413vlow1.9_transp.H Window3d | Transp > 44413vlow1.9.H
	< 9000_mig_vlow.H Window3d min1=1 max1=3 n5=256 | Transp plane=13 > junk.H
	./bin/OW_arealsht.x rec_cig=junk.H sou_cig=junk.H sou=4004data_recip_source_pad.H \
	vel=44413vlow1.9.H rec_out=44413areal_rec_gencw_velow1a.H sou_out=44413areal_sou_gencw_velow1a.H ngap=10 rndcte=1231 nrefvel=1 \
	vsal=6 padtraces=41 encode=1 gold=1 encw encz=1 perfect=1 eps=1 injection=1 nreal=5 nodes= 12 OutCode= Random_goldperfect.H \
	> /dev/null 

44413mig_gencw_velow1a.H: # correct velocity areal shot random encode in w
	< 44413areal_rec_gencw_velow1a.H Transp plane=13 > 44413areal_rec_tempw1.H
	< 44413areal_sou_gencw_velow1a.H Transp plane=13 > 44413areal_sou_tempw1.H
	< 44413areal_rec_tempw1.H ./bin/OWmig2d_areal.x sou=44413areal_sou_tempw1.H \
	vel=44413vlow1.9.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=21 oh=-.20 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 ws= xxx_GS.H wr= xxx_GR.H migration=1 > $@ 
	Rm 44413*tempw1.H

44413mig_gencw_vel.H: # correct velocity areal shot random encode in w
	< 44413areal_rec_gencw_velow1.H Transp plane=13 > 44413areal_rec_tempw1.H
	< 44413areal_sou_gencw_velow1.H Transp plane=13 > 44413areal_sou_tempw1.H
	< 44413areal_rec_tempw1.H InvVel/bin/OWmig2d_areal.x sou=44413areal_sou_tempw1.H \
	vel=4000vel.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=21 oh=-.20 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 ws= xxx_GS.H wr= xxx_GR.H migration=1 > $@ 
	Rm 44413*tempw1.H

###############################################
#Using max.separation
#
44413gencw2: # correct velocity areal shot random encode in w
	< InvVel/44413vlow1.9_transp.H Window3d | Transp > 44413vlow1.9.H
	< 9000_mig_vlow.H Window3d min1=1 max1=3 n5=256 | Transp plane=13 > junk.H
	./bin/OW_arealsht.x rec_cig=junk.H sou_cig=junk.H sou=4004data_recip_source_pad.H \
	vel=44413vlow1.9.H rec_out=44413areal_rec_gencw_velow2.H sou_out=44413areal_sou_gencw_velow2.H ngap=10 rndcte=1231 nrefvel=1 \
	vsal=6 padtraces=41 encode=1 gold=1 encw=1 perfect=1 eps=1 injection=1 nreal=5 node= 12 jshift=1  OutCode= Random_goldperfect.H \
	> /dev/null 

44413mig_gencw_velow2.H: # correct velocity areal shot random encode in w
	< 44413areal_rec_gencw_velow2.H Transp plane=13 > 44413areal_rec_tempw2.H
	< 44413areal_sou_gencw_velow2.H Transp plane=13 > 44413areal_sou_tempw2.H
	< 44413areal_rec_tempw2.H InvVel/bin/OWmig2d_areal.x sou=44413areal_sou_tempw2.H \
	vel=44413vlow1.9.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=21 oh=-.20 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 ws= xxx_GS.H wr= xxx_GR.H migration=1 > $@ 
	Rm 44413*tempw2.H

44413gencw2.1: # correct velocity areal shot random encode in w
	< InvVel/44413vlow1.9_transp.H Window3d | Transp > 44413vlow1.9.H
	< 9000_mig_vlow.H Window3d min1=1 max1=3 n5=256 | Transp plane=13 > junk.H
	./bin/OW_arealsht.x rec_cig=junk.H sou_cig=junk.H sou=4004data_recip_source_pad.H \
	vel=44413vlow1.9.H rec_out=44413areal_rec_gencw_velow2.1.H sou_out=44413areal_sou_gencw_velow2.1.H ngap=10 rndcte=1231 nrefvel=1 \
	vsal=6 padtraces=41 encode=1 gold=1 encw=1 perfect=1 eps=1 injection=1 nreal=1 node= 12 jshift=1  OutCode= Random_goldperfect.H \
	> /dev/null 

44413mig_gencw_velow2.1.H: # correct velocity areal shot random encode in w
	< 44413areal_rec_gencw_velow2.1.H Transp plane=13 > 44413areal_rec_tempw2.H
	< 44413areal_sou_gencw_velow2.1.H Transp plane=13 > 44413areal_sou_tempw2.H
	< 44413areal_rec_tempw2.H InvVel/bin/OWmig2d_areal.x sou=44413areal_sou_tempw2.H \
	vel=44413vlow1.9.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=21 oh=-.20 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 ws= xxx_GS.H wr= xxx_GR.H migration=1 > $@ 
	Rm 44413*tempw2.H

##
44413gencw2.%: # correct velocity areal shot random encode in w
	< InvVel/44413vlow1.9_transp.H Window3d | Transp > 44413vlow1.9.H
	< 9000_mig_vlow.H Window3d min1=1 max1=3 n5=256 | Transp plane=13 > junk.H
	./bin/OW_arealsht.x rec_cig=junk.H sou_cig=junk.H sou=4004data_recip_source_pad.H \
	vel=44413vlow1.9.H rec_out=44413areal_rec_gencw_velow2.$*.H sou_out=44413areal_sou_gencw_velow2.$*.H ngap=10 rndcte=1231 nrefvel=1 \
	vsal=6 padtraces=41 encode=1 gold=1 encw=1 perfect=1 eps=1 injection=1 nreal=1 node=10 jshift=2 fac=$* OutCode= Random_goldperfect.H \
	> /dev/null 

44413mig_gencw_velow2.%.H: # correct velocity areal shot random encode in w
	< 44413areal_rec_gencw_velow2.$*.H Transp plane=13 > 44413areal_rec_tempw2.H
	< 44413areal_sou_gencw_velow2.$*.H Transp plane=13 > 44413areal_sou_tempw2.H
	< 44413areal_rec_tempw2.H InvVel/bin/OWmig2d_areal.x sou=44413areal_sou_tempw2.H \
	vel=44413vlow1.9.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. node=10 nh=21 oh=-.20 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 migration=1 > $@ 
	Rm 44413*tempw2.H
##
44413gencw3: # correct velocity areal shot random encode in w
	< InvVel/44413vlow1.9_transp.H Window3d | Transp > 44413vlow1.9.H
	< 9000_mig_vlow.H Window3d min1=1 max1=3 n5=256 | Transp plane=13 > junk.H
	./bin/OW_arealsht.x rec_cig=junk.H sou_cig=junk.H sou=4004data_recip_source_pad.H \
	vel=44413vlow1.9.H rec_out= 44413areal_rec_gencw_velow3.H sou_out= 44413areal_sou_gencw_velow3.H ngap=10 rndcte=1231 nrefvel=1 \
	vsal=6 padtraces=41 encode=1 gold=1 encw=1 perfect=1 eps=1 injection=1 nreal=5 node=8 jshift=3  OutCode=  Random_goldperfect.H correlations=44413_correlations.H \
	> /dev/null 

44413mig_gencw_velow3.H: # correct velocity areal shot random encode in w
	< 44413areal_rec_gencw_velow3.H Transp plane=13 > 44413areal_rec_tempw3.H
	< 44413areal_sou_gencw_velow3.H Transp plane=13 > 44413areal_sou_tempw3.H
	< 44413areal_rec_tempw3.H InvVel/bin/OWmig2d_areal.x sou=44413areal_sou_tempw3.H \
	vel=44413vlow1.9.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=21 oh=-.20 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 node=10 ws= xxx_GS.H wr= xxx_GR.H migration=1 > $@ 
	Rm 44413*tempw3.H

44413gencw31: # correct velocity areal shot random encode in w
	< InvVel/44413vlow1.9_transp.H Window3d | Transp > 44413vlow1.9.H
	< 9000_mig_vlow.H Window3d min1=1 max1=3 n5=256 | Transp plane=13 > junk.H
	./bin/OW_arealsht.x rec_cig=junk.H sou_cig=junk.H sou=4004data_recip_source_pad.H \
	vel=44413vlow1.9.H rec_out=44413areal_rec_gencw_velow31.H sou_out=44413areal_sou_gencw_velow31.H ngap=10 rndcte=1231 nrefvel=1 \
	vsal=6 padtraces=41 encode=1 gold=1 encw=1 perfect=1 eps=1 injection=1 nreal=1 node=8 jshift=3  OutCode= Random_goldperfect.H correlations= 44413_correlations.H \
	> /dev/null 

44413mig_gencw_velow31.H: # correct velocity areal shot random encode in w
	< 44413areal_rec_gencw_velow31.H Transp plane=13 > 44413areal_rec_tempw3.H
	< 44413areal_sou_gencw_velow31.H Transp plane=13 > 44413areal_sou_tempw3.H
	< 44413areal_rec_tempw3.H InvVel/bin/OWmig2d_areal.x sou=44413areal_sou_tempw3.H \
	vel=44413vlow1.9.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=21 oh=-.20 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 ws= xxx_GS.H wr= xxx_GR.H migration=1 > $@ 
	Rm 44413*tempw3.H

#############################
44413gencw4: # correct velocity areal shot random encode in w
	< InvVel/44413vlow1.9_transp.H Window3d | Transp > 44413vlow1.9.H
	< 9000_mig_vlow.H Window3d min1=1 max1=3 n5=256 | Transp plane=13 > junk.H
	./bin/OW_arealsht.x rec_cig=junk.H sou_cig=junk.H sou=4004data_recip_source_pad.H \
	vel=44413vlow1.9.H rec_out=44413areal_rec_gencw_velow4.H sou_out=44413areal_sou_gencw_velow4.H ngap=10 rndcte=1231 nrefvel=1 \
	vsal=6 padtraces=41 encode=1 gold=1 encw=1 perfect=1 eps=1 injection=1 nreal=1 node=12 jshift=4  OutCode= Random_goldperfect.H correlations= 44413_correlations.H \
	> /dev/null 

44413mig_gencw_velow4.H: # correct velocity areal shot random encode in w
	< 44413areal_rec_gencw_velow4.H Transp plane=13 > 44413areal_rec_tempw3.H
	< 44413areal_sou_gencw_velow4.H Transp plane=13 > 44413areal_sou_tempw3.H
	< 44413areal_rec_tempw3.H InvVel/bin/OWmig2d_areal.x sou=44413areal_sou_tempw3.H \
	vel=44413vlow1.9.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=21 oh=-.20 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=1 encode=1 node=12 migration=1 > $@ 
	Rm 44413*tempw3.H

44413mig_gencw_velow4.%.H: # correct velocity areal shot random encode in w
	< 44413areal_rec_gencw_velow4.H Window3d squeeze=n n4=1 f4=$* | Transp plane=13 > 44413areal_rec_tempw3.H
	< 44413areal_sou_gencw_velow4.H Window3d squeeze=n n4=1 f4=$* | Transp plane=13 > 44413areal_sou_tempw3.H
	< 44413areal_rec_tempw3.H InvVel/bin/OWmig2d_areal.x sou=44413areal_sou_tempw3.H \
	vel=44413vlow1.9.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=21 oh=-.20 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=1 encode=1 node=12 migration=1 > $@ 
	Rm 44413*tempw3.H

#################################
##
44413gencw4.%: # correct velocity areal shot random encode in w
	< InvVel/44413vlow1.9_transp.H Window3d | Transp > 44413vlow1.9.H
	< 9000_mig_vlow.H Window3d min1=1 max1=3 n5=256 | Transp plane=13 > junk.H
	./bin/OW_arealsht.x rec_cig=junk.H sou_cig=junk.H sou=4004data_recip_source_pad.H \
	vel=44413vlow1.9.H rec_out=44413areal_rec_gencw_velow4.fac$*.H sou_out=44413areal_sou_gencw_velow4.fac$*.H ngap=10 rndcte=1231 nrefvel=1 \
	vsal=6 padtraces=41 encode=1 gold=1 encw=1 perfect=1 eps=1 injection=1 nreal=1 node=12 jshift=4 fac=$* OutCode= Random_goldperfect.H correlations= 44413_correlations.H \
	> /dev/null 

44413mig_gencw_velow4.fac%.H: # correct velocity areal shot random encode in w
	< 44413areal_rec_gencw_velow4.fac$*.H Transp plane=13 > 44413areal_rec_tempw3.H
	< 44413areal_sou_gencw_velow4.fac$*.H Transp plane=13 > 44413areal_sou_tempw3.H
	< 44413areal_rec_tempw3.H InvVel/bin/OWmig2d_areal.x sou=44413areal_sou_tempw3.H \
	vel=44413vlow1.9.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=21 oh=-.20 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=1 encode=1 node=12 migration=1 > $@ 
	Rm 44413*tempw3.H

44413mig_gencw_vel2.H: # correct velocity areal shot random encode in w
	< 44413areal_rec_gencw_velow2.H Transp plane=13 > 44413areal_rec_tempw2.H
	< 44413areal_sou_gencw_velow2.H Transp plane=13 > 44413areal_sou_tempw2.H
	< 44413areal_rec_tempw2.H InvVel/bin/OWmig2d_areal.x sou=44413areal_sou_tempw2.H \
	vel=4000vel.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=21 oh=-.20 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 ws= xxx_GS.H wr= xxx_GR.H migration=1 > $@ 
	Rm 44413*tempw2.H

Anneal_Corr_indfreq_wide.H:
	< ~biondo/Code/Anneal/GenAnneal/Test/Claudio/Anneal_Corr_indfreq_wide.H Window3d n1=100 f1=50 n6= 5| Transp > \
	xx.H ; Math file1=xx.H exp="file1*2" > $@

44413opti1: # correct velocity areal shot random encode in w
	< 9000_mig_vlow.H Window3d min1=1 max1=3 n5=256 | Transp plane=13 > junk.H
	< ~biondo/Code/Anneal/GenAnneal/Test/Claudio/Anneal_Corr_indfreq_wide.H Window3d n1=100 f1=50 n6=5| Transp > \
	xx.H ; Math file1=xx.H exp="file1*2" > Anneal_Corr_indfreq_wide.H
	./bin/OW_arealsht.x rec_cig=junk.H sou_cig=junk.H sou=4004data_recip_source_pad.H InCode=Anneal_Corr_indfreq_wide.H \
	vel=44413vlow1.9.H rec_out=44413areal_rec_optiw1_velow1.H sou_out=44413areal_sou_optiw1_velow1.H ngap=10 rndcte=1231 nrefvel=1 \
	vsal=6 padtraces=41 encode=1 encw=1 eps=1 injection=1 nreal=1 node=8 > /dev/null 

44413mig_opti_velow1.H: # correct velocity areal shot random encode in w
	< 44413areal_rec_optiw1_velow1.H Transp plane=13 > 44403areal_rec_tempw.H
	< 44413areal_sou_optiw1_velow1.H Transp plane=13 > 44403areal_sou_tempw.H
	< 44403areal_rec_tempw.H InvVel/bin/OWmig2d_areal.x sou=44403areal_sou_tempw.H \
	vel=44413vlow1.9.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=21 oh=-.20 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 node=8 migration=1> $@ 
	Rm 44403*tempz.H

44413opti5: # correct velocity areal shot random encode in w
	< 9000_mig_vlow.H Window3d min1=1 max1=3 n5=256 | Transp plane=13 > junk.H
	< ~biondo/Code/Anneal/GenAnneal/Test/Claudio/Anneal_Corr_indfreq_wide.H Window3d n1=100 f1=50 n6=5| Transp > \
	xx.H ; Math file1=xx.H exp="file1*2" > Anneal_Corr_indfreq_wide.H
	./bin/OW_arealsht.x rec_cig=junk.H sou_cig=junk.H sou=4004data_recip_source_pad.H InCode=Anneal_Corr_indfreq_wide.H \
	vel=44413vlow1.9.H rec_out=44413areal_rec_optiw5_velow1.H sou_out=44413areal_sou_optiw5_velow1.H ngap=10 nrefvel=1 \
	vsal=6 padtraces=41 encode=1 encw=1 eps=1 node=8 injection=1 nreal=5 > /dev/null 

44413mig_opti_velow5.H: # correct velocity areal shot random encode in w
	< 44413areal_rec_optiw5_velow1.H Transp plane=13 > 44403areal_rec_tempw.H
	< 44413areal_sou_optiw5_velow1.H Transp plane=13 > 44403areal_sou_tempw.H
	< 44403areal_rec_tempw.H InvVel/bin/OWmig2d_areal.x sou=44403areal_sou_tempw.H \
	vel=44413vlow1.9.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=21 oh=-.20 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 node=8 migration=1> $@ 
	Rm 44403*tempz.H

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

## No combination of modeling experiments 
Window3d < 44413mig_enc0_velow1.H min3=1 max3=3 | Transp plane=23 | Grey eout=1 gainpanel=a pclip=100 label1="z(km)" label2="h(km)" label3="x(km)" title="No combination of modeling experiments" > 44413mig_enc0_velow1.T

## No encoding of modeling experiments 
Window3d < 44413mig_enca_velow1.H min3=1 max3=3 | Transp plane=23  | Grey eout=1 gainpanel=a pclip=100 title="No encoding" label1="z(km)" label2="h(km)" label3="x(km)" > 44413mig_enca_velow1.T

## Conv.Random encoding of modeling experiments 
Window3d < 44413mig_encw_velow11.H min3=1 max3=3 | Transp plane=23  | Grey eout=1 gainpanel=a pclip=100 title="1 realization conv.code" label1="z(km)" label2="h(km)" label3="x(km)" > 44413mig_encw_velow11.T

## Conv.Random encoding of modeling experiments 
Window3d < 44413mig_encw_velow1.H min3=1 max3=3 | Transp plane=23  | Grey eout=1 gainpanel=a pclip=100 title="5 realizations conv.code" label1="z(km)" label2="h(km)" label3="x(km)" > 44413mig_encw_velow1.T

##Correl.random codes
#Window3d j1=5 < 44413_rnd_correl.H | Transp plane=12 reshape=1,3 | Graph min2=-600 max2=600 label1=lag label2=correlation title="Conventional random codes"| Tube&

## Correl.Gold codes
#Window3d j1=5 < 44413_correlations.H | Transp plane=12 reshape=1,3 | Graph min2=-600 max2=600 label1=lag label2=correlation title="Gold codes" | Tube&

## Gold sequential selection 
Window3d < 44413mig_gencw_velow4.H min3=1 max3=3 | Transp plane=23  | Grey eout=1 gainpanel=a pclip=100 title="Gold codes sequentially chosen" label1="z(km)" label2="h(km)" label3="x(km)" > 44413mig_gencw_velow4.T

## Gold random selection 1 realiz.
Window3d < 44413mig_gencw_velow31.H min3=1 max3=3 | Transp plane=23  | Grey eout=1 gainpanel=a pclip=100 title="Gold codes randomly chosen - 1 realiz." label1="z(km)" label2="h(km)" label3="x(km)" > 44413mig_gencw_velow31.T

## Gold random selection 5 realiz.
Window3d < 44413mig_gencw_velow3.H min3=1 max3=3 | Transp plane=23  | Grey eout=1 gainpanel=a pclip=100 title="Gold codes randomly chosen - 5 realiz." label1="z(km)" label2="h(km)" label3="x(km)" > 44413mig_gencw_velow3.T

## Gold optimal selection 1 realiz.
#Window3d < 44413mig_gencw_velow2.1.H min3=1 max3=3 | Transp plane=23  | Grey eout=1 gainpanel=a pclip=100 title="Gold codes optimally chosen" label1="z(km)" label2="h(km)" label3="x(km)" > 44413mig_gencw_velow2.1.T


rand.v: # Correlation functions of random codes.
	Window3d < 44413_rnd_correl.H n1=1 f1=0  | Graph label1="lag" label2="correlation" title="a)" wheretitle=t \
	min2=-100 max2=600 screenratio=4 screenht=4 screenwd=4 wantaxis1=y wantaxis2=y labelsz=3 titlesz=3 > /dev/null out=1.v
	Window3d < 44413_rnd_correl.H n1=1 f1=83 | Graph label1="lag" label2="correlation" title="b)" wheretitle=t \
	min2=-100 max2=600 screenratio=4 screenht=4 screenwd=4 wantaxis1=y wantaxis2=n labelsz=3 titlesz=3 > /dev/null out=2.v
	vppen < 1.v xcenter=-3.0 ycenter=-.5 vpstyle=n > junk1.v
	vppen < 2.v xcenter=-6.2 ycenter=-.5 vpstyle=n > junk2.v
	vp_Overlay junk1.v junk2.v > Fig/$@

prefpairs.v: # Correlation functions of preferred pairs of m-sequences.
	Window3d < 44413_prefpair.H n1=1 f1=10 | Graph label1="lag" label2="correlation" title=" " \
	screenratio=4 screenht=4 screenwd=4 wantaxis1=y wantaxis2=y labelsz=3 titlesz=3 > /dev/null out=Fig/$@

gold184.v: # Correlation functions of Gold codes.
	Window3d < 44413_correlations3.H n1=1 f1=0  | Graph label1="lag" label2="correlation" title="a)" wheretitle=t \
	min2=-100 max2=600 screenratio=4 screenht=4 screenwd=4 wantaxis1=y wantaxis2=y labelsz=3 titlesz=3 > /dev/null out=1.v
	Window3d < 44413_correlations3.H n1=1 f1=83 | Graph label1="lag" label2="correlation" title="b)" wheretitle=t \
	min2=-600 max2=100 screenratio=4 screenht=4 screenwd=4 wantaxis1=y wantaxis2=n labelsz=3 titlesz=3 > /dev/null out=2.v
	vppen < 1.v xcenter=-3.0 ycenter=-.5 vpstyle=n > junk1.v
	vppen < 2.v xcenter=-6.2 ycenter=-.5 vpstyle=n > junk2.v
	vp_Overlay junk1.v junk2.v > Fig/$@

perma.v: # Areal shot migration of synthesized data with no phase encoding applied.
	Window3d < 44413mig_enca_velow1.H min3=1 max3=3 | Transp plane=23 | Grey eout=1 gainpanel=a pclip=100 label1="z(km)" \
	label2="h(km)" label3="x(km)" title="" > a0.T
	< a0.T Cubeplot flat=1 frame1=86 frame2=50 frame3=10 wantcoordlabel=-1 point1=.7 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@

perm0.v: # Areal shot migration of synthesized data with no combination of the modeling experiments into super areal data.
	Window3d < 44413mig_enc0_velow1.H min3=1 max3=3 | Transp plane=23 | Grey eout=1 gainpanel=a pclip=100 label1="z(km)" \
	label2="h(km)" label3="x(km)" title="" > a0.T
	< a0.T Cubeplot flat=1 frame1=86 frame2=50 frame3=10 wantcoordlabel=-1 point1=.7 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@

conv1r.v: #Areal shot migration of one realization of synthesized data with conventional random phase encoding.}
	Window3d < 44413mig_encw_velow11.H min3=1 max3=3 | Transp plane=23 | Grey eout=1 gainpanel=a pclip=100 label1="z(km)" \
	label2="h(km)" label3="x(km)" title="" > a0.T
	< a0.T Cubeplot flat=1 frame1=86 frame2=50 frame3=10 wantcoordlabel=-1 point1=.7 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@

conv5r.v: #Areal shot migration of five realizations of synthesized data with conventional random phase encoding.}
	Window3d < 44413mig_encw_velow1.H min3=1 max3=3 | Transp plane=23 | Grey eout=1 gainpanel=a pclip=100 label1="z(km)" \
	label2="h(km)" label3="x(km)" title="" > a0.T
	< a0.T Cubeplot flat=1 frame1=86 frame2=50 frame3=10 wantcoordlabel=-1 point1=.7 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@

opti1.v: #Areal shot migration of five realizations of synthesized data with conventional random phase encoding.}\par
	Window3d < 44413mig_opti_velow1.H min3=1 max3=3 | Transp plane=23 | Grey eout=1 gainpanel=a pclip=100 label1="z(km)" \
	label2="h(km)" label3="x(km)" title="" > a0.T
	< a0.T Cubeplot flat=1 frame1=86 frame2=50 frame3=10 wantcoordlabel=-1 point1=.7 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@

opti5.v: #Areal shot migration of five realizations of synthesized data with conventional random phase encoding.}\par
	Window3d < 44413mig_opti_velow5.H min3=1 max3=3 | Transp plane=23 | Grey eout=1 gainpanel=a pclip=100 label1="z(km)" \
	label2="h(km)" label3="x(km)" title="" > a0.T
	< a0.T Cubeplot flat=1 frame1=86 frame2=50 frame3=10 wantcoordlabel=-1 point1=.7 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@

gold1x.v: #Areal shot migration of one realization of synthesized data with Gold phase encoding. Gold codes selected sequentially such that the cross-correlation peak separation is multiple of 10.}
	Window3d < 44413mig_gencw_velow4.H min3=1 max3=3 | Transp plane=23 | Grey eout=1 gainpanel=a pclip=100 label1="z(km)" \
	label2="h(km)" label3="x(km)" title="" > a0.T
	< a0.T Cubeplot flat=1 frame1=86 frame2=50 frame3=10 wantcoordlabel=-1 point1=.7 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@

gold1r.v: #Areal shot migration of one realization of synthesized data with Gold phase encoding. Gold codes are randomly selected.}
	Window3d < 44413mig_gencw_velow31.H min3=1 max3=3 | Transp plane=23 | Grey eout=1 gainpanel=a pclip=100 label1="z(km)" \
	label2="h(km)" label3="x(km)" title="" > a0.T
	< a0.T Cubeplot flat=1 frame1=86 frame2=50 frame3=10 wantcoordlabel=-1 point1=.7 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@

gold5r.v: #Areal shot migration of five realizations of synthesized data with Gold phase encoding. Gold codes are randomly selected.}
	Window3d < 44413mig_gencw_velow3.H min3=1 max3=3 | Transp plane=23 | Grey eout=1 gainpanel=a pclip=100 label1="z(km)" \
	label2="h(km)" label3="x(km)" title="" > a0.T
	< a0.T Cubeplot flat=1 frame1=86 frame2=50 frame3=10 wantcoordlabel=-1 point1=.7 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@

gold1o.v: #Areal shot migration of one realization of synthesized data with Gold phase encoding. Gold codes are selected such that the crosstalk is shifted out of the migration domain.}
	Window3d < 44413mig_gencw_velow2.1.H min3=1 max3=3 | Transp plane=23 | Grey eout=1 gainpanel=a pclip=100 label1="z(km)" \
	label2="h(km)" label3="x(km)" title="" > a0.T
	< a0.T Cubeplot flat=1 frame1=86 frame2=50 frame3=10 wantcoordlabel=-1 point1=.7 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@