include ${SEPINC}/SEP.top

PERSDEV=../..
HOME=$(PERSDEV)

MYBINDIR = $(PERSDEV)/bin/$(SEP_ARCH)
BINDIR = $(MYBINDIR)

ifndef INPMYHOST
INPMYHOST=$(MYHOST)
endif

SCRNUM=1
SCRNUM_sad=2
SCRNUM_glad=2
SCRNUM_mad=2
SCRNUM_bad=2
SCRDIR=WIT_WAVRES_TRIP_LENS
ifndef HOSTFILE
HOSTFILE= $(HOME)/Lib/makefile/Hostfile
endif
HOST = glad
#include $(HOSTFILE)

MODEL_DIR=../Model
MIG_DIR=../Mig
DATADIR=/data/trip_lens

TEMPDIR=$(MODEL_DIR)/SCR
TEMPPATH=${TEMPDIR}/DATA_

include makefile-figrep


ifeq ($(MYHOST),sad)
#MAINREV=Rev2D_spectral
MAINREV=Rev2D_ho_fd
else
MAINREV=Rev2D_ho_fd
MAINREV-x2-t2=Rev2D_ho_fd
MAINREV-x4-t2=Rev2D_ho_fd
MAINREV-x10-t2=Rev2D_ho_fd
MAINREV-x10-t4=Rev2D_ho_fd
MAINREV-spect=Rev2D_spectral
endif

Model-data-win-$(MYHOST).H: Model-data-rev-time-$(MYHOST).H Win-data.P
	<Model-data-rev-time-$(MYHOST).H Window3d par=Win-data.P \
	datapath=$(TEMPPATH)  \
	>$@

Model-data-VLowFreq-$(MYHOST).H: VLowFreq-data.P \
			Model-data-rev-time-$(MYHOST).H 
	<Model-data-rev-time-$(MYHOST).H \
	Transp plane=12 | \
	Bandpass par=VLowFreq-data.P | \
	Transp plane=12 \
	datapath=$(TEMPPATH)  \
	>$@


Model-data-LowFreq-$(MYHOST).H: LowFreq-data.P #Model-data-rev-time-$(MYHOST).H 
	<Model-data-rev-time-$(MYHOST).H \
	Transp plane=12 | \
	Bandpass par=LowFreq-data.P | \
	Transp plane=12 \
	datapath=$(TEMPPATH)  \
	>$@

Model-data-test-$(MYHOST).H: Test-data.P #Model-data-rev-time-$(MYHOST).H 
	<Model-data-rev-time-$(MYHOST).H Window3d par=Test-data.P | \
	Transp plane=12 | \
	Bandpass par=Test-data.P | \
	Transp plane=12 \
	datapath=$(TEMPPATH)  \
	>$@

Model-data-rev-time-$(MYHOST).H: Model-rev-time.P Rev-time-x10-t4.P \
		Slow-compl-trip-lens-$(MYHOST).H  Wavelet-nz3-$(MYHOST).H 
	$(BINDIR)/$(MAINREV).x \
	par=Model-rev-time.P  par=Rev-time-x10-t4.P \
	Wavelet=Wavelet-nz3-$(MYHOST).H \
	Wavefield_sou=Wave-data-rev-time-$(MYHOST).H \
	Data=Model-data-rev-time-$(MYHOST).H \
	Slow=Slow-compl-trip-lens-$(MYHOST).H  \
	datapath=$(TEMPPATH) filesize=2000  nfiles=1\
	>Model-data-rev-time-$(MYHOST).H.head \
	2>$(TEMPDIR)/Model-data-rev-time-$(MYHOST).H.log
#	echo "z_data_0=0 tau_data_0=0 " >> $@



Lin-Wave-rev-time-sou-%-$(MYHOST).H \
Lin-Wave-rev-time-rec-%-$(MYHOST).H \
Lin-Wave-rev-time-back-%-$(MYHOST).H \
Lin-Model-rev-time-%-$(MYHOST).H: Lin-Model-rev-time-%.P \
		Slow-%.H DSlow-%.H  Wavelet-nz3-%-$(MYHOST).H \
		$(BINDIR)/$(MAINREV).x 
	$(BINDIR)/$(MAINREV).x \
	par=Lin-Model-rev-time-$*.P \
	Wavelet=Wavelet-nz3-$*-$(MYHOST).H \
	Wavefield_back=Lin-Wave-rev-time-back-$*-$(MYHOST).H \
	Wavefield_sou=Lin-Wave-rev-time-sou-$*-$(MYHOST).H \
	Wavefield_rec=Lin-Wave-rev-time-rec-$*-$(MYHOST).H \
	Data=Lin-Model-rev-time-$*-$(MYHOST).H \
	Slow=Slow-$*-data.H \
	Image=DSlow-$*.H \
	datapath=$(TEMPPATH) filesize=2000  nfiles=1\
	>Lin-Model-rev-time-$*-$(MYHOST).H.head \
	2>$(TEMPDIR)/Lin-Model-rev-time-$*-$(MYHOST).H.log
#	echo "z_data_0=0 tau_data_0=0 " >> $@

