PyGRT
PyGRT copied to clipboard
Dengda98
An Efficient and Integrated C/Python Package for Computing Synthetic Seismograms, Strain, Rotation and Stress Tensor in a Layered Half-Space Model (Dynamic & Static Cases)
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...
e.g. `-La20/1e-2/10` Check the output stats can reproduce this tiny bug.
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` 算法的小组合作作业，有想法的同学也可以提意见。 临近毕业，事多，这里做个备忘。