##############################
#
# Figures
#
ZOtime.v:
	Window3d < /net/koko/claudio/Research/Prst_ExpRefl/OneWay/nfs-splitspread-shots.H n2=1 min2=0 min1=2 max1=6 min3=25000 max3=45000| \
	Grey pclip=90 title=" " label1="t(s)" label2="x(ft)" > /dev/null out=temp.v
	vp_annotate < temp.v batch=y text=ZOtime.txt | vppen vpstyle=n > ${R}/$@

sig2a01.v:
	< 00sig2a_mig.T Window3d min2=30500 max2=39600 | Cubeplot wantcoordlabel=-1 frame1=1 frame2=110 frame3=16 point1=1. point2=0.7 title=" " \
	> /dev/null out=temp.v
	vp_annotate < temp.v batch=y text=sig2a01.txt | vppen vpstyle=n > ${R}/$@
	

sig2b01.v:
	< 00sig2b_mig.T Window3d min2=30500 max2=39600 | Cubeplot wantcoordlabel=-1 frame1=1 frame2=110 frame3=16 point1=1. point2=0.7 title=" " \
	> /dev/null out=temp.v
	vp_annotate < temp.v batch=y text=sig2b01.txt | vppen vpstyle=n > ${R}/$@
	

sig2b02.v:
	< mig_Sis_pad_800_strati_sort1.T Cubeplot wantcoordlabel=-1 frame1=1 frame2=78 frame3=16 point1=1. point2=0.7 title=" " > /dev/null \
	out=${R}/$@
	

flt2b01.v:
	< prim_01_subt_newa1.T Cubeplot wantcoordlabel=-1 frame1=1 frame2=78 frame3=16 point1=1. point2=0.7 title=" " > /dev/null \
	out=${R}/$@
	

GREYD= Grey label1="z (ft) " label2="x (ft)" labelsz=4 titlesz=4 geophysics=1  crowd1=.75  crowd2=.75 title=" " min2=30600 max2=39600 > /dev/null 

diag-Sis.v: # LtL_Sis_az41_ax41_wave_marine_strati_t.H 
	Window3d < LtL_Sis_az41_ax41_wave_marine_strati_t.H n3=1 min3=0 n4=1 min4=0 | ${GREYD} newclip=1 bclip=-.5 eclip=2.5 polarity=-1 allpos=y out=${R}/$@ 

refle-new.v: 
	< sigsbee2b_reflection_coefficients_77.H Window3d min1=8000 n1=381 min2=30000 n2=133 j2=3| Transp reshape=2,4 > reflection_coefficients_target_conv.H 
	< reflection_coefficients_target_conv.H Window3d n3=381  > pp.H
	Math file1=pp.H exp="-1.*file1" > 1.H 
	< 1.H ${GREYD} title=" " pclip=95 o1num=9000  out=${R}/$@ 

mig-stack.v: # mig_Sis_pad_800_strati_nointerp_ang_stack.H 
	Cp mig_Sis_pad_800_strati_nointerp_ang_stack.H 1s.H
	echo o1=8000 >> 1s.H
	< 1s.H ${GREYD} title=" " pclip=87 o1num=9000  out=${R}/$@ 

mig-stack-filt.v: # mig_Sis_strati_ang_claudio.H #filter migration in angle
	<  mig_Sis_strati_ang_claudio.H Window3d | Transp | Stack3d  > 1s.H
	echo o1=8000 >> 1s.H
	< 1s.H ${GREYD} title=" " pclip=87 o1num=9000  out=${R}/$@ 
 
migoff1.v: # mig_Sis_pad_800_strati.H 
	< mig_Sis_pad_800_strati.H Scale scale_to=1.| Window3d min1=30000 n1=133 min5=8000 n5=381 min3=0 n3=1 |Transp > 1.H 
	< 1.H Grey min2=30600 max2=39600 label1="z (ft)" label2="x (ft)" title=" " clip=0.03 polarity=-1 screenratio=.69 o1num=9000  out=1.v geophysics=1 > /dev/null
	< 1.H Overlay type=polyline x=10.,10. y=1.3,8.95 out=11.v >/dev/null
	vp_Overlay 1.v 11.v > 111.v
	< mig_Sis_pad_800_strati.H Scale scale_to=1. | Window3d  n1=1 min1=36800 min5=8000 n5=381 min3=-600 n3=17 |Transp > 2.H
	< 2.H Grey label2="h (ft)" title=" " clip=0.03 crowd1=1. polarity=-1 screenratio=2. out=2.v geophysics=1> /dev/null
	< 2.H Overlay type=polyline x=2.55,2.55 y=1.3,8.95 out=22.v >/dev/null
	vp_Overlay 2.v 22.v > 222.v
	vp_SideBySideAniso 111.v 222.v | vppen vpstyle=n hshift=2 vshift=-.5> ${R}/$@ 

