A finite-difference algorithm appropriate for modeling acoustic waves in a fully heterogeneous moving 3D media has been developed. The model is characterized by: acoustic velocity, density, and the three components of the background media velocity. The approach solves a set of coupled 1st order velocity-pressure differential equations appropriate for an adiabatic divergence-free background velocity. The equations are staggered in time and space and the algorithm uses second order temporal and fourth order spatial finite-differences. Since approximations are not adopted in the solution of the equations all arrivals are modeled with fidelity providing the spatial and temporal grids are chosen appropriately. The algorithm can include either a pressure or velocity free surface on the bottom boundary and absorbing boundaries on other model flanks. Designed to run on large scale parallel computational platforms, the algorithm has been validated for four machine architectures. Comparisons are presented to an analytic solution for a constant wind model and fast-field program results for a vertically stratified wind model. Data resulting from simulations through a kinematic turbulence wind profile developed with the quasi-wavelet method are also presented. Sandia National Laboratories is operated by Sandia Corporation, a Lockheed Martin Company, for the USDOE under Contract No. 94-AL85000.