#! /bin/ksh
# shell for smoothing a uniformly sampling velocity  and for generating a
# slowness perturbation
set -v


n1=201 r1=5 dz=0.1 fz=0 labelz="Depth (m)"
n2=801 r2=5 dx=0.1 fx=0 labelx="Midpoint (m)"
d1r=0.5 n1r=41 n2r=161 d2r=0.5 f2r=0.
vsmoo3d<vfile r1=$r1 r2=$r2 d1=$dz d2=$dx n1=$n1 n2=$n2 n3=1 \
        verbose=1 iter=3 slowness=1>data.s
 resamp2<data.s n1_in=$n1 n2_in=$n2 f2_out=0. d1_in=$dz d2_in=$dx \
         n1_out=$n1r n2_out=$n2r d1_out=$d1r d2_out=$d2r>data.r

# calculate slowness perturbation, smooth it and resample it
slowp<vfile vsfile=data.s nz=$n1 nx=$n2|
vsmoo3d r1=$r1 r2=$r2 d1=$dz d2=$dx n1=$n1 n2=$n2 n3=1 \
	verbose=0 iter=3 slowness=0|
 resamp2 n1_in=$n1 n2_in=$n2 d1_in=$dz d2_in=$dx \
	 n1_out=$n1r n2_out=$n2r d1_out=$d1r d2_out=$d2r>vpfile
	 
 # plot smoothed velocity
psimage <data.r  style=seismic \
	n1=$n1r d1=$d1r f1=$fz grid1=dot label1="$labelz" \
	n2=$n2r d2=$d2r f2=$fx grid2=dot label2="$labelx" \
	ybox=4 hbox=4 \
	title="Smoothed Velocity"\
	|rsh dobrin -l mxu images

rm -f data.s
 exit