def goNmoVigGather():
""" Applies NMO correction to Viking Graben gather """
st,sx,f = readVigGather(800)
nt,nx = st.count,sx.count
dt,dx = st.delta,sx.delta
ft,fx = st.first,sx.first
ts=0.010,0.646,0.846,1.150,1.511,1.825,2.490,2.985
vs=1.510,1.575,1.687,1.817,1.938,1.980,2.446,2.735
vnmoP = CubicInterpolator(ts,vs).interpolate(rampfloat(ft,dt,nt))
vnmoM = fillfloat(1.475,nt)
vnmo = vnmoM
texp,tbal,smax = 2.0,100,9.00
tmin,tmax,perc = 0.4,2.4,98.0
#f = tpow(texp,st,f)
if tbal>0:
f = balance(tbal,f)
nmo = NormalMoveout()
nmo.setStretchMax(smax)
g = nmo.apply(st,sx,vnmo,f)
h = nmo.stackAndReplicate(g)
plotGather(st,sx,f,tmin=tmin,tmax=tmax,perc=perc,title="input");
plotGather(st,sx,g,tmin=tmin,tmax=tmax,perc=perc,title="output");
plotGather(st,sx,h,tmin=tmin,tmax=tmax,perc=perc,title="stack");
