# ----------------------------------------------------------------------
#  Makefile - Created by Brad Artman
#        Date:
# ----------------------------------------------------------------------
#       Usage:  
#  Parameters: 
# Description:
#             
# ----------------------------------------------------------------------
include ${SEPINC}/SEP.top
RESDIR=./Fig
R=$(RESDIR)
RESULTSNR=Vactive random first200 rz shortlong zoom1 maps areamap shal200 longact migdmig cmigact manydat manyimg vplane first200.2 manydatN manyimgN noeventdat noevent
RESULTSER=
RESULTSCR=
FIGER=$(addprefix $R/, $(addsuffix .pdf, $(RESULTSER)))
FIGNR=$(addprefix $R/, $(addsuffix .pdf, $(RESULTSNR)))
FIGCR=$(addprefix $R/, $(addsuffix .pdf, $(RESULTSCR)))
FIGS= ${FIGER} ${FIGCR} ${FIGNR} 

#----------------------------------------------------------------
#  Migration Rule:s
PASS_HYD=/mnt/sepdata/Volumes/5b/PassHydrophone
NUM=0 200 400 600 800 1000 1200 1400 1600 1800 2000 2200 2400

Fourth.traces.freq.dec.%.H: $(PASS_HYD)/Fourth.traces.freq.dec.H
	< $< Window3d max1=38 f3=$* maxsize=200 > $@
WINLIST=$(addprefix Fourth.traces.freq.dec.,$(addsuffix .H, $(NUM)))
bigmig:#${WINLIST}
	echo ${WINLIST}
PRG=/net/koko/brad/bin/LINUX/ShotMig
NB=2287
mig.Fourth%.H: Fourth.traces.freq.dec.%.H ValVel.H stations.most.H
	echo ' ' >> nohup.out
	echo ' ' >> nohup.out
	echo $< >> nohup.out
	echo ' ' >> nohup.out
	mkdir /scr1/brad/Fourth.$*.stat
	ln -s /scr1/brad/Fourth.$*.stat stat
	date
	nohup Parallel command=${PRG} antialias=n device=eth0 \
	files="$*.H,ValVel.H,stations.most.H,$@" taper=n xyreverse=y  \
	tags="D,V,HD,R" axis=3,3,2,6 usage="INPUT,INPUT,INPUT,OUTPUT" passive=y    \
	file_type="DISTRIBUTE,COPY,COPY,COPY" pverb=10 nblock=$(NB) decon=n \
	nx=460 ny=176 ox=4825. oy=4800. dx=25. dy=25. nzs=79 nws=252 #restart=1 
	rm -f stat
	date
	echo $< into $@ >> donelog
	xmessage -nearmouse -buttons DONE "Finished $@ " &
	echo "would you like to Rm $< "


#----------------------------------------------------------------
figpdf:${FIGS}
	@echo Done 
dn=/dev/null
GRD=Grey label1='time(s)' label2='x(m)' geophysics=1 title=' ' crowd=.89 > $(dn) labelsz=10 
GRZ=Grey label1='depth(m)' label2='x(m)' geophysics=1 title=' ' crowd=.89 > $(dn) tpow=1.5 labelsz=10 
vp=vppen vpstyle=n 
vpa=vp_annotate batch=y
TR=Transp reshape=4,5,6 <
CB=Cubeplot geophysics=1 title=' ' crowd=.89  > $(dn) labelsz=10 point1=.75 point2=.75 label1='z(m)' label2='x(m)' label3='y(m)'
$R/maps.v: $R/actgeom.v stations.cmplx.H
	$(vp) < $R/actgeom.v xsize=2.3 ysize=2.9 xshift=-.2 yshift=-.1 > a.v 
	< a.v pldb | sed -e 's/y(m)/N(m)/g' | sed -e 's/x(m)/E(m)/g' |plas > a2.v
	Graph crowd=.89 < stations.cmplx.H symbol='*' max1=9050. min1=4900 max2=16200 \
	 min2=4800 symbolsz=4 geophysics=1 title=' ' label1='y(m)' label2='x(m)' \
	 labelsz=13 > /dev/null out=b.v o2num=5000 d2num=2000. d1num=2000.
	$(vp) < b.v xsize=1.7 ysize=3.02 xshift=2.5 yshift=.17 > b2.v 
	< b2.v pldb | sed -e '11779c\m 2022 81' | sed -e '11780c\G 335 0 0 335' |\
	 sed -e '11755c\G 0 -335 335 0' | plas > b3.v
	vp_Overlay a2.v b3.v > $@ 
#	rm -f a.v b.v a2.v b2.v b3.v
$R/vplane.v: img.planes.T V.T
	< img.planes.T  Window3d  \
	   f1=18 f2=34 n1=360 n2=404 n3=176 > a0.T
	< V.T Window3d  Window3d  \
	   f1=18 f2=34 n1=360 n2=404 n3=176 > a1.T
	Merge_field < a0.T vel=a1.T scale=50 > a2.T
	< a2.T $(CB) frame1=316 frame2=279 frame3=65  out=a.v color=j
	$(vp) < a.v  xsize=3. ysize=3. yshift=0.22  xshift=.1 > $@
