FK manual
fk 是 Prof.Lupei Zhu 写的一个用于计算 水平分层模型 下的理论格林函数并合成 理论地
震图的代码包。代码是开源的，可以直接编译使用。代码中包含了如下几个重要的命令 / 脚
本：
• fk ：用于计算格林函数的主程序；
• st_fk ：用于计算静态格林函数的主程序；
• trav ：用于计算 P、S 初至到时的辅助程序；
• sachd ：用于修改 SAC 头段的辅助程序；
• fk.pl ：对 fk 、 st_fk 、 trav 和 sachd 的封装，一般情况下 直接使用 fk.pl
脚本即可；
• syn ：用于将格林函数合成为理论地震图三分量的程序；
因而，实际操作的时候，只需要调用 fk.pl 生成格林函数，再调用 syn 将格林函数合成为
三分量地震图即可。
使用说明
先把用法贴在这里作为参考:
fk.pl -Mmodel/depth[/f_or_k] [-D] [-Hf1/f2] [-Nnt/dt/smth/dk/taper]
[-Ppmin/pmax[/kmax]] [-Rrdep] [-SsrcType] [-Uupdn] [-Xcmd] distances
速度模型
-Mmodel/depth/k_or_f 中的 model 为模型文件名， k 和 f 的作用在下面会解释。
格式
fk 的输入速度模型是一维水平分层速度模型，其格式为:
thickness vs vp/vs [rho Qs Qp]
或:
thickness vs vp [rho Qs Qp]
其中
• 列 1：该层的厚度（km）
• 列 2：S 波速度（km/s）
• 列 3：波速比或 P 波速度
• 列 4：密度（ g/cm3
• ）
• 列 5：S 波的 Q 值
• 列 6：P 波的 Q 值
评论1