RESDIR = ../Figs
include ${SEPINC}/SEP.top
include ../../Src/thesis.include
UF90LIBS=${BOBLIBS}
UF90INCLUDES=${BOBINCS}
RESULTSER =qdome-three qdome-geophysics covar model-covar
RESULTSER = well-logs qdome-lap qdome-reg-cont qdome-prec-cont gp-mis model-covar covar-change geostat-mis gp-continuous patch qdome-refs
UPSTEXFLAGS=invras=y
PLOT_SIZE = min1=0. max1=1.592 min2=0. max2=.99

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



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




default: view
#
#IN HERE
#
NOT=title=" " label1=" " label2=" " wantaxis=n

patch.V:	covar1.P covar2.P covar3.P covar4.P ${BINDIR}/patch_bound.x
	Add add=-1.04 < panel.H >nodc.H
	${BINDIR}/patch_bound.x < nodc.H >a1.V par=covar1.P
	${BINDIR}/patch_bound.x < nodc.H >a2.V par=covar2.P
	${BINDIR}/patch_bound.x < nodc.H >a3.V par=covar3.P
	${BINDIR}/patch_bound.x < nodc.H >a4.V par=covar4.P
	Grey < nodc.H ${NOT} xll=.5 yll=.5 xur=13.12067 yur=9.727999 >b.V crowd=1.
	Vppen erase=o vpstyle=n < b.V a1.V a2.V a3.V a4.V >$@

${RESDIR}/covar-change.v:	panel.H ${BINDIR}/covar.x covar1.P covar2.P covar3.P covar4.P patch.V
	Add add=-1.04 < panel.H >nodc.H
	echo d1=1. d2=1 o1=0 o2=0 >> nodc.H
	${BINDIR}/covar.x <  nodc.H par=covar1.P |Grey title=" " label1=" " label2=" " wantaxis=n >covar1.V crowd=.90 title="1" titlesz=16 crowd2=.8
	${BINDIR}/covar.x <  nodc.H  par=covar2.P |Grey title=" " label1=" " label2=" " wantaxis=n >covar2.V crowd=.90 title="2" titlesz=16 crowd2=.8
	${BINDIR}/covar.x <  nodc.H  par=covar3.P |Grey title=" " label1=" " label2=" " wantaxis=n >covar3.V crowd=.90 title="3" titlesz=16 crowd2=.8
	${BINDIR}/covar.x <  nodc.H  par=covar4.P |Grey title=" " label1=" " label2=" " wantaxis=n >covar4.V crowd=.90 title="4" titlesz=16 crowd2=.8
	Vppen < covar1.V covar2.V covar3.V covar4.V  >covar.V  gridnum=2,2 vpstyle=n sxscale=1 yscale=1
	Vppen  < patch.V  covar.V >change.V out=$@ gridnum=2,1 vpstyle=n yscale=2

${RESDIR}/gp-continuous.v:	${BINDIR}/nfilter_it.x panel.H qdome-lap.H gpcont.P
	${BINDIR}/nfilter_it.x par=gpcont.P < input1.H >gpc.H covar=panel.H
	Grey < gpc.H out=$@ title=" " labl1=" " label2=" " ${PAN_LAB} >gpc.V


fault.H:	fault.pck ${BINDIR}/to_seplib.x ${BINDIR}/create_reflectors.x
	${BINDIR}/to_seplib.x <fault.pck junk.HH 2
	Dd esize=4 < junk.HH  >c.H
	${BINDIR}/create_reflectors.x < c.H >d.H n1out=50
	Real < d.H |Smooth rect1=5 >r.H
	Imag < d.H | Smooth rect1=5 >i.H
	Cmplx i.H r.H >$@

amps2.H:	fault.H panel2.H ${BINDIR}/add_faults.x fault.P
	< panel2.H >$@ ${BINDIR}/add_faults.x par=fault.P faults=fault.H



	

