#####################################################
#
# . . Migration
#

Prep.%.H:
	echo "d1=5 d2=10" >> PS.rays.H
	echo "d1=10" >> PS.data.H
	< PS.rays.H Real | Interp n2out=$* d2out=4.95 type=1 > xx.H
	< PS.rays.H Imag | Interp n2out=$* d2out=4.95 type=1 > zz.H
	Cmplx xx.H zz.H > ray.H
#	Cp PS.rays.H ray.H
	< PS.vel1.H $(B)/CC2RC.x adj=0 rays=ray.H > vel.H 
	< PS.data.H Real | Window | Interp n1out=$* d1out=4.95 > vxx.H
	< PS.data.H Imag | Window | Interp n1out=$* d1out=4.95 > vzz.H
	Cmplx vxx.H vzz.H > dat.H
	echo "d1=4.95" >> dat.H
#	< PS.data.H Window3d > dat.H
TEST1.%.H: Prep.%.H #PS.data.H PS.vel.H PS.rays.H                   
	$(B)/RWE2D_source.x rays="ray.H" Rimage="r$@" \
		vel="vel.H" rwf="dat.H" image="$@" \
		nref=625 verbose=3 min_region_pct=1. del_dist=1. \
		verb=1 niter_lloyd=40 ntap=6 nloop=5 norm=1 \
		xmin=3000 xmax=13000 zmin=0 zmax=5000 \
		dxx=10 dzz=10 nsx=2 nsz=2 > $@

TEST2.%.H: Prep.%.H #PS.data.H PS.vel.H PS.rays.H                   
	$(B)/RWE2D_test.x rays="ray.H" Rimage="r$@" \
		vel="vel.H" rwf="dat.H" image="$@" \
		nref=625 verbose=3 min_region_pct=1. del_dist=1. \
		verb=1 niter_lloyd=40 ntap=6 nloop=5 norm=1 \
		xmin=3000 xmax=13000 zmin=0 zmax=5000 \
		dxx=10 dzz=10 nsx=2 nsz=2 > $@