Introduction - If you have any usage issues, please Google them yourself
1, parameter initialization
2, the formation of boundary clipping data
3, wave field initialization
4, Zmax layer wave field extended to the depth of Zmax-1
5, Zi+1 layer wave field extended to the depth of Zi
6, repeat 5, starting Iz Nz-1 until Iz 1, the frequency of the test line- spatial domain wave field
7, frequency space domain wave field to frequency inverse Fourier transform, time- space field