##################################################
#
# Plane reflector 2km/s - Modeling and preparing
#

mig06.v: # With 201 shots separated 0.02km slower velocity
	 Window3d < 4006_mig_brad.H min1=0 max1=4 max5=1 | Transp plane=12 reshape=2,3 \
	| Grey eout=1 gainpanel=a pclip=99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=41 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=21 title=" " > \
	/dev/null out=./Fig/$@

mig06ang.v:
	 Window3d < 4006_mig_brad_ang.H |Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip=99 title=" " label1="z" label2="x" label3="ang"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=121 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=61 title=" " > \
	/dev/null out=./Fig/$@


mig08.v: # With 201 shots separated 0.02km correct velocity
	 Window3d < 4008_mig_brad.H min1=0 max1=4 max5=1 | Transp plane=12 reshape=2,3 \
	| Grey eout=1 gainpanel=a pclip=99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=41 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=21 title=" " > \
	/dev/null out=./Fig/$@

rot06rec.v:
	 Window3d < 4006_rot_sig_rec.H min2=0 max2=4 max1=1 \
	| Grey eout=1 gainpanel=a pclip=99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=41 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=21 title=" " > \
	/dev/null out=./Fig/$@

rot06sou.v:
	 Window3d < 4006_rot_sig_sou.H min2=0 max2=4 max1=1 \
	| Grey eout=1 gainpanel=a pclip=99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=41 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=21 title=" " > \
	/dev/null out=./Fig/$@

########################
#
# On 04/04/08, back to 
# exploding reflector
# propagation
#
########################
#
# Show crosstalk problem ngap=41
#
#	Areal data 0 - ngap=61
arrec51.v:
	Transp plane=12 reshape=2,4 < 3990areal_rec_comb1_enc0_vel.H | Transf mazsize=1000 \
	| Reverse which=1 opt=ignore | Window3d max1=1 n3=1 f3=18 | \
	Grey pclip=99.9 gainpanel=a title="a)" wheretitle=b label1="t" label2="x" screenratio=4 screenht=4 screenwd=4 \
	> /dev/null out=./Fig/$@

arsou51.v:
	Transp plane=12 reshape=2,4 < 3990areal_sou_comb1_enc0_vel.H | Transf mazsize=1000 \
	| Reverse which=1 opt=no | Window3d min1=-1 n3=1 f3=18 | \
	Grey pclip=99.9 gainpanel=a title="a)" wheretitle=b label1="t" label2="x" screenratio=4 screenht=4 screenwd=4 \
	 > /dev/null out=./Fig/$@


armig51T.v: #	Migration of Areal data 1 with IC in time
	Window3d min3=0 max3=4 < 3990mig_comb1_enc0_vel.H | Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip= 99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=41 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=20 title=" " > \
	/dev/null out=./Fig/$@

armig51Tang.v: #	Angle of areal data 1
	Window3d < 3990zang_comb1_enc0_vel.H |Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip=99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=121 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=61 title=" " > \
	/dev/null out=./Fig/$@

armig51W.v: # Migration of Areal data 1 with IC in freq
	Window3d min3=0 max3=4 < 3991mig_comb1_enc0_vel.H | Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip= 99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=41 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=20 title=" " > \
	/dev/null out=./Fig/$@


armig51Wang.v: #	Angle of Areal data 1 with IC in freq
	Window3d < 3992ang_comb1_enc0_vel.H |Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip=99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=121 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=61 title=" " > \
	/dev/null out=./Fig/$@
	
#	Areal data 1
arrec41.v:
	Transp plane=12 reshape=2,4 < 4100areal_rec_comb1_enc0_vel.H | Transf mazsize=1000 \
	| Reverse which=1 opt=ignore | Window3d max1=1 n3=1 f3=18 | \
	Grey pclip=99.9 gainpanel=a title="b)" wheretitle=b label1="t" label2="x" screenratio=4 screenht=4 screenwd=4 \
	wantaxis1=n wantaxis2=y > /dev/null out=./Fig/$@

arsou41.v:
	Transp plane=12 reshape=2,4 < 4100areal_sou_comb1_enc0_vel.H | Transf mazsize=1000 \
	| Reverse which=1 opt=no | Window3d min1=-1 n3=1 f3=18 | \
	Grey pclip=99.9 gainpanel=a title="b)" wheretitle=b label1="t" label2="x" screenratio=4 screenht=4 screenwd=4 \
	 wantaxis1=n wantaxis2=y > /dev/null out=./Fig/$@


