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

RESULTOLD =     gulfFoltTotHLapTop nseaFoltTotHLapTop traceCor zeroFoltHLap 
RESULTSER = 	zeroFoltCor 
RESULTSCR = 	gulfFoltTotCorr nseaFoltTotCorr 

# horizontal correlation:

ZEROARG = olap=.5f,.5f,0.5f psze=5,5,5 
${RESDIR}/zeroFoltCor.v :  ${DAT}/zeroFoltRaw.HH 
	<$<     ${JX}.TraceCorrelatorToRsfNStat ${ZEROARG}              >junk.H
	<junk.H Add scale=-1                                            >scal.H
	<scal.H Add add=1                                               >addd.H
	<addd.H ${TPA}                                                  >junk.A
	<junk.A ${CPLT} ${PFRAME} title="Horizontal Correlation" out=$@ >junk.V

${RESDIR}/zeroFoltHLap.v : ${DAT}/zeroFoltRaw.HH 
	${JAVAX} rsf.util.Laplacian dim=2                  > lapl.H
	< $<     Window f1=30 n1=1 out=stdout              > wind.H
	< wind.H Cp            > junk.H 
	echo "dim=2 \n"       >> junk.H 
	< junk.H Cp out=stdout > wind.H 
	< wind.H ${JAVAX} rsf.util.Convolver filter=lapl.H > junk.H
	< junk.H ${GMA} title="Horizontal Laplace"  out=$@ > junk.V; 

.PRECIOUS: gulfFoltTotHLap.H nseaFoltTotHLap.H
%TotHLap.H : ${DAT}/%TotRaw.HH 
	${JAVAX} rsf.util.Laplacian dim=2                  > filt.H
	< $<     Window f1=${$*Z} n1=1 out=stdout          > wind.H 
	< wind.H Cp            > junk.H 
	echo "dim=2 \n"       >> junk.H 
	< junk.H Cp out=stdout > wind.H 
	< wind.H ${JAVAX} rsf.util.Convolver filter=filt.H > $@
hlap = 'title="Horizontal Laplace"'
${RESDIR}/%TotHLapTop.v : %TotHLap.H
	< $<     ${GPA} ${TOPGREY} allpos=n ${hlap} out=$@  >junk.V

NSEAARG = olap=.5f,.5f,.5f psze=5,3,3 
.PRECIOUS: gulfFoltTotCorr.H nseaFoltTotCorr.H
%TotCorr.H : ${DAT}/%TotRaw.HH 
	< $< ${JX}.TraceCorrelatorToRsfNStat ${NSEAARG} > corr.H
	< corr.H Add add=1				> $@
corr = 'title="Horizontal correlation"'
${RESDIR}/%TotCorr.v : %TotCorr.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} ${ppcarg} out=$@3 > junk.V

### OLD ###

#${RESDIR}/%TotCorrTop.v : %TotCorr.H
#	<$<     ${TRAN}                                    >tran.H
#	<tran.H ${GPA} ${TOPGREY} allpos=y ${corr} out=$@3 >junk.V
#	<$<     ${ZSID}           allpos=y ${corr} out=$@  >junk.V
#${RESDIR}/%TotCorrSid.v : %TotCorr.H
#	< $< ${XSID} allpos=y      out=xplt.v > junk.V
#	< $< ${YSID} allpos=y      out=yplt.v > junk.V 
#	vppen erase=n vpstyle=n xplt.v yplt.v > $@



 
GPLT = ${GPS} polarity=-1 
${RESDIR}/%TotCor.v : ${DAT}/%TotRaw.HH 
	< $<     ${JX}.TraceCorrelatorToRsfNStat ${NSEAARG}        > corr.H 
	< corr.H ${TRAN}                                           > tran.H
	< tran.H ${GPA} polarity=-1 title="Horizontal Correlation"  out=$@3 \
                                                                   > junk.V
	< corr.H ${ZPLT}                                out=zplt.v > junk.V
	< corr.H ${XPLT}                                out=xplt.v > junk.V
	< corr.H ${YPLT}                                out=yplt.v > junk.V
	vppen erase=n vpstyle=n zplt.v xplt.v yplt.v               > $@

# Horiztonal Laplacian:
# Maybe I should write a small program to compute each slice 
# and output cubes as in all the other cases! 

${RESDIR}/%TotHLap.v : ${DAT}/%TotRaw.HH 
	${JAVAX} rsf.util.Laplacian dim=2                    > lap.H
	${JAVAX} rsf.util.LapGeneral bigFilt=3,3 smlFilt=1,1 > lap.H
	< $<     Window f1=30 n1=1 out=stdout                > wind.H
	< wind.H ${JAVAX} rsf.util.Convolver filter=lap.H    > junk.H
	< corr.H ${GMS}                           out=zplt.v > junk.V

clean: jclean 

include ${SEPINC}/SEP.bottom







