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

RESULTOLD = comparison 
RESULTSER = zeroFoltPP zeroFoltPPC
RESULTSCR = title gulfFoltTotPPC nseaFoltTotPPC 

# PPS 	Project Project Substract
# PPSN 	Project Project Substract Norm 
# PPSLN Project Project Substract PixelWiseNorm 
# PPS 	Project Project Semblance or Correlatioin 
# PP 	Project Project DisplayPlaneWave

lftplt = 'allpos=n title="Gulf Salt dome .." xll=1    yll=1 xur=6.3   yur=9.5'
rgtplt = 'allpos=y title=".. and its faults" xll=7.35 yll=1 xur=12.65 yur=9.5'
stdplt =  Window out=stdout n1=1 f1=15 >junkw.H; <junkw.H Grey titlesz=8 labelsz=4 
${RESDIR}/title.v : gulfFoltTotPPC.H ${DAT}/gulfFoltTotRaw.HH
	< ${DAT}/gulfFoltTotRaw.HH ${stdplt} ${lftplt} > junk.V  out=left.v
	<        gulfFoltTotPPC.H  ${stdplt} ${rgtplt} > junk.V  out=rght.v
	vppen erase=n vpstyle=n rght.v left.v > $@

topplt = 'allpos=n title="Gulf Salt dome .." xll=1 yll=5.9 xur=12.65 yur=9.5'
botplt = 'allpos=y title=".. and its faults" xll=1 yll=1   xur=12.65 yur=4.6'
${RESDIR}/comparison.v : gulfFoltTotPPC.H ${DAT}/gulfFoltTotRaw.HH
	< ${DAT}/gulfFoltTotRaw.HH ${stdplt} ${topplt} > junk.V  out=top.v
	<        gulfFoltTotPPC.H  ${stdplt} ${botplt} > junk.V  out=bot.v
	vppen erase=n vpstyle=n top.v bot.v > $@

PTCHARG = olap=.0f,.0f,0.0f psze=20,10,10  
${RESDIR}/zeroFoltPP.v :   ${DAT}/zeroFoltRaw.HH 
	<$<     ${JX}.ProProNStat ${PTCHARG}                            >junk.H
	<junk.H ${TMA}                                                  >junk.A
	<junk.A ${CPLT} ${PFRAME} titlesz=20 title="Plane wave estimates" out=j1.v >junk.V
	<$< ${TMA}                                                  >junk.A
	<junk.A ${CPLT} ${PFRAME} titlesz=20 title="Synthetic fault"      out=j0.v >junk.V
	vp_SideBySideAniso j0.v j1.v                                    >$@

ZEROARG = olap=.5f,.5f,.5f psze=7,3,3  
${RESDIR}/zeroFoltPPS.v :   ${DAT}/zeroFoltRaw.HH 
	< $<     ${JX}.ProProSubNStat ${ZEROARG}                       > junk.H
	< junk.H ${TMA}                                                > junk.A
	< junk.A ${CPLT} ${PFRAME} title="Wave: Difference" out=$@     > junk.V
${RESDIR}/zeroFoltPPC.v :   ${DAT}/zeroFoltRaw.HH 
	< $<     ${JX}.ProProCorNStat ${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=" " out=$@    > junk.V

#NSEAARG = olap=.5f,.5f,.5f psze=9,9,9
NSEAARG = olap=.5f,.5f,.5f psze=7,3,3
.PRECIOUS: gulfFoltTotPPC.H nseaFoltTotPPC.H 
%TotPPC.H : ${DAT}/%TotRaw.HH 
	< $<     ${JX}.ProProCorNStat ${NSEAARG}        > lapn.H
	< lapn.H Add scale=-1                           > scal.H
	< scal.H Add add=1                              > $@
ppcarg = 'title="Correlated wave estimate"'
${RESDIR}/%TotPPC.v : %TotPPC.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

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

### OLD ###
### PPS: 
### The difference makes for funny singular events. I don't think it is a 
### bug, but believe it has to do with the wrong dip estimate in the case
### of more than a single event. But I don't really know. 
###

.PRECIOUS: gulfFoltTotPPS.H 
%TotPPS.H : ${DAT}/%TotRaw.HH 
	< $<     ${JX}.ProProSubNStat ${NSEAARG} > $@
ppsarg = 'title="Subtracted wave estimate"'
${RESDIR}/%TotPPSTop.v : %TotPPS.H
	< $<     ${TRAN}                               > tran.H
	< tran.H ${GMA} ${TOPGREY} ${ppsarg}   out=$@3 > junk.V
	< $<     ${ZSID}           ${ppsarg}   out=$@  > junk.V
${RESDIR}/%TotPPSSid.v : %TotPPS.H
	< $< ${XSID}  > junk.V  out=xplt.v 
	< $< ${YSID}  > junk.V  out=yplt.v 
	vppen erase=n vpstyle=n xplt.v yplt.v > $@


clean : jclean
include ${SEPINC}/SEP.bottom

