RESDIR = ../Figs
include ${SEPINC}/SEP.top
include ../../Src/thesis.include
UF90LIBS=${BOBLIBS} -lsepmathf90
UF90INCLUDES=${BOBINCS}
RESULTSER =fault-model fault.interp fault.cont fault.amp
UPSTEXFLAGS=invras=y
PLOT_SIZE = min1=0. max1=1.592 min2=0. max2=.99
RESULTSER=sea.init sea.pef sea.movie sea.var well.movie

GREY_OPT = >/dev/null title=' ' label1=' ' label2=' ' \
  bias=-0.03514 allpos=y 


bin.H:	../../Data/seabeam.HH
	Window < ../../Data/seabeam.HH >$@

movie.histo.v: sea.movie.H ${BINDIR}/mk_histo.x
	Window n3=1 < sea.movie.H|\
	${BINDIR}/mk_histo.x  sea.movie.H >a.H map=bin.H known=b.H

pef.histo.v: sea.interp.H ${BINDIR}/mk_histo.x
	${BINDIR}/mk_histo.x < sea.interp.H >a.H map=bin.H known=b.H
	Histogram < a.H min=-1. max=1. >aa.H  ninterv=400
	<aa.H Scale dscale=.0000597407 >a1.H
	Histogram < b.H min=-1. max=1.  >bb.H ninterv=400
	<bb.H Scale dscale=.0001128541 >b1.H
	Cat b1.H a1.H axis=2 |Graph legend=1 curvelabel="Known:Estimated" label1="Value" label2="Frequency" title=" " out=$@ legendloc=tr

sea.movie.H:	bin.H	${BINDIR}/Miss_n.x sea.P ${BINDIR}/Pef.x
	${BINDIR}/Pef.x par=sea_interp.P < bin.H >filt.H
	${BINDIR}/Miss_n.x par=sea.P < bin.H >$@ filt=filt.H

sea.interp.H:	bin.H	${BINDIR}/Miss.x sea.P ${BINDIR}/Pef.x
	${BINDIR}/Pef.x par=sea_interp.P < bin.H >filt.H
	${BINDIR}/Miss.x par=sea.P < bin.H >$@ filt=filt.H



${RESDIR}/sea.init.v:	bin.H
	Grey < bin.H out=$@ >/dev/null transp=n yreverse=n label1=longitude label2=latitude title=" "

${RESDIR}/sea.pef.v:	sea.interp.H
	Grey < sea.interp.H out=$@ >/dev/null transp=n yreverse=n label1=longitude label2=latitude title=" "

sea.var.H:	${BINDIR}/calc_var.x sea.movie.H
	${BINDIR}/calc_var.x < sea.movie.H >$@ 

${RESDIR}/sea.movie.v:	sea.movie.H
	Grey  < sea.movie.H crowd=1.  transp=n yreverse=n | Vppen >mov.V out=$@ gridnum=3,3 vpstyle=n

${RESDIR}/sea.var.v:	sea.var.H
	Window < sea.var.H n3=1|Grey  transp=n yreverse=n label1=longitude label2=latitude >a.V title=" " allpos=y
	Window < sea.var.H f3=1 |Grey transp=n yreverse=n label1=longitude label2=latitude >b.V title=" "
	Vppen < a.V b.V >c.V out=$@ vpstyle=n gridnum=2,1
	

PAN_LAB=o1=0 o2=0 d1=.025 d2=.1 label1="Depth(km)" label2="X Position(km)"
COV_LAB=label1="Z Offset(km)" label2="X Offset(km)" d1=.025 d2=.1
	

GRAPH_OPT = >/dev/null title=' ' transp=y wantaxis=n ${PLOT_SIZE}  symbol="*" wantframe=n crowd=.88

${RESDIR}/fault-model.v:	interp.prob.H model.prob.H
	Grey ${PAN_LAB} < model.prob.H >a.V bias=1.04 title=" "
	Grey ${PAN_LAB} < interp.prob.H >b.V  title=" "
	Vppen < a.V b.V >fault.V out=$@ vpstyle=n gridnum=2,1 yscale=2

${RESDIR}/fault.interp.v:	interp.fault.H
	Grey ${PAN_LAB} < interp.fault.H >b.V  title=" " out=$@

${RESDIR}/fault.cont.v: cont.fault.H
	Grey ${PAN_LAB} < cont.fault.H >b.V  title=" " out=$@

${RESDIR}/fault.amp.v: amp.fault.H
	Grey ${PAN_LAB} < amp.fault.H >b.V  title=" " out=$@ wantscalebar=y allpos=y bias=.1 barlabel="Amplitude"


model.%.H qd.slope.%.H:	${BINDIR}/Qdome.x qdome.%.P
	${BINDIR}/Qdome.x >/dev/null mod=m.H par=qdome.$*.P mod2=qd.slope.$*.H
	Smooth rect1=4 rect2=3 < m.H >model.$*.H par=qdome.$*.P

well.%.H overlay.%.H: ${BINDIR}/make_wells.x model.%.H wells.P
	${BINDIR}/make_wells.x <model.$*.H par=wells.P >well.$*.H locations=overlay.$*.H

amp.%.H:	fault.H model.%.H ${BINDIR}/add_faults.x fault.%.P
	< model.$*.H >$@ ${BINDIR}/add_faults.x par=fault.$*.P faults=fault.H

slopes.%.H sloping.%.H:	model.%.H ${BINDIR}/find_slope.x slope.P
	${BINDIR}/find_slope.x <qd.slope.$*.H  par=slope.P >j1.H slopes.H=slopes.$*.H
	< j1.H Smooth rect1=5 rect2=5 >sloping.$*.H

${RESDIR}/well.movie.v:	interp.nofault.H
	Grey < interp.nofault.H crowd=1. |Vppen yscale=3 >w.V out=$@ vpstyle=n gridnum=3,1

interp2.%.H:	${BINDIR}/Helicon2d.x  well.%.H prec-cont.P sloping.%.H amp.%.H
	${BINDIR}/Helicon2d.x < well.$*.H >$@ par=prec-cont.P angle=sloping.$*.H \
		input=input.$*.H  scales=amp.$*.H

interp.%.H:	${BINDIR}/Helicon_n.x  well.%.H prec-cont.P sloping.%.H amp.%.H
	${BINDIR}/Helicon_n.x < well.$*.H >$@ par=prec-cont.P angle=sloping.$*.H \
		input=input.$*.H  scales=amp.$*.H

cont.%.H:	${BINDIR}/nfilter_it.x gpcont.P model.%.H interp.%.H
	${BINDIR}/nfilter_it.x par=gpcont.P < input.$*.H >$@ covar=model.$*.H

${RESDIR}/cont.%.v:	 cont.%.H
	Grey < cont.$*.H out=$@ title=" " labl1=" " label2=" " ${PAN_LAB} >gpc.V

fault.H:	fault.HH ${BINDIR}/create_reflectors.x
	Dd esize=4 < fault.HH >ff.H
	${BINDIR}/create_reflectors.x  < ff.H >gg.H n1=50
	Real < gg.H >r.H
	Imag < gg.H |Smooth rect1=5 >i.H
	Cmplx r.H i.H >$@

clean: jclean

	


	

include ${SEPINC}/SEP.bottom