#	rm -f a.v; Rm  a0.T  a1.T  a2.T
#	$(vpa) < a.v text=${TXTDIR}/cube_a |\
#	  $(vp) xsize=3. ysize=3. yshift=0.22  xshift=.05 > $@
#	< a2.T $(CB) frame1=228 frame2=254 frame3=27 out=
#check:$R/vplane50.v $R/vplane60.v $R/vplane70.v $R/vplane80.v 
#	tube $^
$R/manydat.v: t.manyplanes.H t.manyevnts.H
	< t.manyplanes.H Window3d n2=300 | $(GRD) label2='trace #' out=a.v
	< t.manyevnts.H Window3d n2=1000 | $(GRD) label2='trace #' out=b.v
	$(vpa) < a.v text=${TXTDIR}/intext_a |\
	 $(vp) xsize=1.5 ysize=2.5 yshift=0.22 > a2.v xshift=.05
	$(vpa) < b.v text=${TXTDIR}/intext_b |\
	  $(vp) xsize=1.5 ysize=2.5 yshift=0.22 > b2.v xshift=1.6
	vp_Overlay a2.v b2.v > $@
	rm -f a.v b.v a2.v b2.v 
$R/manyimg.v: img.manyplanes.T img.manyevents.T
	< img.manyplanes.T Window3d  n1=450 |\
	  $(CB) frame1=100 frame1=316 frame2=300 frame3=65 out=a.v
	< img.manyevents.T Window3d  n1=450 |\
	  $(CB) frame1=100 frame1=316 frame2=300 frame3=65 out=b.v
	$(vpa) < a.v text=${TXTDIR}/cube_a |\
	 $(vp) xsize=3. ysize=3. yshift=0.22 > a2.v xshift=.1
	$(vpa) < b.v text=${TXTDIR}/cube_b |\
	  $(vp) xsize=3. ysize=3. yshift=0.22 > b2.v xshift=3.05
	vp_Overlay a2.v b2.v > $@
	rm -f a.v b.v a2.v b2.v 
$R/manydatN.v: t.noiseplanes.H t.noiseevnts.H
	< t.noiseplanes.H Window3d n2=300 | $(GRD) label2='trace #' out=a.v
	< t.noiseevnts.H Window3d n2=1000 | $(GRD) label2='trace #' out=b.v
	< noisetrace.H Window3d max1=8. |\
	 Graph crowd2=1. crowd1=.89 label1='t(s)' label2='amp' > /dev/null out=c.v \
	 labelsz=60 wantaxis2=n wantframe=n title=' '
	$(vpa) < a.v text=${TXTDIR}/intext_a |\
	 $(vp) xsize=1.5 ysize=2.5 yshift=0.72 > a2.v xshift=.05
	$(vpa) < b.v text=${TXTDIR}/intext_b |\
	  $(vp) xsize=1.5 ysize=2.5 yshift=0.72 > b2.v xshift=1.6
	$(vpa) < c.v text=${TXTDIR}/llgraph_c |\
	  $(vp) xsize=3.5 ysize=.6 yshift=0.26 > c2.v xshift=-.1
	vp_Overlay a2.v b2.v c2.v > $@
	rm -f a.v b.v c.v a2.v b2.v c2.v
$R/manyimgN.v: img.noiseplanes.T img.noiseevents.T
	< img.noiseplanes.T Window3d  n1=450 |\
	  $(CB) frame1=100 frame1=316 frame2=300 frame3=65 out=a.v
	< img.noiseevents.T Window3d  n1=450 |\
	  $(CB) frame1=100 frame1=316 frame2=300 frame3=65 out=b.v
	$(vpa) < a.v text=${TXTDIR}/cube_a |\
	 $(vp) xsize=3. ysize=3. yshift=0.22 > a2.v xshift=.1
	$(vpa) < b.v text=${TXTDIR}/cube_b |\
	  $(vp) xsize=3. ysize=3. yshift=0.22 > b2.v xshift=3.05
	vp_Overlay a2.v b2.v > $@
	rm -f a.v b.v a2.v b2.v 
$R/Vactive.v:img.active.tp.T
	echo o1=60 >> $<
	< $< $(CB) frame1=100 frame2=150 frame3=60  out=a.v 
	< $< $(CB) frame1=325 frame2=300 frame3=120 out=b.v 
	$(vpa) < a.v text=${TXTDIR}/cube_a |\
	  $(vp) xsize=3. ysize=3. yshift=0.22 > a2.v xshift=.1
	$(vpa) < b.v text=${TXTDIR}/cube_b |\
	  $(vp) xsize=3. ysize=3. yshift=0.22 > b2.v xshift=3.05
	vp_Overlay a2.v b2.v > $@
	rm -f a.v b.v a2.v b2.v 
	echo o1=62.5 >> $<
