##################################################
#
# 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=1 point2=0.7 flat=1 frame1=1 frame2=78 frame3=21 title=" " > \
	/dev/null out=./${RESDIR}/$@

#	| Grey eout=1 gainpanel=a pclip= 99 title=" " label1="depth(km)" label2="distance(km)" label3="offset(km)" \
#	labelsz=4 labelfat=n titlesz=4 wantaxis1=y wantaxis2=y wantaxis3=y > 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=1 point2=0.7 flat=1 frame1=1 frame2=78 frame3=21 title=" " \
#	screenratio=1 screenht=3 screenwd=6 labelsz=4 labelfat=n titlesz=4 wantaxis1=y wantaxis2=y wantaxis3=y >/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=1 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=1 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=1 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=1 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="d)" wheretitle=b label1="time(s)" label2="distance(km)" screenratio=4 screenht=4 screenwd=4 \
	wantaxis2=n 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="a) " label1="depth(km)" label2="distance(km)" label3="offset(km)" \
	labelsz=4 labelfat=n titlesz=4 > 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=1 point2=0.7 flat=1 frame1=1 frame2=78 frame3=20 title=" " \
	screenratio=1 screenht=3 screenwd=6 labelsz=4 labelfat=n titlesz=4 title="a) " >/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="angle(deg)"> 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=1 point2=0.7 flat=1 frame1=1 frame2=78 frame3=61 title=" " \
	screenratio=4 screenht=4 screenwd=4 >/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="a) " label1="depth(km)" label2="distance(km)" label3="offset(km)" \
	labelsz=4 labelfat=n titlesz=4 > 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=1 point2=0.7 flat=1 frame1=1 frame2=78 frame3=20 title=" " \
	screenratio=1 screenht=3 screenwd=6 labelsz=4 labelfat=n titlesz=4 title="a) " >/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=1 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="e)" wheretitle=b label1="time(s)" label2="distance(km)" screenratio=4 screenht=4 screenwd=4 \
	labelsz=4 labelfat=n titlesz=4 wantaxis1=n wantaxis2=n > /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="b) " label1="depth(km)" label2="distance(km)" label3="offset(km)" \
	labelsz=4 labelfat=n titlesz=4 wantaxis1=n > 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=1 point2=0.7 flat=1 frame1=1 frame2=78 frame3=20 title=" " wantaxis1=n \
	screenratio=1 screenht=3 screenwd=6 labelsz=4 labelfat=n titlesz=4 title="b) " >/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="angle(deg)"> 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=1 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="b) " label1="depth(km)" label2="distance(km)" label3="offset(km)" \
	labelsz=4 labelfat=n titlesz=4 wantaxis1=n > 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=1 point2=0.7 flat=1 frame1=1 frame2=78 frame3=20 title=" " wantaxis1=n \
	screenratio=1 screenht=3 screenwd=6 labelsz=4 labelfat=n titlesz=4 title="b) " >/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=1 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="angle(deg)"> 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="f)" wheretitle=b label1="time(s)" label2="distance(km)" screenratio=4 screenht=4 screenwd=4 \
	labelsz=4 labelfat=n titlesz=4 wantaxis1=n wantaxis2=n > /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="c) " label1="depth(km)" label2="distance(km)" label3="offset(km)" \
	labelsz=4 labelfat=n titlesz=4 > 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=1 point2=0.7 flat=1 frame1=1 frame2=78 frame3=20 title=" " \
	screenratio=1 screenht=3 screenwd=6 labelsz=4 labelfat=n titlesz=4 title="c) " >/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="angle(deg)"> 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=1 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="c) " label1="depth(km)" label2="distance(km)" label3="offset(km)" \
	labelsz=4 labelfat=n titlesz=4 > 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=1 point2=0.7 flat=1 frame1=1 frame2=78 frame3=20 title=" " \
	screenratio=1 screenht=3 screenwd=6 labelsz=4 labelfat=n titlesz=4 title="c) " >/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=1 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="angle(deg)"> 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 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}/$@

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

arm11phz.v: # Migration of phs-encoding 1 in z
	Window3d min3=0 max3=4 < 4300mig_comb1_enc1z_vel_rnd01_pad.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=1 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_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=1 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_01pad.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=1 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_pad.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=1 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}/$@

#
#
#
#
#  
#	Areal data rand phs-encoding 1 in z
arm01phz.v: # Migration of phs-encoding 1 in z
	Window3d min3=0 max3=4 < 4310mig_comb1_enc1z_vel_rnd0100.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}/$@

