Merge branch 'main' into taup_curve

source/cps/introduction.rst

@@ -3,7 +3,7 @@
 
 CPS，全称 Computer Programs in Seismology，
 着重于理解和解释地震波在地球地壳和上地幔中的传播。
-其可以计算平层介质中的理论地震图、反演面波频散和远震 P 波接收函数、
+其可以计算\ **水平分层模型**\ 下的理论地震图、反演面波频散和远震 P 波接收函数、
 反演震源深度、震源机制和地震矩、处理地震波形数据以及绘图。
 整个程序包很好地集成在一起，所有程序都使用相同的地球速度模型格式、波形格式和图形格式。
 

source/fk/index.rst

@@ -0,0 +1,31 @@
+fk
+==
+
+:软件名称: fk
+:主页: http://www.eas.slu.edu/People/LZhu/home.html
+:源码地址: http://www.eas.slu.edu/People/LZhu/home.html
+:作者: `Lupei Zhu <http://www.eas.slu.edu/People/LZhu/home.html>`__
+:最新版本: v3.3 (2015-02-03)
+:适用平台: Linux、macOS
+:编程语言: Fortran、C、Perl
+:许可协议: MIT
+
+----
+
+:教程贡献者: `田冬冬 <https://me.seisman.info/>`__ (作者)、`姚家园 <https://github.com/core-man>`__ (审稿)
+:教程更新日期: 2021-01-04
+
+----
+
+.. note::
+
+   fk v3.3 发布于 2015 年，其之后程序中的一些小问题被修复，但并未更新版本号。
+   因而 fk v3.3 实际上有多个不同的小版本。
+
+   本教程仅在 Linux 和 macOS 平台下使用 v3.3 （下载于 2021-01-04）测试通过，
+   其余版本或平台可能有所差异。
+
+.. toctree::
+   :hidden:
+
+   introduction

source/fk/introduction.rst

@@ -0,0 +1,63 @@
+简介
+====
+
+fk 可以用于计算\ **水平分层模型**\ 下的理论格林函数并合成理论地震图。
+
+功能
+----
+
+1. 支持计算爆炸源、双力偶和单力源下的格林函数和理论地震图
+2. 支持计算静态格林函数
+3. 理论上可以计算至任意高频
+
+命令及脚本
+----------
+
+fk 中提供了如下几个命令和脚本。通常只需要调用 ``fk.pl`` 生成格林函数，再
+调用 ``syn`` 将格林函数合成为三分量理论地震图：
+
+- ``fk``：用于计算格林函数的主程序，也是 fk 的核心程序
+- ``st_fk``：用于计算静态格林函数的主程序
+- ``fk.pl``：对 ``fk`` 和 ``st_fk`` 的封装，一般情况下直接使用该脚本即可
+- ``syn``：用于将格林函数合成为三分量理论地震图的程序
+- ``fk2mt``：将 fk 生成的格林函数转换为地震矩张量的每个分量所对应的格林函数
+- ``trav``：用于计算 P、S 初至到时的辅助程序
+- ``sachd``：用于修改 SAC 头段的辅助程序
+
+参考文献
+--------
+
+若想要了解 fk 的基本原理，可以阅读如下文章：
+
+- Haskell, N. A. (1964).
+  Radiation pattern of surface waves from point sources in a multi-layered medium.
+  *Bulletin of the Seismological Society of America*, 54(1), 377–393.
+- Takeuchi, H., & Saito, M. (1972).
+  Seismic surface waves.
+  *Methods in computational physics*, 11, 217-295.
+- Wang, C. Y., & Herrmann, R. B. (1980).
+  A numerical study of P-, SV-, and SH-wave generation in a plane layered medium.
+  *Bulletin of the Seismological Society of America*, 70(4), 1015–1036.
+- Zhu, L., & Rivera, L. A. (2002).
+  A note on the dynamic and static displacements from a point source in multilayered media.
+  *Geophysical Journal International*, 148(3), 619–627.
+  https://doi.org/10.1046/j.1365-246X.2002.01610.x
+
+阅读的相关建议：
+
+1. 若想了解如何使用 fk，直接阅读本教程即可
+2. 若想了解 fk 中每个选项的物理含义，需阅读 Zhu and Rivera (2002)
+3. 若想理解代码的实现细节，则需要在 Zhu and Rivera (2002) 的基础上，阅读其余三篇文章，
+   至少要阅读 Haskell (1964)
+4. 需要注意，尽管这几篇文献介绍的是同一种方法，但很多符号的定义是有区别的。
+   在推导代码中的公式时应以 Zhu and Rivera (2002) 为准。
+   Zhu and Rivera (2002) 区别于前面其他文献的地方主要在于，其重新定义了传播矩阵，
+   并将静态解与动态解统一到同一个公式中
+
+引用
+----
+
+- Zhu, L., & Rivera, L. A. (2002).
+  A note on the dynamic and static displacements from a point source in multilayered media.
+  *Geophysical Journal International*, 148(3), 619–627.
+  https://doi.org/10.1046/j.1365-246X.2002.01610.x

source/index.rst

@@ -44,12 +44,20 @@
 
    - `GMT <https://docs.gmt-china.org/latest/>`__
 
+.. rubric:: 合成地震图：
+
+.. hlist::
+   :columns: 1
+
+   - :doc:`fk/index`
+
 .. toctree::
    :maxdepth: 3
    :caption: 地震学软件
    :hidden:
 
    cps/index
+   fk/index
    HinetPy <https://seisman.github.io/HinetPy/zh_CN/index.html>
    GMT <https://docs.gmt-china.org/latest/>
    obspy/index