$R/noeventdat.v:regular.H novent.H
	< regular.H Window3d max1=2. n2=150 | $(GRD) out=a.v label2='trace #'
	< novent.H  Window3d max1=2. n2=150 | $(GRD) out=b.v label2='trace #'
	$(vpa) < a.v text=${TXTDIR}/trtext_a |\
	 $(vp) xsize=1.5 ysize=2.5 yshift=0.21 > a2.v xshift=.05
	$(vpa) < b.v text=${TXTDIR}/trtext_b |\
	  $(vp) xsize=1.5 ysize=2.5 yshift=0.21 > b2.v xshift=1.55
	vp_Overlay a2.v b2.v > $@
	rm -f a.v b.v a2.v b2.v 
#$R/noeventdat.ps:$R/noeventdat.v
#	pstexpen Fig/noeventdat.v Fig/noeventdat.ps \
#	 color=y fat=1 fatmult=1.5 invras=n

$R/noevent.v:img.gathers.T img.noevent.T
	< img.gathers.T Window3d n1=469 > a.T
	< a.T $(CB) frame1=231 frame2=107 frame3=85 out=a.v 
	Rm a.T
	< img.noevent.T $(CB) frame1=231 frame2=107 frame3=85 out=b.v 
	$(vpa) < a.v text=${TXTDIR}/cube_a |\
	  $(vp) xsize=3. ysize=3. yshift=0.22 > a2.v xshift=.1
	$(vpa) < b.v text=${TXTDIR}/cube_b |\
	  $(vp) xsize=3. ysize=3. yshift=0.22 > b2.v xshift=3.05
	vp_Overlay a2.v b2.v > $@
	rm -f a.v b.v a2.v b2.v 
	echo o1=62.5 >> $<

$R/random.v:img.mixed.T img.noise2.T
	< img.noise2.T $(CB) frame1=8 frame2=150 frame3=60 out=a.v 
	< img.mixed.T $(CB) frame1=8 frame2=300 frame3=120 out=b.v 
	$(vpa) <a.v text=${TXTDIR}/cube_a |\
	  $(vp) xsize=3. ysize=3. yshift=0.22 > a2.v xshift=.1
	$(vpa) <b.v text=${TXTDIR}/cube_b |\
	  $(vp) xsize=3. ysize=3. yshift=0.22 > b2.v xshift=3.05
	vp_Overlay a2.v b2.v > $@
	rm -f a.v b.v a2.v b2.v 

$R/rz.v: img.noise2.T #img.active.T
	< img.noise2.T Window3d  max1=600. > trz.T
	< trz.T $(CB) frame1=20 frame2=208 frame3=58 out=a.v 
	$(vp) a.v xsize=3. ysize=2. yshift=.1 > $@ #a2.v
#	< img.active.T Window3d  max1=600. > trzb.T
#	< trzb.T $(CB) frame1=20 frame2=208 frame3=58 out=b.v 
#	$(vp) b.v xsize=3. ysize=2. yshift=.1 > b2.v
	rm -f a.v b.v ;Rm trz.T # trzb.T
