Sound propagation algorithm for a sonar simulator is required to run in real-time and should be able to model the range and depth dependence of the Korean ocean environments. Ray model satisfies these requirements and we developed an algorithm for range-dependent ocean environments. In this algorithm, we considered depth-dependence of sound speed through rays based on a rectangular cell method and layer method. Range-dependence of sound speed was implemented based on a split-step method in the range direction. Eigen-ray is calculated through an interpolation of ray bundles and Gaussian interpolation function was used. The received time signal of sonar was simulated by Fourier transform of eigen-ray solution in the frequency domain. Finally, for the verification of proposed algorithm, we compared the results of transmission loss with other validated models such as BELLHOP, SNUPE, KRAKEN and OASES, for the Pekeris waveguide, wedge, and deep ocean environments. As a result, we obtained satisfactory agreements among them.