${RESDIR}/model-covar.v:	panel.H ${BINDIR}/covar.x model-covar.P
	Add add=-1.04 < panel.H >nodc.H
	echo d1=1. d2=1 o1=0 o2=0 >> nodc.H
	${BINDIR}/covar.x <  nodc.H >covar.H par=model-covar.P
	Grey < covar.H out=$@ >/dev/null title=" " label1=" " label2=" " \
    ${COV_LAB} o1=-1.625 o2=-6.5

${RESDIR}/qdome-prec-cont.v:	panel.H slopes.H sloping.H qdome-prec-cont.H
	Grey < qdome-prec-cont.H  ${GREY_OPT}   allpos=n out=$@  crowd=.75 bias=0.0 ${PAN_LAB}

${RESDIR}/qdome-reg-cont.v:	panel.H slopes.H sloping.H qdome-reg-cont.H
	Grey < qdome-reg-cont.H  ${GREY_OPT}   allpos=n out=$@  ${PAN_LAB} bias=0.0

covar.H: model-covar.P ${BINDIR}/covar.x panel.H
	Add add=-1.04 < panel.H >nodc.H
	echo d1=1. d2=1 o1=0 o2=0 >> nodc.H
	${BINDIR}/covar.x <  nodc.H >covar.H par=model-covar.P

krig.H:	panel.H covar.H ${BINDIR}/krig_it.x krig.P
	${BINDIR}/krig_it.x <  input1.H >$@ par=krig.P covar=covar.H

#${RESDIR}/geostat-mis.v:	panel.H qdome-lap.H model-covar.P krig.P ${BINDIR}/covar.x ${BINDIR}/krig_it.x  covar.H
#	Add add=-1.04 < panel.H >nodc.H
#	echo d1=1. d2=1 o1=0 o2=0 >> nodc.H
#	${BINDIR}/krig_it.x <  input1.H >krig.H par=krig.P covar=covar.H covar_est=est.H
#	Smooth < krig.H |Grey   title=" " label1=" " label2=" " ${PAN_LAB} >pt1.V 
#	 < est.H  | Grey   title=" " label1=" " label2=" " ${COV_LAB} o1=-1.625 d1=-6.5 >pt2.V 
#	Vppen gridnum=2,1 < pt2.V pt1.V >geostat-mis.V out=$@ vpstyle=n yscale=2

${RESDIR}/geostat-mis.v:	panel.H  model-covar.P krig.P ${BINDIR}/covar.x ${BINDIR}/krig2.x  covar.H wells1g.H
	Add add=-1.04 < panel.H >nodc.H
	echo d1=1. d2=1 o1=0 o2=0 >> nodc.H
	${BINDIR}/krig2.x <  wells1g.H >krig.H par=krig.P covar=covar.H covar_est=est.H
	Smooth < krig.H rect1=4 rect2=4|Grey   title=" " label1=" " label2=" " >pt1.V   ${PAN_LAB}
	Smooth < est.H  rect1=4 rect2=4|Grey   title=" " label1=" " label2=" " >pt2.V  ${COV_LAB} o1=-1.625 d1=-6.5
	Vppen gridnum=2,1 < pt2.V pt1.V >geostat-mis.V out=$@ vpstyle=n yscale=2



${RESDIR}/patch.v:	 gp-patch.H geostat-patch.H
	Grey < gp-patch.H   title=" " label1=" " label2=" "  >gp.V  ${PAN_LAB}
	Smooth rect1=4 rect2=3 < geostat-patch.H|Grey   title=" " label1=" " label2=" "  >geostat.V  ${PAN_LAB}
	Vppen < gp.V geostat.V gridnum=2,1  yscale=2 vpstyle=n >patch2.V out=$@