TXTDIR=../Txt
#$R/first200.v:img.active.tpow.462.T img.First.traces.freq.dec.0.tp.T 
#	< img.active.tpow.462.T Window3d  \
#	   f1=18 f2=34 n1=360 n2=404 n3=176 > a0.T
#	< img.First.traces.freq.dec.0.tp.T Window3d  \
#	   f1=18 f2=34 n1=360 n2=404 n3=176   > a1.T
#	< a0.T $(CB) frame1=316 frame2=279 frame3=65 out=a.v 
#	< a1.T $(CB) frame1=316 frame2=279 frame3=65 out=b.v 
#	Rm a0.T a1.T
#	< img.active.tpow.462.T Window3d  \
#	   max1=2000 f1=46 f2=8 n2=441 f3=27 n3=134 > a0.T
#	< img.First.traces.freq.dec.0.tp.T Window3d  \
#	   max1=2000 f1=46 f2=8 n2=441 f3=27 n3=134   > a1.T
#	< a0.T $(CB) frame1=240 frame2=108 frame3=81 out=c.v o1num=1000
#	< a1.T $(CB) frame1=240 frame2=108 frame3=81 out=d.v o1num=1000
#	cat ${TXTDIR}/cube_a ${TXTDIR}/first200oval1 > aa
#	< a.v $(vpa) text=aa  |\
#	  $(vp) xsize=3. ysize=3. > a2.v yshift=3.1 xshift=.1
#	cat ${TXTDIR}/cube_b ${TXTDIR}/first200oval1 > bb
#	< b.v $(vpa) text=bb  |\
#	  $(vp) xsize=3. ysize=3. > b2.v yshift=3.1 xshift=3.1
#	cat ${TXTDIR}/cube_c ${TXTDIR}/first200oval2 ${TXTDIR}/first200oval3 > cc
#	< c.v $(vpa) text=cc  |\
#	  $(vp) xsize=3. ysize=3. yshift=.22 > c2.v xshift=.1
#	cat ${TXTDIR}/cube_d ${TXTDIR}/first200oval2 ${TXTDIR}/first200oval3 > dd
#	< d.v $(vpa) text=dd  |\
#	  $(vp) xsize=3. ysize=3. yshift=.22 > d2.v xshift=3.1 
#	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 a0.T a1.T
$R/first200.v:img.active.tpow.462.T img.First.traces.freq.dec.0.tp.T img.manyevents.T
	< img.active.tpow.462.T Window3d  \
	   f1=18 f2=34 n1=360 n2=404 n3=176 > a0.T
	< img.First.traces.freq.dec.0.tp.T Window3d  \
	   f1=18 f2=34 n1=360 n2=404 n3=176   > a1.T
	< img.manyevents.T Window3d  \
	   f1=18 f2=34 n1=360 n2=404 n3=176   > a2.T
	< a0.T $(CB) frame1=316 frame2=279 frame3=65 out=a.v 
	< a1.T $(CB) frame1=316 frame2=279 frame3=65 out=c.v 
	< a2.T $(CB) frame1=316 frame2=279 frame3=65 out=e.v 
	Rm a0.T a1.T a2.T
	< img.active.tpow.462.T Window3d  \
	   max1=2000 f1=46 f2=8 n2=441 f3=27 n3=134 > a0.T
	< img.First.traces.freq.dec.0.tp.T Window3d  \
	   max1=2000 f1=46 f2=8 n2=441 f3=27 n3=134   > a1.T
	< img.manyevents.T Window3d  \
	   max1=2000 f1=46 f2=8 n2=441 f3=27 n3=134   > a2.T
	< a0.T $(CB) frame1=240 frame2=108 frame3=81 out=b.v o1num=1000
	< a1.T $(CB) frame1=240 frame2=108 frame3=81 out=d.v o1num=1000
	< a2.T $(CB) frame1=240 frame2=108 frame3=81 out=f.v o1num=1000
	Rm a0.T a1.T a2.T
	cat ${TXTDIR}/cube_a ${TXTDIR}/first200oval1 > aa
	cat ${TXTDIR}/cube_b ${TXTDIR}/first200oval2 ${TXTDIR}/first200oval3 > bb
	cat ${TXTDIR}/cube_c ${TXTDIR}/first200oval1 > cc
	cat ${TXTDIR}/cube_d ${TXTDIR}/first200oval2 ${TXTDIR}/first200oval3 > dd
	cat ${TXTDIR}/cube_e ${TXTDIR}/first200oval1 > ee
	cat ${TXTDIR}/cube_f ${TXTDIR}/first200oval2 ${TXTDIR}/first200oval3 > ff
	< a.v $(vpa) text=aa  |\
	  $(vp) xsize=3. ysize=3. > a2.v yshift=6. xshift=.1
	< b.v $(vpa) text=bb  |\
	  $(vp) xsize=3. ysize=3. > b2.v yshift=6. xshift=3.1
	< c.v $(vpa) text=cc  |\
	  $(vp) xsize=3. ysize=3. > c2.v yshift=3.1 xshift=.1
	< d.v $(vpa) text=dd  |\
	  $(vp) xsize=3. ysize=3. > d2.v yshift=3.1 xshift=3.1 
	< e.v $(vpa) text=ee  |\
	  $(vp) xsize=3. ysize=3. > e2.v yshift=.22 xshift=.1 
	< f.v $(vpa) text=ff  |\
	  $(vp) xsize=3. ysize=3. > f2.v yshift=.22 xshift=3.1 
	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 aa bb cc dd ee ff
