##################################################
#
# Plane reflector 2km/s - Modeling and preparing
#
# With 201 shots separated 0.02km correct velocity
mig08.v: # 4008_mig_brad.H
	 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/$@

########################
#
# On 04/04/08, back to 
# exploding reflector
# propagation
#
########################
#
# Show crosstalk problem ngap=41
#
#	Areal data 0 - ngap=61
arrec51.v: # 3990areal_rec_comb1_enc0_vel.H
	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 \
	labelsz=3.5 labelfat=n titlesz=3 titlesz=2> /dev/null out=./Fig/$@

arsou51.v: # 3990areal_sou_comb1_enc0_vel.H
	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 \
	labelsz=3.5 labelfat=n titlesz=3 > /dev/null out=./Fig/$@

#	Migration of Areal data 1 with IC in time
armig51T.v: # 3990mig_comb1_enc0_vel.H
	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/$@

# Migration of Areal data 1 with IC in freq
armig51W.v: # 3991mig_comb1_enc0_vel.H
	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/$@

#	Areal data 1
arrec41.v: # 4100areal_rec_comb1_enc0_vel.H
	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 labelsz=3.5 labelfat=n titlesz=3 > /dev/null out=./Fig/$@

arsou41.v: # 4100areal_sou_comb1_enc0_vel.H
	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 labelsz=3.5 labelfat=n titlesz=3 > /dev/null out=./Fig/$@

#	Migration of Areal data 1 with IC in time
armig41T.v: # 4100mig_comb1_enc0_vel.H
	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/$@

# Migration of Areal data 1 with IC in freq
armig41W.v: # 4101mig_comb1_enc0_vel.H
	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/$@

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

#	Areal data 2
arrec11.v: # 4200areal_rec_comb1_enc0_vel.H
	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 labelsz=3.5 labelfat=n titlesz=3 > /dev/null out=./Fig/$@

arsou11.v: # 4200areal_sou_comb1_enc0_vel.H
	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 labelsz=3.5 labelfat=n titlesz=3 > /dev/null out=./Fig/$@

#	Migration of Areal data 2 with IC in time
armig11T.v: # 4200mig_comb1_enc0_vel.H
	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/$@

# Migration of Areal data 2 with IC in freq
armig11W.v: # 4201mig_comb1_enc0_vel.H
	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/$@

#	
#
# Finished showing crosstalk problem ngap=11
#
############################################################
#
# Encoding from the model space for ngap=11
#
#
#  
#	Areal data rand phs-encoding 1 in z
arr11phz.v: # 4300areal_rec_comb1_enc1z_vel_rnd01_pad.H 
	Transp plane=12 reshape=2,4 < 4300areal_rec_comb1_enc1z_vel_rnd01.H | Transf mazsize=1000 \
	| Reverse which=1 opt=ignore | Window3d max1=1 n3=1 f3=5 | \
	Grey pclip=99.9 gainpanel=a title="b) " label1="t" label2="x" screenratio=4 screenht=4 screenwd=4 \
	wantaxis1=n wantaxis2=y labelsz=3.5 labelfat=n titlesz=3 > /dev/null out=./Fig/$@

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

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

# Summing 4 z-phs-encoding migration
arm11phz4.v: # 4300mig_comb1_enc1z_vel_rnd4real.H
	Window3d min3=0 max3=4 < 4300mig_comb1_enc1z_vel_rnd4real.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/$@

# Areal data rand phs-encoding in w
arr11phw.v: # 4301areal_rec_comb1_enc1w_vel_04pad.H
	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) " label1="t" label2="x" screenratio=4 screenht=4 screenwd=4 \
	wantaxis1=n wantaxis2=y labelsz=3.5 labelfat=n titlesz=3 > /dev/null out=./Fig/$@

ars11phw.v: # 4301areal_sou_comb1_enc1w_vel_04pad.H
	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) " label1="t" label2="x" screenratio=4 screenht=4 screenwd=4 \
	wantaxis1=y wantaxis2=y labelsz=3.5 labelfat=n titlesz=3 > /dev/null out=./Fig/$@

# Migration of phs-encoding in w
arm11phw.v: # 4301mig_comb1_enc1w_vel_04pad.H
	Window3d min3=0 max3=4 < 4301mig_comb1_enc1w_vel_01pad.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/$@

#
#########################################
#
# Arranging Figures together
#
arsource.v: arsou51.v arsou41.v arsou11.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: arrec51.v arrec41.v arrec11.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: arr11phz.v ars11phz.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: ars11phw.v arr11phw.v
	vppen < Fig/ars11phw.v xcenter=-4 ycenter=-.5 vpstyle=n > junk2.v
	vppen < Fig/arr11phw.v xcenter=-7.2 ycenter=-.5 vpstyle=n > junk3.v
	vp_Overlay junk2.v junk3.v > Fig/$@

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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