arm01phz15.v: # Summing 4 z-phs-encoding migration
	Window3d min3=0 max3=4 < 4310mig_comb1_enc1z_vel_rnd15real_pad.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}/$@

arm01phw.v: # Migration of phs-encoding 2 in w
	Window3d min3=0 max3=4 < 4320mig_comb1_enc1w_vel_rnd0100.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}/$@

arm01phw05.v: # Migration of phs-encoding 2 in w 44 areal shots
	Window3d min3=0 max3=4 < 4320mig_comb1_enc1w_vel_rnd5real_pad.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}/$@

#
# 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="angle(deg)"> 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="angle(deg)"> 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="angle(deg)"> 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="angle(deg)"> 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="angle(deg)"> 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}/$@

ardata.v:
	vppen < ${RESDIR}/arsou51.v xcenter=-0.8 ycenter=-4.5 vpstyle=n > junk1.v
	vppen < ${RESDIR}/arsou41.v xcenter=-4 ycenter=-4.5 vpstyle=n > junk2.v
	vppen < ${RESDIR}/arsou11.v xcenter=-7.2 ycenter=-4.5 vpstyle=n > junk3.v
	vppen < ${RESDIR}/arrec51.v xcenter=-0.8 ycenter=-.5 vpstyle=n > junk4.v
	vppen < ${RESDIR}/arrec41.v xcenter=-4 ycenter=-.5 vpstyle=n > junk5.v
	vppen < ${RESDIR}/arrec11.v xcenter=-7.2 ycenter=-.5 vpstyle=n > junk6.v
	vp_Overlay junk1.v junk2.v junk3.v junk4.v junk5.v junk6.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}/armig41W.v xcenter=-6 ycenter=-5 vpstyle=n > junk1.v
	vppen < ${RESDIR}/armig51W.v xcenter=0 ycenter=-5. vpstyle=n > junk2.v
	vppen < ${RESDIR}/armig11W.v xcenter=-3 ycenter=-1 vpstyle=n > junk3.v
	vp_Overlay junk1.v junk2.v junk3.v > ${RESDIR}/$@
	
armigT.v:
	vppen < ${RESDIR}/armig41T.v xcenter=-6 ycenter=-5 vpstyle=n > junk1.v
	vppen < ${RESDIR}/armig51T.v xcenter=0 ycenter=-5. vpstyle=n > junk2.v
	vppen < ${RESDIR}/armig11T.v xcenter=-3 ycenter=-1 vpstyle=n > junk3.v
	vp_Overlay junk1.v junk2.v junk3.v > ${RESDIR}/$@
#########################################
#
#
#
#
migar1scig.v:
	Window3d n3=1 min3=1.75 < 9200mig_comb1_enc0_vel.H | Scale scale_to=1 | Grey \
	title="a)" wheretitle=b label1="depth(m)" label2="offset(km)" screenratio=4 screenht=4 screenwd=4 \
	labelsz=4 labelfat=n titlesz=4  > /dev/null out=${RESDIR}/$@

migsp1scig.v:
	Window3d n1=1 min1=1.75 < 9200_mig_brad.H | Scale scale_to=1 | Transp | Grey \
	title="b)" wheretitle=b label1=" " label2="offset(km)" screenratio=4 screenht=4 screenwd=4 \
	labelsz=4 labelfat=n titlesz=4  > /dev/null out=${RESDIR}/$@ 

migscig1.v:
	vppen < ${RESDIR}/migar1scig.v xcenter=-2.5 ycenter=-.5 vpstyle=n > junk2.v
	vppen < ${RESDIR}/migsp1scig.v xcenter=-6.2 ycenter=-.5 vpstyle=n > junk3.v
	vp_Overlay junk2.v junk3.v > ${RESDIR}/$@

migar1acig.v:
	Window3d n3=1 min3=1.75 < 9200mig_comb1_enc0_vel.H | Scale scale_to=1 | Pad beg2=40 end2=40 end1=100 > xx.H
	< xx.H /net/koko/bob/SEP_BOB/bin/OFF2ANG amax=40 eps=20 maxmem=1000 verb=y > yy.H ; < yy.H Window3d max1=1 | Grey \
	title="a)" wheretitle=b label1="depth(m)" label2="angle(deg)" screenratio=4 screenht=4 screenwd=4 \
	labelsz=4 labelfat=n titlesz=4 > /dev/null out=${RESDIR}/$@