$R/first200.2.v:img.active.tpow.462.T img.First.traces.freq.dec.0.tp.T img.manyevents.T
	< img.active.tpow.462.T Window3d  \
	   f1=18 f2=34 n1=360 n2=404 n3=176 > a0.T
	< img.First.traces.freq.dec.0.tp.T Window3d  \
	   f1=18 f2=34 n1=360 n2=404 n3=176   > a1.T
	< img.manyplanes.T Window3d  \
	   f1=18 f2=34 n1=360 n2=404 n3=176   > a2.T
	< a0.T $(CB) frame1=316 frame2=279 frame3=65 out=a.v 
	< a1.T $(CB) frame1=316 frame2=279 frame3=65 out=c.v 
	< a2.T $(CB) frame1=316 frame2=279 frame3=65 out=e.v 
	Rm a0.T a1.T a2.T
	< img.active.tpow.462.T Window3d  \
	   max1=2000 f1=46 f2=8 n2=441 f3=27 n3=134 > a0.T
	< img.First.traces.freq.dec.0.tp.T Window3d  \
	   max1=2000 f1=46 f2=8 n2=441 f3=27 n3=134   > a1.T
	< img.manyplanes.T Window3d  \
	   max1=2000 f1=46 f2=8 n2=441 f3=27 n3=134   > a2.T
	< a0.T $(CB) frame1=240 frame2=108 frame3=81 out=b.v o1num=1000
	< a1.T $(CB) frame1=240 frame2=108 frame3=81 out=d.v o1num=1000
	< a2.T $(CB) frame1=240 frame2=108 frame3=81 out=f.v o1num=1000
	Rm a0.T a1.T a2.T
	cat ${TXTDIR}/cube_a ${TXTDIR}/first200oval1 > aa
	cat ${TXTDIR}/cube_b ${TXTDIR}/first200oval2 ${TXTDIR}/first200oval3 > bb
	cat ${TXTDIR}/cube_c ${TXTDIR}/first200oval1 > cc
	cat ${TXTDIR}/cube_d ${TXTDIR}/first200oval2 ${TXTDIR}/first200oval3 > dd
	cat ${TXTDIR}/cube_e ${TXTDIR}/first200oval1 > ee
	cat ${TXTDIR}/cube_f ${TXTDIR}/first200oval2 ${TXTDIR}/first200oval3 > ff
	< a.v $(vpa) text=aa  |\
	  $(vp) xsize=3. ysize=3. > a2.v yshift=6. xshift=.1
	< b.v $(vpa) text=bb  |\
	  $(vp) xsize=3. ysize=3. > b2.v yshift=6. xshift=3.1
	< c.v $(vpa) text=cc  |\
	  $(vp) xsize=3. ysize=3. > c2.v yshift=3.1 xshift=.1
	< d.v $(vpa) text=dd  |\
	  $(vp) xsize=3. ysize=3. > d2.v yshift=3.1 xshift=3.1 
	< e.v $(vpa) text=ee  |\
	  $(vp) xsize=3. ysize=3. > e2.v yshift=.22 xshift=.1 
	< f.v $(vpa) text=ff  |\
	  $(vp) xsize=3. ysize=3. > f2.v yshift=.22 xshift=3.1 
	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 aa bb cc dd ee ff


$R/fake200.v:img.manyevents.T img.First.traces.freq.dec.0.tp.T
	< img.manyevents.T Window3d  \
	   f1=18 f2=34 n1=360 n2=404 n3=176 > a0.T
	< img.First.traces.freq.dec.0.tp.T Window3d  \
	   f1=18 f2=34 n1=360 n2=404 n3=176   > a1.T
	< a0.T $(CB) frame1=316 frame2=279 frame3=65 out=a.v 
	< a1.T $(CB) frame1=316 frame2=279 frame3=65 out=b.v 
	Rm a0.T a1.T
	< img.manyevents.T Window3d  \
	   max1=2000 f1=46 f2=8 n2=441 f3=27 n3=134 > a0.T
	< img.First.traces.freq.dec.0.tp.T Window3d  \
	   max1=2000 f1=46 f2=8 n2=441 f3=27 n3=134   > a1.T
	< a0.T $(CB) frame1=240 frame2=108 frame3=81 out=c.v o1num=1000
	< a1.T $(CB) frame1=240 frame2=108 frame3=81 out=d.v o1num=1000
	cat ${TXTDIR}/cube_a ${TXTDIR}/first200oval1 > aa
	< a.v $(vpa) text=aa  |\
	  $(vp) xsize=3. ysize=3. > a2.v yshift=3.1 xshift=.1
	cat ${TXTDIR}/cube_b ${TXTDIR}/first200oval1 > bb
	< b.v $(vpa) text=bb  |\
	  $(vp) xsize=3. ysize=3. > b2.v yshift=3.1 xshift=3.1
	cat ${TXTDIR}/cube_c ${TXTDIR}/first200oval2 ${TXTDIR}/first200oval3 > cc
	< c.v $(vpa) text=cc  |\
	  $(vp) xsize=3. ysize=3. yshift=.22 > c2.v xshift=.1
	cat ${TXTDIR}/cube_d ${TXTDIR}/first200oval2 ${TXTDIR}/first200oval3 > dd
	< d.v $(vpa) text=dd  |\
	  $(vp) xsize=3. ysize=3. yshift=.22 > d2.v xshift=3.1 
	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 a0.T a1.T
