MAKE_DEPEND=yes
include RepMakefile

include /usr/local/SEP/include/SEP.top
#include /opt/SEP/SEP6.4.6/include/SEP.top
# SEP report
LATOPTS  = style=html
# SEG abstract
#LATOPTS = option=cmfonts class=seg97
#LATOPTS = option=web 'style=jfc,makeidx,html,amsmath'
# Reviewer 
#LATOPTS=class=geophysics
SIGDIR = /mnt/sepdata/1a/sigsbee_multiple
OLEVEL=2
UF90FLAGS =-O3 -C 
MAINLIBS=-lsep2df90 -lsep3df90 -lsep3d -lsepf90 -lsep -L${LIBDIR}  -lsepgeef90
UF77LIBS = ${BEILIB}

WHEREFFTW = /net/koko/linux_local/FFTW/FFTW2/lib/lib
UF90LIBS=${MPILIBS} ${MAINLIBS} -lsepparf90 -lseppar_mpi -lsepgeef90 -lsep2df90 -static-libcxa -lguide -lseppar  -lsuperset -lsupersetf90 -lsuperset -lsepauxf90 -lsepmathf90 -lsep3df90 -lsep3d -lsepf90 -lsep $(WHEREFFTW)fftw.a $(WHEREFFTW)rfftw.a #-Bstatic

#COLOR=y

RESULTSER = ZOtime sig2a01 sig2b01 sig2b02 comp-ang invang4 invang-zero invoff4 invoff-zero invz-filt invz migang1-filt migang1 migang mig-gather migoff1-filt migoff1 migoff mig-stack-filt mig-stack refle-new res-ir-off res_ir res-ir diff-filtered filtered not-filtered resinvz-filt resinvz diag-Sis 
RESULTSNR = 

ER_v =  ZOtime.v sig2a01.v sig2b01.v sig2b02.v \
invz.v migang1-filt.v \
migang1.v migoff1-filt.v \
migoff1.v mig-stack-filt.v mig-stack.v \
refle-new.v \
diff-filtered.v filtered.v not-filtered.v \
resinvz-filt.v resinvz.v diag-Sis.v

ER_pdf =  ZOtime.pdf sig2a01.pdf sig2b01.pdf sig2b02.pdf \
invz.pdf migang1-filt.pdf \
migang1.pdf migoff1-filt.pdf \
migoff1.pdf mig-stack-filt.pdf mig-stack.pdf \
refle-new.pdf \
diff-filtered.pdf filtered.pdf not-filtered.pdf \
resinvz-filt.pdf resinvz.pdf diag-Sis.pdf

# ZOtime
# sig2a01 
# sig2b01 
# sig2b02 
# comp-ang Thesis-Fig/Regularization/Makefile
# invang4 Thesis-Fig/Regularization/Makefile
# invang-zero Thesis-Fig/Regularization/Makefile
# invoff4 Thesis-Fig/Regularization/Makefile
# invoff-zero Thesis-Fig/Regularization/Makefile
# invz-filt Thesis-Fig/Regularization/Makefile
# invz Thesis-Fig/Regularization/Makefile
# migang1-filt Thesis-Fig/Regularization/Makefile
# migang1 Thesis-Fig/Regularization/Makefile
# migang Thesis-Fig/Regularization/Makefile
# mig-gather Thesis-Fig/Regularization/Makefile
# migoff1-filt Thesis-Fig/Regularization/Makefile
# migoff1 Thesis-Fig/Regularization/Makefile
# migoff Thesis-Fig/Regularization/Makefile
# mig-stack-filt Thesis-Fig/Regularization/Makefile
# mig-stack Thesis-Fig/Regularization/Makefile
# refle-new Thesis-Fig/Regularization/Makefile
# res-ir-off Thesis-Fig/Regularization/Makefile
# res_ir Thesis-Fig/Regularization/Makefile
# diff-filtered 
# filtered 
# not-filtered 
# resinvz-filt Thesis-Fig/Regularization/Makefile
# resinvz Thesis-Fig/Regularization/Makefile
# diag-Sis 






# er_burn: clean
er_build: ${ER_v}
er_pdf: ${ER_pdf}

BIN=./bin
PAR=./par
RESDIR=./Fig
R = ${RESDIR}
FIG=./Fig
DATA=./data
SRCDIR=./src

#########################################
################################
#
# KxKh filter
#
kxkh: kxkh_flt.x
	rm -f ${BIN}/kxkh_flt.x;rm -f kxkh_flt.x;make kxkh_flt.x
	mv -f kxkh_flt.x ${BIN}/kxkh_flt.x

kfilt: kfilt.x
	rm -f ${BIN}/kfilt.x;rm -f kfilt.x;make kfilt.x
	mv -f kfilt.x ${BIN}/kfilt.x

kxkh2: kxkh_flt2.x
	rm -f ${BIN}/kxkh_flt2.x;rm -f kxkh_flt2.x;make kxkh_flt2.x
	mv -f kxkh_flt2.x ${BIN}/kxkh_flt2.x

########################################################################################################
# These are the rules for the sep134 multiple !!!!!
#########################################################################################################################
###################Best Parameters##################################
khkx10:
	< mig_Sis_pad_800_strati_sort.H Pad beg2=50 end2=51 beg3=50 end3=50 | ${BIN}/kfilt.x Rot_s_cig=9999filt.H apply=y domain=kxkh npadz=00 npadh=00 npadx=00 khflt1=-0.06 kxflt1=-0.05 khflt2=-0.001 kxflt2=-0.005 khflt3=0.001 kxflt3=-0.005 khflt4=0.06 kxflt4=-0.05 khflt5=0.06 kxflt5=0.05 khflt6=0.001 kxflt6=0.005 khflt7=-0.001 kxflt7=0.005 khflt8=-0.06 kxflt8=0.05 > /dev/null

