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(MODEL1D_STATE *mstat, real_t k0, real_t predk, size_t nr, real_t *rs, K_INTEG *K, FILE *ptam_fstatsnr[nr][2], GRT_KernelFunc kerfunc)

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

参数:

• mstat – [inout] MODEL1D_STATE 结构体指针

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

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

• nr – [in] 震中距数量

• rs – [in] 震中距数组

• K – [inout] 用于存储积分的结构体

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

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