#define RM_H /bin/rm
$#define RM_H Rm
#define USRCMLIB -L/usr/local/seplibV2/lib -lsepattr -lsepCM5
#define CMFLAGS -Bstatic
#define XTPANELDIR ../../Xtpanel
#define FIGLIST wnrecicorr 
#include <SEP.defs>
#include "../../Cakerules/SEP.martin.rules"
#define BINDIR ../../Bin/MTYPE
#define FIGDIR ../Fig
#define FIGSUP ../Fig2
#define CUBELABEL label1="Time [sec]" label2="Midpoint [feet]" label3="Offset [feet]"

#define PEMFULL /SDA/martin/PemFull/CDP

#define ATTR attrfile=PEMFULL
#define SCALE scale=1e-15

/* normalize, this is the figure in the thesis */
FIGDIR/wnrecicorr.v: [[ sub -i X wncorX.3-4.T12.H NINECOMP ]] 
	Merge wncorxx.3-4.T12.H wncorxy.3-4.T12.H wncorxz.3-4.T12.H axis=2 space=y >j1.H
	Merge wncoryx.3-4.T12.H wncoryy.3-4.T12.H wncoryz.3-4.T12.H axis=2 space=y >j2.H
	Merge wncorzx.3-4.T12.H wncorzy.3-4.T12.H wncorzz.3-4.T12.H axis=2 space=y >j3.H
	Merge j1.H j2.H j3.H axis=1 space=y >jj.H
	Taplot pclip=90 <jj.H |Ta2vplot wantlabel=n  wantaxis=n wanttitle=n title="Zero Lag Recip. Corr."  out=FIGDIR/wnrecicorr.v NULL


FIGDIR/wnrecicorr.v3&: [[ sub -i X FIGSUP2/wncorX.3-4.T12.v NINECOMP ]]
	
FIGSUP2/wncor%1.3-4.T12.v: /* wncor%1.3-4.T12.H */
	Taplot pclip=90 <wncor%1.3-4.T12.H |Ta2vplot title=" %1 Zero Lag Recip. Corr." out=FIGSUP/wncor%1.3-4.T12.v  NULL

normalized&: [[ sub -i X wncorX.3-4.H NINECOMP ]]  /* NINECOMP ]] */


zerolag&: [[ sub -i X cor2X NINECOMP ]]

windowlag&: [[ sub -i X wncorX.3-4.H xy xz yx yy yz zx zy zz ]]  /* NINECOMP ]] */

window&: [[ sub -i X windX.3-4.H NINECOMP ]]

cor2%&:
	time BINDIR/Corr2 SCALE gpow=0.3 < PEMFULL/%cube.H >cor2%.H

cor%&:
	time BINDIR/Corr scale=1.e-15 gpow=0.3 < PEMFULL/%cube.H >cor%.H

pele&:
	BINDIR/Corr2 SCALE <PEMFULL/xxcube.H >junk2

wind%.3-4.H:
	Window < PEMFULL/%cube.H f1=750 n1=250 >wind%.3-4.H

wcor%.3-4.H:  wind%.3-4.H 
	BINDIR/Corr2 SCALE <wind%.3-4.H >wcor%.3-4.H

wncor%.3-4.H:    wind%.3-4.H   BINDIR/Corr3
	BINDIR/Corr3 SCALE <wind%.3-4.H >wncor%.3-4.H

%.T12.H:    %.H  
	Transp plane=12 <%.H >%.T12.H

cmcorr%1&:
	BINDIR/Corr < PEMFULL/%1cube.H attrfile=PEMFULL/h%1cube.H 



recicorr.v: 
	Merge cor2xx.H cor2xy.H cor2xz.H axis=2 space=y >j1.H
	Merge cor2yx.H cor2yy.H cor2yz.H axis=2 space=y >j2.H
	Merge cor2zx.H cor2zy.H cor2zz.H axis=2 space=y >j3.H
	Merge j1.H j2.H j3.H axis=1 space=y >jj.H
	Window n1=1 <jj.H >jjw.H
	Taplot pclip=90 <jjw.H |Ta2vplot wantlabel=n title="Zero Lag Recip. Corr."  out=recicorr.v NULL

FIGDIR/wrecicorr.v: 
	Merge wcorxx.3-4.H wcorxy.3-4.H wcorxz.3-4.H axis=2 space=y >j1.H
	Merge wcoryx.3-4.H wcoryy.3-4.H wcoryz.3-4.H axis=2 space=y >j2.H
	Merge wcorzx.3-4.H wcorzy.3-4.H wcorzz.3-4.H axis=2 space=y >j3.H
	Merge j1.H j2.H j3.H axis=1 space=y >jj.H
	Taplot pclip=90 <jj.H |Ta2vplot wantlabel=n  wantaxis=n wanttitle=n title="Zero Lag Recip. Corr."  out=FIGDIR/wrecicorr.v NULL


FIGDIR/wnrecicorr2.v: 
	Merge wncorxx.3-4.H wncorxy.3-4.H wncorxz.3-4.H axis=2 space=y >j1.H
	Merge wncoryx.3-4.H wncoryy.3-4.H wncoryz.3-4.H axis=2 space=y >j2.H
	Merge wncorzx.3-4.H wncorzy.3-4.H wncorzz.3-4.H axis=2 space=y >j3.H
	Merge j1.H j2.H j3.H axis=1 space=y >jj.H
	Taplot pclip=90 <jj.H |Ta2vplot wantlabel=n  wantaxis=n wanttitle=n title="Zero Lag Recip. Corr."  out=FIGDIR/wnrecicorr.v NULL

/* do the NMO */

cdp%1cubebig.NMO.%2.H : PEMFULL/%1cube.H
	NMO vel=%2 <PEMFULL/%1cube.H >cdp%1cubebig.NMO.%2.H

BINDIR/Corr3: Corr3.o
	FLINK(Corr3.o,BINDIR/Corr3)

clean&: jclean
	

#include <SEP.obj.rules>
#include <SEP.prog.rules>
#include <SEP.idoc.rules>