migsp1acig.v:
	Window3d n1=1 min1=1.75 < 9200_mig_brad.H | Scale scale_to=1 | Transp | Pad beg2=40 end2=40 end1=100 > xx.H
	< xx.H /net/koko/bob/SEP_BOB/bin/OFF2ANG amax=40 eps=20 maxmem=1000 verb=y > yy.H ; < yy.H Window3d max1=1 | Grey \
	title="b)" wheretitle=b label1=" " label2="angle(deg)" screenratio=4 screenht=4 screenwd=4 \
	labelsz=4 labelfat=n titlesz=4 > /dev/null out=${RESDIR}/$@ 

migacig1.v:
	vppen < ${RESDIR}/migar1acig.v xcenter=-2.5 ycenter=-.5 vpstyle=n > junk2.v
	vppen < ${RESDIR}/migsp1acig.v xcenter=-6.2 ycenter=-.5 vpstyle=n > junk3.v
	vp_Overlay junk2.v junk3.v > ${RESDIR}/$@

migar1rec.v:
	Transp plane=13 < 9201areal_rec_comb1_enc0_vel.H | Transf | Reverse which=1 opt=ignore | Window3d max1=1 | Grey \
	title="b)" wheretitle=b label1=" " label2="distance(km)" screenratio=4 screenht=4 screenwd=4 \
	labelsz=4 labelfat=n titlesz=4 > /dev/null out=${RESDIR}/$@ 

migar1sou.v:
	Transp plane=13 < 9201areal_sou_comb1_enc0_vel.H | Transf | Reverse which=1 opt=no | Window3d min1=-1 | Grey \
	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}/$@ 

migar1.v:
	vppen < ${RESDIR}/migar1sou.v xcenter=-2.5 ycenter=-.5 vpstyle=n > junk2.v
	vppen < ${RESDIR}/migar1rec.v xcenter=-6.2 ycenter=-.5 vpstyle=n > junk3.v
	vp_Overlay junk2.v junk3.v > ${RESDIR}/$@

#
# ... now with low velocity and dip
#
migar1sciglo.v:
	Window3d n3=1 min3=2.5 < 9100mig_comb1_enc0_vel.H | Scale scale_to=1 | Grey \
	title=" " label1="depth(km)" label2="distance(km)" label3="offset(km)" pclip=100 > /dev/null out=${RESDIR}/$@

migsp1sciglo.v:
	Window3d n1=1 min1=2.5 < 9100_mig_brad.H | Scale scale_to=1 | Transp | Grey \
	title=" " label1="depth(km)" label2="distance(km)" label3="offset(km)" pclip=100 > /dev/null out=${RESDIR}/$@ 

migar1aciglo.v:
	Window3d n3=1 min3=2.5 < 9100mig_comb1_enc0_vel.H | Scale scale_to=1 | Pad beg2=40 end2=40 end1=100 > xx.H
	< xx.H /net/koko/bob/SEP_BOB/bin/OFF2ANG amax=40 eps=20 maxmem=1000 verb=y > yy.H ; < yy.H Window3d max1=1 | Grey \
	title=" " label1="depth(km)" label2="distance(km)" label3="angle(deg)" pclip=100 > /dev/null out=${RESDIR}/$@

migar2aciglo.v:
	Window3d n3=1 min3=2.5 < 9101mig_comb1_enc0_vel.H | Scale scale_to=1 | Pad beg2=40 end2=40 end1=100 > xx.H
	< xx.H /net/koko/bob/SEP_BOB/bin/OFF2ANG amax=40 eps=20 maxmem=1000 verb=y > yy.H ; < yy.H Window3d max1=1 | Grey \
	title=" " label1="depth(km)" label2="distance(km)" label3="angle(deg)" pclip=100 > /dev/null out=${RESDIR}/$@

migsp1aciglo.v:
	Window3d n1=1 min1=2.5 < 9100_mig_brad.H | Scale scale_to=1 | Transp | Pad beg2=40 end2=40 end1=100 > xx.H
	< xx.H /net/koko/bob/SEP_BOB/bin/OFF2ANG amax=40 eps=20 maxmem=1000 verb=y > yy.H ; < yy.H Window3d max1=1 | Grey \
	title=" " label1="depth(km)" label2="distance(km)" label3="angle(deg)" pclip=100 > /dev/null out=${RESDIR}/$@ 

