ptam.h

以下代码实现的是 峰谷平均法 ，参考：

    1. 张海明. 2021. 地震学中的Lamb问题（上）. 科学出版社
    2. Zhang, H. M., Chen, X. F., & Chang, S. (2003). 
          An efficient numerical method for computing synthetic seismograms 
          for a layered half-space with sources and receivers at close or same depths. 
          Seismic motion, lithospheric structures, earthquake and volcanic sources: 
          The Keiiti Aki volume, 467-486.

Author

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

Date

2024-07-24

Functions

void grt_PTA_method(GRT_MODEL1D *mod1d, real_t k0, real_t predk, size_t nr, real_t *rs, cplxIntegGrid sum_J0[nr], bool calc_upar, cplxIntegGrid sum_uiz_J0[nr], cplxIntegGrid sum_uir_J0[nr], FILE *ptam_fstatsnr[nr][2], GRT_KernelFunc kerfunc)

峰谷平均法 Peak-Trough Averaging Method，最后收敛的积分结果以三维数组的形式返回，

参数:

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

• k0 – [in] 先前的积分已经进行到了波数k0

• predk – [in] 先前的积分使用的积分间隔dk，因为峰谷平均法使用的 积分间隔会和之前的不一致，这里传入该系数以做预先调整

• nr – [in] 震中距数量

• rs – [in] 震中距数组

• sum_J0 – [out] 积分值

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

• sum_uiz_J0 – [out] uiz的积分值

• sum_uir_J0 – [out] uir的积分值

• ptam_fstatsnr – [out] 峰谷平均法过程文件指针数组

• kerfunc – [in] 计算核函数的函数指针