44413vlow1.9.H:
	Spike n1=1000 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:
	< 44413vlow1.9.H Transp plane=12 reshape=1,3 > $@

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:
	@date
	/net/koko/claudio/Research/Others/ShotMig/ShotMig \
	D=9000data_prep1.H D_d=9000data_source1.H V=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

44413enc0:
	< 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 
	Rm junk.H

44413mig_enc0_velow1.H:
	< 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 ./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

44413encw11:
	< 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 
	Rm junk.H

44413mig_encw_velow11.H:
	< 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:
	< 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=1 correlations= 44413_rnd_correl.H > /dev/null 
	Rm junk.H

44413mig_encw_velow1.H:
	< 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

44413gencw2:
	< 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 
	Rm junk.H

44413mig_gencw_velow2.1.H:
	< 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

44413gencw31: # correct velocity areal shot random encode in w
	< 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 
	Rm junk.H

44413mig_gencw_velow31.H:
	< 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

44413gencw3:
	< 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=1 jshift=3 \
	OutCode=  Random_goldperfect.H prefpair=44413_prefpair.H correlations=44413_correlations3.H \
	> /dev/null 
	Rm junk.H

44413mig_gencw_velow3.H:
	< 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

44413gencw4:
	< 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 
	Rm junk.H

44413mig_gencw_velow4.H:
	< 44413areal_rec_gencw_velow4.H Transp plane=13 > 44413areal_rec_tempw4.H
	< 44413areal_sou_gencw_velow4.H Transp plane=13 > 44413areal_sou_tempw4.H
	< 44413areal_rec_tempw4.H InvVel/bin/OWmig2d_areal.x sou=44413areal_sou_tempw4.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*tempw4.H

44413opti1:
	< 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 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 nrefvel=1 \
	vsal=6 padtraces=41 encode=1 encw=1 eps=1 injection=1 nreal=1 node=12 > /dev/null 
	Rm junk.H

44413mig_opti_velow1.H:
	< 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*tempw.H

44413opti5:
	< 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 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=12 injection=1 nreal=5 > /dev/null 
	Rm junk.H

44413mig_opti_velow5.H:
	< 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*tempw.H

###############################################
#
# Scaling phases and migrating up to 3km
#
44413encw1.eps5:
	< 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.eps5.H sou_out=44413areal_sou_encw_velow1.eps5.H \
	OutCode= Random_conventional.H ngap=10 rndcte=1231 nrefvel=1 \
	vsal=6 padtraces=41 encode=1 random=1 encw=1 eps=5 injection=1 nreal=1 correlations= 44413_rnd_correl.H > /dev/null 

44413mig_encw_velow1.eps5.3km.H:
	< 44413areal_rec_encw_velow1.eps5.H Transp plane=13 > 44413areal_rec_tempw1.H
	< 44413areal_sou_encw_velow1.eps5.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=3. 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.eps3:
	< 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.eps3.H sou_out=44413areal_sou_encw_velow1.eps3.H \
	ngap=10 rndcte=1231 nrefvel=1 \
	vsal=6 padtraces=41 encode=1 random=1 encw=1 eps=3 injection=1 nreal=1 > /dev/null 

44413mig_encw_velow1.eps3.3km.H:
	< 44413areal_rec_encw_velow1.eps3.H Transp plane=13 > 44413areal_rec_tempw1.H
	< 44413areal_sou_encw_velow1.eps3.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=3. 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.eps55:
	< 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.eps55.H sou_out=44413areal_sou_encw_velow1.eps55.H \
	OutCode= Random_conventional.H ngap=10 rndcte=1231 nrefvel=1 \
	vsal=6 padtraces=41 encode=1 random=1 encw=1 eps=5 injection=1 nreal=5 correlations= 44413_rnd_correl.H > /dev/null 

44413mig_encw_velow1.eps55.3km.H:
	< 44413areal_rec_encw_velow1.eps55.H Transp plane=13 > 44413areal_rec_tempw1.H
	< 44413areal_sou_encw_velow1.eps55.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=3. 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_gencw_velow2.1.3km.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=3. 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

############
#
# Figures


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

####################
#
# PDF's

rand.pdf:
	pstexpen ${R}/rand.v ${R}/rand.ps color=n fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/rand.ps

prefpairs.pdf:
	pstexpen ${R}/prefpairs.v ${R}/prefpairs.ps color=n fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/prefpairs.ps

gold184.pdf:
	pstexpen ${R}/gold184.v ${R}/gold184.ps color=n fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/gold184.ps

perma.pdf:
	pstexpen ${R}/perma.v ${R}/perma.ps color=n fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/perma.ps

perm0.pdf:
	pstexpen ${R}/perm0.v ${R}/perm0.ps color=n fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/perm0.ps

conv1r.pdf:
	pstexpen ${R}/conv1r.v ${R}/conv1r.ps color=n fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/conv1r.ps

conv5r.pdf:
	pstexpen ${R}/conv5r.v ${R}/conv5r.ps color=n fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/conv5r.ps

opti1.pdf:
	pstexpen ${R}/opti1.v ${R}/opti1.ps color=n fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/opti1.ps

opti5.pdf:
	pstexpen ${R}/opti5.v ${R}/opti5.ps color=n fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/opti5.ps

gold1x.pdf:
	pstexpen ${R}/gold1x.v ${R}/gold1x.ps color=n fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/gold1x.ps

gold1r.pdf:
	pstexpen ${R}/gold1r.v ${R}/gold1r.ps color=n fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/gold1r.ps

gold5r.pdf:
	pstexpen ${R}/gold5r.v ${R}/gold5r.ps color=n fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/gold5r.ps

gold1o.pdf:
	pstexpen ${R}/gold1o.v ${R}/gold1o.ps color=n fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/gold1o.ps