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

RESDIR    = ../Fig
RESULTSER = transOrig intrnOrig tcai icai icaf

default: view 

INTNL=250
TRANS=265

iput = ../Data/wz.24.HH
iput.H: ${iput}
	< ${iput} Window f2=36 n2=1 f1=180 j1=2 n1=250 > junk.H
	echo "title=Seismogram"            >> junk.H
	Cp < junk.H out=stdout              > $@

DOTARGS = title=\"\" labels=trace:wavelet:reflectivity labelsz=14 
${RESDIR}/transOrig.v : SynTrace.java waveimp.p 
	Wavelet par=waveimp.p out=stdout                            > wavlt.H 
	${JAVAX} ${PACKAGE}.SynTrace conv=transient wavelet=wavlt.H > trace.H 
	< wavlt.H Pad n1out=${TRANS}                                > pw.H
	< rflct.H Pad n1out=${TRANS}                                > pr.H 
	< trace.H Pad n1out=${TRANS}                                > pt.H
	Merge space=n axis=2 < pt.H pw.H pr.H                       > junk.H
	< junk.H Dots ${DOTARGS} out=$@                             > /dev/null

${RESDIR}/intrnOrig.v : SynTrace.java waveimp.p
	Wavelet par=waveimp.p out=stdout                            > wavlt.H 
	${JAVAX} ${PACKAGE}.SynTrace                wavelet=wavlt.H > trace.H 
	< wavlt.H Pad n1out=${INTNL}                                > pw.H
	< rflct.H Pad n1out=${INTNL}                                > pr.H 
	< trace.H Pad n1out=${INTNL}                                > pt.H
	Merge space=n axis=2 < pt.H pw.H pr.H                       > junk.H
	< junk.H Dots ${DOTARGS} out=$@                             > /dev/null

TCAIARGS = $DOTARGS title=\"\" labels=original:estimation:trace:wavelet
${RESDIR}/tcai.v : TcaiDecon.java SynTrace.java waveimp.p iput.H 
	Wavelet par=waveimp.p out=stdout                            > wavlt.H 
	${JAVAX} ${PACKAGE}.SynTrace \
			in=iput.H conv=transient wavelet=wavlt.H   > trace.H 
	${JAVAX} ${PACKAGE}.TcaiDecon                               > estim.H 
	< wavlt.H Pad n1out=${TRANS}                                > pw.H
	< trace.H Pad n1out=${TRANS}                                > pt.H
	< estim.H Pad n1out=${TRANS}                                > pp.H 
	< rflct.H Pad n1out=${TRANS}                                > pr.H 
	Merge space=n axis=2 < pr.H pp.H pt.H pw.H                  > junk.H
	< junk.H Dots ${TCAIARGS} out=$@                            > /dev/null

ICAIARGS = $DOTARGS title=\"\" labels=original:estimation:trace:wavelet
${RESDIR}/icai.v : IcaiDecon.java SynTrace.java waveimp.p
	Wavelet par=waveimp.p out=stdout                            > wavlt.H 
	${JAVAX} ${PACKAGE}.SynTrace                wavelet=wavlt.H > trace.H 
	${JAVAX} ${PACKAGE}.IcaiDecon                               > estim.H 
	< wavlt.H Pad n1out=${INTNL}                                > pw.H
	< trace.H Pad n1out=${INTNL}                                > pt.H
	< estim.H Pad n1out=${INTNL}                                > pp.H 
	< rflct.H Pad n1out=${INTNL}                                > pr.H 
	Merge space=n axis=2 < pr.H pp.H pt.H pw.H                  > junk.H
	< junk.H Dots ${ICAIARGS} out=$@                            > /dev/null

ICAFARGS = $DOTARGS title=\"\" labels=original:estimation:trace:reflectivity
${RESDIR}/icaf.v : IcafDecon.java SynTrace.java waveimp.p
	Wavelet par=waveimp.p out=stdout                            > wavlt.H 
	${JAVAX} ${PACKAGE}.SynTrace                wavelet=wavlt.H > trace.H 
	${JAVAX} ${PACKAGE}.IcafDecon                               > estim.H 
	< rflct.H Pad n1out=${INTNL}                                > pr.H 
	< trace.H Pad n1out=${INTNL}                                > pt.H
	< estim.H Pad n1out=${INTNL}                                > pp.H 
	< wavlt.H Pad n1out=${INTNL}                                > pw.H
	Merge space=n axis=2 < pw.H pp.H pt.H pr.H                  > junk.H
	< junk.H Dots ${ICAFARGS} out=$@                            > /dev/null




${RESDIR}/wavelet.v : trace.H reflect1.H wavelet1.H 
	${JAVAX} ${PACKAGE}.IcafDecon < trace.H reflectivity=reflect1.H > wavelet2.H 
	< wavelet1.H Dots out=junky1.v > /dev/null 
	echo title="Known Ricker Wvlt"      >>wavelet1.H 
	< wavelet2.H Dots out=junky2.v > /dev/null 
	echo title="Estimated Wvlt"      >>wavelet2.H 
	vp_OverUnderAniso junky1.v junky2.v > ${RESDIR}/wavelet.v 

${RESDIR}/reflect.v : trace.H reflect1.H wavelet1.H 
	${JAVAX} ${PACKAGE}.IcaiDecon < trace.H wavelet=wavelet1.H > reflect2.H
	< reflect1.H Dots ${DOTARGS} out=junky1.v > /dev/null 
	< reflect2.H Dots ${DOTARGS} out=junky2.v > /dev/null 
	vp_OverUnderAniso junky1.v junky2.v > ${RESDIR}/reflect.v 

${RESDIR}/decon.v : ${RESDIR}/trace.v ${RESDIR}/wavelet.v 
	vp_SideBySideIso ${RESDIR}/trace.v ${RESDIR}/wavelet.v > ${RESDIR}/decon.v

clean : jclean

include ${SEPINC}/SEP.bottom