gp-patch.H:	panel.H sloping.H qdome-lap.H  ${BINDIR}/Pef.x pef.P ${BINDIR}/Miss.x ${BINDIR}/Helicon.x covar1.P covar2.P covar3.P covar4.P ${BINDIR}/patch2.x ${BINDIR}/patch1.x
	cp panel.H dumb.H ; echo esize=4 >> dumb.H
	cp input1.H dumb2.H ; echo esize=4 >> dumb2.H
	Window < dumb.H par=covar1.P | ${BINDIR}/Pef.x >pef.H  par=pef.P
	Window < dumb2.H par=covar1.P | ${BINDIR}/Miss.x filt=pef.H >t1.H prec=1 niter=360
	Window < dumb.H par=covar2.P | ${BINDIR}/Pef.x >pef.H  par=pef.P
	Window < dumb2.H par=covar2.P | ${BINDIR}/Miss.x filt=pef.H >t2.H prec=1 niter=360
	Window < dumb.H par=covar3.P | ${BINDIR}/Pef.x >pef.H  par=pef.P
	Window < dumb2.H par=covar3.P | ${BINDIR}/Miss.x filt=pef.H >t3.H prec=1 niter=360
	Window < dumb.H par=covar4.P | ${BINDIR}/Pef.x >pef.H  par=pef.P
	Window < dumb2.H par=covar4.P | ${BINDIR}/Miss.x filt=pef.H >t4.H prec=1 niter=360
	Window < dumb2.H >a1.H par=covar1.P
	Window < dumb2.H >a2.H par=covar2.P
	Window < dumb2.H >a3.H par=covar3.P
	Window < dumb2.H >a4.H par=covar4.P
	${BINDIR}/patch2.x  < t2.H file2=t3.H >tmp.H
	${BINDIR}/patch1.x  < t1.H file2=tmp.H >tmp2.H
	${BINDIR}/patch1.x  < tmp2.H file2=t4.H >gp-patch.H

krig-patch%.H:	covar%.P ${BINDIR}/covar.x ${BINDIR}/krig2.x wells1g.H panel.H krig.P
	Add add=-1.04 < panel.H >nodc.H
	echo o1=0 o2=0 d1=1 d2=1 >> nodc.H
	${BINDIR}/covar.x <  nodc.H >covar$*.H par=covar$*.P
	<wells1g.H  ${BINDIR}/krig2.x >$@ par=krig.P covar=covar$*.H  covar_est=cv$*.H par=covar$*.P
	

	
geostat-patch.H: krig-patch1.H krig-patch2.H krig-patch3.H krig-patch4.H 	${BINDIR}/patch2.x ${BINDIR}/patch1.x
	${BINDIR}/patch2.x  < krig-patch2.H file2=krig-patch3.H >tmp.H
	${BINDIR}/patch1.x  < krig-patch1.H file2=tmp.H >tmp2.H
	${BINDIR}/patch1.x  < tmp2.H file2=krig-patch4.H >$@
	

	

${RESDIR}/gp-mis.v:	panel.H sloping.H qdome-lap.H  ${BINDIR}/Pef.x pef.P ${BINDIR}/Miss.x ${BINDIR}/Helicon.x
	cp panel.H dumb.H ; echo esize=4 >> dumb.H
	${BINDIR}/Pef.x < dumb.H >pef.H  par=pef.P
	cp input1.H dumb.H ; echo esize=4 >> dumb.H
	${BINDIR}/Miss.x   < dumb.H filt=pef.H >temp.H prec=0 niter=220
	Spike n1=200 n2=200 k1=100 k2=100 >sp.H
	${BINDIR}/Helicon.x <sp.H  div=1 adj=0 filt=pef.H |${BINDIR}/Helicon.x  filt=pef.H div=1 adj=1 >imp1.H
	Grey < temp.H   title=" " label1=" " label2=" "  >pt1.V ${PAN_LAB}
	Grey < imp1.H    title=" " label1=" " label2=" "  >pt2.V ${COV_LAB} o2=-10 o1=-2.5
	Vppen gridnum=2,1 < pt2.V pt1.V >gp-mis.V out=$@ vpstyle=n yscale=2



