include ${SEPINC}/SEP.top
RESDIR=./Fig
R=$(RESDIR)
RESULTSNR=cost1 cost2
RESULTSER=
RESULTSCR=dat1200 mig.shot1200 mig.all mig.shot1200.3500 mig.shotsmall1 mig.shotsmall2 dmig1200 dmigall dmig1200 dmig1200.3500 dmigsmall dmigsmallsmall smallcomp smallsmallcomp directsum allcomp disccomp migsums

FIGER=$(addprefix $R/, $(addsuffix .pdf, $(RESULTSER)))
FIGNR=$(addprefix $R/, $(addsuffix .pdf, $(RESULTSNR)))
FIGCR=$(addprefix $R/, $(addsuffix .pdf, $(RESULTSCR)))
FIGS= ${FIGER} ${FIGCR} ${FIGNR} 
figpdf:${FIGS}
	@echo Done 
#LATOPTS= option=thesis #class=brad 
dn=/dev/null
GRD=Grey label1='time(s)' label2='x(m)' geophysics=1 title=' ' crowd=.89 > $(dn) labelsz=10 wantframe=0
GRP=Graph label1='time(s)' label2='x(m)' geophysics=1 title=' ' crowd=.89 > $(dn) wantframe=0 labelsz=10
GRZ=Grey label1='depth(m)' label2='x(m)' geophysics=1 title=' ' crowd=.89 > $(dn) tpow=1.5 labelsz=10 wantframe=0
vp=vppen vpstyle=n <
vpn=vppen vpstyle=n 
TR=Transp reshape=4,5,6 <
TXTDIR=../Txt
T=$(TXTDIR)
vpa=vp_annotate batch=y text=${TXTDIR}/ 
ysh=yshift=.12 
$R/allcomp.v: dmig.all.H mig.corrs.H isum.dmig.all.H mig.corrs.syn.shift.H
	$(TR) isum.dmig.all.H | $(GRZ) out=a.v 
	$(TR) mig.corrs.H     |Add scale=-1|$(GRZ) out=b.v
	$(TR) dmig.all.H      | $(GRZ) out=c.v
	$(TR) mig.corrs.syn.shift.H|$(GRZ) out=d.v
	$(vpa) text=$T/intext_a < a.v |\
	 $(vpn) > a2.v xsize=2. ysize=3. yshift=3.05 xshift=.05
	$(vpa) text=$T/intext_b < b.v |\
	 $(vpn) > b2.v xsize=2. ysize=3. yshift=3.05 xshift=2.05
	$(vpa) text=$T/intext_c < c.v |\
	 $(vpn) > c2.v xsize=2. ysize=3. $(ysh) xshift=.05
	$(vpa) text=$T/intext_d < d.v |\
	 $(vpn) > d2.v xsize=2. ysize=3. $(ysh) xshift=2.05
	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
$R/directsum.v: da.3500.H dasmallsmallrand.H dasmallrand.H
	< da.3500.H ${GRD} out=a.v
	< dasmallsmallrand.H ${GRD} out=b.v
	< dasmallrand.H ${GRD} out=c.v
	$(vpa) text=${TXTDIR}/intext_a < a.v |\
	 $(vpn)> a2.v xsize=2. ysize=3. $(ysh) xshift=.05
	$(vpa) text=${TXTDIR}/intext_b < b.v |\
	 $(vpn)> b2.v xsize=2. ysize=3. $(ysh) xshift=2.05
	$(vpa) text=${TXTDIR}/intext_c < c.v |\
	 $(vpn)> c2.v xsize=2. ysize=3. $(ysh) xshift=4.05
	vp_Overlay a2.v b2.v c2.v > $@
	rm -f a.v b.v c.v a2.v b2.v c2.v a1.v b1.v s1.v s2.v 