migoff1-filt.v: # ~claudio/Research/Prst_ExpRefl/RTM/prim_01_subt_sort.H#filter migration offset needs rule
	< prim_01_subt_sort.H Scale scale_to=1. | Window3d min1=30000 n1=133 min5=8000 n5=381 min3=0 n3=1 |Transp > 1.H 
	< 1.H Grey   min2=30600 max2=39600 label1="z (ft)" label2="x (ft)" title=" " clip=0.03 polarity=-1 screenratio=.69 o1num=9000  out=1.v geophysics=1 > /dev/null
	< 1.H Overlay type=polyline x=10.,10. y=1.3,8.95 out=11.v >/dev/null
	vp_Overlay 1.v 11.v > 111.v
	< prim_01_subt_sort.H Scale scale_to=1. | Window3d  n1=1 min1=36800 min5=8000 n5=381 min3=-600 n3=17 |Transp > 2.H
	< 2.H Grey label2="h (ft)" title=" " clip=0.03 crowd1=1. polarity=-1 screenratio=2. out=2.v geophysics=1> /dev/null
	< 2.H Overlay type=polyline x=2.55,2.55 y=1.3,8.95 out=22.v >/dev/null
	vp_Overlay 2.v 22.v > 222.v
	vp_SideBySideAniso 111.v 222.v | vppen vpstyle=n hshift=2 vshift=-.5> ${R}/$@ 

fragreya1=eout=1  gainpanel=a  label1="z (ft)" label2="x (ft)" label3="angle (degree)" title=" "  
frames1a= frame1=240 frame2=30 frame3=35 title=" "  crowd=.65 screenratio=1 point1=.75 point2=.75 polarity=-1 wantcoordlabel=-1 labelsz=3 titlesz=3

migang1.v: # mig_Sis_pad_800_strati_nointerp_ang_cow.H 
	< mig_Sis_pad_800_strati_nointerp_ang_cow.H Scale scale_to=1.| Window3d f1=9 n1=1  > 1.H
	< 1.H Grey min2=30600 max2=39600 label1="z (ft)" label2="x (ft)" title=" " clip=0.05    polarity=-1 screenratio=.69 o1num=9000  out=1.v geophysics=1 > /dev/null
	< 1.H Overlay type=polyline x=10.,10. y=1.3,8.95 out=11.v >/dev/null
	vp_Overlay 1.v 11.v > 111.v
	< mig_Sis_pad_800_strati_nointerp_ang_cow.H Scale scale_to=1. | Window3d  n4=1 min4=36800 |Transp > 2.H
	< 2.H Grey label2="angle (degree)" title=" " clip=0.05 crowd1=1. polarity=-1 screenratio=2. out=2.v  geophysics=1> /dev/null
	< 2.H Overlay type=polyline x=.8,.8 y=1.3,8.95 out=22.v >/dev/null
	vp_Overlay 2.v 22.v > 222.v
	vp_SideBySideAniso 111.v 222.v | vppen vpstyle=n hshift=2 vshift=-.5> ${R}/$@ 

migang1-filt.v: # mig_Sis_strati_ang_claudio.H 
	< mig_Sis_strati_ang_claudio.H Scale scale_to=1 | Window3d f1=9 n1=1  > 1.H
	< 1.H Grey  min2=30600 max2=39600 label1="z (ft)" label2="x (ft)" title=" " clip=0.05   polarity=-1 screenratio=.69 o1num=9000  out=1.v geophysics=1 > /dev/null
	< 1.H Overlay type=polyline x=10.,10. y=1.3,8.95 out=11.v >/dev/null
	vp_Overlay 1.v 11.v > 111.v
	< mig_Sis_strati_ang_claudio.H Scale scale_to=1 |Window3d n4=1 min4=36800 |Transp > 2.H
	< 2.H Grey label2="angle (degree)" title=" " clip=0.05 crowd1=1. polarity=-1 screenratio=2. out=2.v  geophysics=1> /dev/null
	< 2.H Overlay type=polyline x=.8,.8 y=1.3,8.95 out=22.v >/dev/null
	vp_Overlay 2.v 22.v > 222.v
	vp_SideBySideAniso 111.v 222.v | vppen vpstyle=n hshift=2 vshift=-.5> ${R}/$@

GREYK=Grey clip=10 allpos=y label1="kx (1/ft)" label2="kh (1/ft)" title=" " geophysics=1 min1=-0.03 max1=0.03 min2=-0.03 max2=0.03> /dev/null
filtered.v: # need rules
	Window3d min1=14150 n1=1 < filtered.H | Transp plane=12  >1.H
	< 1.H ${GREYK} out=${R}/$@

not-filtered.v: #need rules
	Window3d min1=14150 n1=1 < not_filtered.H | Transp plane=12 >1.H
	< 1.H ${GREYK} out=${R}/$@

diff-filtered.v: # need rules
	Window3d min1=14150 n1=1 < not_filtered.H | Transp plane=12 >1.H
	Window3d min1=14150 n1=1 < filtered.H | Transp plane=12 >2.H
	Math file1=1.H file2=2.H exp="file1-file2" > 3.H
	< 3.H Grey ${GREYK} out=${R}/$@
	
invz.v: 
	Window3d <mmov_new_sym.H f6=6 n6=1> 1s.H
	echo o1=8000 >> 1s.H
	< 1s.H ${GREYD} title=" " pclip=85 o1num=9000  out=${R}/$@ 

invz-filt.v: 
	Window3d <inv_Sis_az41ax41_newsolver_sym_strati_prec_claudio2.0.H> 1s.H
	echo o1=8000 >> 1s.H
	< 1s.H ${GREYD} title=" " pclip=87 o1num=9000  out=${R}/$@ 

resinvz.v: 
	Window3d <rmov0.H f6=9 n6=1> 1s.H
	echo o1=8000 >> 1s.H
	< 1s.H ${GREYD} title=" " pclip=100 o1num=9000  out=${R}/$@ 

resinvz-filt.v: 
	Window3d <rmov_claudio0.0.H f6=9 n6=1> 1s.H
	echo o1=8000 >> 1s.H
	< 1s.H ${GREYD} title=" " pclip=100 o1num=9000  out=${R}/$@

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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