Numerical simulations of line-profile variation beyond a single-surface approximation for oscillations in roAp stars