include ${SEPINC}/SEP.top
BINDIR=./
SRCDIR=../Src
UF90INCLUDES = -I.
UF90LIBS=-lsupersetf90 -lsuperset -lgeef90 -lsepauxf90
HASK=../data/hask.HH
RESDIR=../Figs
RESULTSER=huber mobil compare
UPSTEXFLAGS=invras=y

hask_90.H: ${HASK}
	Window < ${HASK} >hask.H
	<hask.H Create3d keyname1=offset_x >$@


hask_nl.H:  hask_90.H  ${BINDIR}/Demultiple.x demult_nl.P hask_pick.H
	${BINDIR}/Demultiple.x < hask_90.H >$@  multiples=hask_nl_mult.H \
    par=demult_nl.P  vscan.H=nl_v.H cut=hask_pick.H full_both=full.H\
		full_sig=sig.H

hask_lin.H:  hask_90.H  ${BINDIR}/Demultiple.x demult_lin.P hask_pick.H
	${BINDIR}/Demultiple.x < hask_90.H >$@  multiples=hask_lin_mult.H \
    par=demult_lin.P  vscan.H=lin_v.H cut=hask_pick.H

mobil_nl.H:  ../data/mobil_90.HH  ${BINDIR}/Demultiple.x demult_nl.P mobil_pick.H
	${BINDIR}/Demultiple.x < ../data/mobil_90.HH >$@  multiples=mobil_nl_mult.H \
    par=demult_nl.P  vscan.H=mobil_nl_v.H cut=mobil_pick.H full_both=full.H\
		full_sig=sig.H

mobil_lin.H:  ../data/mobil_90.HH  ${BINDIR}/Demultiple.x demult_lin.P mobil_pick.H
	${BINDIR}/Demultiple.x < ../data/mobil_90.HH >$@  multiples=mobil_lin_mult.H \
    par=demult_lin.P  vscan.H=mobil_lin_v.H cut=mobil_pick.H full_both=/dev/null full_sig=/dev/null

${RESDIR}/mobil.v:	../data/mobil_90.HH
	Grey out=$@ >/dev/null title=" " label1="Time(s)" label2="Offset (m)"\
  <../data/mobil_90.HH o2=275 d2=50
  

${RESDIR}/huber.v:	${BINDIR}/hub_graph.x
	 ${BINDIR}/hub_graph.x  o1=-10 cut=2. d1=.1 n1=201 |Graph\
	 out=$@ title=" " label1="Residual" label2="Functional value" >/dev/null\
   legend=1 curvelabel="Huber:L2:L1" legendloc=tr max2=15. dash=0,1,2\
	 plotfat=6 legendtitle="Functional"

${RESDIR}/compare.v:	mobil_nl.H mobil_lin.H
	Envelope < mobil_lin_v.H tc1=5 tc2=2 |Window n3=1| Grey allpos=y out=b.v\
    title="L2" label="Tau (s)"    >/dev/null label2="Velocity (m/s)" max1=3
	Envelope < mobil_nl_v.H tc1=5 tc2=2 | Window n3=1|Grey allpos=y out=a.v \
    title="Huber" label="Tau (s)"    >/dev/null label2="Velocity (m/s)" max1=3
	vp_SideBySideAniso b.v a.v >$@

mobil_pick.H:	hask_pick.H
	cp hask_pick.H $@

hask_pick.H:  pick2.txt ${BINDIR}/create_sem.x ${BINDIR}/to_seplib.x
	${BINDIR}/to_seplib.x < pick2.txt a.H 2
	Dd < a.H esize=4 >b.H
	${BINDIR}/create_sem.x < b.H  o1=0. d1=.004 n1=1024 |Smooth rect1=20>$@

test.view: hask_pick.H hask_nl.H
	Graph < hask_pick.H  d1=.004 transp=y out=a.v yreverse=y >/dev/null \
  max2=3.15 min2=1.2
	Window < l_v.H n3=1 |Envelope |Grey out=b.v >/dev/null  max2=3.15 min2=1.2
	vp_Overlay b.v a.v |tube

${BINDIR}/to_seplib.x:	../Src/to_seplib.x
	cp ../Src/to_seplib.x $@

cut.H:	cut.HH
	Dd < cut.HH  >$@ eout=4

clean:  jclean
	${RM} *.H@*



include ${SEPINC}/SEP.bottom