MAKE_DEPEND=yes
include ${SEPINC}/SEP.top
include ../../Src/thesis.include
RESDIR = ../Figs
RESULTSER= reflectors-cor reflectors-est  compare-small
RESULTSNR=
UF90INCLUDES=${BOBINCS}
UF90LIBS = -lgemops -lsepauxf90 -lsepaux ${MATHLIBF90}
UF90LIBS=${BOBLIBS} ${MATHLIBF90} ${SEPMATHF90}

GRAPH= Graph min2=14000 max2=18000 min1=0 max1=2.0  >/dev/null  wantframe=n title=" " yreverse=y wantaxis=n transp=y xll=2. yll=1.2 xur=11.2 yur=9
GREY = Grey min1=0 max1=2.0 min2=14.000 max2=18.000 >/dev/null title=" " \
       label1="Time(s)" label2="X_Position(km)" label3="velocity(km/s)" wantscalebar=y allpos=y \
			bias=1.5 


###########
##FIGURES##
###########

${RESDIR}/compare-small.v:	real.H ${BINDIR}/calc_var.x
	${BINDIR}/calc_var.x < real.H >gg.H bi=175 ei=275
	Window < gg.H f3=2 f1=175 n1=100 n3=1 |Smooth rect1=3 rect2=3 >a1.H
	Window < gg.H f3=2 f1=175 n1=100H f3=3|Smooth rect1=3 rect2=3 >a2.H
	Add scale=1,-1 a1.H a2.H  >a3.H
	<a1.H Scale dscale=1000 |Grey  bias=3700 clip=1300 label1="Time(s)" label2="X Position(m)" title=" " wantscalebar=y allpos=y barlabel="Velocity (m/s)" >a1.V
	<a2.H Scale dscale=1000 |Grey  bias=3700 clip=1300 label1="Time(s)" label2="X Position(m)" title=" " wantscalebar=y allpos=y barlabel="Velocity (m/s)" >a2.V
	<a3.H Scale dscale=1000 |Grey  label1="Time(s)" label2="X Position(m)" title=" " wantscalebar=y  barlabel="Velocity (m/s)" >a3.V
	Vppen < a1.V a2.V a3.V out=$@ >cp.V gridnum=3,1 vpstyle=n yscale=3
	
			

reflectors.v: layer1.H layer2.H layer3.H layer4.H layer5.H layer6.H 
	Window<layer1.H n1=2>j.H; echo esize=8 n2=1 n1=19 >>j.H; <j.H ${GRAPH} out=1.v
	Window<layer2.H n1=2>j.H; echo esize=8 n2=1 n1=19 >>j.H; <j.H ${GRAPH} out=2.v
	Window<layer3.H n1=2>j.H; echo esize=8 n2=1 n1=15 >>j.H; <j.H ${GRAPH} out=3.v
	Window<layer4.H n1=2>j.H; echo esize=8 n2=1 n1=16 >>j.H; <j.H ${GRAPH} out=4.v
	Window<layer5.H n1=2>j.H; echo esize=8 n2=1 n1=18 >>j.H; <j.H ${GRAPH} out=5.v
	Window<layer6.H n1=2>j.H; echo esize=8 n2=1 n1=7 >>j.H; <j.H ${GRAPH} out=6.v
	vp_Overlay  1.v 2.v 3.v 4.v 5.v 6.v >$@

	
${RESDIR}/reflectors-cor.v: time_int.H reflectors.v
	< time_int.H Add scale=.001 |${GREY} out=a.v o2=13.97 d2=.009375
	vp_Overlay a.v reflectors.v >$@

${RESDIR}/reflectors-est.v: int_angle.H reflectors.v
	${GREY} < int_angle.H out=a.v o2=13.59 d2=.040
	vp_Overlay a.v reflectors.v >$@

############
#INVERSION##
############

interval.H:	pick.H  ${BINDIR}/Vrms2int.x int.P
	${BINDIR}/Vrms2int.x < pick.H >$@ weight=weights.H par=int.P

