include $(SEPINC)/SEP.top
RESULTSNR=events water waterslice modldat modldat2 evnttest evnttest2 Fourth.mute synNMO oneside
RESULTSER=
RESULTSCR=
#LATOPTS= option=thesis
figs:
	$(MAKE) $(addprefix $R/, $(addsuffix .ps, $(RESULTSNR)))
RESDIR=Fig
R=$(RESDIR)
FIGER=$(addprefix $R/, $(addsuffix .pdf, $(RESULTSER)))
FIGNR=$(addprefix $R/, $(addsuffix .pdf, $(RESULTSNR)))
FIGCR=$(addprefix $R/, $(addsuffix .pdf, $(RESULTSCR)))
FIGS= ${FIGER} ${FIGCR} ${FIGNR}
FIGERV=$(addprefix $R/, $(addsuffix .v, $(RESULTSER)))
FIGNRV=$(addprefix $R/, $(addsuffix .v, $(RESULTSNR)))
FIGCRV=$(addprefix $R/, $(addsuffix .v, $(RESULTSCR)))
FIGSV= ${FIGERV} ${FIGCRV} ${FIGNRV}
tubev: #$(FIGSV)
	tube $(FIGSV)
figpdf:${FIGS}
	@echo Done
dn=/dev/null
vp=vppen vpstyle=n <
VP=Vppen vpstyle=n > $(dn) <
GRD=Grey label1='time(s)' label2='x(m)' geophysics=1 title=' ' crowd=.89 > $(dn)
GRZ=Grey label1='depth(m)' label2='x(m)' geophysics=1 title=' ' crowd=.89 > $(dn) tpow=0.5
TR=Transp reshape=4,5,6 <
TXTDIR=../Txt
T=$(TXTDIR)
vpa=vp_annotate batch=y
vp=vppen vpstyle=n < 
vpn=vppen vpstyle=n 

#$R/First.nmostack.v  $R/Fourth.nmostack.v 
RESDIR=./Fig
R=$(RESDIR)

#----------------------------------------------------------------
#   Only correlations going toward Armore production noise
nmo.H:Fourth.time.dec.cmpoffkeys.H vrms.H
	Nmo3d < $< velfile=vrms.H maxsize=1000 verb=1 \
	 off_keys=7,8 data_keys=5,6 vel_axes=2,3   >$@
Fourth.nmo.decon.sort.H:nmo.H
	Sort3d maxsize=300 verb=10 max_memory=300 nkeys=4 \#grid_fold=1 \
	key1=g_x  key2=g_y  key3=cmp_x  key4=cmp_y   \
	og1=5040. og2=5013. og3=5040. og4=5013.  \
	dg1=50.   dg2=300.  dg3=25.   dg4=150.   \
	ng1=221   ng2=14    ng3=442   ng4=27     \
	< $< > $@  

Fourth.onesided.nmo.H:  # only positive y
	Window_key < nmo.H maxsize=1000 synch=1 verb=1\
	 key1="offset_y" maxk1=10. > Fourth.onesided.nmo.H
Fourth.sort.onesided.nmo.H:Fourth.onesided.nmo.H
	Sort3d maxsize=300 verb=10 max_memory=300 nkeys=4 \
	 key1=g_x key2=g_y key3=cmp_x key4=cmp_y \
	 og1=5040. og2=5013. og3=5040. og4=5013. \
	 dg1=50. dg2=300. dg3=25. dg4=150. \
	 ng1=221 ng2=14 ng3=442 ng4=27 	< $< > $@
Fourth.stack.onesided.nmo.H:Fourth.sort.onesided.nmo.H
	Stack3d < $< > $@ maxsize=500 axes=4
Fourth.onesided-x.nmo.H:  # only positive x
	Window_key < nmo.H maxsize=1000 synch=1 verb=1\
	 key1="offset_x" maxk1=0. > Fourth.onesided-x.nmo.H