#
# Cubes low velocity and dip
#
sourot.v:
	< 9100_rot_sig_sou.H Transp plane=23 | Window3d min3=1 max3=3 min1=.5 max1=1 min2=-.4 n2=81| Grey gainpanel=a pclip=100 \
	title=" " label1="depth(km)" label3="distance(km)" label2="offset(km)" eout=1 | Transp plane=23 > x.T
	< x.T Window3d f1=0 f2=0 f3=0 f4=0 f5=0 n1=125 n2=200 n3=81 n4=1 n5=1  >a0.T
	< a0.T Cubeplot flat=1 frame1=53 frame2=100 frame3=40 wantcoordlabel=0 point1=.7 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@

recrot.v:
	< 9100_rot_sig_rec.H Transp plane=23 | Window3d min3=1 max3=3 min1=.5 max1=1 min2=-.4 n2=81| Grey gainpanel=a pclip=100 \
	title=" " label1="depth(km)" label3="distance(km)" label2="offset(km)" eout=1 | Transp plane=23 > x.T
	< x.T Window3d f1=0 f2=0 f3=0 f4=0 f5=0 n1=125 n2=200 n3=81 n4=1 n5=1  >a0.T
	< a0.T Cubeplot flat=1 frame1=53 frame2=100 frame3=40 wantcoordlabel=0 point1=.7 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@

sourot1.v:
	< 9100_rot_sig_sou.H Transp plane=23 | Window3d min3=1 max3=3 min1=.5 max1=1 min2=-.4 n2=81| Grey gainpanel=a pclip=100 \
	title="a) " label1="depth(km)" label3="distance(km)" label2="offset(km)" eout=1 | Transp plane=23 > x.T
	< x.T Window3d f1=0 f2=0 f3=0 f4=0 f5=0 n1=125 n2=200 n3=81 n4=1 n5=1  >a0.T
	< a0.T Cubeplot flat=1 frame1=53 frame2=100 frame3=40 wheretitle=t wantcoordlabel=-1 point1=1 point2=.7 title="a) " screenratio=4 screenht=4 screenwd=6 \
	labelsz=3 labelfat=n titlesz=4 > \
	/dev/null out=./${RESDIR}/$@

recrot1.v:
	< 9100_rot_sig_rec.H Transp plane=23 | Window3d min3=1 max3=3 min1=.5 max1=1 min2=-.4 n2=81| Grey gainpanel=a pclip=100 \
	title="b) " label1="depth(km)" label3="distance(km)" label2="offset(km)" eout=1 | Transp plane=23 > x.T
	< x.T Window3d f1=0 f2=0 f3=0 f4=0 f5=0 n1=125 n2=200 n3=81 n4=1 n5=1  >a0.T
	< a0.T Cubeplot flat=1 frame1=53 frame2=100 frame3=40 wheretitle=t wantcoordlabel=-1 point1=1 point2=.7 title="b) " screenratio=4 screenht=4 screenwd=6 \
	labelsz=3 labelfat=n titlesz=4 > \
	/dev/null out=./${RESDIR}/$@

datarot.v:
	vppen < ${RESDIR}/sourot1.v xcenter=-2.5 ycenter=-5.5 vpstyle=n > junk2.v
	vppen < ${RESDIR}/recrot1.v xcenter=-2.5 ycenter=-.5 vpstyle=n > junk3.v
	vp_Overlay junk2.v junk3.v > ${RESDIR}/$@

diprot.v:
	< 9102mig_comb1_enc0_vel.H Window3d min3=1 max3=3 min1=.5 max1=1 min2=-.4 n2=81| Grey gainpanel=a pclip=100 eout=1 \
	title=" " label1="depth(km)" label3="distance(km)" label2="offset(km)" | Transp plane=23 > 9102mig_comb1_enc0_vel.T
	< 9102mig_comb1_enc0_vel.T Window3d f1=0 f2=0 f3=0 f4=0 f5=0 n1=125 n2=200 n3=81 n4=1 n5=1  >a0.T
	< a0.T Cubeplot flat=1 frame1=70 frame2=91 frame3=40 wantcoordlabel=0 point1=.7 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@

