- loop over shots {
- calculate accelerations {
- calculate strain
- stress = b(x) strain
- apply stress boundary conditions
- (add in stress sources)
- constrain stress (eg. free surface)

- calculate divergence of stress
- accel = a(x) div of stress
- }

- add in acceleration sources
- time update the wave field
- record seismograms

- }

- Strain calculation
- Stress calculation
- Acceleration calculation
- Convolution operator
- The wave operator
- Time Update

11/17/1997