include ${SEPINC}/SEP.top

RES_PATH = /net/santorin/scr2/louis/6planes-2nd
RES_PATH2 = /net/santorin/scr4/louis/6planes-2nd

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
P = ./Par

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

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

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

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

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

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

#6planes-pad.H: 6planes.HH
#	Cp 6planes.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=32 > $@ \
#	datapath=$(RES_PATH)/
#	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_PATH2)/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_PATH)/DATA_
#	echo label1="cmp-x" label2="cmp-y" >> $@
#	echo label3="offset-x" label4="offset-y" >> $@
#	echo label5="freq" >> $@

#--------------------------------
6planes-NApad.H: 6planes.HH
	Cp 6planes.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" >> $@

#---------------------------------
#data-freq-hy-pad-straight.H: data-freq-transp.H
#	Cp data-freq-transp.H junk.H
#	< junk.H /paul/plib/bin/SGI64/FakeAxis \
#	axis=4 nn=1 no=0. nd=25. >> 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=32 > $@ \
#	datapath=$(RES_PATH)/
#	Rm junk.H

#---------------------------------
spike.H:
	Spike n1=512 n2=48 n3=1024 d1=0.004 d2=25 d3=25 \
	k1=250 k2=24 k3=528 > $@
	echo "n3=32 n4=32 d4=25" >> $@

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

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

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

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

# idem = no difference
#CA-mig.H: data-freq-transp.H slowness.H $(P)/CA-mig.P
CA-mig.H: slowness.H $(P)/CA-mig.P
	/bin/time $(BINDIR)/CAM.x par=$(P)/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 $(P)/full-mig.P
	/bin/time $(BINDIR)/Full3D_salt.x par=$(P)/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-transp-hy-pad.H slowness.H $(P)/full-hy-mig.P
	/bin/time $(BINDIR)/Full3D_salt.x par=$(P)/full-hy-mig.P \
	Slow=slowness.H \
	Data=data-freq-transp-hy-pad.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 $(P)/fulltst-hy-mig.P
#	/bin/time $(BINDIR)/Full3Dtst_salt.x par=$(P)/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 $(P)/NA-mig.P
NA-mig.H: NA-mig.P
	/bin/time $(BINDIR)/NAM.x par=$(P)/NA-mig.P \
	Slow=slowness.H \
	Data=data-freq-hy-NApad-transp.H \
	>NA-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-2.H: $(P)/NA-mig.P
	/bin/time $(BINDIR)/NAM-2.x par=$(P)/NA-mig.P \
	Slow=slowness.H \
	Data=data-freq-hy-NApad-transp.H \
	>NA-mig-head-2.H \
	Image=$@ \
	datapath=$(RES_PATH)/DATA_ \
	2>$@.log

# NAM DSR on common-azimuth narrow-azimuth padded data (5-D)
# CIGs are px projections
NA-dat-PS-mig.H: $(P)/NA-mig.P
	/bin/time $(BINDIR)/Full3D_salt.x par=$(P)/NA-mig.P \
	Slow=slowness.H \
	Data=data-freq-hy-NApad-transp.H \
	>NA-dat-PS-mig-head.H \
	Image=$@ \
	datapath=$(RES_PATH)/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

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

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

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

.PRECIOUS: 

clean:jclean	

include ${SEPINC}/SEP.bottom