# linear moveout
#
subroutine lmo( conj,add, slow, tau0, t0,dt, x0,dx, modl,nt,nx,  data       )
integer         conj,add,                                nt,nx,       it,ix,iu
real t, x, tau,           slow, tau0, t0,dt, x0,dx, modl(nt,nx), data(nt,nx)
call conjnull( conj, add,                           modl,nt*nx,  data,nt*nx)
do ix= 1, nx {	x= x0 + dx * (ix-1)
do it= 1, nt {	t= t0 + dt * (it-1)
		tau = tau0 + t - x * slow
		 iu = 1.5 + tau/dt
		if( 0 < iu  &  iu <= nt)
			if( conj == 0 )
				data(it,ix) = data(it,ix) + modl(iu,ix)
			else
				modl(iu,ix) = modl(iu,ix) + data(it,ix)
		}}
return; end