include ${SEPINC}/SEP.top

RES_PATH  = /net/santorin/scr2/louis/6planes-4th
RES_PATH2 = /net/santorin/scr4/louis/6planes-4th
RES_PATHL = /scr1/louis

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

T12 = Transp core=50000 plane=12
T23 = Transp core=50000 plane=23

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

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

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

debugdata-win.H: ../5planes-4th.HH
	Window3d < ../5planes-4th.HH > $@ \
	n1=1750 nfiles=5 filesize=2000 synch=1

debugdata-freq.H: debugdata-win.H
	/net/kana/book/3DI/Bin/SGI64/Transf f_min=0 f_max=70 dim=4 \
	< debugdata-win.H > $@ \
	out=/net/kana/scr1/louis/DATA_$(@F)@ \
	2>$@.log

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

debugdata-sub-freq-transp.H: debugdata-win.H
	Window3d < debugdata-win.H j1=2 > j1.H datapath=$(RES_PATH)/
	$(3DSI-BINDIR)/Transf f_min=0 f_max=60 dim=4 < j1.H > j2.H \
	out=$(RES_PATH)/DATA_$(@F)@ \
	2>j2.H.log
	$(T12) < j2.H | $(T23) > $@ datapath=$(RES_PATH)/
	Rm j1.H
	Rm j2.H

debugdata-freq-transp-hy-pad.H: debugdata-sub-freq-transp.H
	Cp debugdata-sub-freq-transp.H junk.H
	< junk.H /net/kana/paul/plib/bin/SGI64/FakeAxis \
	axis=4 nn=1 no=0. nd=50. >> junk.H
	echo label1="cmp-x" label2="cmp-y" >> junk.H
	echo label3="offset-x" label4="offset-y" >> junk.H
	echo label5="freq" >> junk.H
	$(BINDIR)/Zero_pad.x < junk.H n4out=16 > $@ \
	datapath=$(RES_PATH2)/ nfiles=16 filesize=2000
#	datapath=$(RES_PATH)/ nfiles=18 filesize=2000
	Rm junk.H



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

# idem = no difference
CA-mig.H: debugdata-freq-transp.H slowness.H $P/CA-mig.P
	time $(BINDIR)/CAM.x par=$P/CA-mig.P \
	Slow=slowness.H \
	Data=debugdata-freq-transp.H \
	>CA-mig-head.H \
	Image=$@ \
	datapath=$(RES_PATHL)/DATA_ \
	2>$@.log

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

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

# NAM DSR on common-azimuth narrow-azimuth padded data (5-D)
# CIGs are px projections 

#NA-mig.H: debugdata-freq-hy-NApad-transp.H slowness.H $P/NA-mig.P
NA-mig.H: debugdata-freq-transp-hy-pad.H slowness.H $P/NA-mig.P
	time $(BINDIR)/NAM.x par=$P/NA-mig.P \
	Slow=slowness.H \
	Data=debugdata-freq-transp-hy-pad.H \
	>NA-mig-head.H \
	Image=$@ \
	datapath=$(RES_PATHL)/DATA_ \
	2>$@.log

# impulse response (test)
spike-mig.H: spike-freq-transp.H slowness.H $P/spike-mig.P
	time $(3DSI-BINDIR)/ComAz_salt.x par=$P/spike-mig.P \
	Slow=slowness.H \
	Data=spike-freq-transp.H \
	>spike-mig-head.H \
	Image=$@ \
	datapath=$(RES_PATH)/DATA_ \
	2>$@.log

#--------------------------------------
.PRECIOUS: 

clean:jclean

include ${SEPINC}/SEP.bottom