1、参数初始化；

2、形成边界削波数据；

3、波场初始化；

4、Zmax层波场延拓到深度Zmax-1；

5、Zi+1层波场延拓到深度Zi；

6、重复5，从Iz Nz-1开始，直到Iz 1，得测线上的频率—空间域波场；

7、频率-空间域波场对频率做反傅里叶变换，得时间-空间波场；

-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