include ${SEPINC}/SEP.top
-include ../paper/graphics.mf 

RESULTOLD =     nseaFoltTotXPefLN 
RESULTSER = 	zeroFoltXPef zeroFolt1Pef3d zeroFoltXPefLN zeroFoltXPefAA 

RESULTSCR = 	gulfFoltTot1Pef3d nseaFoltTot1Pef3d \
		gulfFoltTotXPefAA nseaFoltTotXPefAA \
		gulfFoltTotXPefLN 

pef = pef.pef

zero3darg = olap=.5f,.5f,0.5f psze=12,12,12 fsze=4,2,2  
${RESDIR}/zeroFolt1Pef3d.v : ${DAT}/zeroFoltRaw.HH 
	< $<     ${JAVAX} ${pef}.IcaiPefNStat ${zero3darg}           > junk.H 
	< junk.H ${TMA}                                              > junk.A
	< junk.A ${CPLT} ${PFRAME} title="PEF: one 3D filter" out=$@ > junk.V

zerolnargs = olap=.5f,.5f,0.5f psze=12,12,12 pefi=4,4,4
${RESDIR}/zeroFoltXPefLN.v : ${DAT}/zeroFoltRaw.HH 
	<$<     ${JX}.XPefLocalNormNStat ${zerolnargs}                  >junk.H
	<junk.H ${TPA}                                                  >junk.A
	<junk.A ${CPLT} ${PFRAME} title="XPEF: Three 2D filters" out=$@ >junk.V

${RESDIR}/zeroFoltXPefAA.v : ${DAT}/zeroFoltRaw.HH 
	<$<     ${JX}.XPefAndAdjNStat ${zerolnargs}                     >junk.H
	<junk.H ${TPA}                                                  >junk.A
	<junk.A ${CPLT} ${PFRAME} title="XPEF: Three 2D filters" out=$@ >junk.V

lplt = 'xll=1    yll=1 xur=4.58  yur=9.5 title="XPef: Z"'
mplt = 'xll=4.98 yll=1 xur=8.56  yur=9.5 title="XPef: Y"'
rplt = 'xll=8.96 yll=1 xur=12.65 yur=9.5 title="XPef: X"'
xpefargs = olap=.5f,.5f,0.5f psze=12,12,12
${RESDIR}/zeroFoltXPef.v :   ${DAT}/zeroFoltRaw.HH 
	< $<      ${JAVAX} ${pef}.IcaiPefNStat ${xpefargs} fsze=4,2,1 > junk0.H 
	< junk0.H ${TMA}                                              > junk0.A;
	< junk0.A ${CPLT} ${PFRAME} ${rplt}         out=junk0.v       > junk0.V
	< $<      ${JAVAX} ${pef}.IcaiPefNStat ${xpefargs} fsze=4,1,2 > junk1.H 
	< junk1.H ${TMA}                                              > junk1.A;
	< junk1.A ${CPLT} ${PFRAME} ${mplt}         out=junk1.v       > junk1.V
	< $<      ${JAVAX} ${pef}.IcaiPefNStat ${xpefargs} fsze=1,4,2 > junk2.H 
	< junk2.H ${TMA}                                              > junk2.A; 
	< junk2.A ${CPLT} ${PFRAME} ${lplt}         out=junk2.v       > junk2.V
	vppen erase=n vpstyle=n  junk2.v junk1.v junk0.v                     > $@

NSEAARG = olap=.5f,.5f,0.5f psze=20,12,12 fsze=4,2,2
.PRECIOUS: gulfFoltTot1Pef3d.H nseaFoltTot1Pef3d.H 
%Tot1Pef3d.H : ${DAT}/%TotRaw.HH 
	< $<    ${JAVAX} ${pef}.IcaiPefNStat  ${NSEAARG}   > pef3.H
	${JAVAX} rsf.util.Tent nsamples=5,1,1              > filt.H
	< pef3.H ${JAVAX} rsf.util.Convolver filter=filt.H > $@
pef3d = title="PEF: one 3D filter"
${RESDIR}/%Tot1Pef3d.v : %Tot1Pef3d.H
	< $< ${ZPLT} allpos=n > junk.V  out=zplt.v 
	< $< ${XPLT} allpos=n > junk.V  out=xplt.v 
	< $< ${YPLT} allpos=n > junk.V  out=yplt.v 
	vppen erase=n vpstyle=n zplt.v xplt.v yplt.v > $@
	< $< ${TRAN}                                 > tran.H
	< tran.H ${GMA} ${TOPGREY} ${pef3d}  out=$@3 > junk.V

${RESDIR}/%TotXPefAAMrk.v : ${RESDIR}/%TotXPefAA.v
	< $< vp_annotate text=../paper/$*TotMrk.txt batch=y > $@

.PRECIOUS: gulfFoltTotXPefLN.H nseaFoltTotXPefLN.H 
%TotXPefLN.H : ${DAT}/%TotRaw.HH 
	< $<     ${JX}.XPefLocalNormNStat ${NSEAARG}       > xpln.H
	${JAVAX} rsf.util.Tent nsamples=5,1,1              > filt.H
	< xpln.H ${JAVAX} rsf.util.Convolver filter=filt.H > $@

xlnarg = 'title="XPEF: Three 2D filters"'
${RESDIR}/%TotXPefLN.v : %TotXPefLN.H
	< $< ${ZPLT} allpos=y > junk.V  out=zplt.v 
	< $< ${XPLT} allpos=y > junk.V  out=xplt.v 
	< $< ${YPLT} allpos=y > junk.V  out=yplt.v 
	vppen erase=n vpstyle=n zplt.v xplt.v yplt.v > $@
	#< $< ${TRAN}                                 > tran.H
	#< tran.H ${GPA} ${TOPGREY} ${xlnarg} out=$@3 > junk.V

.PRECIOUS: gulfFoltTotXPefAA.H nseaFoltTotXPefAA.H 
%TotXPefAA.H : ${DAT}/%TotRaw.HH 
	< $<     ${JX}.XPefAndAdjNStat ${NSEAARG}          > xpaa.H
	${JAVAX} rsf.util.Tent nsamples=5,1,1              > filt.H
	< xpaa.H ${JAVAX} rsf.util.Convolver filter=filt.H > $@
xaaarg = 'title="Correlated wave estimate"'
${RESDIR}/%TotXPefAA.v : %TotXPefAA.H
	< $< ${ZPLT} allpos=n > junk.V  out=zplt.v 
	< $< ${XPLT} allpos=n > junk.V  out=xplt.v 
	< $< ${YPLT} allpos=n > junk.V  out=yplt.v 
	vppen erase=n vpstyle=n zplt.v xplt.v yplt.v > $@
	#< $< ${TRAN}                                 > tran.H
	#< tran.H ${GPA} ${TOPGREY} ${xaaarg} out=$@3 > junk.V

clean : jclean

include ${SEPINC}/SEP.bottom