diprot1.v:
	< 9102mig_comb1_enc0_vel.H Window3d min3=1 max3=3 min1=.5 max1=1 min2=-.4 n2=81| Grey gainpanel=a pclip=100 eout=1 \
	title=" " label1="depth(km)" label3="distance(km)" label2="offset(km)" | Transp plane=23 > 9102mig_comb1_enc0_vel.T
	< 9102mig_comb1_enc0_vel.T Window3d f1=0 f2=0 f3=0 f4=0 f5=0 n1=125 n2=200 n3=81 n4=1 n5=1  >a0.T
	< a0.T Cubeplot flat=1 frame1=70 frame2=91 frame3=40 wantcoordlabel=-1 point1=1 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@

dipnorot.v:
	< 9103mig_comb1_enc0_vel.H Window3d min3=1 max3=3 min1=.5 max1=1 min2=-.4 n2=81| Grey gainpanel=a pclip=100 eout=1 \
	title=" " label1="depth(km)" label3="distance(km)" label2="offset(km)" | Transp plane=23 > 9103mig_comb1_enc0_vel.T
	< 9103mig_comb1_enc0_vel.T Window3d f1=0 f2=0 f3=0 f4=0 f5=0 n1=125 n2=200 n3=81 n4=1 n5=1  >a0.T
	< a0.T Cubeplot flat=1 frame1=74 frame2=100 frame3=40 wantcoordlabel=0 point1=.7 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@

dipnorot1.v:
	< 9103mig_comb1_enc0_vel.H Window3d min3=1 max3=3 min1=.5 max1=1 min2=-.4 n2=81| Grey gainpanel=a pclip=100 eout=1 \
	title=" " label1="depth(km)" label3="distance(km)" label2="offset(km)" | Transp plane=23 > 9103mig_comb1_enc0_vel.T
	< 9103mig_comb1_enc0_vel.T Window3d f1=0 f2=0 f3=0 f4=0 f5=0 n1=125 n2=200 n3=81 n4=1 n5=1  >a0.T
	< a0.T Cubeplot flat=1 frame1=74 frame2=100 frame3=40 wantcoordlabel=-1 point1=1 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@

dipshtpr.v:
	< 9100_mig_brad.H Window3d min1=1 n1=201 min5=0.5 max1=1 min3=-.4 n3=81 | Transp plane=13 | Grey gainpanel=a pclip=100 \
	title=" " label1="depth(km)" label3="distance(km)" label2="offset(km)" eout=1 | Transp plane=23 > 9100_mig_brad.T
	< 9100_mig_brad.T  Window3d f1=0 f2=0 f3=0 f4=0 f5=0 n1=125 n2=200 n3=81 n4=1 n5=1  >a0.T
	< a0.T Cubeplot flat=1 frame1=53 frame2=100 frame3=40 wantcoordlabel=0 point1=.7 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@

dipshtpr1.v:
	< 9100_mig_brad.H Window3d min1=1 n1=201 min5=0.5 max1=1 min3=-.4 n3=81 | Transp plane=13 | Grey gainpanel=a pclip=100 \
	title=" " label1="depth(km)" label3="distance(km)" label2="offset(km)" eout=1 | Transp plane=23 > 9100_mig_brad.T
	< 9100_mig_brad.T  Window3d f1=0 f2=0 f3=0 f4=0 f5=0 n1=125 n2=200 n3=81 n4=1 n5=1  >a0.T
	< a0.T Cubeplot flat=1 frame1=53 frame2=100 frame3=40 wantcoordlabel=-1 point1=1 point2=0.7 title=" " screenratio=4 screenht=4 screenwd=6 \
	labelsz=3 labelfat=n titlesz=4> \
	/dev/null out=./${RESDIR}/$@

dipshtprV.v:
	< 9110_mig_brad.H Window3d min5=0.5 max5= 1 min1=1 max1=3 | Transp plane=13 | Grey eout=1 gainpanel=a pclip=100 \
	title=" " label1="depth(km)" label3="distance(km)" label2="offset(km)" | Transp plane=23 > 9110_mig_brad.T
	< 9110_mig_brad.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=124 n2=200 n3=81 n4=1 n5=1  >a0.T
	< a0.T Cubeplot flat=1 frame1=74 frame2=100 frame3=40 wantcoordlabel=0 point1=.7 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@

nodipshtpr.v:
	< 9200_mig_brad.H Window3d min1=1 n1=151 max1=1 min3=-.4 n3=81 | Transp plane=13 | Grey gainpanel=a pclip=100 \
	title=" " label1="depth(km)" label3="distance(km)" label2="offset(km)" eout=1 | Transp plane=23 > 9200_mig_brad.T
	< 9200_mig_brad.T  Window3d f1=0 f2=0 f3=0 f4=0 f5=0 n3=81 n4=1 n5=1  >a0.T
	< a0.T Cubeplot flat=1 frame1=178 frame2=100 frame3=40 wantcoordlabel=-1 point1=1 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@
 
