shotloc1 = zshot=6000 zsou=6000 yshot=14000 xsou=14000
shotloc2 = zshot=3990 zsou=3990 yshot=4250.0 xsou=4250.0
shotloc3 = zshot=3990 zsou=3990 yshot=6612.5 xsou=6612.5
shotloc4 = zshot=4390 zsou=4390 yshot=6987.5 xsou=6987.5

fmeminmax = min1=0.000 max1=2048. min2=1. max2=3.
hwtminmax = min1=1000. max1=1024. min2=1. max2=3.

vel.v: 
	Cp vel.H jvel.H
	echo "d1=.01 d2=0.0125" >>jvel.H
	< jvel.H $(grey) $(ratio) bias=3.3 title=' ' out=$@

#*********************************
# 2-D Fast marching eikonal
#*********************************
fme%.H:
	< vel2.H FMeikonal $(shotloc$*) >$@
fme%.v: fme%.H
	< fme$*.H Window f1=3 |\
	Contour c0=.123 dc=0.2 nc=30 $(dn) $(ratio) $(noframe) plotcol=7 out=$@
FME%.v: fme%.v 
	vp_Overlay vel2.v fme$*.v >$@

$R/FMpicks%.v: fme%.H
	< fme$*.H Window n1=1 f1=0 | Add add=0.015               > jnk1.H
	Spike n1=515 o1=0 d1=1    | Add add=3.    | Tpow tpow=1 > jnk2.H
	Cmplx jnk2.H jnk1.H | $(graph) $(fmeminmax) $(noframe) ${dn} out=junkfmtt.v
	vp_Overlay  $R/fdseis$*.v junkfmtt.v >$@
	vp_Movie $@ $R/fdseis$*.v            >$@3
	Rm jnk2.H jnk1.H
	\rm junkfmtt.v

#*********************************
# 2-D HWT  par=$P/wf.P 
#*********************************
Hwt1.H: 
	< vel2.H Pad extend=1 beg2=1000 end2=1000 end1=1001 > vv.H
	< vv.H Hwt2d zsou=15000 xsou=14000 par=$P/wf.P > $@



hwt1.v:
	< Hwt1.H Window n1=1 f1=123 | Transp | ${graph} $(noframe) plotcol=1 symbol="." plotfat=3 ${dn} out=ray.v
	< Hwt1.H Window j2=15 |${graph} $(noframe) plotcol=2 ${dn} out=ray2.v
	vp_Overlay ray.v ray2.v >$@

HWT1.v:
	< vv.H Grey newclip=1 color=j ${dn} out=vv.v
	vp_Overlay vv.v hwt1.v >$@

#*********************************
# 2-D Comparison
#*********************************
$R/Shot%.v:      $R/fdsnap%.v  $R/FME%.v  $R/HWT%.v
	vp_Movie $R/fdsnap$*.v $R/FME$*.v $R/HWT$*.v >$@

#*********************************
# 3-D velo salt
#*********************************
#segsalt = /net/kana/wrk1/velocity/SEG-EAEG/velocity/salt_slow_transp_cmp_desp.h
#segsalt = $(SEGSALT_DATA)/slowness.HH
segsalt = /data/data_syn/seg_eage_salt/slowness.HH
TAPLOT = Taplot pclip=100 gainpanel=all allpos=y
CUBEPLOT = Cubeplot frame1=10 frame2=40 frame3=45 color=j wheretitle=l title=" " \
	label2="X location" label3="Y location" label1="depth" point1=.3 wantcoordlabel=-1

slowsalt.H: ${segsalt}
	cp ${segsalt} .
	echo "in=${segsalt}@" >> slowness.HH
	<   slowness.HH Window j1=5 j2=10 j3=10 |\
	Smooth rect1=3 rect2=3 rect3=3 >$@
velosalt.H: slowsalt.H
	<   slowsalt.H Vconvert invelopt=1 outvelopt=2 >$@

# 3-D FME on salt
3dfme%.H: velo%.H $P/wf%.P
	< velo$*.H Transp plane=23 | FMeikonal par=$P/wf$*.P |\
	Transp plane=23 >$@
$R/3dfme%.v: 3dfme%.H
	< 3dfme$*.H ${TAPLOT} > jt$*.T
	< jt$*.T ${CUBEPLOT} >${dn} out=$@
	Rm jt$*.T

# 3-D HWT on salt
3dhwt%.H: velo%.H $P/wf%.P
	< velo$*.H Hwt3d par=$P/wf$*.P interp=1 verb=1 fast=1 >$@
$R/3dhwt%.v: 3dhwt%.H
	< 3dhwt$*.H ${TAPLOT} > jt$*.T
	< jt$*.T ${CUBEPLOT} >${dn} out=$@
	Rm jt$*.T