$R/2fake200.v:img.manyplanes.T img.First.traces.freq.dec.0.tp.T
	< img.manyplanes.T Window3d  \
	   f1=18 f2=34 n1=360 n2=404 n3=176 > a0.T
	< img.First.traces.freq.dec.0.tp.T Window3d  \
	   f1=18 f2=34 n1=360 n2=404 n3=176   > a1.T
	< a0.T $(CB) frame1=316 frame2=279 frame3=65 out=a.v 
	< a1.T $(CB) frame1=316 frame2=279 frame3=65 out=b.v 
	Rm a0.T a1.T
	< img.manyplanes.T Window3d  \
	   max1=2000 f1=46 f2=8 n2=441 f3=27 n3=134 > a0.T
	< img.First.traces.freq.dec.0.tp.T Window3d  \
	   max1=2000 f1=46 f2=8 n2=441 f3=27 n3=134   > a1.T
	< a0.T $(CB) frame1=240 frame2=108 frame3=81 out=c.v o1num=1000
	< a1.T $(CB) frame1=240 frame2=108 frame3=81 out=d.v o1num=1000
	cat ${TXTDIR}/cube_a ${TXTDIR}/first200oval1 > aa
	< a.v $(vpa) text=aa  |\
	  $(vp) xsize=3. ysize=3. > a2.v yshift=3.1 xshift=.1
	cat ${TXTDIR}/cube_b ${TXTDIR}/first200oval1 > bb
	< b.v $(vpa) text=bb  |\
	  $(vp) xsize=3. ysize=3. > b2.v yshift=3.1 xshift=3.1
	cat ${TXTDIR}/cube_c ${TXTDIR}/first200oval2 ${TXTDIR}/first200oval3 > cc
	< c.v $(vpa) text=cc  |\
	  $(vp) xsize=3. ysize=3. yshift=.22 > c2.v xshift=.1
	cat ${TXTDIR}/cube_d ${TXTDIR}/first200oval2 ${TXTDIR}/first200oval3 > dd
	< d.v $(vpa) text=dd  |\
	  $(vp) xsize=3. ysize=3. yshift=.22 > d2.v xshift=3.1 
	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 a0.T a1.T


$R/shortlong.v:img.First.traces.freq.dec.0.tp.T img.sum.First.T2 #img.sum.First.bp.01.T
	< img.sum.First.T2 Window3d  \
	   f1=18 f2=34 n1=360 n2=404 n3=176 > a1.T
	< img.First.traces.freq.dec.0.tp.T Window3d  \
	   f1=18 f2=34 n1=360 n2=404 n3=176   > a0.T
	< a0.T $(CB) frame1=316 frame2=279 frame3=65 out=a.v 
	< a1.T $(CB) frame1=316 frame2=279 frame3=65 out=b.v 
	Rm a0.T a1.T
	<  img.sum.First.T2 Window3d  \
	   max1=2000 f1=46 f2=8 n2=441 f3=27 n3=134 > a1.T
	< img.First.traces.freq.dec.0.tp.T Window3d  \
	   max1=2000 f1=46 f2=8 n2=441 f3=27 n3=134   > a0.T
	< a0.T $(CB) frame1=240 frame2=108 frame3=81 out=c.v o1num=1000
	< a1.T $(CB) frame1=240 frame2=108 frame3=81 out=d.v o1num=1000
	cat ${TXTDIR}/cube_a ${TXTDIR}/first200oval2 ${TXTDIR}/first200oval3 > aa
	< c.v $(vpa) text=aa  |\
	  $(vp) xsize=3. ysize=3. > a2.v yshift=3.1 xshift=.1
	cat ${TXTDIR}/cube_b ${TXTDIR}/first200oval2 ${TXTDIR}/first200oval3 > bb
	< d.v $(vpa) text=bb  |\
	  $(vp) xsize=3. ysize=3. > b2.v yshift=3.1 xshift=3.1
	cat ${TXTDIR}/cube_c ${TXTDIR}/first200oval1 > cc
	< a.v $(vpa) text=cc  |\
	  $(vp) xsize=3. ysize=3. yshift=.22 > c2.v xshift=.1
	cat ${TXTDIR}/cube_d ${TXTDIR}/first200oval1 > dd
	< b.v $(vpa) text=dd  |\
	  $(vp) xsize=3. ysize=3. yshift=.22 > d2.v xshift=3.1 
	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 a0.T a1.T
