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

RESULTOLD =     gulfFoltZomRaw    nseaFoltCmbRaw 
RESULTSER = 	gulfFoltTotRawVol nseaFoltTotRawVol \
		gulfFoltTotRaw    nseaFoltTotRaw \
		gulfFoltFltRaw    nseaFoltFltRaw \
		                  nseaFoltZomRaw \
		gulfFoltTotMrk    nseaFoltTotMrk

DATDIR = ../Data
DAT    = /matt/Data/mig3d

cplt = Cubeplot movie=0 wheretitle="b" titlesz=12 titlefat=1 flat=n wantaxes=n wantcoordlabel=y point1=.5 point2=.66 frame1=${nseaFoltZ} frame2=${nseaFoltX} frame3=${nseaFoltY} title=" " 
${RESDIR}/%TotRawVol.v : %TotRaw.p ${DAT}/%.HH plot.p 
	<${DAT}/$*.HH          ${WIND} par=$*TotRaw.p           > junk0.H
	<junk0.H               ${JX}.Normalisor                 > junk1.H
	<junk1.H               ${TMA}                           > junk1.A
	<junk1.A               ${cplt} wantcoordlabel=-1 out=$@ > junk.V    

.PRECIOUS: ${DATDIR}/gulfFoltTotRaw.HH ${DATDIR}/nseaFoltTotRaw.HH
${DATDIR}/%TotRaw.HH : ${DAT}/%.HH %TotRaw.p 
	<${DAT}/$*.HH          ${WIND} par=$*TotRaw.p > junk.H
	<junk.H                ${JX}.Normalisor       > ${DATDIR}/$*TotRaw.HH 

