For reverse time migration one step of the algorithm loop on page
is carried out to simulate a shot.
Then that synthetic wave field
is backward propagated, simultaneously with the real recorded data,
which are fed in time-reverse at the geophone locations. At each time step the
cross correlation between the strains of the two wavefields is
summed for each spatial position.