$R/synex.v: e3500.H dmig.e2000.H  dmig.e3500.H
	$(TR) dmig.e2000.H | $(GRZ) out=b.v
	$(TR) dmig.e3500.H | $(GRZ) out=c.v
	< e3500.H $(GRD) max1=3. out=a.v
	$(vp) a.v > a2.v xsize=2. ysize=3. $(ysh) xshift=.02
	$(vp) b.v > b2.v xsize=2. ysize=3. $(ysh) xshift=2.02
	$(vp) c.v > c2.v xsize=2. ysize=3. $(ysh) xshift=4.02
	vp_Overlay a2.v b2.v c2.v > $@
	rm -f a.v b.v c.v a2.v b2.v c2.v

$R/mig.all.v:corr.all.csum.H mig.corrs.H
	echo "(3200.,0.) (3200.,4.1)" > t.txt
	echo esize=0 n1=2 in=t.txt >> t.txt
	Sdd < t.txt esize=8 > s1.H
	< s1.H $(GRP) label1=' ' label2=' ' wantaxis=n out=s1.v min1=1200.\
	   max1=6780. min2=0. max2=4.1 plotfat=3 
	echo "(5580.,0.) (5580.,4.1)" > t.txt
	echo esize=0 n1=2 in=t.txt >> t.txt
	Sdd < t.txt esize=8 > s2.H
	rm -f t.txt
	< s2.H $(GRP) label1=' ' label2=' ' wantaxis=n out=s2.v min1=1200. \
	  max1=6780. min2=0. max2=4.1 plotfat=3 
	Rm s1.H s2.H
	< corr.all.csum.H Window3d n3=1 f3=100 | $(GRD) out=a1.v
	< corr.all.csum.H Window3d n3=1 f3=220 | $(GRD) out=b1.v
	vp_Overlay a1.v s1.v > a.v 
	vp_Overlay b1.v s2.v > b.v
	$(TR) mig.corrs.H  |Add scale=-1|$(GRZ) out=c1.v
	vp_Overlay c1.v s1.v s2.v > c.v
	$(vpa) text=${TXTDIR}/intext_a < a.v |\
	 $(vpn)> a2.v xsize=2. ysize=3. $(ysh) xshift=.05
	$(vpa) text=${TXTDIR}/intext_b < b.v |\
	 $(vpn)> b2.v xsize=2. ysize=3. $(ysh) xshift=2.05
	$(vpa) text=${TXTDIR}/intext_c < c.v |\
	 $(vpn)> c2.v xsize=2. ysize=3. $(ysh) xshift=4.05
	vp_Overlay a2.v b2.v c2.v > $@
	rm -f a.v b.v c.v a2.v b2.v c2.v a1.v b1.v s1.v s2.v 
$R/dat1200.v:shots.pad.H corr.1200.H
	< corr.1200.H Window3d n3=1        | $(GRD) out=a.v
	< corr.1200.H Window3d n3=1 f3=140 | $(GRD) out=c.v
	< corr.1200.H Window3d n3=1 f3=220 | $(GRD) out=e.v
	< shots.pad.H Window3d n3=1        | $(GRD) out=b.v
	< shots.pad.H Window3d n3=1 f3=140 | $(GRD) out=d.v
	< shots.pad.H Window3d n3=1 f3=220 | $(GRD) out=f.v
	vp_annotate < a.v batch=y text=${TXTDIR}/lrtext_a |\
	 $(vpn) > a2.v xsize=2. ysize=2.66 yshift=5.07 xshift=.05
	vp_annotate < b.v batch=y text=${TXTDIR}/lltext_b |\
	 $(vpn) > b2.v xsize=2. ysize=2.66 yshift=5.07 xshift=2.05
	vp_annotate < c.v batch=y text=${TXTDIR}/lrtext_c |\
	 $(vpn) > c2.v xsize=2. ysize=2.66 yshift=2.6 xshift=.05
	vp_annotate < d.v batch=y text=${TXTDIR}/lltext_d |\
	 $(vpn) > d2.v xsize=2. ysize=2.66 yshift=2.6 xshift=2.05
	vp_annotate < e.v batch=y text=${TXTDIR}/lrtext_e |\
	 $(vpn) > e2.v xsize=2. ysize=2.66 $(ysh) xshift=.05
	vp_annotate < f.v batch=y text=${TXTDIR}/lltext_f |\
	 $(vpn) > f2.v xsize=2. ysize=2.66 $(ysh) xshift=2.05
	vp_Overlay a2.v b2.v c2.v d2.v e2.v f2.v > $@
	rm -f a.v b.v c.v d.v e.v f.v a2.v b2.v c2.v d2.v e2.v f2.v
