##################################################
#
# 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="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@

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=./${RESDIR}/$@


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="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@

rot06rec.v:
	 Window3d < 4006_rot_sig_rec.H min2=0 max2=4 max1=1 \
	| Grey eout=1 gainpanel=a pclip=99 title=" " label1="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@

rot06sou.v:
	 Window3d < 4006_rot_sig_sou.H min2=0 max2=4 max1=1 \
	| Grey eout=1 gainpanel=a pclip=99 title=" " label1="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@

########################
#
# On 04/04/08, back to 
# exploding reflector
# propagation
#
########################
#
# Only one CIG
#
arrec01.v:
	Transp plane=12 reshape=2,4 < 3890areal_rec_comb1_enc0_vel.H | Transf mazsize=1000 \
	| Reverse which=1 opt=ignore | Window3d max1=1 n3=1| \
	Grey pclip=99.9 gainpanel=a title="a)" wheretitle=b label1="time(s)" label2="distance(km)" screenratio=4 screenht=4 screenwd=4 \
	labelsz=4 labelfat=n titlesz=4 > /dev/null out=./${RESDIR}/$@

arsou01.v:
	Transp plane=12 reshape=2,4 < 3890areal_sou_comb1_enc0_vel.H | Transf mazsize=1000 \
	| Reverse which=1 opt=no | Window3d min1=-1 n3=1 | \
	Grey pclip=99.9 gainpanel=a title="a)" wheretitle=b label1="time(s)" label2="distance(km)" screenratio=4 screenht=4 screenwd=4 \
	 labelsz=4 labelfat=n titlesz=4 > /dev/null out=./${RESDIR}/$@


armig01T.v: #	Migration of Areal data 1 with IC in time
	Window3d < 3891mig_comb1_enc0_vel.H min3=2 max3=4| Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip=99.9 title=" " label1="depth(km)" label2="distance(km)" label3="offset(km)"> xxx.T
	 < xxx.T  Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=100 n3=41 n4=1 n5=1 > a0.T
	< a0.T Cubeplot  frame1=128 frame2=49 frame3=20 wantcoordlabel=-1 point1=.7 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@

armig01W.v: # Migration of Areal data 1 with IC in freq
	Window3d < 3890mig_comb1_enc0_vel.H min3=2 max3=4| Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip=99.9 title=" " label1="depth(km)" label2="distance(km)" label3="offset(km)"> xxx.T
	 < xxx.T  Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=100 n3=41 n4=1 n5=1 > a0.T
	< a0.T Cubeplot  frame1=128 frame2=49 frame3=20 wantcoordlabel=-1 point1=.7 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@

arrec03.v:
	Transp plane=12 reshape=2,4 < 3893areal_rec_comb1_enc0_vel.H | Transf mazsize=1000 \
	| Reverse which=1 opt=ignore | Window3d max1=1 n3=1| \
	Grey pclip=99.9 gainpanel=a title="a)" wheretitle=b label1="time(s)" label2="distance(km)" screenratio=4 screenht=4 screenwd=4 \
	labelsz=4 labelfat=n titlesz=4 > /dev/null out=./${RESDIR}/$@

arsou03.v:
	Transp plane=12 reshape=2,4 < 3893areal_sou_comb1_enc0_vel.H | Transf mazsize=1000 \
	| Reverse which=1 opt=no | Window3d min1=-1 n3=1 | \
	Grey pclip=99.9 gainpanel=a title="a)" wheretitle=b label1="time(s)" label2="distance(km)" screenratio=4 screenht=4 screenwd=4 \
	 labelsz=4 labelfat=n titlesz=4 > /dev/null out=./${RESDIR}/$@

armig03T.v: #	Migration of Areal data 1 with IC in time
	Window3d < 3894mig_comb1_enc0_vel.H min3=2 max3= 4| Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip=99.9 title=" " label1="depth(km)" label2="distance(km)" label3="offset(km)"> xxx.T
	 < xxx.T  Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=100 n3=41 n4=1 n5=1 > a0.T
	< a0.T Cubeplot  frame1=128 frame2=49 frame3=20 wantcoordlabel=-1 point1=.7 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@