armig41T.v: #	Migration of Areal data 1 with IC in time
	Window3d min3=0 max3=4 < 4100mig_comb1_enc0_vel.H | Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip= 99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=41 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=20 title=" " > \
	/dev/null out=./Fig/$@

armig41Tang.v: #	Angle of areal data 1
	Window3d < 4100zang_comb1_enc0_vel.H |Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip=99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=121 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=61 title=" " > \
	/dev/null out=./Fig/$@

armig41W.v: # Migration of Areal data 1 with IC in freq
	Window3d min3=0 max3=4 < 4101mig_comb1_enc0_vel.H | Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip= 99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=41 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=20 title=" " > \
	/dev/null out=./Fig/$@


armig41Wang.v: #	Angle of Areal data 1 with IC in freq
	Window3d < 4102ang_comb1_enc0_vel.H |Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip=99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=121 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=61 title=" " > \
	/dev/null out=./Fig/$@
	

armig1T.v: #	Migration of separate Areal data 1 with IC in time
	Window3d min3=0 max3=4 < 4105mig_comb0_enc0_vel.H | Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip= 99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=41 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=20 title=" " > \
	/dev/null out=./Fig/$@

armig1Tang.v: # Angle of separate Areal data 1 with IC in time
	Window3d < 4106ang_comb0_enc0_vel.H |Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip=99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=121 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=61 title=" " > \
	/dev/null out=./Fig/$@

#
# Finished showing crosstalk problem ngap=41
#
############################################################
#
# Show crosstalk problem ngap=11 
# Make the adcigs (04/06/08)
#
#

#	Areal data 2
arrec11.v:
	Transp plane=12 reshape=2,4 < 4200areal_rec_comb1_enc0_vel.H | Transf mazsize=1000 \
	| Reverse which=1 opt=ignore | Window3d max1=1 n3=1 f3=5 | \
	Grey pclip=99.9 gainpanel=a title="c)" wheretitle=b label1="t" label2="x" screenratio=4 screenht=4 screenwd=4 \
	wantaxis1=n wantaxis2=y > /dev/null out=./Fig/$@

arsou11.v:
	Transp plane=12 reshape=2,4 < 4200areal_sou_comb1_enc0_vel.H | Transf mazsize=1000 \
	| Reverse which=1 opt=no | Window3d min1=-1 n3=1 f3=5 | \
	Grey pclip=99.9 gainpanel=a title="c)" wheretitle=b label1="t" label2="x" screenratio=4 screenht=4 screenwd=4 \
	wantaxis1=n wantaxis2=y > /dev/null out=./Fig/$@

armig11T.v: #	Migration of Areal data 2 with IC in time
	Window3d min3=0 max3=4 < 4200mig_comb1_enc0_vel.H | Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip= 99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=41 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=20 title=" " > \
	/dev/null out=./Fig/$@

armig11Tang.v: #	Angle of areal data 2
	Window3d < 4200zang_comb1_enc0_vel.H |Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip=99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=121 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=61 title=" " > \
	/dev/null out=./Fig/$@

armig11W.v: # Migration of Areal data 2 with IC in freq
	Window3d min3=0 max3=4 < 4201mig_comb1_enc0_vel.H | Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip= 99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=41 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=20 title=" " > \
	/dev/null out=./Fig/$@


armig11Wang.v: #	Angle of Areal data 2 with IC in freq
	Window3d < 4202ang_comb1_enc0_vel.H |Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip=99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=121 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=61 title=" " > \
	/dev/null out=./Fig/$@
	

armig2T.v: # Migration of separate Areal data 2 with IC in time
	Window3d min3=0 max3=4 < 4205mig_comb0_enc0_vel.H | Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip= 99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=41 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=20 title=" " > \
	/dev/null out=./Fig/$@

armig2Tang.v: # Angle of separate Areal data 2 with IC in time
	Window3d < 4206ang_comb0_enc0_vel.H |Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip=99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=121 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=61 title=" " > \
	/dev/null out=./Fig/$@
#	
#
# Finished showing crosstalk problem ngap=11
#
############################################################
#
# Encoding from the model space for ngap=11
#
#
#  
#	Areal data rand phs-encoding 1 in z
arr11phz.v:
	Transp plane=12 reshape=2,4 < 43003areal_rec_comb1_enc1z_vel.H | Transf mazsize=1000 \
	| Reverse which=1 opt=ignore | Window3d max1=1 n3=1 f3=5 | \
	Grey pclip=99.9 gainpanel=a title=" " label1="t" label2="x" screenratio=4 screenht=4 screenwd=4 \
	wantaxis1=n wantaxis2=y > /dev/null out=./Fig/$@

