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.