armig03W.v: # Migration of Areal data 1 with IC in freq
	Window3d < 3893mig_comb1_enc0_vel.H min3=2 max3= 4| Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip=99.9 title=" " label1="depth(km)" label2="distance(km)" label3="offset(km)"> xxx.T
	 < xxx.T  Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=100 n3=41 n4=1 n5=1 > a0.T
	< a0.T Cubeplot  frame1=128 frame2=49 frame3=20 wantcoordlabel=-1 point1=.7 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@

#	< stdin Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=250 n2=150 n3=41 n4=1 n5=1  >a0.T
#	< a0.T Cubeplot flat=1 frame1=125 frame2=99 frame3=20


#
# 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="time(s)" label2="distance(km)" screenratio=4 screenht=4 screenwd=4 \
	labelsz=4 labelfat=n titlesz=4 > /dev/null out=./${RESDIR}/$@

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="time(s)" label2="distance(km)" screenratio=4 screenht=4 screenwd=4 \
	 labelsz=4 labelfat=n titlesz=4 > /dev/null out=./${RESDIR}/$@


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="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@

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="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@

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="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@


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="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@
	
#	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="time(s)" label2="distance(km)" screenratio=4 screenht=4 screenwd=4 \
	labelsz=4 labelfat=n titlesz=4 wantaxis1=n wantaxis2=y > /dev/null out=./${RESDIR}/$@

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="time(s)" label2="distance(km)" screenratio=4 screenht=4 screenwd=4 \
	 labelsz=4 labelfat=n titlesz=4 wantaxis1=n wantaxis2=y > /dev/null out=./${RESDIR}/$@


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="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@

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="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@

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="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@


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="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@
	

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="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@

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="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@

#
# Finished showing crosstalk problem ngap=41
#
############################################################
#
# Show crosstalk problem ngap=11 
# Make the adcigs (04/06/08)
#
# label1="depth(km)" label2="distance(km)" label3="offset(km)"

#	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="time(s)" label2="distance(km)" screenratio=4 screenht=4 screenwd=4 \
	labelsz=4 labelfat=n titlesz=4 wantaxis1=n wantaxis2=y > /dev/null out=./${RESDIR}/$@

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="time(s)" label2="distance(km)"  screenratio=4 screenht=4 screenwd=4 \
	labelsz=4 labelfat=n titlesz=4 wantaxis1=n wantaxis2=y > /dev/null out=./${RESDIR}/$@

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="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@

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="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@

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="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@


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="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@
	

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="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@

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="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@
#	
#
# 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 < 4300areal_rec_comb1_enc1z_vel_rnd05.H | Transf mazsize=1000 \
	| Reverse which=1 opt=ignore | Window3d n3=1 f3=5 | \
	Grey pclip=99.9 gainpanel=a title="b)" wheretitle=b label1="time(s)" label2="distance(km)"  screenratio=8 screenht=8 screenwd=8 \
	labelsz=4 labelfat=n titlesz=4 wantaxis1=y wantaxis2=y > /dev/null out=./${RESDIR}/$@

ars11phz.v:
	Transp plane=12 reshape=2,4 < 4300areal_sou_comb1_enc1z_vel_rnd05.H | Transf mazsize=1000 \
	| Reverse which=1 opt=no | Window3d min1=-3 n3=1 f3=5 | \
	Grey pclip=99.9 gainpanel=a title="a)" wheretitle=b label1="time(s)" label2="distance(km)"  screenratio=8 screenht=8 screenwd=8 \
	labelsz=4 labelfat=n titlesz=4 wantaxis1=y wantaxis2=y  > /dev/null out=./${RESDIR}/$@

arm11phz.v: # Migration of phs-encoding 1 in z
	Window3d min3=0 max3=4 < 4300mig_comb1_enc1z_vel_rnd05a.H | Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip= 99 title=" " label1="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@

arm11phz4.v: # Summing 4 z-phs-encoding migration
	Window3d min3=0 max3=4 < 4300mig_comb1_enc1z_vel_rnd4reala.H | Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip= 99 title=" " label1="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@