$R/longact.v:img.active.tp.T img.sum.First.T2 #img.sum.First.bp.01.T
	< img.sum.First.T2 Window3d  \
	   f1=18 f2=34 n1=360 n2=404 n3=176 > a1.T
	< img.active.tp.T  Window3d  \
	   f1=18 f2=34 n1=360 n2=404 n3=176   > a0.T
	< a0.T $(CB) frame1=316 frame2=279 frame3=65 out=e.v o1num=1000
	< a1.T $(CB) frame1=316 frame2=279 frame3=65 out=f.v o1num=1000
	< a0.T $(CB) frame1=316 frame2=108 frame3=65 out=c.v o1num=1000
	< a1.T $(CB) frame1=316 frame2=108 frame3=65 out=d.v o1num=1000
	Rm a0.T a1.T
	<  img.sum.First.T2 Window3d  \
	   max1=2000 f1=46 f2=8 n2=441 f3=27 n3=134 > a1.T
	< img.active.tp.T  Window3d  \
	   max1=2000 f1=46 f2=8 n2=441 f3=27 n3=134   > a0.T
	< a0.T $(CB) frame1=240 frame2=108 frame3=81 out=a.v o1num=1000
	< a1.T $(CB) frame1=240 frame2=108 frame3=81 out=b.v o1num=1000
	cat ${TXTDIR}/cube_a ${TXTDIR}/first200oval2 ${TXTDIR}/first200oval3 > aa
	< a.v $(vpa) text=aa  |\
	  $(vp) xsize=3. ysize=3. > a2.v yshift=6.1 xshift=.1
	cat ${TXTDIR}/cube_b ${TXTDIR}/first200oval2 ${TXTDIR}/first200oval3 > bb
	< b.v $(vpa) text=bb  |\
	  $(vp) xsize=3. ysize=3. > b2.v yshift=6.1 xshift=3.1
	cat ${TXTDIR}/cube_c ${TXTDIR}/first200oval2 > cc
	< c.v $(vpa) text=cc  |\
	  $(vp) xsize=3. ysize=3. yshift=3.1 > c2.v xshift=.1
	cat ${TXTDIR}/cube_d ${TXTDIR}/first200oval2 > dd
	< d.v $(vpa) text=dd  |\
	  $(vp) xsize=3. ysize=3. yshift=3.1 > d2.v xshift=3.1 
	cat ${TXTDIR}/cube_e ${TXTDIR}/first200oval1 > ee
	< e.v $(vpa) text=ee  |\
	  $(vp) xsize=3. ysize=3. yshift=.22 > e2.v xshift=.1
	cat ${TXTDIR}/cube_f ${TXTDIR}/first200oval1 > ff
	< f.v $(vpa) text=ff  |\
	  $(vp) xsize=3. ysize=3. yshift=.22 > f2.v xshift=3.1 
	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 aa bb cc dd ee ff 
	Rm a0.T a1.T
$R/shal200.v:img.noise2.T img.First.traces.freq.dec.0.tp.T #img.active.shal.T
	echo o1=60 >> img.noise2.T
	echo o1=60 >> img.First.traces.freq.dec.0.tp.T
	< img.noise2.T Window3d  max1=1000. > trz.T
	< trz.T $(CB) frame1=20 frame2=100 frame3=58 out=b.v o1num=260 d1num=200 
	< img.First.traces.freq.dec.0.tp.T Window3d  max1=1000. > trzb.T
	< trzb.T $(CB) frame1=20 frame2=100 frame3=58 out=a.v o1num=260 d1num=200 
#	< img.active.shal.T $(CB) frame1=20 frame2=100 frame3=58 out=d.v o1num=260 d1num=200 
	$(vpa) < b.v text=${TXTDIR}/cube_b | $(vp) xsize=3. ysize=2. yshift=3.1 > b2.v xshift=3.05 
	$(vpa) < a.v  text=${TXTDIR}/cube_a | $(vp) xsize=3. ysize=2. yshift=3.1 > a2.v xshift=.1
#	$(vpa) < a.v  text=${TXTDIR}/cube_c | $(vp) xsize=3. ysize=2. yshift=.2 > c2.v xshift=.1
#	$(vpa) < d.v  text=${TXTDIR}/cube_d | $(vp) xsize=3. ysize=2. yshift=.2 > d2.v xshift=3.05
#	vp_Overlay a2.v b2.v c2.v d2.v > $@
	vp_Overlay a2.v b2.v > $@
	rm -f a.v b.v c.v d.v a2.v b2.v c2.v d2.v ;Rm trz.T trzb.T
	echo o1=62.5 >> img.noise2.T
	echo o1=62.5 >> img.First.traces.freq.dec.0.tp.T
$R/zoom1.v:img.active.tpow.462.T img.First.traces.freq.dec.0.tp.T ${TXTDIR}/zoom1oval
	< img.active.tpow.462.T Window3d  \
	    f1=222 f2=14 f3=0 f4=0 f5=0 n1=380 n2=436 n3=176 n4=1 n5=1  >a0.T
	< img.First.traces.freq.dec.0.tp.T Window3d  \
	    f1=222 f2=14 f3=0 f4=0 f5=0 n1=380 n2=436 n3=176 n4=1 n5=1  >a1.T
	< a0.T $(CB) frame1=118 frame2=266 frame3=69 out=a.v 
	< a1.T $(CB) frame1=118 frame2=266 frame3=69 out=b.v 
	cat  ${TXTDIR}/cube_a ${TXTDIR}/zoom1oval > aa 
	cat  ${TXTDIR}/cube_b ${TXTDIR}/zoom1oval > bb
	< a.v $(vpa) text=aa  |\
	  $(vp) xsize=3. ysize=3. yshift=.22 > a2.v xshift=.1
	< b.v $(vpa) text=bb  |\
	  $(vp) xsize=3. ysize=3. yshift=.22 > b2.v xshift=3.1 
	vp_Overlay a2.v b2.v > $@
	rm -f a.v b.v a2.v b2.v aa bb ;Rm a0.T a1.T