#
# Additional figures
#
arm112eve.v: # Summing 4 z-phs-encoding migration
	Window3d min3=0 max3=4 < 1003mig_combine.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}/$@

xtalk:
	< 3.H Window3d j3=5 min2=2 max2=4 max3=0.235 | Window3d n3=1 f3=0 | Grey gainpanel=e pclip=100 \
	title=" " label1="depth(km)" label2="distance(km)" > /dev/null out=./${RESDIR}/1.v
	< 3.H Window3d j3=5 min2=2 max2=4 max3=0.235 | Window3d n3=1 f3=1 | Grey gainpanel=e pclip=100 \
	title=" " label1="depth(km)" label2="distance(km)" > /dev/null out=./${RESDIR}/2.v
	< 3.H Window3d j3=5 min2=2 max2=4 max3=0.235 | Window3d n3=1 f3=2 | Grey gainpanel=e pclip=100 \
	title=" " label1="depth(km)" label2="distance(km)" > /dev/null out=./${RESDIR}/3.v
	< 3.H Window3d j3=5 min2=2 max2=4 max3=0.235 | Window3d n3=1 f3=3 | Grey gainpanel=e pclip=100 \
	title=" " label1="depth(km)" label2="distance(km)" > /dev/null out=./${RESDIR}/4.v
	< 3.H Window3d j3=5 min2=2 max2=4 max3=0.235 | Window3d n3=1 f3=4 | Grey gainpanel=e pclip=100 \
	title=" " label1="depth(km)" label2="distance(km)" > /dev/null out=./${RESDIR}/5.v
	< 3.H Window3d j3=5 min2=2 max2=4 max3=0.235 | Window3d n3=1 f3=5 | Grey gainpanel=e pclip=100 \
	title=" " label1="depth(km)" label2="distance(km)" > /dev/null out=./${RESDIR}/6.v
	< 3.H Window3d j3=5 min2=2 max2=4 max3=0.235 | Window3d n3=1 f3=6 | Grey gainpanel=e pclip=100 \
	title=" " label1="depth(km)" label2="distance(km)" > /dev/null out=./${RESDIR}/7.v
	< 3.H Window3d j3=5 min2=2 max2=4 max3=0.235 | Window3d n3=1 f3=7 | Grey gainpanel=e pclip=100 \
	title=" " label1="depth(km)" label2="distance(km)" > /dev/null out=./${RESDIR}/8.v

# echo esize=8 >> 3893_GR.H ; < 3893_GR.H Transf maxsize=1000 | Reverse which=1 opt=ignore | Transp plane=12 reshape=1,3 | Pad beg2=20 | Window3d n2=200 > 2.H
# echo esize=8 >> 3893_GS.H ; < 3893_GS.H Transf maxsize=1000 | Reverse which=1 opt=ignore | Transp plane=12 reshape=1,3  > 1.H
# Math file1=1.H file2=2.H exp=file1+file2 > 4.H

xtalk2:
	< 4.H Window3d j3=5 min2=2 max2=4 max3=0.235 | Window3d n3=1 f3=0 | Grey gainpanel=e pclip=100 \
	title=" " label1="depth(km)" label2="distance(km)" > /dev/null out=./${RESDIR}/1a.v
	< 4.H Window3d j3=5 min2=2 max2=4 max3=0.235 | Window3d n3=1 f3=1 | Grey gainpanel=e pclip=100 \
	title=" " label1="depth(km)" label2="distance(km)" > /dev/null out=./${RESDIR}/2a.v
	< 4.H Window3d j3=5 min2=2 max2=4 max3=0.235 | Window3d n3=1 f3=2 | Grey gainpanel=e pclip=100 \
	title=" " label1="depth(km)" label2="distance(km)" > /dev/null out=./${RESDIR}/3a.v
	< 4.H Window3d j3=5 min2=2 max2=4 max3=0.235 | Window3d n3=1 f3=3 | Grey gainpanel=e pclip=100 \
	title=" " label1="depth(km)" label2="distance(km)" > /dev/null out=./${RESDIR}/4a.v
	< 4.H Window3d j3=5 min2=2 max2=4 max3=0.235 | Window3d n3=1 f3=4 | Grey gainpanel=e pclip=100 \
	title=" " label1="depth(km)" label2="distance(km)" > /dev/null out=./${RESDIR}/5a.v