Lin-Mig-rev-time-sou-%-$(MYHOST).H \
Lin-Mig-rev-time-rec-%-$(MYHOST).H \
Lin-Mig-rev-time-back-%-$(MYHOST).H \
Lin-Mig-rev-time-%-$(MYHOST).H: Lin-Mig-rev-time-%.P \
		Slow-%.H Wavelet-nz3-%-$(MYHOST).H \
		$(BINDIR)/$(MAINREV).x \
		Lin-Model-rev-time-%-$(MYHOST).H 
	$(BINDIR)/$(MAINREV).x \
	par=Lin-Mig-rev-time-$*.P \
	Wavelet=Wavelet-nz3-$*-$(MYHOST).H \
	Wavefield_back=Lin-Mig-rev-time-back-$*-$(MYHOST).H \
	Wavefield_sou=Lin-Mig-rev-time-sou-$*-$(MYHOST).H \
	Wavefield_rec=Lin-Mig-rev-time-rec-$*-$(MYHOST).H \
	Data=Lin-Model-rev-time-$*-$(MYHOST).H \
	Slow=Slow-$*-data.H \
	Image=Lin-Mig-rev-time-$*-$(MYHOST).H \
	datapath=$(TEMPPATH) filesize=2000  nfiles=1 \
	>Lin-Mig-rev-time-$*-$(MYHOST).H.head \
	2>$(TEMPDIR)/Lin-Mig-rev-time-$*-$(MYHOST).H.log


Wavelet-nz1-$(MYHOST).H : Wavelet.P
	Wavelet par=Wavelet.P >tmp.H \
	datapath=$(TEMPPATH) 
	Scale <tmp.H rscale=.5 >tmp1.H \
	datapath=$(TEMPPATH) 
	Merge axis=3 space=n tmp1.H tmp.H tmp1.H >$@ \
	datapath=$(TEMPPATH) 
	echo "d2=5 o2=0" >> $@
	echo "d3=5 o3=-5" >> $@

Wavelet-nz3-$(MYHOST).H : Wavelet.P 
	rm -f tmp.H
	rm -f tmp1.H
	rm -f tmp2.H
	rm -f tmp3.H
#	sleep 2
	 /opt/SEP/SEP7.0/bin/Wavelet par=Wavelet.P >tmp.H \
	datapath=$(TEMPPATH) 
#	sleep 2
	Scale <tmp.H rscale=.5 >tmp1.H \
	datapath=$(TEMPPATH) 
#	sleep 2
	Merge axis=2 space=n tmp1.H tmp.H tmp1.H >tmp2.H \
	datapath=$(TEMPPATH) 
#	sleep 2
	Scale <tmp2.H rscale=.5 >tmp3.H \
	datapath=$(TEMPPATH) 
#	sleep 2
	Merge axis=3 space=n tmp3.H tmp2.H tmp3.H >$@ \
	datapath=$(TEMPPATH) 
	echo "d2=5 o2=-5" >> $@
	echo "d3=5 o3=-5" >> $@

# velocity and slowness


Slow-data.H  :Vel-data.P
	Vel par=Vel-data.P  \
	datapath=$(TEMPPATH) | Transp plane=12 |\
	Vconvert invelopt=1 outvelopt=2 \
	datapath=$(TEMPPATH)  \
	>$@

#Slow-trip-lens-$(MYHOST).H	: /data/trip_lens/velocity_xdr.H
Slow-trip-lens-$(MYHOST).H: $(DATADIR)/velocity_xdr.H
	<$< \
	Vconvert invelopt=1 outvelopt=2  \
	datapath=$(TEMPPATH) \
	>$@
	echo "d3=5 " >> $@

Slow-0-trip-lens-$(MYHOST).H: Vel-0-trip-lens.P
	Vel par=$< \
	datapath=$(TEMPPATH) \
	>$@

Slow-00-trip-lens-$(MYHOST).H: Vel-0-trip-lens.P Vel-00-trip-lens.P
	Vel par=$< par=Vel-00-trip-lens.P \
	datapath=$(TEMPPATH) \
	>$@

Slow-top-trip-lens-$(MYHOST).H: Vel-top-trip-lens.P
	Vel par=$< \
	datapath=$(TEMPPATH) \
	>$@

Slow-diff-%-trip-lens-$(MYHOST).H: \
		 Slow-trip-lens-$(MYHOST).H Slow-0-trip-lens-$(MYHOST).H
	Add scale=1,-1 \
	Slow-trip-lens-$(MYHOST).H Slow-0-trip-lens-$(MYHOST).H | \
	Scale rscale=$* \
	datapath=$(TEMPPATH) \
	>$@

Slow-sc-%-trip-lens-$(MYHOST).H: \
		Slow-diff-%-trip-lens-$(MYHOST).H Slow-0-trip-lens-$(MYHOST).H \
		Slow-00-trip-lens-$(MYHOST).H
	Add scale=1,1,1 \
	Slow-diff-$*-trip-lens-$(MYHOST).H Slow-0-trip-lens-$(MYHOST).H  \
	Slow-00-trip-lens-$(MYHOST).H \
	datapath=$(TEMPPATH) \
	>$@

Slow-compl-trip-lens-$(MYHOST).H: Slow-top-trip-lens-$(MYHOST).H \
			Slow-sc-.1-trip-lens-$(MYHOST).H
	Merge axis=1 space=n \
	Slow-top-trip-lens-$(MYHOST).H \
	Slow-sc-.1-trip-lens-$(MYHOST).H \
	datapath=$(TEMPPATH) \
	>$@

Slow-0-compl-trip-lens-$(MYHOST).H: Slow-top-trip-lens-$(MYHOST).H \
			Slow-0-trip-lens-$(MYHOST).H
	Merge axis=1 space=n \
	Slow-top-trip-lens-$(MYHOST).H \
	Slow-0-trip-lens-$(MYHOST).H \
	datapath=$(TEMPPATH) \
	>$@


clean	:	
	@-${RM}  .make.dependencies*
	@-${RM}  log*
	@-${RM}  *.H
	@-${RM}  *.Vp
	@-${RM}  *.h
	@-${RM}  *.h@@
	@-${RM}  *.h@@@@
	@-${RM}  *.H@@
	@-${RM}  *.H@@@@
	@-${RM}  *.head
	@-${RM}  $(TEMPDIR)/*


include ${SEPINC}/SEP.bottom
