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 movie.histo

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


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


dist%.V: sea.distir.H
	Window <sea.distir.H n3=1 f3=$* |${BINDIR}/mk_histo.x  >a.H map=bin.H known=b.H
	Histogram < a.H ninterv=251 dinterv=.024>aa.H  min=-3. max=3.
	<aa.H Scale  dscale=.00005974429 >bb.H
	<bb.H Window j1=10 |Graph  label1="Value" label2="Frequency" title=" " >a1.V max2=.037   dash="$*" plotcol=5 symbolsz=10  wantaxis=n wantframe=n symbol="$*"
	<bb.H Window j1=10 |Graph  label1="Value" label2="Frequency" title=" " >a2.V max2=.037   dash="$*" plotcol=5 symbolsz=10  wantaxis=n wantframe=n 
	Vppen erase=o <a1.V a2.V >$@


${RESDIR}/distrib.v:	sea.distir.H
	Grey < sea.distir.H gainpanel=e label1=longitude label2=latitude title=" " |Vppen >d.V out=$@ gridnum=2,2 vpstyle=n

frame%.V: sea.movie.H
	Window <sea.movie.H n3=1 f3=$* |${BINDIR}/mk_histo.x  >a.H map=bin.H known=b.H
	Histogram < a.H ninterv=251 dinterv=.008>aa.H  min=-1. max=1.
	<aa.H Scale  dscale=.00005974429 >bb.H
	<bb.H Window j1=10 |Graph  label1="Value" label2="Frequency" title=" " >$@ max2=.017   dash="$*" plotcol=5 symbolsz=12 symbol="$*" 


${RESDIR}/movie.distir.v: sea.distir.H ${BINDIR}/mk_histo.x dist0.V dist1.V dist2.V dist3.V # frame6.V frame7.V
	Window <sea.movie.H n3=1 f3=0 |${BINDIR}/mk_histo.x >a.H map=bin.H known=b.H
	<b.H Histogram ninterv=251 dinterv=.024 min=-3. max=3. | Scale dscale=.0001128541 |Graph  >a1.V label1="Value" label2="Frequency" title=" " max2=.037 
	Vppen < a1.V dist0.V dist1.V dist2.V dist3.V   >sto.V out=$@ erase=n  vpstyle=n 

movie.histo.V: sea.movie.H ${BINDIR}/mk_histo.x frame0.V frame1.V frame2.V frame3.V frame4.V frame5.V # frame6.V frame7.V
	Window <sea.movie.H n3=1 f3=0 |${BINDIR}/mk_histo.x >a.H map=bin.H known=b.H
	<b.H Histogram ninterv=251 dinterv=.008 min=-1. max=1. | Scale dscale=.0001128541 |Graph  >a1.V label1="Value" label2="Frequency" title=" " max2=.017 wantaxis=n
	Vppen < a1.V frame1.V frame2.V frame3.V frame4.V frame5.V >$@ erase=n  vpstyle=n 

${RESDIR}/non-stat.v:	${BINDIR}/Miss_d.x sea.P bin.H 
	${BINDIR}/Pef.x par=sea_interp.P < bin.H >filt.H
	${BINDIR}/Miss_d.x par=sea.P < bin.H  do_ns=y filt=filt.H |Grey out=$@ >dd.V title=" "

${RESDIR}/check.v:	${BINDIR}/Miss_d.x check.P bin.H 
	${BINDIR}/Pef.x par=sea_interp.P < bin.H >filt.H
	${BINDIR}/Miss_d.x par=sea.P < bin.H  do_ns=n filt=filt.H par=check.P myres=myres.H >ff.H 
	echo n1=160 n2=160 n3=2 >> myres.H	
	Cat ff.H myres.H| Grey out=$@ >dd.V gainpanel=e wantscalebar=y
	

${RESDIR}/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=251 dinterv=.008
	<aa.H Scale dscale=.0000597407 >a1.H
	Histogram < b.H min=-1. max=1.  >bb.H ninterv=251 dinterv=.008
	<bb.H Scale dscale=.0001128541 >b1.H
	< b1.H Graph label1="Value" label2="Probability Density" title=" "> a1.V max2=.05 
	< a1.H Window j1=2 | Graph label1="Value" label2="Frequency" title=" " > a2.V max2=.05   symbol="*" plotcol=5 symbolsz=12
	Vppen < a1.V a2.V >ff.V erase=n  vpstyle=n out=$@

sea.distir.H:	bin.H	${BINDIR}/Miss_d.x sea.P ${BINDIR}/Pef.x
	${BINDIR}/Pef.x par=sea_interp.P < bin.H >filt.H
	${BINDIR}/Miss_d.x par=sea.P < bin.H >a1.H filt=filt.H scale=.02
	${BINDIR}/Miss_d.x par=sea.P < bin.H >a2.H filt=filt.H scale=.04
	${BINDIR}/Miss_d.x par=sea.P < bin.H >a3.H filt=filt.H scale=.08
	${BINDIR}/Miss_d.x par=sea.P < bin.H >a4.H filt=filt.H scale=.16
	Cat a1.H a2.H a3.H a4.H >$@

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 movie.histo.V
	Grey  < sea.movie.H crowd=1.  transp=n yreverse=n n3=8  >a1.V  
	Vppen < a1.V movie.histo.V  out=$@ gridnum=3,3 vpstyle=n >dd.V

${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

	
morgan-pef.H:
	${BINDIR}/Morgan_n.x < m-bin.H >$@ variance=var.H par=morgan.P filt=lap.H scale=0.

morgan.movie.H:	m-bin.H ${BINDIR}/Morgan_n.x morgan.P
	${BINDIR}/Morgan_n.x < m-bin.H >$@ variance=var.H par=morgan.P filt=lap.H


	

include ${SEPINC}/SEP.bottom