${RESDIR}/qdome-lap.v:	panel.H slopes.H sloping.H qdome-lap.H
	Grey < qdome-lap.H  ${GREY_OPT}   allpos=n out=$@  crowd=.75 bias=0.0 ${PAN_LAB}

${RESDIR}/well-logs.v:	panel.H slopes.H sloping.H qdome-lap.H
	Grey < panel.H  ${GREY_OPT} bias=1.05   allpos=n out=junk1.v  ${PAN_LAB}
	Grey < input1.H  ${GREY_OPT} out=junk6.v  ${PAN_LAB}
	vp_SideBySideIso junk1.v junk6.v >$@

${RESDIR}/qdome-refs.v:	panel.H slopes.H sloping.H qdome-lap.H
	Grey < panel.H  ${GREY_OPT} bias=1.05   allpos=n  crowd=.75 >a.V\
    wanttitle=n ${PAN_LAB}
	Real < slopes.H |Smooth rect1=5 >r.H
	Imag < slopes.H >i.H
	Cmplx r.H i.H |Graph yreverse=y min1=0 max2=.99 min2=0 max1=1.592 crowd=.75\
   >b.V symbol="*"  symbolsz=16  plotcol=0 transp=1\
   wantaxis=n wanttitle=n
	Cmplx r.H i.H |Graph yreverse=y min1=0 max2=.99 min2=0 max1=1.592 crowd=.75\
   >c.V symbol="*"  symbolsz=8  plotcol=7 transp=1\
   wantaxis=n wanttitle=n 
	Grey < sloping.H  ${GREY_OPT}  allpos=n  crowd=.75 >d.V ${PAN_LAB} wantscalebar=y barlabel="Slope(degrees)"
	Vppen vpstyle=n erase=o < a.V b.V c.V >refs.V 
	Vppen vpstyle=n  < refs.V d.V > qdome-refs.V out=$@ gridnum=2,1 yscale=2


qdome-lap.H:	${BINDIR}/Helicon2d.x  wells1.H lap.P sloping.H
	${BINDIR}/Helicon2d.x < wells1.H >$@ par=lap.P angle=sloping.H \
		input=input1.H

qdome-reg-cont.H:	${BINDIR}/Helicon2d.x  wells1.H reg-cont.P sloping.H
	${BINDIR}/Helicon2d.x < wells1.H >$@ par=reg-cont.P angle=sloping.H \
		input=input1.H

qdome-prec-cont.H:	${BINDIR}/Helicon2d.x  wells1.H prec-cont.P sloping.H
	${BINDIR}/Helicon2d.x < wells1.H >$@ par=prec-cont.P angle=sloping.H \
		input=input1.H

qdome2-prec-cont.H:	${BINDIR}/Helicon2d.x  well21.H prec-cont.P sloping2.H amps2.H
	${BINDIR}/Helicon2d.x < well21.H >$@ par=prec-cont.P angle=sloping2.H \
		input=input2.H  scales=amps2.H


${RESDIR}/qdome-geophysics.v:	panel.H slopes.H sloping.H qdome-inverse1.H
	Grey < panel.H  ${GREY_OPT} bias=1.05   allpos=n out=junk1.v  crowd=.75
	Window < slopes.H  j1=2| Graph ${GRAPH_OPT}  out=junk2.v d1=.008 o1=0 o2=0 d2=.01 yreverse=y  plotcol=2 symbolsz=20 crowd=.75
	Window < slopes.H j1=2 |Graph  ${GRAPH_OPT}  out=junk3.v d1=.008 o1=0 o2=0 d2=.01 yreverse=y plotcol=7 symbolsz=7 crowd=.75
	vp_Overlay junk1.v junk2.v junk3.v > junk4.v
	Grey < sloping.H  ${GREY_OPT} bias=-16.69 wantscalebar=y barlabel="dip (degrees)"\
    out=junk5.v crowd=.75 labelsz=16
	Grey < qdome-input1.H  ${GREY_OPT} out=junk6.v crowd=.75 labelsz=16
	vp_SideBySideAniso junk4.v junk5.v  junk6.v>$@
	-rm junk1.v junk2.v junk3.v junk4.v junk5.v junk6.v