xtalk3.v: # Migration of phs-encoding 2 in w 44 areal shots
	Window3d min3=0 max3=4 < 1003mig_combine.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}/$@

tmshv.v:




##############################################################
#
# Sigsbee examples
#
sigshp.v:
	Window3d < ../2Wpar/9903sigsb2a_mig.H min1=15000 max1=30000 min5=9000 max5=16025| Transp plane=13 | Grey eout=1 gainpanel=a pclip=99 \
	title=" " label1="depth(km)" label3="distance(km)" label2="offset(km)" screenratio=4 screenht=4 screenwd=4 \
	labelsz=2 labelfat=n titlesz=4 | Transp plane=23  > a0.T
	< a0.T Cubeplot flat=1 frame1=0 frame2=72 frame3=17 wantcoordlabel=-1 point1=1 point2=0.7 title=" " > crowd2=1 \
	/dev/null out=./${RESDIR}/$@
	tube ./${RESDIR}/$@ &
# no horizons
signicT.v:
	Window3d < 2000sig_mig_comb1_enc00_vel00.H n4=1 f4=14 min3=15000 max3=30000 min1=9000 | Grey eout=1 gainpanel=a pclip=99 \
	title=" " label1="depth(km)" label3="distance(km)" label2="offset(km)" | Transp plane=23  > a0.T
	< a0.T Cubeplot flat=1 frame1=0 frame2=72 frame3=17 wantcoordlabel=-1 point1=1 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@
	tube ./${RESDIR}/$@ &

signicW.v:
	Window3d < 2000sig_mig_comb1_enc01_vel00.H n4=1 f4=14 min3=15000 max3=30000 min1=9000 | Grey eout=1 gainpanel=a pclip=99 \
	title=" " label1="depth(km)" label3="distance(km)" label2="offset(km)" | Transp plane=23  > a0.T
	< a0.T Cubeplot flat=1 frame1=0 frame2=72 frame3=17 wantcoordlabel=-1 point1=1 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@
	tube ./${RESDIR}/$@ &

signphz.v:
	Window3d < 2001sig_mig_comb1_enc1z_vel00.H n4=1 f4=14 min3=15000 max3=30000 min1=9000 | Grey eout=1 gainpanel=a pclip=99 \
	title=" " label1="depth(km)" label3="distance(km)" label2="offset(km)" | Transp plane=23  > a0.T
	< a0.T Cubeplot flat=1 frame1=0 frame2=72 frame3=17 wantcoordlabel=-1 point1=1 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@
	tube ./${RESDIR}/$@ &

signphz4.v:
	Window3d < 2001sig_mig_comb1_enc1z_vel4sum.H min3=15000 max3= 30000 min1=9000 | Grey eout=1 gainpanel=a pclip=99 \
	title=" " label1="depth(km)" label3="distance(km)" label2="offset(km)" | Transp plane=23  > a0.T
	< a0.T Cubeplot flat=1 frame1=0 frame2=72 frame3=17 wantcoordlabel=-1 point1=1 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@
	tube ./${RESDIR}/$@ &

signphw.v:
	Window3d < 2002sig_mig_comb1_enc1w_vel00.H n4=1 f4=14 min3=15000 max3=30000 min1=9000 | Grey eout=1 gainpanel=a pclip=99 \
	title=" " label1="depth(km)" label3="distance(km)" label2="offset(km)" | Transp plane=23  > a0.T
	< a0.T Cubeplot flat=1 frame1=0 frame2=72 frame3=17 wantcoordlabel=-1 point1=1 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@
	tube ./${RESDIR}/$@ &

signphw4.v:
	Window3d < 2002sig_mig_comb1_enc1w_vel4sum.H min3=15000 max3= 30000 min1=9000 | Grey eout=1 gainpanel=a pclip=99 \
	title=" " label1="depth(km)" label3="distance(km)" label2="offset(km)" | Transp plane=23  > a0.T
	< a0.T Cubeplot flat=1 frame1=0 frame2=72 frame3=17 wantcoordlabel=-1 point1=1 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@
	tube ./${RESDIR}/$@ &