$R/mig.shot1200.v: mig.1200.H mig.shots1200.win4000.H mig.corr.1200win4000.1.H  mig.corr.1200win4000.20.H #mig.active.1200.H 
	$(TR) mig.1200.H              | $(GRZ) out=a.v pclip=96 
	$(TR) mig.shots1200.win4000.H | $(GRZ) out=b.v 
	$(TR) mig.corr.1200win4000.1.H| $(GRZ) out=c.v 
	$(TR)  mig.corr.1200win4000.20.H | $(GRZ) out=d.v 
	$(vpa) text=$T/intext_a < a.v |\
	$(vpa) text=$T/intext_a < a.v |\
	 $(vpn) > a2.v xsize=2. ysize=3. yshift=3.05 xshift=.05
	$(vpa) text=$T/intext_b < b.v |\
	 $(vpn) > b2.v xsize=2. ysize=3. yshift=3.05 xshift=2.05
	$(vpa) text=$T/intext_c < c.v |\
	 $(vpn) > c2.v xsize=2. ysize=3. $(ysh) xshift=.05
	$(vpa) text=$T/intext_d < d.v |\
	 $(vpn) > d2.v xsize=2. ysize=3. $(ysh) xshift=2.05
	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
$R/mig.shot1200.3500.v: mig.corr.1200.3500.rand.H  mig.corr.1200.3500.shift.H mig.corr.1200csum3500.H  mig.corr.1200csum3500rand.H
	$(TR) mig.corr.1200.3500.rand.H   | $(GRZ) out=b.v
	$(TR) mig.corr.1200.3500.shift.H  | $(GRZ) out=a.v
	$(TR) mig.corr.1200csum3500.H     | $(GRZ) out=c.v
	$(TR) mig.corr.1200csum3500rand.H | $(GRZ) out=d.v
	$(vpa) text=$T/intext_a < a.v |\
	 $(vpn) > a2.v xsize=2. ysize=3. yshift=3.05 xshift=.05
	$(vpa) text=$T/intext_b < b.v |\
	 $(vpn) > b2.v xsize=2. ysize=3. yshift=3.05 xshift=2.05
	$(vpa) text=$T/intext_c < c.v |\
	 $(vpn) > c2.v xsize=2. ysize=3. $(ysh) xshift=.05
	$(vpa) text=$T/intext_d < d.v |\
	 $(vpn) > d2.v xsize=2. ysize=3. $(ysh) xshift=2.05
	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
$R/mig.shotsmall1.v:mig.corr.small.H  mig.corr.smallrand.H 
	$(TR) mig.corr.small.H    | $(GRZ) out=a.v
	$(TR) mig.corr.smallrand.H   | $(GRZ) out=b.v
	$(vpa) text=$T/intext_a < a.v |\
	 $(vpn) > a2.v xsize=2. ysize=3. $(ysh) xshift=.05
	$(vpa) text=$T/intext_b < b.v |\
	 $(vpn) > b2.v xsize=2. ysize=3. $(ysh) xshift=2.05
	vp_Overlay a2.v b2.v > $@
	rm -f a.v b.v  a2.v b2.v 
$R/mig.shotsmall2.v:  mig.corr.smallsmall.H  mig.corr.smallsmallrand.H
	$(TR) mig.corr.smallsmall.H     | $(GRZ) out=a.v
	$(TR) mig.corr.smallsmallrand.H   | $(GRZ) out=b.v
	$(vpa) text=$T/intext_a < a.v |\
	 $(vpn) > a2.v xsize=2. ysize=3. $(ysh) xshift=.05
	$(vpa) text=$T/intext_b < b.v |\
	 $(vpn) > b2.v xsize=2. ysize=3. $(ysh) xshift=2.05
	vp_Overlay a2.v b2.v > $@
	rm -f a.v b.v  a2.v b2.v 
