###########################################
#
# Figures in the report
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=./fig/$@

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

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

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

arr11phz.v:
	Transp plane=12 reshape=2,4 < 43051areal_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="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=./fig/$@

ars11phz.v:
	Transp plane=12 reshape=2,4 < 43051areal_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="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=./fig/$@

arm11phz.v: # Migration of phs-encoding 1 in z
	Window3d min3=0 max3=4 < 43051mig_comb1_enc1z_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=./fig/$@

arm11phz4.v: # Summing 4 z-phs-encoding migration
	Window3d min3=0 max3=4 < 4305mig_comb1_enc1z_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=./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="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=./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="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=./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="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=./fig/$@

t51.v:
	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=" " \
	screenht=3.5 screenwd=6 > /dev/null out=./fig/$@
t41.v:
	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=" " \
	screenht=3.5 screenwd=6 > /dev/null out=./fig/$@

t11.v:
	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=" " \
	screenht=3.5 screenwd=6 > /dev/null out=./fig/$@


#########################################
#
# Arranging figures together
#
arsource.v:
	vppen < fig/arsou51.v xcenter=-0.8 ycenter=-1.5 vpstyle=n > junk1.v
	vppen < fig/arsou41.v xcenter=-4 ycenter=-1.5 vpstyle=n > junk2.v
	vppen < fig/arsou11.v xcenter=-7.2 ycenter=-1.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=.5 ycenter=-.5 vpstyle=n > junk2.v
	vppen < fig/arr11phz.v xcenter=-6.2 ycenter=-.5 vpstyle=n > junk3.v
	vp_Overlay junk2.v junk3.v > fig/$@

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

armigT.v:
	vppen < fig/t11.v xcenter=-0 ycenter=1. vpstyle=n > junk1.v
	vppen < fig/t41.v xcenter=-0 ycenter=-3. vpstyle=n > junk2.v
	vppen < fig/t51.v xcenter=-0 ycenter=-7. vpstyle=n > junk3.v
	vp_Overlay junk1.v junk2.v junk3.v > fig/$@