include ${SEPINC}/SEP.top

VEL_PATH = /net/kana/book/3DI/Data/SEG/Salt/Vel
STACK_FREQ_PATH = /net/kana/book/3DI/Tmp/split/SEG/Salt

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

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: grad.P
	$(BINDIR)/slow_grad.x par=grad.P > $@ datapath=$(RES_PATH)/

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

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

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

data-sub-freq-transp.H: 5planes-4th.HH
	Window3d < 5planes-4th.HH j1=2 > j1.H
	$(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

data-freq-transp-hy-pad.H: data-sub-freq-transp.H
	Cp data-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_PATH)/
#	datapath=$(RES_PATH)/ nfiles=18 filesize=2000
	Rm junk.H

#--------------------------------
#5planes-NApad.H: 5planes-4th.HH
#	Cp 5planes-4th.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: 5planes-NApad.H
#	$(3DSI-BINDIR)/Transf f_min=0 f_max=55 dim=5 \
#	< 5planes-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" >> $@
#
#---------------------------------
#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 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

# idem = no difference
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

# idem = no difference
CA-mig-linux.H: data-freq-transp.H slowness.H CA-mig-linux.P
	$(BINDIR)/CAM.x par=CA-mig-linux.P \
	Slow=slowness.H \
	Data=data-freq-transp.H \
	>CA-mig-linux-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: data-freq-transp.H slowness.H OPW-mig.P
	/bin/time $(BINDIR)/Full3D_salt.x par=OPW-mig.P \
	Slow=slowness.H \
	Data=data-freq-transp.H \
	>OPW-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
OPW-mig-linux.H: data-freq-transp.H slowness.H OPW-mig.P
	time $(BINDIR)/Full3D_salt.x par=OPW-mig.P \
	Slow=slowness.H \
	Data=data-freq-transp.H \
	>OPW-mig-linux-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-linux.H: data-freq-transp-hy-pad.H slowness.H full-hy-mig.P
	time $(BINDIR)/Full3D_salt.x par=full-hy-mig.P \
	Slow=slowness.H \
	Data=data-freq-transp-hy-pad.H \
	>full-hy-mig-linux-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: data-freq-transp-hy-pad.H slowness.H full-hy-mig.P
	/bin/time $(BINDIR)/Full3D_salt.x par=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

# 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

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

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

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

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

ZO-mig.H: data-zo-freq-transp.H slowness.H ZO-mig.P
	/bin/time $(BINDIR)/CAM.x par=ZO-mig.P \
	Slow=slowness.H \
	Data=data-zo-freq-transp.H \
	>ZO-mig-head.H \
	Image=$@ \
	datapath=$(RES_PATH)/DATA_ \
	2>$@.log

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


#--------------------------------------
# delayed crap

do-it-all-test:
	sleep 10
	cp ../../avant.v test.v
	gmake slowness.H

do-it-all:
	sleep 115200
	cp ../../4th-data.HH 5planes-4th.HH
	echo 'o1=0. gff=-1 hff=-1' >> 5planes-4th.HH
	gmake CA-mig.H
	gmake OPW-mig.H
	gmake full-hy-mig.H

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

clean:jclean	

include ${SEPINC}/SEP.bottom