int_angle.H:	pick.H  ${BINDIR}/Vrms2int.x int.P nfilt.vel0.H
	${BINDIR}/Extend_it.x beg1=185 beg2=185 end1=85 end2=35 <  pick.H >p.H
	${BINDIR}/Extend_it.x beg1=185 beg2=185 end1=85 end2=35 <  weights.H >w.H
	${BINDIR}/Vrms2int.x < p.H >it.H weight=w.H par=int.P ang=1 \
   filt_file=nfilt.vel0.H
	Window < it.H f1=185 f2=185 n1=650 n2=175 >$@

real.H:	pick.H  ${BINDIR}/Vrms2int.x real.P nfilt.vel0.H
	${BINDIR}/Extend_it.x beg1=185 beg2=185 end1=85 end2=35 <  pick.H >p.H
	${BINDIR}/Extend_it.x beg1=185 beg2=185 end1=85 end2=35 <  weights.H >w.H
	${BINDIR}/Vrms2int.x < p.H >it.H weight=w.H par=real.P ang=1 \
   filt_file=nfilt.vel0.H
	Window < it.H f1=185 f2=185 n1=650 n2=175 >$@

${RESDIR}/dix-real.v:	real.H
	Scale dscale=1000 < real.H |Window n3=4 |Grey allpos=y bias=1800. clip=3400  label1="Time(s)" label2="X Position(m)" max1=2.1 title=" " wantframenum=0 |Vppen gridnum=2,2 vpstyle=n >real.V out=$@


########
#ANGLES#
########

#%.H:  %.txt
#	Dd < $*.txt >$*.H ein=0 esize=4

#####
#RMS#
#####

pick.H:	pick.P ${BINDIR}/Pick2.x velan.H limits.H
	${BINDIR}/Pick2.x < velan.H >junk.H 	 weight=weights.H par=pick.P \
    limits=limits.H
	Smooth  < junk.H >$@ par=pick.P


angles.H: layer1.H layer2.H layer3.H layer4.H pick.H ${BINDIR}/temp_refs.x layer5.H layer6.H
	${BINDIR}/temp_refs.x < pick.H file1=layer1.H file2=layer2.H  \
   file3=layer3.H file4=layer4.H file5=layer5.H  nlayers=6 file6=layer6.H >j.H
	<j.H Smooth rect1=30 rect2=10 >$@


########
#LIMITS#
########

limits.H:	min_max.H 
	Window < min_max.H   n1=640 >limits.H dtout=.004

min_max.H:	${BINDIR}/min_max.x vel.H
	${BINDIR}/min_max.x	 <vel.H >$@


######
#DATA#
######

SEP97= /net/kana/wrk1/sep97/bob1/Data

vel.H:	${SEP97}/vel.HH
	Window < ${SEP97}/vel.HH >$@

velan.H:	${SEP97}/velan.HH
	Window < ${SEP97}/velan.HH >$@

time_int.H:	${SEP97}/time_int.HH
	Window < ${SEP97}/time_int.HH >$@

%.zero.H: %.txt ${BINDIR}/create_reflectors.x ${BINDIR}/to_seplib.x
	${BINDIR}/to_seplib.x < $*.txt junk.HH 3
	Dd esize=4 < junk.HH  >c.H
	<c.H Window n1=2 | ${BINDIR}/create_reflectors.x  n1=300  >a.H
	Imag < a.H >r.H
	Real <a.H |Smooth rect1=15 >i.H
	Cmplx r.H i.H >$@

refs.H:	layer1.zero.H layer2.zero.H layer3.zero.H layer4.zero.H layer5.zero.H layer6.zero.H
	Cat layer1.zero.H layer2.zero.H layer3.zero.H layer4.zero.H layer5.zero.H layer6.zero.H axis=2 >$@

nfilt.vel0.H:  nfilt.P refs.H# ${BINDIR}/create_nfilter.x
	Spike n1=640 n2=175 o1=0 o2=13593 d2=40 d1=.004 >aa.H
	${BINDIR}/create_nfilter.x < aa.H par=nfilt.P  >$@ \
    reflectors=refs.H



clean: jclean
	${RM} *log*


include ${SEPINC}/SEP.bottom