$R/dmigall.v:dmig.all.H smig.all.csum.H isum.dmig.all.H dmig.all.noshift.H dmiglong.all.H 
	$(TR) dmig.all.noshift.H | $(GRZ) out=a.v 
	$(TR) dmig.all.H | $(GRZ) out=b.v 
	$(TR) smig.all.csum.H | $(GRZ) out=d.v 
	$(TR) isum.dmig.all.H | $(GRZ) out=c.v 
	$(TR) dmiglong.all.H  | $(GRZ) out=e.v
	$(vpa) text=$T/intext_a < a.v |\
	 $(vpn) > a2.v xsize=2. ysize=3. yshift=3.05 xshift=.05
	$(vpa) text=$T/intext_b < b.v |\
	 $(vpn) > b2.v xsize=2. ysize=3. yshift=3.05 xshift=2.05
	$(vpa) text=$T/intext_c < c.v |\
	 $(vpn) > c2.v xsize=2. ysize=3. $(ysh) xshift=.05
	$(vpa) text=$T/intext_d < d.v |\
	 $(vpn) > d2.v xsize=2. ysize=3. $(ysh) xshift=2.05
	$(vpa) text=$T/intext_e < e.v |\
	 $(vpn) > e2.v xsize=2. ysize=6.1 yshift=.23 xshift=4.
	vp_Overlay a2.v b2.v c2.v d2.v e2.v > $@
	rm -f a.v b.v c.v d.v e.v a2.v b2.v c2.v d2.v e2.v


$R/dmig1200.v:dmig.1200.H mig.1200.H img.w.1500.H  
	$(TR)  dmig.1200.H | $(GRZ) out=a.v pclip=96 
	$(TR)  mig.1200.H    | $(GRZ) out=b.v 
	Window <img.w.1500.H | $(GRZ) out=c.v pclip=96 
	$(vpa) text=${TXTDIR}/trtext_a < a.v |\
	 $(vpn)> a2.v xsize=2. ysize=3. $(ysh) xshift=.05
	$(vpa) text=${TXTDIR}/trtext_b < b.v |\
	 $(vpn)> b2.v xsize=2. ysize=3. $(ysh) xshift=2.05
	$(vpa) text=${TXTDIR}/trtext_c < c.v |\
	 $(vpn)> c2.v xsize=2. ysize=3. $(ysh) xshift=4.05
	vp_Overlay a2.v b2.v c2.v > $@
	rm -f a.v b.v c.v a2.v b2.v c2.v

$R/dmig1200.3500.v: dmig.t.1200.3500.rand.H dmig.t.1200.3500.tsum.H dmig.t.1200.3500.isum.H mig.corr.1200.3500.rand.H #smig.1200csum3500rand.H
	$(TR) dmig.t.1200.3500.rand.H | $(GRZ) out=b.v 
	$(TR) dmig.t.1200.3500.tsum.H | $(GRZ) out=a.v 
	$(TR) dmig.t.1200.3500.isum.H | $(GRZ) out=c.v 
	$(TR) mig.corr.1200.3500.rand.H | Add scale=-1| $(GRZ)  out=d.v 
	$(vpa) text=$T/intext_a < a.v |\
	 $(vpn) > a2.v xsize=2. ysize=3. yshift=3.05 xshift=.05
	$(vpa) text=$T/intext_b < b.v |\
	 $(vpn) > b2.v xsize=2. ysize=3. yshift=3.05 xshift=2.05
	$(vpa) text=$T/intext_c < c.v |\
	 $(vpn) > c2.v xsize=2. ysize=3. $(ysh) xshift=.05
	$(vpa) text=$T/intext_d < d.v |\
	 $(vpn) > d2.v xsize=2. ysize=3. $(ysh) xshift=2.05
	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

