modal_util.h

Author

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

Date

2025-08

频散相关的辅助函数

Functions

void grt_output_cdisp(const char *filepath, const char *full_command, const char *modelpath, EIGENV_INFO *eigmet)

输出相速度频散结果

参数:

• filepath – [in] 输出路径

• full_command – [in] 计算频散的完整命令

• modelpath – [in] 模型路径

• eigmet – [in] 频散数据结构体指针

void grt_output_udisp(const char *filepath, EIGENFN_INFO *eigfnmet)

输出群速度频散结果

参数:

• filepath – [in] 输出路径

• eigfnmet – [in] 本征函数数据结构体指针

void grt_output_eigenfns(const char *filepath, const int ncols, EIGENFN_INFO *eigfnmet)

输出本征函数结果

参数:

• filepath – [in] 输出路径

• ncols – [in] 4 for Rayleigh, 2 for Love

• eigfnmet – [in] 本征函数数据结构体指针

void grt_output_energy_integrals(const char *filepath, EIGENFN_INFO *eigfnmet)

输出能量积分结果

参数:

• filepath – [in] 输出路径

• eigfnmet – [in] 本征函数数据结构体指针

void grt_output_sensitivity(const char *filepath, const char *char_uc, EIGENFN_INFO *eigfnmet)

输出相速度/群速度敏感核

参数:

• filepath – [in] 输出路径

• char_uc – [in] ‘c’ for phase velocity, ‘u’ for group velocity

• eigfnmet – [in] 本征函数数据结构体指针

void grt_group_sensitivity(EIGENFN_INFO *eigfnmet)

根据计算好的相速度、群速度和相速度敏感核，计算群速度敏感核

根据群速度与相速度的关系： U = c^2 / ( c - w*dc/dw ) 推导得到敏感核（其中 X 表示 alpha/beta/rho ） dU/dX = U/c * (2 - U/c) * dc/dX + (U/c)^2 * w * d(dc/dX)/dw 后一项存在对相速度敏感核在频率上做差分， 因此输入的频率间隔大小会直接影响群速度敏感核精度

参数:

• eigfnmet – [inout] 本征函数数据结构体指针

void grt_read_dispersion(const char *filepath, EIGENV_INFO *eigmet, char **pt_modelpath)

读取相/群速度频散结果

参数:

• filepath – [in] 读入频散文件路径

• eigmet – [out] 频散数据结构体指针

• pt_modelpath – [out] 模型路径