Fourth.sort.onesided-x.nmo.H:Fourth.onesided-x.nmo.H
	Sort3d maxsize=300 verb=10 max_memory=300 nkeys=4 \
	 key1=g_x key2=g_y key3=cmp_x key4=cmp_y \
	 og1=5040. og2=5013. og3=5040. og4=5013. \
	 dg1=50. dg2=300. dg3=25. dg4=150. \
	 ng1=221 ng2=14 ng3=442 ng4=27 	< $< > $@
Fourth.stack.onesided-x.nmo.H:Fourth.sort.onesided-x.nmo.H
	Stack3d < $< > $@ maxsize=500 axes=4

Fourth.onesided-xy.nmo.H:# only positive x&y
	Window_key < nmo.H maxsize=1000 synch=1 verb=1 \
	 key1="offset_y" maxk1=10. key2="offset_x" maxk2=0. \
	  > Fourth.onesided-xy.nmo.H
Fourth.sort.onesided-xy.nmo.H:Fourth.onesided-xy.nmo.H
	Sort3d maxsize=300 verb=10 max_memory=300 nkeys=4 \
	 key1=g_x key2=g_y key3=cmp_x key4=cmp_y \
	 og1=5040. og2=5013. og3=5040. og4=5013. \
	 dg1=50. dg2=300. dg3=25. dg4=150. \
	 ng1=221 ng2=14 ng3=442 ng4=27 < $< > $@
Fourth.stack.onesided-xy.nmo.H:Fourth.sort.onesided-xy.nmo.H
	Stack3d < $< > $@ maxsize=500 axes=4
Fourth.bp.stack.onesided-xy.nmo.T: Fourth.bp.stack.onesided-xy.nmo.H
	< $< Grey eout=1 gainpanel=a > $@
$R/oneside.v:# Fourth.bp.stack.onesided-xy.nmo.T
	< Fourth.bp.stack.onesided-xy.nmo.T Window3d f1=0 f2=0 f3=0 f4=0 \
	  f5=0 n1=720 n2=442 n3=27 n4=1 n5=1  >a0.T
	< a0.T Cubeplot flat=1 frame1=158 frame2=289 frame3=15 geophysics=1 \
	 title=' ' crowd=.89 wantcoordlabel=-1 point1=.75 point2=.75 \
	  label1='t(s)' label2='x(m)' label3='shot(m)' out=a.v > $(dn)
	 < a.v $(vpn) xsize=3. ysize=3. yshift=0.12 > $@
	rm -f a.v
#----------------------------------------------------------------
$R/synNMO.v: e3500.H stk.e1800.H stk.e3500.H
	< stk.e1800.H  $(GRD) out=b.v
	< stk.e3500.H  $(GRD) out=c.v
	< e3500.H $(GRD) max1=3. out=a.v
	$(vpa) < a.v  text=${TXTDIR}/lrtext_a |\
	  $(vpn) xsize=2. ysize=3. yshift=.1           > a2.v
	$(vpa) < b.v text=${TXTDIR}/lrtext_b |\
	  $(vpn) xsize=2. ysize=3. yshift=.1 xshift=2. > b2.v
	$(vpa) < c.v text=${TXTDIR}/lrtext_c |\
	  $(vpn) xsize=2. ysize=3. yshift=.1 xshift=4. > c2.v 
	vp_Overlay a2.v b2.v c2.v > $@
	rm -f a.v b.v c.v a2.v b2.v c2.v
$R/events.v:  $R/First.bp.nmostack.v  $R/Fourth.bp.nmostack.v $R/First.bp.yplus1.nmostack.v  $R/Fourth.bp.yplus1.nmostack.v
	$(vpa) < $R/First.bp.nmostack.v text=${TXTDIR}/lrcmp_a |\
	  $(vpn) xsize=3. ysize=2.6 yshift=2.7 > c.v
	$(vpa) < $R/First.bp.yplus1.nmostack.v text=${TXTDIR}/llcmp_b |\
	  $(vpn) xsize=3. ysize=2.6 yshift=2.7 > d.v xshift=3.
	$(vpa) < $R/Fourth.bp.nmostack.v text=${TXTDIR}/lrcmp_c |\
	  $(vpn) xsize=3. ysize=2.6 yshift=.1 > e.v
	$(vpa) < $R/Fourth.bp.yplus1.nmostack.v text=${TXTDIR}/llcmp_d |\
	 $(vpn) xsize=3. ysize=2.6 yshift=.1 > f.v xshift=3.
	vp_Overlay c.v d.v e.v f.v > $@
	rm -f  c.v d.v e.v f.v

