By using the elastodynamic equations rather than wave equations for seismogram synthesis, multi-attribute wave fields can be computed in a single execution of one program. In the present implementation, for 2-D models, the wave fields and seismograms produced include any or all of the following: two components of each of particle acceleration, particle velocity, and particle displacement; two components of normal strain; shear strain; two components of normal stress; shear stress; and the dilatation and curl of the particle displacement. If flexible source and receiver configurations are also included, a single piece of software can be used for most seismic modeling applications. This significantly reduces the need for development and maintenance of separate programs. The algorithm is illustrated using waves synthesized for an explosive source in a layered medium. Snapshots and seismograms at both surface and borehole arrays are presented.