说明：  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 )