$R/water.v:$R/First.bp.water.v $R/First.bp.yplus1.water.v
	$(vpa) < $R/First.bp.water.v text=${TXTDIR}/lrcmp_a |\
	  $(vpn) xsize=3. ysize=2.6 yshift=.1 > a.v
	$(vpa) < $R/First.bp.yplus1.water.v text=${TXTDIR}/llcmp_b |\
	  $(vpn) xsize=3. ysize=2.6 yshift=.1 > b.v xshift=3.
	vp_Overlay a.v b.v> $@
	rm -f  a.v b.v
$R/waterslice.v:$R/First.water.f1.88.v
	$(vp) $< xsize=2. ysize=4. yshift=.1 > $@
$R/shotmodel.v:shotmodel.H
	< $< Grey label1='time(s)' label2='trace#' geophysics=1 title=' ' crowd=.89 > $(dn) out=$@ max1=3.4
$R/modldat.v: $R/First.evntmute.data.v $R/First.evntwind.data.v $R/shotmodel.v
	$(vpa) < $R/shotmodel.v   text=${TXTDIR}/llshot_a |\
	  $(vpn) xsize=6. ysize=2.6 yshift=5.3 > a.v
	$(vpa) < $R/First.evntmute.data.v text=${TXTDIR}/llshot_b |\
	  $(vpn) xsize=6. ysize=2.6 yshift=2.7 > b.v
	$(vpa) < $R/First.evntwind.data.v text=${TXTDIR}/llshot_c |\
	  $(vpn) xsize=6. ysize=2.6 yshift=0.1 > c.v
	vp_Overlay a.v b.v c.v > $@
	rm -f  a.v b.v c.v 
$R/modldat2.v:$R/Fourth.evntmute.data.v $R/First.evntmute.data.v
	$(vpa) < $R/First.evntmute.data.v  text=${TXTDIR}/llshot_a |\
	  $(vpn) xsize=6. ysize=2.6 yshift=2.7> a.v
	$(vpa) < $R/Fourth.evntmute.data.v text=${TXTDIR}/llshot_b |\
	  $(vpn) xsize=6. ysize=2.6 yshift=.1 > b.v
	vp_Overlay a.v b.v > $@
	rm -f  a.v b.v
$R/evnttest.v:$R/First.eventwind.v $R/First.modl.v $R/First.evntmute.v $R/First.bp.nmostack.v
	$(vpa) < $R/First.bp.nmostack.v text=${TXTDIR}/lrcmp_a |\
	  $(vpn) xsize=3. ysize=2.6 yshift=2.7 > a.v	
	$(vpa) < $R/First.evntmute.v text=${TXTDIR}/llcmp_b |\
	  $(vpn) xsize=3. ysize=2.6 yshift=2.7 > b.v xshift=3.
	$(vpa) < $R/First.eventwind.v text=${TXTDIR}/lrcmp_c |\
	  $(vpn) xsize=3. ysize=2.6 yshift=0.1 > c.v 
	$(vpa) < $R/First.modl.v      text=${TXTDIR}/llcmp_d |\
	 $(vpn) xsize=3. ysize=2.6 yshift=0.1 > d.v xshift=3.
	vp_Overlay a.v b.v c.v d.v > $@
	rm -f  a.v b.v c.v d.v
$R/evnttest2.v:$R/Fourth.nmostack.v $R/Fourth.evntmute.v
	$(vpa) < $R/Fourth.nmostack.v text=${TXTDIR}/lrcmp_a |\
	  $(vpn) xsize=3. ysize=2.6 yshift=0.1 > a.v 
	$(vpa) < $R/Fourth.evntmute.v text=${TXTDIR}/llcmp_b |\
	  $(vpn) xsize=3. ysize=2.6 yshift=0.1 > b.v xshift=3.
	vp_Overlay a.v b.v> $@
	rm -f  a.v b.v
Fourth.bp.mute.stack.T : Fourth.bp.mute.stack.H
	Grey < $< eout=1 gainpanel=e pclip=100 > $@
pan.%.H:Fourth.bp.mute.stack.H
	<$< Window3d n3=8 f3=$*| Bandpass fhi=10 > $@
#|Pad n1out=1024 n2out=512  > $@

fxd.%.H: pan.%.H
	Fx2d < $< > $@
fk.%.H: pan.%.H
	< $< Rtoc | Ft3d sign1=1 sign2=-1 center2=1 > $@
$R/Fmute.slice.v:Fourth.bp.mute.stack.H #1669
	Window3d n1=1 min1=0.0917 < $< | \
	Grey geophysics=1 title=' ' crowd=.89 label1='x(m)' label2='y(m)'\
	  yreverse=n transp=n out=a.v > /dev/null
	$(vp) a.v xsize=3. ysize=1.3 yshift=.05 > $@
	rm -f a.v

$R/Fourth.mute.v:Fourth.bp.mute.stack.H
	< $< Grey eout=1 pclip=99 gainpanel=e > a0.T
	< $< Grey eout=1 pclip=95 tpow=1. gainpanel=e > a1.T
	Cubeplot < a0.T frame1=13 frame2=220 frame3=13 out=b.v > $(dn) \
	  geophysics=1 title=' ' crowd=.89 wantcoordlabel=-1 point1=.75 point2=.75 \
	  label1='t(s)' label2='x(m)' label3='shot(m)'
	Cubeplot < a0.T frame1=125 frame2=220   frame3=15   out=a.v > $(dn) \
	  geophysics=1 title=' ' crowd=.89 wantcoordlabel=-1 point1=.75 point2=.75 \
	  label1='t(s)' label2='x(m)' label3='shot(m)'
	cat ${TXTDIR}/cmpmute_a ${TXTDIR}/cube_a ${TXTDIR}/cmpevent > aa
	< a.v $(vpa) text=aa |\
	  $(vpn) xsize=3. ysize=3. yshift=3.1 > a2.v
	cat ${TXTDIR}/cmpmute_b ${TXTDIR}/cube_b ${TXTDIR}/cmpevent > bb
	< b.v $(vpa) text=bb |\
	  $(vpn) xsize=3. ysize=3. yshift=3.1 > b2.v xshift=3.
	Cubeplot < a1.T frame1=13 frame2=220 frame3=13 out=d.v > $(dn) \
	  geophysics=1 title=' ' crowd=.89 wantcoordlabel=-1 point1=.75 point2=.75 \
	  label1='t(s)' label2='x(m)' label3='shot(m)'
	Cubeplot < a1.T frame1=125 frame2=220   frame3=15   out=c.v > $(dn) \
	  geophysics=1 title=' ' crowd=.89 wantcoordlabel=-1 point1=.75 point2=.75 \
	  label1='t(s)' label2='x(m)' label3='shot(m)'
	cat ${TXTDIR}/cmpmute_a ${TXTDIR}/cube_c ${TXTDIR}/cmpevent > cc
	< c.v $(vpa) text=cc |\
	  $(vpn) xsize=3. ysize=3. yshift=0.12 > c2.v
	cat ${TXTDIR}/cmpmute_b ${TXTDIR}/cube_d ${TXTDIR}/cmpevent > dd
	< d.v $(vpa) text=dd |\
	  $(vpn) xsize=3. ysize=3. yshift=0.12 > d2.v xshift=3.
	vp_Overlay a2.v b2.v c2.v d2.v > $@
#	rm -f a.v b.v c.v d.v a2.v b2.v c2.v d2.v aa bb cc dd 
#	Rm a1.T a0.T

include $(SEPINC)/SEP.bottom