static_layer.h

以下代码实现的是静态反射透射系数矩阵，参考：

    1. 姚振兴, 谢小碧. 2022/03. 理论地震图及其应用（初稿）.  
    2. 谢小碧, 姚振兴, 1989. 计算分层介质中位错点源静态位移场的广义反射、
         透射系数矩阵和离散波数方法[J]. 地球物理学报(3): 270-280.

Functions

void grt_static_topfree_RU(GRT_MODEL1D *mod1d)

计算自由表面的静态反射系数，公式(6.3.12)

void grt_static_wave2qwv_REV_PSV(GRT_MODEL1D *mod1d)

计算接收点位置的 P-SV 静态接收矩阵，将波场转为位移，公式(6.3.35)

void grt_static_wave2qwv_REV_SH(GRT_MODEL1D *mod1d)

计算接收点位置的 SH 静态接收矩阵，将波场转为位移，公式(6.3.37)

void grt_static_wave2qwv_z_REV_PSV(GRT_MODEL1D *mod1d)

计算接收点位置的ui_z的 P-SV 静态接收矩阵，即将波场转为ui_z。公式本质是推导ui_z关于q_m, w_m, v_m的连接矩阵（就是应力推导过程的一部分）

void grt_static_wave2qwv_z_REV_SH(GRT_MODEL1D *mod1d)

计算接收点位置的ui_z的 SH 静态接收矩阵，即将波场转为ui_z。公式本质是推导ui_z关于q_m, w_m, v_m的连接矩阵（就是应力推导过程的一部分）

void grt_static_RT_matrix_PSV(const GRT_MODEL1D *mod1d, const size_t iy, RT_MATRIX *M)

计算界面的 P-SV 波静态反射透射系数RD/RU/TD/TU 根据公式(6.3.18)

参数:

void grt_static_RT_matrix_SH(const GRT_MODEL1D *mod1d, const size_t iy, RT_MATRIX *M)

计算界面的 SH 波静态反射透射系数RDL/RUL/TDL/TUL 根据公式(6.3.18)

参数:

void grt_static_delay_RT_matrix(const GRT_MODEL1D *mod1d, const size_t iy, RT_MATRIX *M)

为静态 R/T 矩阵添加时间延迟因子

参数: