include ${SEPINC}/SEP.top
BINDIR=./
UF90INCLUDES=-I/net/kana/bob/include/${MTYPE} 
UF90FLAGS=
UF90LIBDIR=/net/kana/bob/lib/${MTYPE} 
UF90INCLUDES=

UF90LIBS=  -lsupersetf90 -lsuperset -lsep2df90



vel_new.H: vel_cor.HH ${BINDIR}/modify_vel.x modify.P ref.vel_cor.HH
	${BINDIR}/modify_vel.x < vel_cor.HH >$@ par=modify.P ref_in=ref.vel_cor.HH ref_out=ref_new.vel_cor.H


moduli_new.H:	vel_new.H ${BINDIR}/mk_mod.x 
	Interp <vel_new.H  d1out=.02 d2out=.02 |Extend_it.x beg2=325 beg1=10 end1=10 end2=325 >v.H
	${BINDIR}/mk_mod.x  < v.H >$@ ratio=.02 refs=ref_new.vel_cor.H

test_new.H:	moduli_new.H wavelet.H iso.P
	Iso2d <	wavelet.H moduli=moduli_new.H par=iso.P >$@



#CREATE A CMP GATHER 

synth3.3d.H:	tot.H sort.P
	echo o2=0. d3=.02 o3=0.0  >> tot.H
	Bandpass < tot.H flo=3 >ff.H
	Window < ff.H >wind.H j1=4
	echo o2=0. d3=.04 o3=0.0  >> wind.H
	Create3d keyname1=offset_x keyname2=s_x < wind.H >tot.h.H
	Headermath < tot.h.H >tot.m.H par=math.P
	Sort3d < tot.m.H > $@ par=sort.P  datapath=/net/kana/data/2d_synthetic/anticline/

synth.3d.H:	tot.H sort.P
	echo o2=0. d3=.04 o3=0. d2  >> tot.H
	Window < tot.H >wind.H j1=4
	Create3d keyname1=offset_x keyname2=s_x < wind.H >tot.h.H
	Headermath < tot.h.H >tot.m.H par=math.P
	Sort3d < tot.m.H > $@ par=sort.P  datapath=/net/kana/bob/data/ANTI_


synth3.H:	synth3.3d.H wind.P
	Infill3d < synth3.3d.H > big.H
	Window < big.H  par=wind.P |Window  min2=2 max2=18  >wind2.H
	Transp plane=23 < wind2.H  |Mute tmute=.2 vmute=1.9 |Transp plane=23>$@

synth2.H:	synth2.3d.H wind.P
	Stack3d axes=1 < synth2.3d.H > big.H
	Window < big.H  par=wind.P |Window  min2=2 max2=18  >wind2.H
	Transp plane=23 < wind2.H  |Mute tmute=.2 vmute=1.9 |Transp plane=23>$@

synth.H:	synth.3d.H wind.P
#	Stack3d axes=1 < synth.3d.H > big.H
#	Window < big.H  par=wind.P |Window min2=2 max2=18 >wind2.H
	Transp plane=23 < wind2.H  |Mute tmute=.2 vmute=1.9 |Transp plane=23>$@

synth.freq.H:	synth2.H freq.P
	Window < synth2.H f3=1 |Reverse which=4 >aa.H
	Cat  aa.H  synth2.H axis=3 |Window f1=3 >bb.H
	echo o1=0. >> bb.H
	Transf < bb.H >synth.freq.H par=freq.P
	echo o2=0 d2=1 n2=1 o3=-2.528 d3=.032 n3=159 n4=128 d4=.262605  o4=0 >>$@
	echo z_data_0=0.  >> $@

mig.%.H:	synth.freq.H mig.P %.mig.H
	Window <  synth.freq.H >in.H squeeze=n  n1=500 n3=159
	echo d3=.032 >> in.H
	/net/kana/biondo/bin/${MTYPE}/Pre2D_salt_ssf.x par=mig.P>log.mig2 Data=in.H Image=image.H Slow=$*.mig.H
	<image.H Window | Transp plane=23 |Transp plane=12 |Transp plane=23 >$@

mig2.%.H:	synth.freq.H mig.P %.mig.H
	Window <  synth.freq.H >in2.H squeeze=n  n1=500 n3=159
	echo d3=.032 >> in2.H
#	/net/kana/biondo/bin/${MTYPE}/Pre2D_salt_ssf.x par=mig.P>log.mig Data=in2.H Image=image2.H Slow=$*.mig.H
	../../Bin/Pre2D_salt_ssf.x par=mig.P>log.mig Data=in2.H Image=image2.H Slow=$*.mig.H
	<image2.H Window | Transp plane=23 |Transp plane=12 |Transp plane=23 >$@

#synth2.3d.H:	synth.H sort2.P math2.P
#	Create3d keyname1=cmp_x keyname2=offset_x < synth.H >ss.H
#	Headermath < ss.H >tt.H par=math2.P
#	Sort3d < tt.H >$@ par=sort2.P
	
kirch.H: synth2.3d.H vel_cor.z.H depth.P
	cp synth2.3d.H in.H
#	Headermath < synth2.3d.H >nn.H eqn1=offset_x*1. key1=offset_x
#	Headermath < nn.H >mn.H key1=s_x eqn1=cmp_x-offset_x/2. key2=g_x eqn2=cmp_x+offset_x/2.
#	Window < nn.H f1=0 n3=1 >in.H
	Window < vel_cor.z.HH  j1=4 j2=8  j1=4 j2=8  >vvv.H
	echo n3=1  o1=0 o1=0.00 >> in.H 
	/net/kana/bob/bin/SGI64/Kirchoff_2d.x travel=times3.synth.H velocity=vvv.H \
   < in.H > $@ par=depth.P

	
#mig2.%.H:	synth.freq.H mig.P %.mig.H
## n3=54 f3=10
##	echo d3=.032 n2=59 n3=128 n4=1 >> inb.H 
##	Reverse < inb.H which=2 >inb1.H
##	Spike n1=1000    n2=10   n3=128 >ina.H mag=0.
##	echo n1=500 esize=8 >> ina.H
##	Cat inb1.H ina.H  axis=2|Reverse which=2 >inb2.H
##	Spike n1=1000    n2=63   n3=128 >in2.H mag=0.
##	echo n1=500 esize=8 >> in2.H
##	Cat  inb2.H in2.H >in3.H axis=2
##	echo n2=1 n3=127 n4=128 >> in3.H
#	Window <  synth.freq.H >inb.H squeeze=n  n1=500 n3=110
#	Scale dscale=1.0 < $*.mig.H >vel.H
#	/net/kana/biondo/bin/${MTYPE}/Pre2D_salt_ssf.x par=mig.P>log.mig2 Data=inb.H Image=image2.H Slow=vel.H
#	<image2.H Window | Transp plane=23 |Transp plane=12 |Transp plane=23 >$@
#


coher.%.H:	mig.%.H ${BINDIR}/semblance.x coher.P ${BINDIR}/rmute.x
	<mig.$*.H ${BINDIR}/rmute.x >smute.H par=coher.P
	< smute.H ${BINDIR}/semblance.x  >$@ par=coher.P

vel0.mig.H:	vel0.z.H
	< vel0.z.H Window  j1=2 | Transp plane=12 |Transp plane=23 > $@

vel_cor.mig.H:	vel0.z.H
	< vel_cor.z.H Window  j1=2 |Extend_it.x  beg2=100 end2=100 end1=10 | Transp plane=12 |Transp plane=23 > $@

	




moduli2.H:	vel_cor.HH ${BINDIR}/mk_mod.x  ref.vel_cor.HH
	Spike n1=167 n2=626 o1=0 o2=0 d1=.032 d2=.032 mag=2. >vv.H
	< vv.H  Window  |Extend_it.x beg2=10 end2=20 end=10  >v.H
	${BINDIR}/mk_mod.x  < v.H >$@ ratio=.4 refs=ref.vel_cor.HH
	
	
moduli3.H:	vel_cor.HH ${BINDIR}/mk_mod.x  ref.vel_cor.HH
	< vel_cor.HH Window |  Transp >a.H
	< a.H Interp d1out=.016 |Transp >b.H
	<b.H  Extend_it.x beg2=220 end2=40 end1=20 |Window max2=14  >v.H
	${BINDIR}/mk_mod.x  < v.H >$@ ratio=.02 refs=ref.vel_cor.HH
	

wavelet.H:
	 Wavelet n1=40 d1=0.002 domain=time wavelet=ricker1 tdelay=0.025 fund=35. fhigh=40.  > $@

test2.H:	moduli.H wavelet.H iso.P
	Iso2d <	wavelet.H moduli=moduli.H par=iso.P >$@

part1.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651>mod1.H 
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=52 moduli=mod1.H verbose=1

part2.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651 f2=104>mod2.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=52 moduli=mod2.H

part3.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651 f2=208>mod3.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=12 moduli=mod3.H

part4.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651 f2=232>mod4.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=12 moduli=mod4.H

part5.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651 f2=256>mod5.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=12 moduli=mod5.H

part6.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651 f2=280 >mod6.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=12 moduli=mod6.H verbose=1


part7.H:	part7a.H part7b.H part7c.H part7d.H part7e.H part7f.H part7g.H
	Cat part7a.H part7b.H part7c.H part7d.H part7e.H part7f.H part7g.H >$@

section1.H:	part1.H part2.H part3.H part4.H part5.H part6.H part7.H part8.H part9.H
	Cat part1.H part2.H part3.H part4.H part5.H part6.H part7.H part8.H part9.H >$@

tot.H: section1.H section2.H section3.H
	Cat section1.H section2.H section3.H >$@

part7a.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651 f2=304 >mod7a.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=2 moduli=mod7a.H verbose=1


part7b.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651 f2=308 >mod7b.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=2 moduli=mod7b.H verbose=1


part7c.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651 f2=312 >mod7c.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=2 moduli=mod7c.H verbose=1


part7d.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651 f2=316 >mod7d.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=2 moduli=mod7d.H verbose=1


part7e.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651 f2=320 >mod7e.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=2 moduli=mod7e.H verbose=1


part7f.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651 f2=324 >mod7f.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=2 moduli=mod7f.H verbose=1


part7g.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651 f2=328 >mod7g.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=2 moduli=mod7g.H verbose=1


#part7.H: moduli_new.H iso.P
#	Window < moduli_new.H n2=651 f2=304 >mod7.H
#	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=12 moduli=mod6.H verbose=1
#

part8.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651 f2=332 >mod8.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=30 moduli=mod8.H 


#NEEDS TO BE STARTED ON KANA
part9.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651 f2=392 >mod9.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=30 moduli=mod9.H

part10.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651 f2=452 >mod10.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=20 moduli=mod10.H

part11.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651 f2=492 >mod11.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=20 moduli=mod11.H

part12.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651 f2=532 >mod12.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=20 moduli=mod12.H

part13.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651 f2=572 >mod13.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=20 moduli=mod13.H

part14.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651 f2=612 >mod14.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=20 moduli=mod14.H

part15.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651 f2=652 >mod15.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=20 moduli=mod15.H

part16.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651 f2=692 >mod16.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=20 moduli=mod16.H

part17.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651 f2=732 >mod17.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=20 moduli=mod17.H

part18.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651 f2=772 >mod18.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=20 moduli=mod18.H


part19.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651 f2=812 >mod19.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=10 moduli=mod19.H


part20.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651 f2=832 >mod20.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=10 moduli=mod20.H

part21.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651 f2=852 >mod21.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=10 moduli=mod21.H

part22.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651 f2=872 >mod22.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=10 moduli=mod22.H

part23.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651 f2=892 >mod23.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=10 moduli=mod23.H

part24.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651 f2=912 >mod24.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=10 moduli=mod24.H


part25.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651 f2=932 >mod25.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=20 moduli=mod25.H

part26.H: moduli_new.H iso.P
	Window < moduli_new.H n2=651 f2=972 >mod26.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=20 moduli=mod26.H

part27.H: moduli_new.H iso.P
	Window < moduli_new.H f2=1012 >mod27.H
	/net/kana/bob/bin/${MTYPE}/Iso2d < wavelet.H >$@  par=iso.P src_1=326 nsrc=20 moduli=mod27.H

include ${SEPINC}/SEP.bottom