ars11phz.v:
	Transp plane=12 reshape=2,4 < 43003areal_sou_comb1_enc1z_vel.H | Transf mazsize=1000 \
	| Reverse which=1 opt=no | Window3d min1=-1 n3=1 f3=5 | \
	Grey pclip=99.9 gainpanel=a title=" " label1="t" label2="x" screenratio=4 screenht=4 screenwd=4 \
	wantaxis1=y wantaxis2=y > /dev/null out=./Fig/$@

arm11phz.v: # Migration of phs-encoding 1 in z
	Window3d min3=0 max3=4 < 4300mig_comb1_enc1z_vel.H | Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip= 99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=41 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=20 title=" " > \
	/dev/null out=./Fig/$@

arm11phz4.v: # Summing 4 z-phs-encoding migration
	Window3d min3=0 max3=4 < 43004mig_sum_comb1_encz1_vel.H | Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip= 99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=41 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=20 title=" " > \
	/dev/null out=./Fig/$@

arm11z4ang.v: # Angle of 4 z-phs-encoding migration 
	Window3d < 43005mig_sum_comb1_encz1_vel_ang.H |Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip=99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=121 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=61 title=" " > \
	/dev/null out=./Fig/$@

arr11phw.v: # Areal data rand phs-encoding 2 in w
	Transp plane=12 reshape=2,4 < 4301areal_rec_comb1_enc1w_vel.H | Transf mazsize=1000 \
	| Reverse which=1 opt=ignore | Window3d max1=1 n3=1 f3=5 | \
	Grey pclip=99.9 gainpanel=a title=" " label1="t" label2="x" screenratio=4 screenht=4 screenwd=4 \
	wantaxis1=n wantaxis2=y > /dev/null out=./Fig/$@

ars11phw.v:
	Transp plane=12 reshape=2,4 < 4301areal_sou_comb1_enc1w_vel.H | Transf mazsize=1000 \
	| Reverse which=1 opt=no | Window3d min1=-1 n3=1 f3=5 | \
	Grey pclip=99.9 gainpanel=a title=" " label1="t" label2="x" screenratio=4 screenht=4 screenwd=4 \
	wantaxis1=y wantaxis2=y > /dev/null out=./Fig/$@

arm11phw.v: # Migration of phs-encoding 2 in w
	Window3d min3=0 max3=4 < 4301mig_comb1_enc1w_vel.H | Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip= 99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=41 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=20 title=" " > \
	/dev/null out=./Fig/$@

arm11wang.v: # Angle of phs-encoding 2 in w
	Window3d < 4302mig_comb1_encw1_vel.H |Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip=99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=121 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=61 title=" " > \
	/dev/null out=./Fig/$@

#
# Encoding in z reduces cross-talks but destroys the signal. 
# So, I used an epsilon to control the amount of randomness.
# The results are better. I'll try later linear encoding.
# Edit the flows.
#
###################################
#
# Slower velocity
#
# Make:
# 

#	Areal data slower velocity
arr11lo.v:
	Transp plane=12 reshape=2,4 < 4500areal_rec_comb1_enc0_velow.H | Transf mazsize=1000 \
	| Reverse which=1 opt=ignore | Window3d max1=1 n3=1 f3=5 | \
	Grey pclip=99.9 gainpanel=a title=" " label1="t" label2="x" > /dev/null out=./Fig/$@

ars11lo.v:
	Transp plane=12 reshape=2,4 < 4500areal_sou_comb1_enc0_velow.H | Transf mazsize=1000 \
	| Reverse which=1 opt=no | Window3d min1=-1 n3=1 f3=5 | \
	Grey pclip=99.9 gainpanel=a title=" " label1="t" label2="x" > /dev/null out=./Fig/$@

arm11Tlo.v: #	Migration of slower vel.data IC in time
	Window3d min3=0 max3=4 < 4500mig_comb1_enc0_velow.H | Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip= 99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=41 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=20 title=" " > \
	/dev/null out=./Fig/$@

arm11Tanglo.v: # Angle of slower velocity data IC time
	Window3d < 4500zang_comb1_enc0_velow.H |Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip=99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=121 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=61 title=" " > \
	/dev/null out=./Fig/$@

arm11Wlo.v: # Migration of Areal data 2 with IC in freq
	Window3d min3=0 max3=4 < 4501mig_comb1_enc0_velow.H | Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip= 99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=41 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=20 title=" " > \
	/dev/null out=./Fig/$@


