A method for calculating seismic velocity for migration purposes as a function of subsurface spatial position that gives the seismic processing analyst direct control of the resulting migration velocity model, and includes a means to ensure that the model is consistent with RMS velocity sweeps computed from available surface seismic or other data. The method involves generating average velocities from the model to compare to the measured RMS velocity data, and comparing the model predictions to actual data on a 3-D interactive visual display allowing quick and easy model adjustment. Other available measured data besides surface seismic data maybe used to develop the model.