static_grn.h

以下代码实现的是 广义反射透射系数矩阵+离散波数法 计算静态格林函数，参考：

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

Author

Zhu Dengda (zhudengda@mail.iggcas.ac.cn)

Date

2025-04-03

Functions

void grt_integ_static_grn(GRT_MODEL1D *mod1d, size_t nr, real_t *rs, K_INTEG_METHOD *Kmet, bool calc_upar, realChnlGrid grn[nr], realChnlGrid grn_uiz[nr], realChnlGrid grn_uir[nr], const char *statsstr)

积分计算Z, R, T三个分量静态格林函数的核心函数

参数:

• mod1d – [inout] GRT_MODEL1D 结构体指针

• nr – [in] 震中距数量

• rs – [in] 震中距数组

• Kmet – [inout] 波数积分相关参数的结构体指针

• calc_upar – [in] 是否计算位移u的空间导数

• grn – [out] 浮点数数组，不同震源不同阶数的静态格林函数的Z、R、T分量

• grn_uiz – [out] 浮点数数组，不同震源不同阶数的ui_z的Z、R、T分量

• grn_uir – [out] 浮点数数组，不同震源不同阶数的ui_r的Z、R、T分量

• statsstr – [in] 积分过程输出目录