# . .  Coordinate System and 
#
# Test0.rays.HH
#
# . . and velocity
#
# Test1.vel.HH
#

########
#
# . .  Source Programs
haha: 
	@-cd Src/; ${MAKE} default ; cd ..

$(XBIN)/%.x: haha

weipar0= min5=0.1 max5=42.5 j5=1 squeeze=n
rwepar0 = nref=1 verbose=1 forward=1 kinematic=1. min_region_pct=.1\
		 del_dist=.05 verb=0 niter_lloyd=20 ntap=20 nloop=3
nx0=400
ox0=0
dx0=10
ny0=1
oy0=0
dy0=20
nz0=400
dz0=5
oz0=0
CCout0  =azn=$(nz0) azo=$(oz0) azd=$(dz0)\
	 axn=$(nx0) axo=$(ox0) axd=$(dx0)\
	 ayn=$(ny0) ayo=$(oy0) ayd=$(dy0)\
	 nsx=1 nsz=1 nsy=0 adj=1
HWTgraph = Graph crowd=${crowd} wantaxis=n title=" " label1=" " label2=" "
XBIN=./Bin/
Test0.vel.HH Test0.rays.HH: 
	matlab < ./matlab/Test0.m 
	Cp Vel0.HH  Test0.vel.HH
Test0.data.H:
	Wavelet n1=2048 d1=0.004 wavelet=ricker2 o1=0. fund=25. tdelay=0.2 > t1.H
	Wavelet n1=2048 d1=0.004 wavelet=ricker2 o1=0. fund=25. tdelay=0.4 > t2.H
	Wavelet n1=2048 d1=0.004 wavelet=ricker2 o1=0. fund=25. tdelay=0.6 > t3.H
	Wavelet n1=2048 d1=0.004 wavelet=ricker2 o1=0. fund=25. tdelay=0.8 > t4.H
	Add t1.H t2.H t3.H t4.H > t5.H
	Pad < t5.H end2=255 extend=1 > $@

ww.Test0.H: Test0.data.H
	<   Test0.data.H Transf wei=y verb=y >t$@
	<  t$@ Window3d $(weipar0) | Transp plane=35 > $@
	echo "d1=10 d2=10" >> $@
	\rm t$@ 
rwepar0 = nref=1 verbose=1 forward=1 kinematic=1. min_region_pct=.1\
	del_dist=.05 verb=0 niter_lloyd=20 ntap=20 nloop=3

Test0.mig.H: ww.Test0.H Test0.rays.HH Test0.vel.HH  $(XBIN)/RWE3D_source.x 
	echo "n4=3 n5=1" >> Test0.rays.HH           
	< ww.Test0.H $(B)/RWE3D_source.x rays=Test0.rays.HH \
		vel=Test0.vel.HH swf=ww.Test0.H $(rwepar0)> $@

#	< Test0.rays.HH Window3d n4=2 j4=2 | Window3d |\
#		Interp n1out=300 d1out=5 o1out=0 n2out=512 d2out=5 o2out=0 type=1> ray.H
#	echo "d1=0.5 d2=0.5" >> ray.H	
#	echo "d1=0.5 " >> wf.H	
#	< Test0.vel.HH Pad end1=150 end2=256 extend=1  > vel.H
#	$(XBIN)/RWE2D_source.x rays="ray.H" \
#		vel="vel.H" rwf="wf.H" swf="wf.H"\
#		image="$@" nref=1625 verbose=1  \
#		min_region_pct=5 del_dist=5 verb=1 \
#		niter_lloyd=0 ntap=10 nloop=3 migapp=512 \
#		norm=1 zmin=0 zmax=1600 xmin=0 xmax=3200 \
#		dxx=2.5 dzz=2.5 nsx=3 nsz=3 Rimage="Rimage0.H"> $@

CCtest0.H: Test0.mig.H $(XBIN)/C2Rsinc3D.x                Test0.rays.HH
	<  Test0.mig.H $(XBIN)/C2Rsinc3D.x $(CCout0) rays=Test0.rays.HH > $@ 

Rays0.v: Test0.rays.HH CCtest0.H
	< Test0.rays.HH Window3d  n4=2 f4=0 j4=2 > tt.H
	< tt.H Window3d n3=1 f3=0 > xx.H
	< tt.H Window3d n3=1 f3=1 > zz.H
	Cmplx xx.H zz.H > r.H
	<  r.H Window3d j1=1 j2=30 > kk1.H
	< kk1.H ${HWTgraph} $(dn) out=jr.v plotcol=2 min2=0
	<  r.H Window3d j1=10 j2=1 | Transp > kk2.H
	< kk2.H ${HWTgraph} label2="Depth (m)" min2=0 \
	label1="Distance (m)" labelrot=n wantaxis=n $(dn) out=jw.v plotcol=1
	vp_Overlay jr.v jw.v >t1.v