arm11z4ang.v: # Angle of 4 z-phs-encoding migration 
	Window3d < 4300ang_comb1_enc1z_vel_rndall.H |Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip=99 title=" " label1="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@

arr11phw.v: # Areal data rand phs-encoding 2 in w
	Transp plane=12 reshape=2,4 < 4301areal_rec_comb1_enc1w_vel_04.H | Transf mazsize=1000 \
	| Reverse which=1 opt=ignore | Window3d max1= 1 n3=1 f3=5 | \
	Grey pclip=99.9 gainpanel=a title="b)" wheretitle=b label1="time(s)" label2="distance(km)"  screenratio=8 screenht=8 screenwd=8 \
	labelsz=4 labelfat=n titlesz=4 wantaxis1=y wantaxis2=y  > /dev/null out=./${RESDIR}/$@

ars11phw.v:
	Transp plane=12 reshape=2,4 < 4301areal_sou_comb1_enc1w_vel_04.H | Transf mazsize=1000 \
	| Reverse which=1 opt=no | Window3d min1= -1 n3=1 f3=5 | \
	Grey pclip=99.9 gainpanel=a title="a)" wheretitle=b label1="time(s)" label2="distance(km)"  screenratio=8 screenht=8 screenwd=8 \
	labelsz=4 labelfat=n titlesz=4 wantaxis1=y wantaxis2=y  > /dev/null out=./${RESDIR}/$@

arm11phw.v: # Migration of phs-encoding 2 in w
	Window3d min3=0 max3=4 < 4301mig_comb1_enc1w_vel_04.H | Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip= 99 title=" " label1="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@

arm11phw4.v: # Migration of phs-encoding 2 in w 44 areal shots
	Window3d min3=0 max3=4 < 4301mig_comb1_enc1w_vel_rnd4real.H | Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip= 99 title=" " label1="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@

arm11wang.v: # Angle of phs-encoding 2 in w 44 areal shots
	Window3d < 4301ang_comb1_encw1_vel.H |Transp plane=23 \
	| Grey eout=1 gainpanel=a pclip=99 title=" " label1="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@

#
# 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="time(s)" label2="distance(km)"  > /dev/null out=./${RESDIR}/$@

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="time(s)" label2="distance(km)"  > /dev/null out=./${RESDIR}/$@

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="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@

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="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@

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="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@


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="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@
	

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="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@

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="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@

############################################################
#
# 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="time(s)" label2="distance(km)"  > /dev/null out=./${RESDIR}/$@

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="time(s)" label2="distance(km)"  > /dev/null out=./${RESDIR}/$@

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="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@

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="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@

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="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@

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="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@

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="depth(km)" label2="distance(km)" label3="offset(km)"> 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=./${RESDIR}/$@

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

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

arrphz.v:
	vppen < ${RESDIR}/ars11phz.v xcenter=.5 ycenter=-.5 vpstyle=n > junk2.v
	vppen < ${RESDIR}/arr11phz.v xcenter=-6.2 ycenter=-.5 vpstyle=n > junk3.v
	vp_Overlay junk2.v junk3.v > ${RESDIR}/$@

arrphw.v:
	vppen < ${RESDIR}/ars11phw.v xcenter=.5 ycenter=-.5 vpstyle=n > junk2.v
	vppen < ${RESDIR}/arr11phw.v xcenter=-6.2 ycenter=-.5 vpstyle=n > junk3.v
	vp_Overlay junk2.v junk3.v > ${RESDIR}/$@

armigW.v:
	vppen < ${RESDIR}/armig51W.v xcenter=-0 ycenter=-.8 vpstyle=n > junk1.v
	vppen < ${RESDIR}/armig51W.v xcenter=-0 ycenter=-4. vpstyle=n > junk2.v
	vppen < ${RESDIR}/armig51W.v xcenter=-0 ycenter=-7.2 vpstyle=n > junk3.v
	vp_Overlay junk1.v junk2.v junk3.v > ${RESDIR}/$@