# Horizons
sigicT.v:
	Window3d < 2050sig_mig_comb1_enc00_vel00.H n4=1 f4=14 min3=15000 max3=30000 min1=9000 | Grey eout=1 gainpanel=a pclip=99.95 \
	title=" " label1="depth(km)" label3="distance(km)" label2="offset(km)" | Transp plane=23  > a0.T
	< a0.T Cubeplot flat=1 frame1=0 frame2=72 frame3=16 wantcoordlabel=-1 point1=1 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@
	tube ./${RESDIR}/$@ &

sigicW.v:
	Window3d < 2050sig_mig_comb1_enc01_vel00.H n4=1 f4=14 min3=15000 max3=30000 min1=9000 | Grey eout=1 gainpanel=a pclip=99.95 \
	title=" " label1="depth(km)" label3="distance(km)" label2="offset(km)" | Transp plane=23  > a0.T
	< a0.T Cubeplot flat=1 frame1=0 frame2=72 frame3=16 wantcoordlabel=-1 point1=1 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@
	tube ./${RESDIR}/$@ &

sigphz.v:
	Window3d < 2051sig_mig_comb1_enc1z_vel4.H n4=1 f4=14 min3=15000 max3=30000 min1=9000 | Grey eout=1 gainpanel=a pclip=99.95 \
	title=" " label1="depth(km)" label3="distance(km)" label2="offset(km)" | Transp plane=23  > a0.T
	< a0.T Cubeplot flat=1 frame1=0 frame2=72 frame3=16 wantcoordlabel=-1 point1=1 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@
	tube ./${RESDIR}/$@ &

sigphz4.v:
	Window3d < 2051sig_mig_comb1_enc1z_vel4sum.H n4=1 f4=14 min3=15000 max3=30000 min1=9000 | Grey eout=1 gainpanel=a pclip=99.95 \
	title=" " label1="depth(km)" label3="distance(km)" label2="offset(km)" | Transp plane=23  > a0.T
	< a0.T Cubeplot flat=1 frame1=0 frame2=72 frame3=16 wantcoordlabel=-1 point1=1 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@
	tube ./${RESDIR}/$@ &

sigphw.v:
	Window3d < 2052sig_mig_comb1_enc1w_vel00.H n4=1 f4=14 min3=15000 max3=30000 min1=9000 | Grey eout=1 gainpanel=a pclip=99.95 \
	title=" " label1="depth(km)" label3="distance(km)" label2="offset(km)" | Transp plane=23  > a0.T
	< a0.T Cubeplot flat=1 frame1=0 frame2=72 frame3=16 wantcoordlabel=-1 point1=1 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@
	tube ./${RESDIR}/$@ &

sigphw4.v:
	Window3d < 2002sig_mig_comb1_enc1w_vel4sum.H min3=15000 max3= 30000 min1=9000 | Grey eout=1 gainpanel=a pclip=99.95 \
	title=" " label1="depth(km)" label3="distance(km)" label2="offset(km)" | Transp plane=23  > a0.T
	< a0.T Cubeplot flat=1 frame1=0 frame2=72 frame3=16 wantcoordlabel=-1 point1=1 point2=0.7 title=" " > \
	/dev/null out=./${RESDIR}/$@
	tube ./${RESDIR}/$@ &

sigphz4ang.v: 
	Window3d < 2051sig_mig_comb1_enc1z_ang4sum.H min1=9000 | Grey eout=1 gainpanel=a pclip=99.95 \
	title=" " label1="depth(km)" label3="distance(km)" label2="angle(deg)" | Transp plane=23  > a0.T
	< a0.T Cubeplot flat=1 frame1=0 frame2=72 frame3=92 wantcoordlabel=-1 point1=1 point2=0.7 title=" " > \
	/dev/null out=temp.v
	vp_annotate < temp.v batch=y text=ang.txt | vppen vpstyle=n > ${RESDIR}/$@
	tube ./${RESDIR}/$@ &

sigshtang.v: 
	Window3d < 2000sig_angsht.H min1=9000 | Grey eout=1 gainpanel=a pclip=99.95 \
	title=" " label1="depth(km)" label3="distance(km)" label2="angle(deg)" | Transp plane=23  > a0.T
	< a0.T Cubeplot flat=1 frame1=0 frame2=72 frame3=92 wantcoordlabel=-1 point1=1 point2=0.7 title=" " > \
	/dev/null out=temp.v
	vp_annotate < temp.v batch=y text=ang.txt | vppen vpstyle=n > ${RESDIR}/$@
	tube ./${RESDIR}/$@ &

