SPICE logo


The program is designed for computation of seismic wavefields in 3D heterogeneous surface geological structures with planar free surface due to surface and near-surface point doublecouple sources

Authors: Jozef Kristek, Peter Moczo
General: Wave propagation
Code level: Production Code
Language: Fortran 90/95 or later
Associated groups: Local Scale
Supercomputing: Shared memory (e.g., openMP)
Grid dimensions: 3D
Coordinate system: Cartesian
Method: Finite differences
Tested operating system: Linux, Unix (e.g. SunOS, HP Unix), Windows
Rheology: Viscoelastic isotropic
Boundary conditions: Free surface (planar), Absorbing boundaries, Reflecting boundaries
Grid type: Rectangular, regular, Other
Solution type: Numerical
Total hits: 9095

Additional resources

Latest Releases

Project Description

Program 3DFD_DVS is designed for the finite-difference simulation of seismic wave propagation and seismic ground motion in a 3D surface heterogeneous viscoelastic structure with a planar free surface. The computational algorithm is based on the explicit heterogeneous finite-difference scheme solving equations of motion in the heterogeneous viscoelastic medium with material discontinuities. The scheme is 4th-order accurate in space and 2nd-order accurate in time. The displacement-velocity-stress scheme is constructed on a staggered finite-difference grid. The computational region is represented by a volume of a parallelepiped with the top side representing a planar free surface, and bottom, rear, front, left and right sides representing either non-reflecting boundaries or planes of symmetry. Different types of non-reflecting boundaries can be chosen on different sides of the computational region. The discontinuous spatial grid is used to cover the computational region. The upper part of the grid has three times smaller grid spacing than the lower part. Each part itself is a uniform rectangular grid. The rheology of the medium corresponds to the generalized Maxwell body. This makes it possible to account both for spatially varying quality factors of the P and S waves and for arbitrary Q-omega law. The wavefield is excited by a set of point double-couple sources. The core memory optimization is applied in order to significantly reduce requirements of the computer’s core memory.