crs3s/capjoint

CAPjoint

CAPjoint 是中科大的 Weiwen Chen 和 Zhe Jia 基于朱露培教授的 CAP 发展而来，其优势在于可以利用远震数据。 CAPjoint 的官方发布地址是 http://home.ustc.edu.cn/~vincentc/CAPjoint

这个 repo 是我在学习 CAPjoint 后整理并用 perl 改写其主要脚本后的脚本集。对 CAPjoint 方法本身并无任何改进。 如果按照下面的流程，原则上可以把官方包内的例子跑出来。详细的笔记请看我的博客

编译与安装

地震学科研一般都要用的工具在此处不指导如何安装，只列出。

1. SAC
2. Tau-P
3. fk
4. GMT4
5. PSSAC2

下载

git clone [email protected]:wangliang1989/capjoint.git

编译 crust2

进入 src/crust2,编译：

gfortran getCN2point.f -o getCN2point

编译 teles

进入 src/tel3,编译：

make

编译 CAP

进入 src/cap,编译：

make

在编译完 teles 和 cap 之后，会分别生成可执行文件 teles 和 cap，请务必把他们拷贝到 example/cmds 下

运行例子

例子放在 example 路径下。

数据预处理

perl process.pl 20140310

计算格林函数

perl grn.pl 20140310

生成权重文件

perl weight.pl 20140310

cap 反演

perl inversion.pl 20140310