We have developed 2D and 3D finite-difference codes for modeling seismic wave propagation. These are based on the velocity-stress formulation of the elastodynamic equation, which is discretized by a staggered-grid scheme second-order accurate in time and second- or fourth-order accurate in space. In this article we describe the outline of the codes and boundary conditions we employed. We developed a free surface condition for 3D finite-difference scheme, which is an extension of a free surface condition for 2.5D finite-difference scheme developed by Okamoto(1994). We also present some numerical examples for demonstrating the accuracy of our codes.
|Number of pages
|Science Reports of the Kyushu University, Department of Earth and Planetary Sciences
|Published - Dec 1 1998
ASJC Scopus subject areas
- Space and Planetary Science
- General Earth and Planetary Sciences