MODELING STRUCTURE

 Following is an outline of the modeling algorithm for a general wave equation (1). Compare this to subroutine wave() .

loop over shots {

calculate accelerations {

calculate strain

stress = b(x) $\times$ strain

apply stress boundary conditions

(add in stress sources)

constrain stress (eg. free surface)

calculate divergence of stress

accel = a(x) $\times$ div of stress

}

add in acceleration sources

time update the wave field

record seismograms

}