weight.H:
	${BIN}/mute1.x cig=mig_Sis_pad_800_strati_sort.H vel=vel_sigsbee_pad.H | Smooth rect3=5 rect1=5 > $@

multiple_model1_new.H:
	Window3d n3=200 min3=29975 n2=33 min2=-1200 < sig_flt_new.H > sig_flt_cut_new.H
	Math file1=mig_Sis_pad_800_strati_sort.H file2=sig_flt_cut_new.H exp="file1-file2" >xdiff_new.H
	Window3d min1=10150 max3=34700 < xdiff_new.H | Pad beg1=166 end3=137 > xdiff_cut_pad_new.H
	Math file1=xdiff_cut_pad_new.H file2=sig_flt_cut_new.H exp="file1+file2" > sig_flt_cutbk_new.H
	Math file1=mig_Sis_pad_800_strati_sort.H file2=sig_flt_cutbk_new.H exp="file1-file2"  >multiple_model_new.H
	< multiple_model_new.H ${BIN}/mute.x > $@

mult_est_01_new.H:
	Transp plane=23 <mig_Sis_pad_800_strati_sort.H >junk1_new.H
	Transp plane=23 <multiple_model.H >junk2_new.H
	Transp plane=23 <weight.H | Smooth rect1=5 rect2=5 rect3=1 tridiag=1 >junk3_new.H
	Math file1=junk3_new.H file2=junk2_new.H exp="file1*file2" >junk4_new.H
	time ~claudio/Research/ResMult/Bin/Ultimatch-new-all.x < junk1_new.H \
		primaries_in=  multiples_in=junk4_new.H mu=0.001 \
		primaries_out=prim_est_01_new.H \
		mode=2 nsmode=2 wmode=1 nfilt=5,5 npatch=101,101 \
		op=1 niter=10 2D=1 eps=120. > $@

prim_01_subt_new.H:
	Math file1=junk1_new.H file2=mult_est_01_new.H exp="file1-file2" > $@

mult_est_01_newa.H:
	Transp plane=23 <mig_Sis_pad_800_strati_sort.H >junk1_new.H
	Transp plane=23 <multiple_model1_new.H >junk2_new.H
	Transp plane=23 <weight.H | Smooth rect1=5 rect2=5 rect3=1 tridiag=1 >junk3_new.H
	Math file1=junk3_new.H file2=junk2_new.H exp="file1*file2" >junk4_new.H
	time ~claudio/research.ln/ResMult/Bin/Ultimatch-new-all.x < junk1_new.H \
		primaries_in=  multiples_in=junk4_new.H mu=0.001 \
		primaries_out=prim_est_01_new.H \
		mode=2 nsmode=2 wmode=1 nfilt=5,5 npatch=101,101 \
		op=1 niter=20 2D=1 eps=120. > $@

prim_01_subt_newa.H:
	Math file1=junk1_new.H file2=mult_est_01_newa.H exp="file1-file2" > $@
####################################################################################################
#
#
9999mult_est_01_newa.H:
	Transp plane=23 <mig_Sis_pad_800_sort.H >junk1_new.H
	Transp plane=23 <multiple_model1_new.H >junk2_new.H
	Transp plane=23 <weight.H | Smooth rect1=5 rect2=5 rect3=1 tridiag=1 >junk3_new.H
	Math file1=junk3_new.H file2=junk2_new.H exp="file1*file2" >junk4_new.H
	time ~claudio/research.ln/ResMult/Bin/Ultimatch-new-all.x < junk1_new.H \
		primaries_in=  multiples_in=junk4_new.H mu=0.001 \
		primaries_out=prim_est_01_new.H \
		mode=2 nsmode=2 wmode=1 nfilt=5,5 npatch=101,101 \
		op=1 niter=20 2D=1 eps=120. > $@

9999prim_01_subt_newa.H:
	Math file1=junk1_new.H file2=9999mult_est_01_newa.H exp="file1-file2" > $@
#
#
#######################################################################################################
mig_Sis_pad_800_strati_sort1.T:
	< ../../../Prst_ExpRefl/RTM/mig_Sis_pad_800_strati_sort.H Window3d min3=30500 max3=39600 min2=-675 max2=675\
	| Grey eout=1 gainpanel=a pclip=90 label1="z(ft)" label3="x(ft)" label2="h(ft)" title=" " o2num=-675 d2num=675. \
	| Transp plane=23 > $@

prim_01_subt_newa1.T:
	< ../../../Prst_ExpRefl/RTM/prim_01_subt_newa.H Window3d min2=30500 max2=39600 \
	| Grey eout=1 gainpanel=a pclip=85 label1="z(ft)" label3="x(ft)" label2="h(ft)" title=" " o2num=-1200 d2num=600. > $@




##########################################

regclean_%:
	@-touch $*; Rm $*; touch $*; \rm $*

clean:
	rm -f *.o *.mod
	rm -f ${RESDIR}/*.pdf ${RESDIR}/*.ps ${RESDIR}/*.v

#include /opt/SEP/SEP6.4.6/include/SEP.bottom
include /usr/local/SEP/include/SEP.bottom