${RESDIR}/%TotRaw.v : ${DATDIR}/%TotRaw.HH
	< $< ${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

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

.PRECIOUS: ${DATDIR}/gulfFoltZomRaw.HH ${DATDIR}/nseaFoltZomRaw.HH
${DATDIR}/%TotRaw.HH: ${DAT}/%.HH %ZomRaw.p
	<${DAT}/%.HH ${WIND} par=$*ZomRaw.p >junk.H
	<junk.H      ${JX}.Normalisor       >$@
zomTit = title=" "
${RESDIR}/%ZomRaw.v: ${DATDIR}/%ZomRaw.HH %ZomRaw.p  plot.p
	<${DATDIR}/$*ZomRaw.HH ${TRAN}                   >tran.H
	<tran.H                ${GMA}  ${zomTit} out=$@3 >junk.V 
	<${DATDIR}/$*ZomRaw.HH ${TMA}                    >junk.A
	<junk.A ${CPLT} par=$*ZomRaw.p ${zomTit} out=$@  >junk.V 

.PRECIOUS: ${DATDIR}/gulfFoltFltRaw.HH ${DATDIR}/nseaFoltFltRaw.HH
${DATDIR}/%FltRaw.HH: %FltRaw.p ${DAT}/%.HH
	<${DAT}/$*.HH ${WIND} par=$*FltRaw.p    >junk.H
	<junk.H       ${JX}.Normalisor          >$@

nseaFoltTit = "North Sea fault"
gulfFoltTit = "Gulf Salt fault"
fltTit = title=" "
cube = Cubeplot movie=0 titlesz=10 titlefat=2 flat=n wantaxes=n wantcoordlabel=-2 wantoutline=y point1=.66 point2=.66 frame1=${nseaFoltZ} frame2=${nseaFoltX} frame3=${nseaFoltY} title=${$*Tit}
lplt = xll=1    yll=1 xur=6.1   yur=9.5
lplt = xll=1    yll=1 xur=6.1   yur=8.75
rplt = xll=7.55 yll=1 xur=12.65 yur=9.5
${RESDIR}/%FltRaw.v : ${DATDIR}/%FltRaw.HH %FltRaw.p  plot.p
	<${DATDIR}/$*FltRaw.HH ${TRAN}                   >tran.H
	<tran.H                ${GMA}  ${fltTit} out=$@3 >junk.V 
	<${DATDIR}/$*FltRaw.HH ${TMA}                    >junk.A
	<junk.A                ${cube} ${lplt} out=cube.v >junk.V 
	<${DATDIR}/$*FltRaw.HH Cp  > junk.H
	echo o1=0 o2=0 label1=sec label2=" " >> junk.H 
	<junk.H                Spectrum o1=0              >spec.H
	<spec.H Graph labelsz=8 titlesz=10 titlefat=2 ${rplt} out=ampl.v >junk.V
	vppen erase=n vpstyle=n ampl.v cube.v > $@

cuby = Cubeplot movie=0 wheretitle=b titlesz=17 titlefat=2 flat=n wantaxes=n wantcoordlabel=-2 wantoutline=y point1=.66 point2=.66 frame1=${nseaFoltZ} frame2=${nseaFoltX} frame3=${nseaFoltY} par=$*FltRaw.p out=cube.v
${RESDIR}/%CmbRaw.v : ${DATDIR}/%FltRaw.HH ${DATDIR}/%ZomRaw.HH 
	<${DATDIR}/$*FltRaw.HH ${TMA} >junk.A
	<junk.A ${cuby} titlesz=20 title="Slump fault" out=flt.v > junk.V 
	<${DATDIR}/$*ZomRaw.HH ${TMA} >junk.A
	<junk.A ${cuby} titlesz=20 title="Flat layers" out=zom.v > junk.V 
	vp_SideBySideAniso zom.v flt.v       > $@

WIND = Window out=stdout 

### TEMP ###

junk:
	#< ${DATDIR}/nseaFoltFltRaw.HH Pad n1out=64 n2out=32 n3out=32 > padd.H
	Spike n1=16 n2=4 n3=8 > padd.H 
	< padd.H Rtoc - > cmpx.H 
	< cmpx.H Ft3d sign1=1 center1=1 > fft1.H
	< fft1.H Ft3d sign2=1 center2=1 > fft2.H
	< fft2.H Ft3d sign3=1 center3=1 > fft3.H
	< fft3.H Cabs > cabs.H
	< cabs.H Window n3=1 f3=16 > wind.H
	< wind.H ${TPA} > junk.V
	Tube < junk.V	

shift: 
	#< /matt/Data/mig3d/gulfFolt.HH Cp > /scrka1/matt/junkhim.H 
	#< /scrka1/matt/junkhim.H Transp plane=23 > /scrka1/matt/junkme.H
	#echo '## next line added by hand. Matthias. ##' >> /scrka1/matt/junkme.H
	#echo 'label1="Time" label2="North-South" label3="East-West"' >> /scrka1/matt/junkme.H
	#< /scrka1/matt/junkme.H Cp out=stdout   > /scrka1/matt/transpose.H
	< /scrka1/matt/transpose.H       Window f1=4 n1=1 > /scrka1/matt/time.H
	< /scrka1/matt/time.H Grey > /scrka1/matt/time.V
	< /scrka1/matt/time.V Tube 

labelchange:
	< /matt/Data/mig3d/gulfFolt.HH Cp                             > /scrka1/matt/junkme.H
	echo '## next line added by hand. Matthias. ##'              >> /scrka1/matt/junkme.H
	echo 'label1="Time" label2="In-line" label3="Cross-line"'    >> /scrka1/matt/junkme.H
	echo 'd1=0.004'                                              >> /scrka1/matt/junkme.H
	< /scrka1/matt/junkme.H Cp out=stdout   > /scrka1/matt/transpose.H
	< /scrka1/matt/transpose.H       Window f1=4 n1=1 > /scrka1/matt/time.H
	< /scrka1/matt/time.H Grey > /scrka1/matt/time.V
	< /scrka1/matt/time.V Tube 


clean: jclean
	${RM} spike.T stdin_picks

distclean: clean
	${TOUCH} junk.HH
	${RM}       *.HH 

include ${SEPINC}/SEP.bottom






