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

RESULTOLD = dip cmp coh  zeroFoltDXGPN \
		nseaFoltTotDXGPNTop nseaFoltTotDXGPNSid nseaFoltTotDXGLN 
RESULTSER =     zeroFoltDXG zeroFoltDXGAA zeroFoltDXGLN             
RESULTSCR = 	gulfFoltTotDXGLN \
		gulfFoltTotDXGAA nseaFoltTotDXGAA 

# DXG 	DipXGrad
# DXGPN DipXGradPatchNorm (not used) 
# DXGLN	DipXGradLocalNorm
# DXGAA DipXGradAndAdj

#ZEROARG = olap=.5f,.5f,0.5f psze=7,3,3  
ZEROARG = olap=.5f,.5f,.5f psze=12,12,12  
${RESDIR}/zeroFoltDXG.v :   ${DAT}/zeroFoltRaw.HH 
	< $<     ${JX}.DipXGradNStat comp=0   ${ZEROARG}        > junk.H 
	< junk.H ${TMA}                                         > junk.A; 
	< junk.A ${CPLT} ${PFRAME} title="Xprod: X" out=junk0.v > junk.V
	< $<     ${JX}.DipXGradNStat comp=1   ${ZEROARG}        > junk.H 
	< junk.H ${TMA}                                         > junk.A; 
	< junk.A ${CPLT} ${PFRAME} title="Xprod: Y" out=junk1.v > junk.V
	< $<     ${JX}.DipXGradNStat comp=2   ${ZEROARG}        > junk.H 
	< junk.H ${TMA}                                         > junk.A; 
	< junk.A ${CPLT} ${PFRAME} title="Xprod: Z" out=junk2.v > junk.V
	vp_SideBySideIso junk0.v junk1.v junk2.v > $@

${RESDIR}/zeroFoltDXGAA.v :  ${DAT}/zeroFoltRaw.HH 
	<$<     ${JX}.DipXGradAndAdjNStat    ${ZEROARG}                >junk.H 
	<junk.H ${TMA}                                                 >junk.A
	<junk.A ${CPLT} ${PFRAME} title="Xprod: Backprojection" out=$@ >junk.V
${RESDIR}/zeroFoltDXGPN.v :  ${DAT}/zeroFoltRaw.HH 
	<$<     ${JX}.DipXGradNormNStat      ${ZEROARG}                 >junk.H
	<junk.H ${TPA}                                                  >junk.A
	<junk.A ${CPLT} ${PFRAME} title="Xprod: Patch magnitude" out=$@ >junk.V
${RESDIR}/zeroFoltDXGLN.v : ${DAT}/zeroFoltRaw.HH 
	<$<     ${JX}.DipXGradLocalNormNStat ${ZEROARG}                 >junk.H
	<junk.H ${TPA}                                                  >junk.A
	<junk.A ${CPLT} ${PFRAME} title="Xprod: Pixel magnitude" out=$@ >junk.V

NSEAARG = olap=.5f,.5f,0.5f psze=12,8,8
.PRECIOUS: gulfFoltTotDXGLN.H nseaFoltTotDXGLN.H
%TotDXGLN.H : ${DAT}/%TotRaw.HH 
	< $<    ${JX}.DipXGradLocalNormNStat  ${NSEAARG}   > $@
dxgln = 'title="Local dip estimate norm"'
${RESDIR}/%TotDXGLN.v : %TotDXGLN.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} ${dxgln}  out=$@3 > junk.V

.PRECIOUS: gulfFoltTotDXGAA.H nseaFoltTotDXGAA.H 
%TotDXGAA.H : ${DAT}/%TotRaw.HH 
	< $<    ${JX}.DipXGradAndAdjNStat  ${NSEAARG}   > $@
dxgaa = 'title="Back projection"'
${RESDIR}/%TotDXGAA.v : %TotDXGAA.H
	< $< ${ZPLT} > junk.V  out=zplt.v 
	< $< ${XPLT} > junk.V  out=xplt.v 
	< $< ${YPLT} > junk.V  out=yplt.v 
	vppen erase=n vpstyle=n zplt.v xplt.v yplt.v > $@
	< $< ${TRAN}                                 > tran.H
	< tran.H ${GMA} ${TOPGREY} ${ppcarg} out=$@3 > junk.V

.PRECIOUS: gulfFoltTotDXGPN.H nseaFoltTotDXGPN.H 
%TotDXGPN.H : ${DAT}/%TotRaw.HH 
	< $<    ${JX}.DipXGradNormNStat  ${NSEAARG}   > $@
${RESDIR}/%TotDXGPNTop.v : %TotDXGPN.H
	<$<     ${TRAN}                                                 >tran.H
	<tran.H ${GMA} ${TOPGREY} title="XProd: Backprojection" out=$@3 >junk.V
	<$<     ${ZSID} allpos=n  title="XProd: Backprojection" out=$@  >junk.V
${RESDIR}/%TotDXGPNSid.v : %TotDXGPN.H
	< $< ${XSID} allpos=n > junk.V  out=xplt.v 
	< $< ${YSID} allpos=n > junk.V  out=yplt.v 
	vppen erase=n vpstyle=n xplt.v yplt.v > $@


### 2-D Shotgather Stuff ### 

${RESDIR}/dip.v : DipXGradNStat.java 
	${JX}.CmpSyn                               > iput.H
	< iput.H ${JX}.DipXGradNStat               > oput.H
	< iput.H ${JX}.DipFillerNStat              > dipy.H
	< iput.H ${TPLT}  out=junk2.v              > junk.V
	< oput.H ${TPLT}  out=junk1.v              > junk.V
	< dipy.H ${TPLT}  out=junk0.v              > junk.V
	vp_SideBySideAniso junk2.v junk1.v junk0.v  > $@ 

${RESDIR}/coh.v : ${DAT}/wz.24.HH DipXGradNStat.java 
	< $<     Window f2=24 n1=1500 out=stdout   > iput.H 	
	< iput.H ${JX}.DipXGradNStat               > oput.H
	< iput.H ${JX}.DipFillerNStat              > dipy.H
	< iput.H ${TPLT}  out=junk2.v              > junk.V
	< oput.H ${TPLT}  out=junk1.v              > junk.V
	< dipy.H ${TPLT}  out=junk0.v              > junk.V
	vp_SideBySideAniso junk2.v junk1.v junk0.v > $@ 

TPLT = Taplot > junk.A; <junk.A Ta2vplot wantframe=n wantaxis=n lines=0 transp=n  title=" " crowd=0.9

clean : jclean

include ${SEPINC}/SEP.bottom
