USE_JAG   = yes
#JDK11    = yes
include ${SEPINC}/SEP.top

RESULTSER = qdome dips seis data modl comp vdip
RESULTSER = sigmoid qdome 
RESULTSER = qdome 

${RESDIR}/qdome.v: Qdome.java
	${JAVAX} ${PACKAGE}.Qdome slip=0.1                         > junk0.H
	< junk0.H Smooth rect1=3 diff1=1 | Smooth rect1=3 > qdome.H 
	< qdome.H Taplot gainpanel=all | Ta2vplot out=$@ > junk.V 

test:
	${MAKE} Qdome.jc  
	${JAVAX} ${PACKAGE}.Qdome slip=0.0                         > junk0.H
	${JAVAX} ${PACKAGE}.Qdome slip=1.0                         > junk1.H
	Add scale=1.,-1. < junk0.H junk1.H > junky.H
	< junky.H  Taplot gainpanel=all pclip=100 | Ta2vplot > junk.V 
	< junk.V   Tube &
	< junky.H  Window n1=1 f1=60                         > junk.H
	< junk.H   Taplot gainpanel=all pclip=100 | Ta2vplot > junk2.V
	< junk2.V  Tube &

${RESDIR}/sigmoid.v: Sigmoid.java
	${JAVAX} ${PACKAGE}.Sigmoid d1=.008 n1=200        > junk.H
	< junk.H Smooth rect1=3 diff1=1 | Smooth rect1=3 > sigm.H 
	< sigm.H Taplot | Ta2vplot out=$@ > junk.V 

### old stuff below ### 

JUNK = /dev/null

default: tube

TARG = f2=49 f3=0
dipm.H: Qdome.x  qdome.p3
	 Qdome.x par=qdome.p3 ${TARG} n1=2500 trace=dipm.H > ${JUNK}

qdome.H dips.H: Qdome.x qdome.p3
	 Qdome.x par=qdome.p3 ${TARG}         trace=dips.H > qdome.H

BYTE = Byte gainpanel=all Xpclip=100 

${RESDIR}/%.A : qdome.H
	<    qdome.H ${BYTE} > ${RE${SDI}R}/qdome.A

${RESDIRR}/qdome.v3: qdome.H 
	<        qdome.H ${BYTE} >junk.H; <junk.H Ta2vplot out=${RE${SDI}R}/qdome.v3 > ${JUNK}

WIGGLE = Wiggle transp=1 yreverse=1 screenratio=1.5
${REDIRR}/dips.v3: dips.H 
	echo title=Dipmeter >> dips.H
	< dips.H ${WIGGLE} out=${RE${SDI}R}/dips.v3 > ${JUNK}

${RESDIRR}/seis.v3 seis.H: qdome.H 
	Window <qdome.H ${TARG} n2=2 > seis.H
	echo title='Seismic' >> seis.H
	< seis.H ${WIGGLE} out=${RE${SDI}R}/seis.v3 > ${JUNK}

VEL = velocity=2. offset=0.
${RESDIR}/modl.v3 ${RE${SDI}R}/data.v3 data.H modl.H: qdome.H Kirchhoff.x 
	< qdome.H Kirchhoff.x adj=0 ${VEL} > data.H
	< data.H  Kirchhoff.x adj=1 ${VEL} > modl.H
	< data.H Taplot >junk.H; <junk.H Ta2vplot out=${RE${SDI}R}/data.v3 > ${JUNK}
	< modl.H Taplot >junk.H; <junk.H  Ta2vplot out=${RE${SDI}R}/modl.v3 > ${JUNK}

${RESDIR}/comp.v3: qdome.H modl.H
	< qdome.H Window n3=1 > junk1.H
	< modl.H   Window n3=1 > junk2.H
	Merge space=no < junk1.H junk2.H > comp.H
	< comp.H Taplot gainpanel=each >junk.H; <junk.H  Ta2vplot out=${RE${SDI}R}/comp.v3 > ${JUNK}

${RESDIR}/vdip.v: Vdip.x data.H dipm.H
	Zero chck.H
	Vdip.x < data.H dipm=dipm.H ${TARG} seis=chck.H nvel=30 minv=1. maxv=4.> vdip.H
	< vdip.H Taplot >junk.H; <junk.H  Ta2vplot out=${RE${SDI}R}/vdip.v > ${JUNK}

VDIP = Vdip.o kirchtarg.o dipfit.o norm2.o
Vdip.x :      ${VDIP}

clean: jclean

include ${SEPINC}/SEP.bottom