#	$(TR) smig.1200csum3500rand.H | $(GRZ) out=d.v 

$R/smallcomp.v: dmig.small.H dmig.smallrand.H mig.corr.small.H  mig.corr.smallrand.H
	$(TR) mig.corr.small.H   | Add scale=-1 | $(GRZ) out=a.v
	$(TR) mig.corr.smallrand.H  | Add scale=-1 | $(GRZ) out=b.v
	$(TR) dmig.small.H |  $(GRZ) out=c.v 
	$(TR) dmig.smallrand.H | $(GRZ) out=d.v 
	$(vpa) text=$T/intext_a < a.v |\
	 $(vpn) > a2.v xsize=2. ysize=3. yshift=3.05 xshift=.05
	$(vpa) text=$T/intext_b < b.v |\
	 $(vpn) > b2.v xsize=2. ysize=3. yshift=3.05 xshift=2.05
	$(vpa) text=$T/intext_c < c.v |\
	 $(vpn) > c2.v xsize=2. ysize=3. $(ysh) xshift=.05
	$(vpa) text=$T/intext_d < d.v |\
	 $(vpn) > d2.v xsize=2. ysize=3. $(ysh) xshift=2.05
	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
$R/smallsmallcomp.v: dmig.smallsmall.H dmig.smallsmallrand.H mig.corr.smallsmall.H  mig.corr.smallsmallrand.H
	$(TR) mig.corr.smallsmall.H | Add scale=-1   | $(GRZ) out=a.v
	$(TR) mig.corr.smallsmallrand.H | Add scale=-1  | $(GRZ) out=b.v
	$(TR) dmig.smallsmall.H  | $(GRZ) out=c.v 
	$(TR) dmig.smallsmallrand.H | $(GRZ) out=d.v 
	$(vpa) text=$T/intext_a < a.v |\
	 $(vpn) > a2.v xsize=2. ysize=3. yshift=3.05 xshift=.05
	$(vpa) text=$T/intext_b < b.v |\
	 $(vpn) > b2.v xsize=2. ysize=3. yshift=3.05 xshift=2.05
	$(vpa) text=$T/intext_c < c.v |\
	 $(vpn) > c2.v xsize=2. ysize=3. $(ysh) xshift=.05
	$(vpa) text=$T/intext_d < d.v |\
	 $(vpn) > d2.v xsize=2. ysize=3. $(ysh) xshift=2.05
	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

$R/dmigsmall.v:dmig.small.H dmig.smallrand.H wmig.small.H wmig.smallrand.H
	$(TR) dmig.small.H  | $(GRZ) out=a.v 
	$(TR) dmig.smallrand.H | $(GRZ) out=b.v 
	$(TR) wmig.small.H | $(GRZ) out=c.v 
	$(TR) wmig.smallrand.H | $(GRZ) out=d.v 
	$(vpa) text=$T/intext_a < a.v |\
	 $(vpn) > a2.v xsize=2. ysize=3. yshift=3.05 xshift=.05
	$(vpa) text=$T/intext_b < b.v |\
	 $(vpn) > b2.v xsize=2. ysize=3. yshift=3.05 xshift=2.05
	$(vpa) text=$T/intext_c < c.v |\
	 $(vpn) > c2.v xsize=2. ysize=3. $(ysh) xshift=.05
	$(vpa) text=$T/intext_d < d.v |\
	 $(vpn) > d2.v xsize=2. ysize=3. $(ysh) xshift=2.05
	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
