include ${SEPINC}/SEP.top

RES_PATH = /net/santorin/scr2/louis/6planes-3rd
RES_PATH2 = /net/santorin/scr3/louis/6planes-3rd
RES_PATH3 = /net/santorin/scr4/louis/6planes-3rd

BINDIR = /net/kana/louis/Bin/$(SEP_ARCH)
3DSI-BINDIR = /net/kana/book/3DI/Bin/$(SEP_ARCH)

TT  = Grey eout=1 pclip=95 gainpanel=a
T12 = Transp core=50000 plane=12
T23 = Transp core=50000 plane=23

######################## Prepare slowness ############################

slowness.H: grad.P
	$(BINDIR)/slow_grad.x par=grad.P > $@ datapath=$(RES_PATH)/

######################### Prepare data ################################

data-freq.H : 6planes-3rd.HH
	/net/kana/book/3DI/Bin/SGI64/Transf f_min=0 f_max=55 dim=4 \
	< 6planes-3rd.HH > $@ \
	out=$(RES_PATH)/DATA_$(@F)@ \
	2>$@.log

data-freq-transp.H: data-freq.H
	$(T12) < data-freq.H | $(T23) > $@ datapath=$(RES_PATH)/

6planes-pad.H: 6planes-3rd.HH
	Cp 6planes-3rd.HH junk.H
	< junk.H /net/kana/paul/plib/bin/SGI64/FakeAxis \
	axis=3 nn=1 no=0. nd=30. >> junk.H
	echo label1="time" label2="offset-x" >> junk.H
	echo label3="offset-y" label4="cmp-x" >> junk.H
	echo label5="cmp-y" >> junk.H
	$(BINDIR)/Zero_pad.x < junk.H n3out=16 > $@ \
	datapath=$(RES_PATH2)/
	Rm junk.H

data-freq-hy-pad.H: 6planes-pad.H
	$(3DSI-BINDIR)/Transf f_min=0 f_max=55 dim=5 \
	< 6planes-pad.H > $@ \
	out=$(RES_PATH3)/DATA_$(@F)@ \
	2>$@.log

data-freq-hy-pad-transp.H: data-freq-hy-pad.H
	/net/kana/paul/plib/bin/SGI64/F12T34 data-freq-hy-pad.H $@ \
	datapath=$(RES_PATH2)/DATA_
	echo label1="cmp-x" label2="cmp-y" >> $@
	echo label3="offset-x" label4="offset-y" >> $@
	echo label5="freq" >> $@

#-------narrow-azimuth data------------------
6planes-NApad.H: 6planes-3rd.HH
	Cp 6planes-3rd.HH junk.H
	< junk.H /net/kana/paul/plib/bin/SGI64/FakeAxis \
	axis=3 nn=1 no=0. nd=25. >> junk.H
	echo label1="time" label2="offset-x" >> junk.H
	echo label3="offset-y" label4="cmp-x" >> junk.H
	echo label5="cmp-y" >> junk.H
	$(BINDIR)/Zero_pad.x < junk.H n3out=5 > $@ \
	datapath=$(RES_PATH)/
	Rm junk.H

data-freq-hy-NApad.H: 6planes-NApad.H
	$(3DSI-BINDIR)/Transf f_min=0 f_max=55 dim=5 \
	< 6planes-NApad.H > $@ \
	out=$(RES_PATH2)/DATA_$(@F)@ \
	2>$@.log

data-freq-hy-NApad-transp.H: data-freq-hy-NApad.H
	/net/kana/paul/plib/bin/SGI64/F12T34 data-freq-hy-NApad.H $@ \
	datapath=$(RES_PATH)/DATA_
	echo label1="cmp-x" label2="cmp-y" >> $@
	echo label3="offset-x" label4="offset-y" >> $@
	echo label5="freq" >> $@

#########################  Beginning computation #####################

# common-azimuth DSR - CIG are px
comaz-mig.H: data-freq-transp.H slowness.H comaz-mig.P
	/bin/time $(3DSI-BINDIR)/ComAz_salt.x par=comaz-mig.P \
	Slow=slowness.H \
	Data=data-freq-transp.H \
	>comaz-mig-head.H \
	Image=$@ \
	datapath=$(RES_PATH)/DATA_ \
	2>$@.log

# common-azimuth DSR - CIG are px
CA-mig.H: data-freq-transp.H slowness.H CA-mig.P
	/bin/time $(BINDIR)/CAM.x par=CA-mig.P \
	Slow=slowness.H \
	Data=data-freq-transp.H \
	>CA-mig-head.H \
	Image=$@ \
	datapath=$(RES_PATH)/DATA_ \
	2>$@.log

# full 3-D DSR on common-azimuth data (4-D,  no padding)
# CIGs are px
full-mig.H: data-freq-transp.H slowness.H full-mig.P
	/bin/time $(BINDIR)/Full3D_salt.x par=full-mig.P \
	Slow=slowness.H \
	Data=data-freq-transp.H \
	>full-mig-head.H \
	Image=$@ \
	datapath=$(RES_PATH)/DATA_ \
	2>$@.log

# full 3-D DSR on common-azimuth padded data (5-D)
# CIGs are px projections 
full-hy-mig.H: data-freq-hy-pad-transp.H slowness.H full-hy-mig.P
	/bin/time $(BINDIR)/Full3D_salt.x par=full-hy-mig.P \
	Slow=slowness.H \
	Data=data-freq-hy-pad-transp.H \
	>full-hy-mig-head.H \
	Image=$@ \
	datapath=$(RES_PATH)/DATA_ \
	2>$@.log

# full 3-D DSR on common-azimuth padded data (5-D)
# CIGs are |p| projections on px axis
fulltst-hy-mig.H: data-freq-hy-pad-transp.H slowness.H fulltst-hy-mig.P
	/bin/time $(BINDIR)/Full3Dtst_salt.x par=fulltst-hy-mig.P \
	Slow=slowness.H \
	Data=data-freq-hy-pad-transp.H \
	>fulltst-hy-mig-head.H \
	Image=$@ \
	datapath=$(RES_PATH)/DATA_ \
	2>$@.log

# NAM DSR on common-azimuth narrow-azimuth padded data (5-D)
# CIGs are px projections 
NA-mig.H: data-freq-hy-NApad-transp.H slowness.H NA-mig.P
	/bin/time $(BINDIR)/NAM.x par=NA-mig.P \
	Slow=slowness.H \
	Data=data-freq-hy-NApad-transp.H \
	>NA-mig-head.H \
	Image=$@ \
	datapath=$(RES_PATH)/DATA_ \
	2>$@.log

#############################  Beginning graphics ############################

%.T : %.H
	$(TT) < $*.H > $@ datapath=$(RES_PATH)/DATA_

#############################  Beginning misc ############################

.PRECIOUS: 

clean:jclean	

include ${SEPINC}/SEP.bottom