gath=
#gath=/scr1/brad/img.gathers.Fourth.compare.T
gath=img.gathers.First.T
sum=img.sum.First.T2 
#sum=/scr1/brad/img.sum.Fourth.T
$R/migdmig.v: $(sum) $(gath)
	< $(sum) Window3d  \
	   f1=18 f2=34 n1=360 n2=404 n3=176 > a1.T
	< $(gath) Window3d  \
	   f1=18 f2=34 n1=360 n2=404 n3=176   > a0.T
	< a0.T $(CB) frame1=316 frame2=279 frame3=65 out=a.v 
	< a1.T $(CB) frame1=316 frame2=279 frame3=65 out=b.v 
	Rm a0.T a1.T
	<  $(sum) Window3d  \
	   max1=2000 f1=46 f2=8 n2=441 f3=27 n3=134 > a1.T
	< $(gath) Window3d  \
	   max1=2000 f1=46 f2=8 n2=441 f3=27 n3=134   > a0.T
	< a0.T $(CB) frame1=240 frame2=108 frame3=81 out=c.v o1num=1000
	< a1.T $(CB) frame1=240 frame2=108 frame3=81 out=d.v o1num=1000
	cat ${TXTDIR}/cube_a ${TXTDIR}/first200oval1 > aa
	< a.v $(vpa) text=aa  |\
	  $(vp) xsize=3. ysize=3. > a2.v yshift=3.1 xshift=.1
	cat ${TXTDIR}/cube_b ${TXTDIR}/first200oval1 > bb
	< b.v $(vpa) text=bb  |\
	  $(vp) xsize=3. ysize=3. > b2.v yshift=3.1 xshift=3.1
	cat ${TXTDIR}/cube_c ${TXTDIR}/first200oval2 ${TXTDIR}/first200oval3 > cc
	< c.v $(vpa) text=cc  |\
	  $(vp) xsize=3. ysize=3. yshift=.22 > c2.v xshift=.1
	cat ${TXTDIR}/cube_d ${TXTDIR}/first200oval2 ${TXTDIR}/first200oval3 > dd
	< d.v $(vpa) text=dd  |\
	  $(vp) xsize=3. ysize=3. yshift=.22 > d2.v xshift=3.1 
	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 a0.T a1.T
$R/cmigact.v: img.gathers.First.T img.active.tp.T 
	< img.active.tp.T Window3d  \
	   f1=18 f2=34 n1=360 n2=404 n3=176 > a1.T
	< img.gathers.First.T Window3d  \
	   f1=18 f2=34 n1=360 n2=404 n3=176   > a0.T
	< a0.T $(CB) frame1=116 frame2=279 frame3=65 out=a.v 
	< a1.T $(CB) frame1=116 frame2=279 frame3=65 out=b.v 
	Rm a0.T a1.T
	< img.active.tp.T Window3d  \
	   max1=2000 f1=46 f2=8 n2=441 f3=27 n3=134 > a1.T
	< img.gathers.First.T Window3d  \
	   max1=2000 f1=46 f2=8 n2=441 f3=27 n3=134   > a0.T
	< a0.T $(CB) frame1=200 frame2=108 frame3=81 out=c.v o1num=1000
	< a1.T $(CB) frame1=200 frame2=108 frame3=81 out=d.v o1num=1000
	cat ${TXTDIR}/cube_a ${TXTDIR}/first200oval1 > aa
	< a.v $(vpa) text=aa  |\
	  $(vp) xsize=3. ysize=3. > a2.v yshift=3.1 xshift=.1
	cat ${TXTDIR}/cube_b ${TXTDIR}/first200oval1 > bb
	< b.v $(vpa) text=bb  |\
	  $(vp) xsize=3. ysize=3. > b2.v yshift=3.1 xshift=3.1
	cat ${TXTDIR}/cube_c ${TXTDIR}/first200oval2 ${TXTDIR}/first200oval3 > cc
	< c.v $(vpa) text=cc  |\
	  $(vp) xsize=3. ysize=3. yshift=.22 > c2.v xshift=.1
	cat ${TXTDIR}/cube_d ${TXTDIR}/first200oval2 ${TXTDIR}/first200oval3 > dd
	< d.v $(vpa) text=dd  |\
	  $(vp) xsize=3. ysize=3. yshift=.22 > d2.v xshift=3.1 
	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 a0.T a1.T
include ${SEPINC}/SEP.bottom