$R/dmigsmallsmall.v:dmig.smallsmall.H dmig.smallsmallrand.H wmig.smallsmall.H wmig.smallsmallrand.H
	$(TR) dmig.smallsmall.H | $(GRZ) out=a.v 
	$(TR) dmig.smallsmallrand.H | $(GRZ) out=b.v 
	$(TR) wmig.smallsmall.H | $(GRZ) out=c.v 
	$(TR) wmig.smallsmallrand.H | $(GRZ) out=d.v 
	$(vpa) text=$T/intext_a < a.v |\
	 $(vpn) > a2.v xsize=2. ysize=3. yshift=3.05 xshift=.05
	$(vpa) text=$T/intext_b < b.v |\
	 $(vpn) > b2.v xsize=2. ysize=3. yshift=3.05 xshift=2.05
	$(vpa) text=$T/intext_c < c.v |\
	 $(vpn) > c2.v xsize=2. ysize=3. $(ysh) xshift=.05
	$(vpa) text=$T/intext_d < d.v |\
	 $(vpn) > d2.v xsize=2. ysize=3. $(ysh) xshift=2.05
	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
cb=geophysics=1 title=' ' crowd=.89 wantcoordlabel=-1 label1='t(s)' label2='x(m)' label3='virtual shot(m)'

#	Grey eout=1 < corr.smallrand.H gainpanel=a | \
#	  Cubeplot frame1=80 frame2=140 frame3=140 out=a.v > $(dn) $(cb) 
#	Grey eout=1 < corr.smallsmallrand.H gainpanel=a | \
#	  Cubeplot frame1=80 frame2=140 frame3=140 out=c.v > $(dn) $(cb) 
$R/disccomp.v:dmig.smallsmallrand.H dmig.smallrand.H corr.smallrand.H corr.smallsmallrand.H
	echo "(4000.,0.) (4000.,4.1)" > t.txt
	echo esize=0 n1=2 in=t.txt >> t.txt
	Sdd < t.txt esize=8 > s.H
	rm -f t.txt
	< s.H $(GRP) label1=' ' label2=' ' wantaxis=n out=s.v min1=1200. \
	  max1=6780. min2=0. max2=4.1 plotfat=3 
	Rm s.H
	Window3d < corr.smallsmallrand.H n3=1 f3=140 |$(GRD) out=a1.v
	Window3d < corr.smallrand.H      n3=1 f3=140 |$(GRD) out=c1.v
	vp_Overlay a1.v s.v > a.v
	vp_Overlay c1.v s.v > c.v
	rm -f s.v a1.v c1.v
	$(TR) dmig.smallsmallrand.H | $(GRZ) out=b.v 
	$(TR) dmig.smallrand.H | $(GRZ) out=d.v 
	$(vpa) < a.v text=../Txt/intext_a | \
	 $(vpn) xsize=2. ysize=3. yshift=3. > a2.v xshift=.05
	$(vpa) text=$T/intext_b < b.v |\
	 $(vpn) > b2.v xsize=2. ysize=3. yshift=3. xshift=2.05
	$(vpa) < c.v text=../Txt/intext_c |\
	 $(vpn) xsize=2. ysize=3. yshift=0.13 > c2.v xshift=.05
	$(vpa) text=$T/intext_d < d.v |\
	 $(vpn) > d2.v xsize=2. ysize=3. yshift=0.13 xshift=2.05
	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 
$R/migsums.v: 5summedTWs.H dmig.isum.rand.H isum.dmig.fivesum.rand.H
	< 5summedTWs.H                  $(GRD) out=a.v 
	$(TR) isum.dmig.fivesum.rand.H |$(GRZ) out=b.v
	$(TR) dmig.isum.rand.H         |$(GRZ) out=c.v
	$(vpa) text=${TXTDIR}/intext_a < a.v |\
	 $(vpn)> a2.v xsize=2. ysize=3. $(ysh) xshift=.05
	$(vpa) text=${TXTDIR}/intext_b < b.v |\
	 $(vpn)> b2.v xsize=2. ysize=3. $(ysh) xshift=2.05
	$(vpa) text=${TXTDIR}/intext_c < c.v |\
	 $(vpn)> c2.v xsize=2. ysize=3. $(ysh) xshift=4.05
	vp_Overlay a2.v b2.v c2.v > $@
	rm -f a.v b.v c.v a2.v b2.v c2.v
include ${SEPINC}/SEP.bottom