Dengda Zhu

At present, **PyGRT** support several wavenumber integration methods. However, no matter C or Python, the way settting the method is not apparent and elegant. Given that the Direct Convergence Method...

repeated k when mix `dwm` and `safim`

1

e.g. `-La20/1e-2/10` Check the output stats can reproduce this tiny bug.

Improvement for different source depths in the same media layer

1

It provides a significant speed improvement for the finite fault inversion. The steps are simple, see [Yao and Zheng (1984)](http://www.geophy.cn/article/id/cjg_5105).

《理论地震图及其应用》预计在2026年出版，相比于初稿版本，公式索引会有变化，因此代码注释中的索引也需要更改，以方便初学者学习对照。

Currently, in `kmax` formula, there's a static offset $\dfrac{5\pi}{\Delta h}$ , which may increase unnecessary calculation when the sampling frequency is low (e.g. at a large distance).

`CAP` 和 `gCAP` 是反演震源机制的经典方法，从 `GRT` 程序包的求解方法、代码框架来看，很适合以模块方式加入这两个算法。 算法原理上并不复杂，前人已经做了很多尝试、验证和代码开发，例如王亮老师的 [oh-my-cap](https://github.com/wangliang1989/oh-my-cap)，可以参考学习。 中科院《理论地震学》课程也有实现 `CAP` 算法的小组合作作业，有想法的同学也可以提意见。 临近毕业，事多，这里做个备忘。

Does it make sense?

1. `+g` behaves strange with `+w` , but single modifer works fine. ``` bash gmt begin panels png gmt subplot begin 2x2 -Ff6c+w1p,-+gyellow -R0/10/0/10 gmt basemap -c gmt subplot end...