UF77LIBS = ${TDFLIB}
RESULTSER = raw nmo
RESULTSER = dune3D gravel2D

#RESDIR = ../Fig

include ${SEPINC}/SEP.top
VPLOT = Ta2vplot crowd=.99 > /dev/null
WINDOW = Window n3=1 n1=600
DATA = ../../Data

default : tube

${RESDIR}/raw.v ${RESDIR}/nmo.v : ${DATA}/dune3D.H  Nmo3.x makefile 
	< ${DATA}/dune3D.H Window xn1=300 | Tpow tpow=-2 > junk.H
####### 						Next line is guesswork.
	echo o2=-4.8 d2=.1  o3=4.8 d3=-.8 >> junk.H
	echo o2=-4.8 d2=.1  o3=2.4 d3=-.4 >> junk.H
	echo o2=-4.8 d2=.1  o3=3.6 d3=-.6 >> junk.H
	echo o2=-4.8 d2=.1  o3=3.6 d3=-.6 Xo1=-.15 >> junk.H
	echo o2=-4.8 d2=.1  o3=2.4 d3=-.4 o1=-.15 >> junk.H
	< junk.H Window Xf3=5 Xn3=1 > junk2.H
	< junk2.H Byte gainpanel=every pclip=85 > junk.A
	< junk.A Ta2vplot transp=n crowd1=.87 >/dev/null out=${RESDIR}/raw.v
	< junk.H Nmo3.x    vel=7    > nmo.H
	< nmo.H Byte gainpanel=every pclip=85 > junk.A
	< junk.A Ta2vplot transp=n crowd1=.87 >/dev/null out=${RESDIR}/nmo.v


${RESDIR}/gravel2D.v : ${DATA}/gravel2D.H     
	< ${DATA}/gravel2D.H Byte tpow=-1 pclip=92|Window n1=500 n2=160> gravel2D.A
	< gravel2D.A Ta2vplot transp=n crowd1=.87 >/dev/null out=${RESDIR}/gravel2D.v

${RESDIR}/dune3D.v ${RESDIR}/dune3D.v3 : ${DATA}/dune3D.H  
	< ${DATA}/dune3D.H Byte gainpanel=every tpow=-2 pclip=85 > dune3D.A
	< dune3D.A Ta2vplot transp=n crowd=.99 > /dev/null   out=${RESDIR}/dune3D.v3
	< dune3D.A ${WINDOW} f3=0  | ${VPLOT} out=frame0.v
	< dune3D.A ${WINDOW} f3=1  | ${VPLOT} out=frame1.v
	< dune3D.A ${WINDOW} f3=2  | ${VPLOT} out=frame2.v
	< dune3D.A ${WINDOW} f3=3  | ${VPLOT} out=frame3.v
	< dune3D.A ${WINDOW} f3=4  | ${VPLOT} out=frame4.v
	< dune3D.A ${WINDOW} f3=5  | ${VPLOT} out=frame5.v
	vp_SideBySideAniso frame5.v frame4.v frame3.v > top.v
	vp_SideBySideAniso frame2.v frame1.v frame0.v > bot.v
	vp_OverUnderAniso  top.v bot.v                        > ${RESDIR}/dune3D.v

clean:  jclean

include ${SEPINC}/SEP.bottom