#	Window3d < CCtest0.H min1=0 max1=1500 min2=0 max2=3200  |\
	Window3d < Test0.mig.H min1=0 max1=1500 min2=0 max2=3200  |\
		Grey pclip=99 out=t2.v $(dn) grid=y g1num=300 g2num=3200  min1=0\
		title=" " label1=" " label2=" " crowd=0.8 wantaxis=y
	vp_Overlay t2.v t1.v > $@
	\rm -f tt.H xx.H zz.H r.H kk1.H kk2.H t1.v t2.v jw.v jr.v

################
#
# . . Figure 2

NewRT.v $(R)/Fig2.pdf: Test0.rays.HH Test0.mig.H CCtest0.H
	< Test0.rays.HH Window3d  n4=2 f4=0 j4=2 > tt.H
	< tt.H Window3d n3=1 f3=0 > xx.H
	< tt.H Window3d n3=1 f3=1 > zz.H
	Cmplx xx.H zz.H > r.H
	<  r.H Window3d j1=1 j2=15 f2=0> kk1.H
	< kk1.H ${HWTgraph} $(dn) out=jr.v plotcol=7 min2=0 max2=1350 max1=3200
	<  r.H Window3d j1=5 j2=1 f1=0| Transp > kk2.H
	< kk2.H ${HWTgraph} label2="Depth (m)" min2=0 max2=1350 max1=3200 \
	label1="Distance (m)" labelrot=y wantaxis=y $(dn) out=jw.v plotcol=7 \
		wherexlabel=t labelsz=10
	vp_Overlay jr.v jw.v >t1.v
	Window3d < CCtest0.H min1=0 max1=1350 min2=0 max2=3200  |\
		Grey pclip=100. out=t2.v $(dn) grid=y g1num=300 g2num=3200  min1=0\
		title=" " label2="Distance (m) " label1="Depth (m) " crowd=0.8 \
		wantaxis=y  labelsz=10
	vp_SideBySideAniso t1.v t2.v > $@
	\rm -f tt.H xx.H zz.H r.H kk1.H kk2.H t1.v t2.v jw.v jr.v
#	pstexpen $@ ./Fig/Fig2.ps fat=1 fatmult=1.5
#	epstopdf ./Fig/Fig2.ps


RTpar = nT=251 dT=0.005 oT=0. oG=-70 nG=131 dG=1 xsou=4000 zsou=0

MARM_DATA = /data/marmousi/
toporays.v: #Test0.rays.HH CCtest0.H
	< toporays.HH Window3d  j2=10 > tt.H
	< tt.H Real > xx.H
	< tt.H Imag > zz.H
	Cmplx xx.H zz.H > r.H
	<  r.H Window3d j1=1 j2=2 > kk1.H
	< kk1.H ${HWTgraph} $(dn) out=t2.v plotcol=7 \
		 min1=0 max1=24000 min2=-2000 max2=8000
	<  r.H Window3d j1=2 j2=1 | Transp > kk2.H
	< kk2.H ${HWTgraph} label2="Depth (m)" min2=0 \
		 min1=0 max1=24000 min2=-2000 max2=8000\
	label1="Distance (m)" labelrot=n wantaxis=y $(dn) out=t1.v plotcol=7
	vp_Overlay t1.v t2.v >$@

Problem.v:
	< $(MARM_DATA)/marmvel.H Smooth rect1=12 rect2=12 repeat=20 tridiag=y > vv.H
	< vv.H Hwt2d $(RTpar) > Mrays.H
	< Mrays.H Window3d j1=1 j2=1 > kk1.H
	< kk1.H ${HWTgraph} $(dn) wantaxis=n wantlabel=n\
		out=jr.v plotcol=7 min2=0 max2=3000 min1=0 max1=8000 
	< Mrays.H Window3d j1=15 j2=1 f1=10 | Transp > kk2.H
	< kk2.H ${HWTgraph} label2="Depth (m)" min2=0 max2=3000 min1=0\
		 max1=8000 \
	label1="Distance (m)" labelrot=n wantaxis=y $(dn) out=jw.v plotcol=7
	vp_Overlay jr.v jw.v >t1.v
# velocity
	 < vv.H Grey newclip=1 out=t2.v $(dn) min1=0 max1=3000\
		min2=0 max2=8000 wantaxis=n\
		title=" " label1="Depth (m) " label2="Distance (m) " \
		crowd=0.8 wantaxis=n wherexlabel=b
	vp_Overlay t2.v t1.v > v1.v 
# Foothills coords
	vp_OverUnderAniso v1.v toporays.v > $@