panel2.H: qdome2-mod.H
	Smooth rect1=4 rect2=3 <qdome2-mod.H > panel2.H

panel.H: qdome-mod.H
	Smooth rect1=4 rect2=3 <qdome-mod.H > panel.H

${RESDIR}/qdome-three.v:	qdome-lap.H qdome-reg.H qdome-prec.H
	Grey < qdome-lap.H out=junk1.v ${GREY_OPT} title="Symetric, 35 itterations "
	Grey < qdome-reg.H out=junk2.v ${GREY_OPT} title="Steering reg, 35 itterations "
	Grey < qdome-prec.H out=junk3.v ${GREY_OPT} title="Steering prec, 2 itterations "
	vp_SideBySideAniso junk1.v junk2.v junk3.v >$@


qdome-reg.H: ${BINDIR}/Helicon2d.x  wells1.H reg.P  sloping.H
	${BINDIR}/Helicon2d.x < wells1.H >$@ par=reg.P angle=sloping.H
	
qdome-prec.H: ${BINDIR}/Helicon2d.x  wells1.H prec.P  sloping.H
	${BINDIR}/Helicon2d.x < wells1.H >$@ par=prec.P angle=sloping.H
	
qdome2-prec.H: ${BINDIR}/Helicon2d.x  well21.H prec.P  sloping2.H
	${BINDIR}/Helicon2d.x < well21.H >$@ par=prec.P angle=sloping2.H
	


qdome2-slope.H qdome2-mod.H qdome2-refl.H qdome2-mod2.H:	${BINDIR}/Qdome.x qdome2.P
	${BINDIR}/Qdome.x >/dev/null mod=qdome2-mod.H par=qdome2.P slope=qdome2-slope.H \
		testing=qdome2-refl.H mod2=qdome2-mod2.H

qdome-slope.H qdome-mod.H qdome-refl.H qdome-mod2.H:	${BINDIR}/Qdome.x qdome.P
	${BINDIR}/Qdome.x >/dev/null mod=qdome-mod.H par=qdome.P slope=qdome-slope.H \
		testing=qdome-refl.H mod2=qdome-mod2.H

slopes2.H sloping2.H:	qdome2-mod2.H ${BINDIR}/find_slope.x slope.P
	${BINDIR}/find_slope.x <qdome2-mod2.H  par=slope.P >j1.H
#	${BINDIR}/find_slopes.x <qdome-mod2.H  par=slope.P  >j1.H
	< j1.H Smooth rect1=15 rect2=15 >sloping2.H

slopes.H sloping.H:	qdome-mod2.H ${BINDIR}/find_slope.x slope.P
	${BINDIR}/find_slope.x <qdome-mod2.H  par=slope.P >j1.H
#	${BINDIR}/find_slopes.x <qdome-mod2.H  par=slope.P  >j1.H
	< j1.H Smooth rect1=15 rect2=15 >sloping.H

qdome-inverse%.H: ${BINDIR}/Helicon2d.x  wells%.H interp.P sloping.H
	${BINDIR}/Helicon2d.x < wells$*.H >qdome-inverse$*.H par=interp.P angle=sloping.H \
		input=qdome-input$*.H
	
wells%.H overlay%.H: ${BINDIR}/make_wells.x panel.H wells%.P
	${BINDIR}/make_wells.x <panel.H par=wells$*.P >wells$*.H locations=overlay$*.H

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

wells%g.H overlay%.H: ${BINDIR}/make_wells.x panel.H wells%.P
	${BINDIR}/make_wells.x <panel.H par=wells$*.P >wells$*g.H locations=overlay$*.H geostat=1


clean: jclean

	

#END OF RELAVANT
	