arm11Wanglo.v: #	Angle of slower velocity data IC freq
	Window3d < 4502ang_comb1_enc0_velow.H |Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip=99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=121 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=61 title=" " > \
	/dev/null out=./Fig/$@
	

arm2Tlo.v: # Migration of separate slower vel.data IC in time
	Window3d min3=0 max3=4 < 4505mig_comb0_enc0_velow.H | Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip= 99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=41 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=20 title=" " > \
	/dev/null out=./Fig/$@

arm2Tanglo.v: # Angle of slower velocity data IC time
	Window3d < 4506ang_comb0_enc0_velow.H |Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip=99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=121 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=61 title=" " > \
	/dev/null out=./Fig/$@

############################################################
#
# Encoding from the model space for ngap=11
#
# Make:
#
#
arr11lophz.v: # Areal data rnd-phs encoding in z
	Transp plane=12 reshape=2,4 < 46003areal_rec_comb1_enc1z_velow.H | Transf mazsize=1000 \
	| Reverse which=1 opt=ignore | Window3d max1=1 n3=1 f3=5 | \
	Grey pclip=99.9 gainpanel=a title=" " label1="t" label2="x" > /dev/null out=./Fig/$@

ars11lophz.v:
	Transp plane=12 reshape=2,4 < 46003areal_sou_comb1_enc1z_velow.H | Transf mazsize=1000 \
	| Reverse which=1 opt=no | Window3d min1=-1 n3=1 f3=5 | \
	Grey pclip=99.9 gainpanel=a title=" " label1="t" label2="x" > /dev/null out=./Fig/$@

arm11phzlo.v: # Migration of phs-encoding 1 in z
	Window3d min3=0 max3=4 < 4600mig_comb1_enc1z_velow.H | Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip= 99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=41 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=20 title=" " > \
	/dev/null out=./Fig/$@

arm11phz4lo.v: # Summing 4 z-phs-encoding migration
	Window3d min3=0 max3=4 < 46004mig_sum_comb1_encz1_velow.H | Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip= 99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=41 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=20 title=" " > \
	/dev/null out=./Fig/$@

arm11z4anglo.v: # Angle of 4 z-phs-encoding migration 
	Window3d < 46005mig_sum_comb1_encz1_velow_ang.H |Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip=99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=121 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=61 title=" " > \
	/dev/null out=./Fig/$@

arm11phwlo.v: # Migration of phs-encoding 2 in w
	Window3d min3=0 max3=4 < 4601mig_comb1_enc1w_velow.H | Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip= 99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=41 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=20 title=" " > \
	/dev/null out=./Fig/$@

arm11wanglo.v: # Angle of phs-encoding 2 in w
	Window3d < 4602ang_comb1_encw1_velow.H |Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip=99 title=" " label1="z" label2="x" label3="h"> xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=200 n3=121 n4=1 n5=1  > a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=1 frame2=78 frame3=61 title=" " > \
	/dev/null out=./Fig/$@

#########################################
#
# Arranging Figures together
#
arsource.v:
	vppen < Fig/arsou51.v xcenter=-0.8 ycenter=-.5 vpstyle=n > junk1.v
	vppen < Fig/arsou41.v xcenter=-4 ycenter=-.5 vpstyle=n > junk2.v
	vppen < Fig/arsou11.v xcenter=-7.2 ycenter=-.5 vpstyle=n > junk3.v
	vp_Overlay junk1.v junk2.v junk3.v > Fig/$@

arrecev.v:
	vppen < Fig/arrec51.v xcenter=-0.8 ycenter=-.5 vpstyle=n > junk1.v
	vppen < Fig/arrec41.v xcenter=-4 ycenter=-.5 vpstyle=n > junk2.v
	vppen < Fig/arrec11.v xcenter=-7.2 ycenter=-.5 vpstyle=n > junk3.v
	vp_Overlay junk1.v junk2.v junk3.v > Fig/$@

arrphz.v:
	vppen < Fig/ars11phz.v xcenter=-4 ycenter=-.5 vpstyle=n > junk2.v
	vppen < Fig/arr11phz.v xcenter=-7.2 ycenter=-.5 vpstyle=n > junk3.v
	vp_Overlay junk2.v junk3.v > Fig/$@

arrphw.v:
	vppen < Fig/ars11phz.v xcenter=-4 ycenter=-.5 vpstyle=n > junk2.v
	vppen < Fig/arr11phz.v xcenter=-7.2 ycenter=-.5 vpstyle=n > junk3.v
	vp_Overlay junk2.v junk3.v > Fig/$@


