At a free surface the stress free boundary condition must be satisfied.
As discussed in the appendix this results in matrix operators
that convert applied tractions to particle velocities and downgoing wave amplitudes,
and operators that convert upgoing wave amplitudes into velocities and
downgoing amplitudes.
![]() |
||
![]() |
||
These four operators are required in performing steps one and six of the algorithm.