${RESDIR}/qdome-input%.v:	qdome-inverse%.H
	Grey < qdome-input$*.H  ${GREY_OPT} out=${RESDIR}/qdome-input$*.v


${RESDIR}/qdome-inverse%.v: qdome-inverse%.H overlay%.H
	Grey < panel.H  ${GREY_OPT} bias=1.00486  allpos=y out=junk1.v title="Correct model"
	Grey < qdome-inverse$*.H out=junk2.v ${GREY_OPT} title=" J 1/A x + (I-J) m_k "
	vp_SideBySideAniso junk1.v junk2.v >${RESDIR}/qdome-inverse$*.v 

${RESDIR}/qdome-inverse%.v3: ${RESDIR}/qdome-inverse%.v panel.H qdome-inverse%.H
	Grey < panel.H  ${GREY_OPT} bias=1.00486  allpos=y out=junk1.v title="Goal"
	Grey < qdome-input$*.H  ${GREY_OPT} out=junk2.v title="Input"
	Grey < qdome-inverse$*.H out=junk4.v ${GREY_OPT} title="Inverse"
	vp_Movie junk1.v junk2.v ${RESDIR}/qdome-inverse$*.v \
    junk4.v >${RESDIR}/qdome-inverse$*.v3
#	-rm junk1.v junk2.v junk4.v
	

${RESDIR}/qdome-dip-field.v:	sloping.H
	Grey < sloping.H  ${GREY_OPT} bias=-16.69 wantscalebar=y \
    out=${RESDIR}/qdome-dip-field.v

${RESDIR}/qdome-reflectors.v:	panel.H slopes.H sloping.H
	Grey < panel.H  ${GREY_OPT} bias=1.05   allpos=n out=junk1.v  crowd=.75
	Window < slopes.H  j1=2| Graph ${GRAPH_OPT}  out=junk2.v d1=.008 o1=0 o2=0 d2=.01 yreverse=y  plotcol=2 symbolsz=20 crowd=.75
	Window < slopes.H j1=2 |Graph  ${GRAPH_OPT}  out=junk3.v d1=.008 o1=0 o2=0 d2=.01 yreverse=y plotcol=7 symbolsz=7 crowd=.75
	vp_Overlay junk1.v junk2.v junk3.v > junk4.v
	Grey < sloping.H  ${GREY_OPT} bias=-16.69 wantscalebar=y barlabel="dip (degrees)"\
    out=junk5.v crowd=.75 labelsz=16
	vp_SideBySideAniso junk4.v junk5.v >$@
#	-rm junk1.v junk2.v junk3.v junk4.v junk5.v


${RESDIR}/qdome-combo1.v:	${RESDIR}/qdome-inverse1.v qdome-inverse1.H
	Grey < qdome-input1.H  ${GREY_OPT} out=junk1.v  title="Well Logs" titlesz=16
	Grey < qdome-inverse1.H out=junk2.v ${GREY_OPT} title="Output" labelsz=16 titlesz=16
	Grey < panel.H  ${GREY_OPT} bias=1.00486  allpos=y out=junk3.v  title="Correct Model" titlesz=16
	vp_SideBySideAniso junk3.v junk1.v junk2.v >$@
#	-rm junk1.v junk2.v junk3.v

${RESDIR}/qdome-combo4.v:	${RESDIR}/qdome-inverse4.v qdome-inverse4.H
	Grey < qdome-input4.H  ${GREY_OPT} out=junk1.v  title="Well Logs" 
	Grey < qdome-inverse4.H out=junk2.v ${GREY_OPT} title="J 1/A x + (I-J) m_k"
	Grey < panel.H  ${GREY_OPT} bias=1.00486  allpos=y out=junk3.v  title="Correct Model"
	vp_SideBySideAniso junk3.v junk1.v junk2.v >${RESDIR}/qdome-combo4.v
#	-rm junk1.v junk2.v junk3.v
	


	

include ${SEPINC}/SEP.bottom