DOC: update command description

Dengda98 · Dengda98 · commit 46e1e691da99 · 2025-08-05T11:42:58.000+08:00
diff --git a/docs/source/API/c_api.rst b/docs/source/API/c_api.rst
@@ -11,5 +11,7 @@ C extension API
 
    C_extension/common
    C_extension/dynamic
+   C_extension/grt.h
+   C_extension/grt_error.h
    C_extension/static
    C_extension/travt
diff --git a/docs/source/Advanced/filon/linear_filon.rst b/docs/source/Advanced/filon/linear_filon.rst
@@ -1,5 +1,3 @@
-.. _linear_filon:
-
 固定间隔的Filon积分法
 =========================
 
@@ -10,7 +8,7 @@
 
 .. warning:: 
 
-    谨慎使用固定间隔的Filon积分法，除非你很清楚原理以及误差来源。推荐使用 :ref:`self_adaptive_filon` 以减少固定积分间隔的影响。 
+    谨慎使用固定间隔的Filon积分法，除非你很清楚原理以及误差来源。推荐使用 :doc:`/Advanced/filon/self_adaptive_filon` 以减少固定积分间隔的影响。 
 
 
 以下介绍程序中使用的 **基于两点线性插值的Filon积分** :ref:`(纪晨, 姚振兴, 1995) <jichen_1995>`  :ref:`(初稿) <yao_init_manuscripts>` ，主要介绍思路，具体公式推导详见对应论文。
@@ -27,7 +25,7 @@
 
     P_m(\omega,r) = \int_0^{k^*} F_m(k, \omega)J_m(kr)kdk + \int_{k^*}^{k_{\text{max}}} F_m(k, \omega)J_m(kr)kdk 
 
-其中积分第一部分仍然使用离散波数积分求解 （见 :ref:`k_integ_rst` 部分），积分第二部分将Bessel函数取渐近表达式，可将其转为以下形式
+其中积分第一部分仍然使用离散波数积分求解 （见 :doc:`/Advanced/k_integ` 部分），积分第二部分将Bessel函数取渐近表达式，可将其转为以下形式
 
 .. math:: 
 
@@ -62,11 +60,11 @@
 
     .. group-tab:: C 
 
-        :command:`grt` 和 :command:`stgrt` 程序支持以下可选参数来使用Filon积分，具体说明详见 :command:`grt -h` 或 :command:`stgrt -h`。
+        :command:`greenfn` 和 :command:`static_greenfn` 模块支持以下可选参数来使用Filon积分，具体说明详见 :command:`grt greenfn -h` 或 :command:`grt static greenfn -h`。
 
         + ``-L<length>[/<Flength>/<Fcut>]``
          
-          + ``<length>``  定义离散波数积分的积分间隔 （见 :ref:`k_integ_rst` 部分）
+          + ``<length>``  定义离散波数积分的积分间隔 （见 :doc:`/Advanced/k_integ` 部分）
           + ``<Flength>`` 定义Filon积分间隔，公式和离散波数积分使用的一致。
           + ``<Fcut>`` 定义了两个积分的分割点， :math:`k^*=` ``<Fcut>`` :math:`/r_{\text{max}}`
          
diff --git a/docs/source/Advanced/filon/self_adaptive_filon.rst b/docs/source/Advanced/filon/self_adaptive_filon.rst
@@ -1,6 +1,4 @@
 
-.. _self_adaptive_filon:
-
 自适应Filon积分法
 ===================
 
@@ -9,7 +7,7 @@
 
 -----------------------------------------------------------
 
-在 :ref:`linear_filon` 部分简单介绍了Filon积分的过程与核心原理，即 **分段计算积分的解析解**，而以下介绍的 **自适应Filon积分法（SAFIM）** ，便是 **动态分段**，在核函数变化剧烈的部分多采样，在核函数变化平缓的部分少采样。具体原理详见 :ref:`(Chen and Zhang, 2001) <chen_2001>` :ref:`(张海明, 2021) <zhang_book_2021>` 。
+在 :doc:`/Advanced/filon/linear_filon` 部分简单介绍了Filon积分的过程与核心原理，即 **分段计算积分的解析解**，而以下介绍的 **自适应Filon积分法（SAFIM）** ，便是 **动态分段**，在核函数变化剧烈的部分多采样，在核函数变化平缓的部分少采样。具体原理详见 :ref:`(Chen and Zhang, 2001) <chen_2001>` :ref:`(张海明, 2021) <zhang_book_2021>` 。
 
 
 参数介绍
@@ -21,13 +19,13 @@
 
     .. group-tab:: C 
 
-        :command:`grt` 和 :command:`stgrt` 程序支持以下可选参数来使用自适应Filon积分，具体说明详见 :command:`grt -h` 或 :command:`stgrt -h`。
+        :command:`greenfn` 和 :command:`static_greenfn` 模块支持以下可选参数来使用自适应Filon积分，具体说明详见 :command:`grt greenfn -h` 或 :command:`grt static greenfn -h`。
 
         + ``-La<length>[/<Ftol>/<Fcut>]``
          
-          + ``<length>``  定义离散波数积分的积分间隔 （见 :ref:`k_integ_rst` 部分, :ref:`linear_filon` 部分）
+          + ``<length>``  定义离散波数积分的积分间隔 （见 :doc:`/Advanced/k_integ` 部分, :doc:`/Advanced/filon/linear_filon` 部分）
           + ``<Ftol>`` 定义自适应采样精度，见 :ref:`(Chen and Zhang, 2001) <chen_2001>`  :ref:`(张海明, 2021) <zhang_book_2021>`，通常1e-2即可。
-          + ``<Fcut>`` 定义了两个积分的分割点， :math:`k^*=` ``<Fcut>`` :math:`/r_{\text{max}}` （见 :ref:`linear_filon` 部分）
+          + ``<Fcut>`` 定义了两个积分的分割点， :math:`k^*=` ``<Fcut>`` :math:`/r_{\text{max}}` （见 :doc:`/Advanced/filon/linear_filon` 部分）
          
     .. group-tab:: Python
 
diff --git a/docs/source/Advanced/integ_converg/integ_converg.rst b/docs/source/Advanced/integ_converg/integ_converg.rst
@@ -1,5 +1,3 @@
-.. _integ_converg_rst:
-
 积分收敛性与峰谷平均法
 ===============================
 
@@ -8,7 +6,7 @@
 
 -----------------------------------------------------------
 
-通过输出核函数文件，观察源点和场点深度接近时，积分收敛性的变化，以及峰谷平均法的作用。 **具体积分表达式以及分类详见** :ref:`gfunc_rst`。
+通过输出核函数文件，观察源点和场点深度接近时，积分收敛性的变化，以及峰谷平均法的作用。 **具体积分表达式以及分类详见** :doc:`/Tutorial/dynamic/gfunc`。
 
 核函数文件
 ---------------
@@ -41,7 +39,7 @@ C和Python导出的核函数文件是一致的，底层调用的是相同的函
 
     .. group-tab:: C 
 
-        :command:`grt.k2a` 程序可将单个核函数文件转为文本格式。
+        :command:`k2a` 模块可将单个核函数文件转为文本格式。
 
         .. literalinclude:: run/run.sh
             :language: bash
@@ -62,15 +60,15 @@ C和Python导出的核函数文件是一致的，底层调用的是相同的函
             :start-after: BEGIN read statsfile
             :end-before: END read statsfile
 
-其中除了波数 ``k`` 外，每条结果的命名格式均为 ``{srcType}_{q/w/v}``，与 :ref:`gfunc_rst` 部分介绍的积分公式中的核函数 :math:`q_m, w_m, v_m` 保持一致。
+其中除了波数 ``k`` 外，每条结果的命名格式均为 ``{srcType}_{q/w/v}``，与 :doc:`/Tutorial/dynamic/gfunc` 部分介绍的积分公式中的核函数 :math:`q_m, w_m, v_m` 保持一致。
 
 .. note:: 
 
     核函数文件中记录的值非最终核函数值。对于动态解，还需乘 :math:`\left(-\dfrac{\Delta k}{4\pi\rho\omega^2}\right)`。
 
 可视化
 -------------
-以下将使用Python进行图件绘制。 **在Python函数中指定震源类型、阶数、积分类型，可自动绘制核函数、被积函数和积分值随波数的变化**，其中积分类型对应 :ref:`gfunc_rst` 部分介绍的4种类型。
+以下将使用Python进行图件绘制。 **在Python函数中指定震源类型、阶数、积分类型，可自动绘制核函数、被积函数和积分值随波数的变化**，其中积分类型对应 :doc:`/Tutorial/dynamic/gfunc` 部分介绍的4种类型。
 
 .. literalinclude:: run/run.py
     :language: python
@@ -145,7 +143,7 @@ C和Python导出的核函数文件是一致的，底层调用的是相同的函
 
 .. note:: 
 
-    :command:`grt.k2a` 程序也支持将 ``PTAM_{ir}_{dist}/PTAM_{iw}_{freq}`` 文件转为文本格式，
+    :command:`k2a` 模块也支持将 ``PTAM_{ir}_{dist}/PTAM_{iw}_{freq}`` 文件转为文本格式，
 
     .. literalinclude:: run/run.sh
         :language: bash
diff --git a/docs/source/Advanced/k_integ.rst b/docs/source/Advanced/k_integ.rst
@@ -1,5 +1,3 @@
-.. _k_integ_rst:
-
 控制波数积分
 ===================
 
@@ -8,13 +6,13 @@
 
 -----------------------------------------------------------
 
-在 :ref:`gfunc_rst` 部分介绍了程序将按以下形式从数值上计算波数积分，
+在 :doc:`/Tutorial/dynamic/gfunc` 部分介绍了程序将按以下形式从数值上计算波数积分，
 
 .. math:: 
 
    P_m(\omega) = \Delta k \sum_{j=0}^{\infty} F_m(k_j,\omega)J_m(k_j r)k_j
 
-其中 :math:`\Delta k = 2\pi/L, k_j=j\Delta k`，:math:`L` 为特征长度，即 :math:`L` 控制波数积分的积分间隔，默认根据 :ref:`gfunc_rst` 部分介绍的约束条件自动确定。
+其中 :math:`\Delta k = 2\pi/L, k_j=j\Delta k`，:math:`L` 为特征长度，即 :math:`L` 控制波数积分的积分间隔，默认根据 :doc:`/Tutorial/dynamic/gfunc` 部分介绍的约束条件自动确定。
 
 程序会取一个较大值 :math:`k_{\text{max}}` 作为波数积分的上限，经验公式为
 
@@ -46,7 +44,7 @@
 
     .. group-tab:: C 
 
-        :command:`grt` 程序支持以下可选参数来控制波数积分，具体说明详见 :command:`grt -h`。  
+        :command:`greenfn` 模块支持以下可选参数来控制波数积分，具体说明详见 :command:`grt greenfn -h`。  
 
         + ``-K<k0>[/<ampk>/<keps>]``  
          
@@ -58,7 +56,7 @@
 
         + ``-L<length>``， 对应公式中的 :math:`L` 的系数, :math:`L=` ``<length>`` :math:`\cdot r_{\text{max}}`
         
-        :command:`stgrt` 程序支持以下可选参数来控制波数积分，参数与上面对应，具体说明详见 :command:`stgrt -h`。
+        :command:`static_greenfn` 模块支持以下可选参数来控制波数积分，参数与上面对应，具体说明详见 :command:`grt static greenfn -h`。
 
         + ``-K<k0>[/<keps>]``  
         + ``-L<length>``
diff --git a/docs/source/Advanced/kernel/kernel.rst b/docs/source/Advanced/kernel/kernel.rst
@@ -7,7 +7,7 @@
 
 -----------------------------------------------------------
 
-在 :ref:`integ_converg_rst` 部分介绍了程序可以导出积分过程中不同频率的核函数值，这使得我们可以绘制核函数的 :math:`f-v` 谱图 [#]_ ，以观察其中的频散特征 。
+在 :doc:`/Advanced/integ_converg/integ_converg` 部分介绍了程序可以导出积分过程中不同频率的核函数值，这使得我们可以绘制核函数的 :math:`f-v` 谱图 [#]_ ，以观察其中的频散特征 。
 
 以薄层模型为例，
 
diff --git a/docs/source/Formula/static_source.rst b/docs/source/Formula/static_source.rst
@@ -95,7 +95,7 @@
 
 其中三个势函数均吸收了 :math:`k` 因子，即 :math:`\phi_m \leftarrow k\phi_m, \psi_m \leftarrow k\psi_m, \chi_m \leftarrow k\chi_m` ，这会体现在后续的震源系数中。
 
-为了求解震源系数，需将 :ref:`static_uniform` 部分的最终表达式展开成 :eq:`u_qwv` 的形式，对应的垂直波函数的系数 :math:`(B_m, D_m, F_m)` 即为震源系数。具体而言，我们可将位移同一分量的不同表达式进行对比得到震源系数。将 :eq:`qwv_pot` 式代入 :eq:`u_qwv` 式，柱坐标系下的位移三分量表达式为
+为了求解震源系数，需将 :doc:`static_uniform` 部分的最终表达式展开成 :eq:`u_qwv` 的形式，对应的垂直波函数的系数 :math:`(B_m, D_m, F_m)` 即为震源系数。具体而言，我们可将位移同一分量的不同表达式进行对比得到震源系数。将 :eq:`qwv_pot` 式代入 :eq:`u_qwv` 式，柱坐标系下的位移三分量表达式为
 
 .. math::
     :label:
diff --git a/docs/source/Formula/static_uniform.rst b/docs/source/Formula/static_uniform.rst
@@ -1,5 +1,3 @@
-.. _static_uniform:
-
 无限均匀介质中集中脉冲力产生的静态位移场
 ===========================================
 
diff --git a/docs/source/Tutorial/dynamic/gfunc.rst b/docs/source/Tutorial/dynamic/gfunc.rst
@@ -1,5 +1,3 @@
-.. _gfunc_rst:
-
 计算动态格林函数
 =================
 
@@ -9,9 +7,9 @@
 -----------------------------------------------------------
 
 
-Python中计算动态格林函数的主函数为 :func:`compute_grn() <pygrt.pymod.PyModel1D.compute_grn>` ，C程序为 :command:`grt`。
+Python中计算动态格林函数的主函数为 :func:`compute_grn() <pygrt.pymod.PyModel1D.compute_grn>` ，C模块为 :command:`greenfn`。
 
-核心计算逻辑来自  :ref:`初稿 <yao_init_manuscripts>` ，具体代码可见与C API中对应同名 :file:`*.c` 文件，其中计算格林函数频谱的主函数为 :file:`grt.c` 里的 :c:func:`integ_grn_spec`。输出结果的坐标系见下图。
+核心计算逻辑来自  :ref:`初稿 <yao_init_manuscripts>` ，具体代码可见与C API中对应同名 :file:`*.c` 文件，其中计算格林函数频谱的主函数为 :file:`grn.c` 里的 :c:func:`integ_grn_spec`。输出结果的坐标系见下图。
 
 .. image:: ../coord.svg
    :align: center
@@ -43,7 +41,7 @@ Python中计算动态格林函数的主函数为 :func:`compute_grn() <pygrt.pym
 
         一些基本信息（包括源点和场点的物性参数）保存在SAC头段变量中，其中 :c:var:`t0` 和 :c:var:`t1` 分别代表初至P波和初至S波的到时。在不设定其它参数时，程序中使用0（发震时刻）作为参考时间，故其等价于走时。
 
-        + :command:`grt.travt` 命令可以显式地再计算初至波走时，
+        + :command:`travt` 模块可以显式地再计算初至波走时，
 
          .. literalinclude:: run/run.sh
             :language: bash
@@ -55,7 +53,7 @@ Python中计算动态格林函数的主函数为 :func:`compute_grn() <pygrt.pym
          .. literalinclude:: run/travt
             :language: text 
 
-        + 如果你没有安装SAC软件，可以使用Python的ObsPy库读取生成的SAC数据，或者使用 :command:`grt.b2a` 工具临时将SAC格式文件转为如下的文本文件：
+        + 如果你没有安装SAC软件，可以使用Python的ObsPy库读取生成的SAC数据，或者使用 :command:`b2a` 模块临时将SAC格式文件转为如下的文本文件：
 
          .. literalinclude:: run/run.sh
             :language: bash
diff --git a/docs/source/Tutorial/dynamic/strain_stress.rst b/docs/source/Tutorial/dynamic/strain_stress.rst
@@ -1,5 +1,3 @@
-.. _strain_stress_rst:
-
 计算动态应变、旋转、应力张量
 =======================================
 
@@ -9,7 +7,7 @@
 -----------------------------------------------------------
 
 
-计算这些张量需要计算位移的空间导数 :math:`\partial z,\partial r,\partial \theta` [#partial]_ 。 从前文关于格林函数的计算 （ :ref:`gfunc_rst` ） 和位移合成的公式 （ :ref:`syn_rst` ） 出发，位移空间导数分两个步骤计算，再根据几何方程和本构方程合成应变、旋转和应力张量。
+计算这些张量需要计算位移的空间导数 :math:`\partial z,\partial r,\partial \theta` [#partial]_ 。 从前文关于格林函数的计算 （ :doc:`/Tutorial/dynamic/gfunc` ） 和位移合成的公式 （ :doc:`/Tutorial/dynamic/syn` ） 出发，位移空间导数分两个步骤计算，再根据几何方程和本构方程合成应变、旋转和应力张量。
 
 .. [#partial] :math:`\partial z,\partial r,\partial \theta` 为 :math:`\dfrac{\partial}{\partial z},\dfrac{\partial}{\partial r},\dfrac{\partial}{\partial \theta}` 的简写。
 
diff --git a/docs/source/Tutorial/dynamic/syn.rst b/docs/source/Tutorial/dynamic/syn.rst
@@ -1,5 +1,3 @@
-.. _syn_rst:
-
 合成动态位移
 =================
 
@@ -13,7 +11,7 @@
     **震源机制参数中（如单力源、矩张量源）Z轴取向下为正。**
 
 
-Python中合成动态位移的主函数为 :func:`gen_syn_from_gf_*() <pygrt.utils.gen_syn_from_gf_DC>` （\*表示对不同震源）  ，C程序为 :command:`grt.syn`。
+Python中合成动态位移的主函数为 :func:`gen_syn_from_gf_*() <pygrt.utils.gen_syn_from_gf_DC>` （\*表示对不同震源）  ，C模块为 :command:`syn`。
 
 使用上节计算的格林函数，合成动态位移（理论地震图）。方便起见，这里统一使用milrow模型，震源深度2km，场点位于地表，震中距10km的格林函数，方位角30°。
 
@@ -195,7 +193,7 @@ Python中合成动态位移的主函数为 :func:`gen_syn_from_gf_*() <pygrt.uti
             :start-after: BEGIN TIME FUNC
             :end-before: END TIME FUNC
 
-        生成的时间函数会以SAC格式保存在对应路径中，文件名为 :file:`sig.sac`。 其它时间函数以及具体参数用法可使用 :command:`grt.syn -h` 查看说明。
+        生成的时间函数会以SAC格式保存在对应路径中，文件名为 :file:`sig.sac`。 其它时间函数以及具体参数用法可使用 :command:`grt syn -h` 查看说明。
 
     .. group-tab:: Python 
 
diff --git a/docs/source/Tutorial/index.rst b/docs/source/Tutorial/index.rst
@@ -6,9 +6,18 @@
 
 -----------------------------------------------------------
 
-**PyGRT** 程序包由C和Python两个编程语言的代码组成，目的是兼并高效性和便捷性。底层复杂运算由C语言编写，编译链接成动态库 ``libgrt.so`` 供Python调用。Python通过 ``ctypes`` 库导入动态库以使用外部函数，以此兼并了C语言的高效和Python语言的便捷。
+**PyGRT** 程序包由C和Python两个编程语言的代码组成，目的是兼并高效性和便捷性。
+底层复杂运算由C语言编写，编译链接成动态库 ``libgrt.so`` 供Python调用。Python通过 ``ctypes`` 库导入动态库以使用外部函数，以此兼并了C语言的高效和Python语言的便捷。
 
-除了Python脚本式运行， **PyGRT** 保留传统命令行式运行C程序，即在建立动态库过程中也编译了一些可执行文件。 **C程序的运行独立于Python，不需要Python环境，以满足更多计算场景。** 每个C程序可使用 ``-h`` 查看帮助。
+除了Python脚本式运行， **PyGRT** 保留传统命令行式运行C程序 :command:`grt` 。
+受 `GMT <https://www.generic-mapping-tools.org/>`_ 的启发， :command:`grt` 程序对多个计算功能使用模块化管理，
+可使用类似于以下格式来运行程序。每个模块可使用 ``-h`` 查看帮助。
+
+.. code-block:: C
+
+    grt <module-name> [<module-options>]
+
+**C程序的运行独立于Python，不需要Python环境，从而满足了更多计算场景。** 
 
 计算 **动态解** 的主要计算流程如下：
 
@@ -42,18 +51,18 @@
         classDef cmdcls fill:#f9f2d9,stroke:#e8d174,stroke-width:2px,color:#333;
         class GG,SS,EE,RR,TT cmdcls
 
-+ C
++ C (module name)
 
 .. mermaid::
     :zoom:
 
     flowchart TB 
 
-        GG(["grt"])
-        SS(["grt.syn"])
-        EE(["grt.strain"])
-        RR(["grt.rotation"])
-        TT(["grt.stress"])
+        GG(["greenfn"])
+        SS(["syn"])
+        EE(["strain"])
+        RR(["rotation"])
+        TT(["stress"])
 
         G["Compute Green's Functions
         (and its Spatial Derivatives)"]
@@ -72,8 +81,18 @@
         classDef cmdcls fill:#f9f2d9,stroke:#e8d174,stroke-width:2px,color:#333;
         class GG,SS,EE,RR,TT cmdcls
 
-计算 **静态解** 的过程相同，只是函数名称和程序名称有变化。
+计算 **静态解** 的过程相同，只是在名称上有如下区别：
+
++ Python 函数名称增加了 "static_" ，例如 :func:`compute_static_grn() <pygrt.pymod.PyModel1D.compute_static_grn>` 。
++ C 程序模块名增加 "static_" 前缀，例如 :command:`static_greenfn` ；或者在 :command:`grt` 与无 "static_" 前缀的模块名之间增加 "static" 命令，
+  即支持以下两种方式（以 :command:`greenfn` 模块为例）：
+
+  .. code-block:: bash
+
+    grt static_greenfn [<module-options>]
+    grt static greenfn [<module-options>]
 
+-----------------
 
 以下是一些简易教程，可快速上手。Github主页的 :rst:dir:`example/` 文件夹中有更多示例，可在 `Github Releases <https://github.com/Dengda98/PyGRT/releases>`_ 中下载。
 
diff --git a/docs/source/Tutorial/prepare.rst b/docs/source/Tutorial/prepare.rst
@@ -10,7 +10,7 @@
 安装 **PyGRT**
 --------------------
 
-详见 :ref:`install_section` 。
+详见 :doc:`/install` 。
 
 
 建立模型文件
@@ -19,7 +19,7 @@
 .. image:: mod.svg
    :align: center
 
-这里的选项卡 **C** 不代表内部是C语言代码，而是 **使用C语言编译出的可执行程序** （与之相关）。后续的选项卡沿用此设定，不再解释。
+这里的选项卡 **C** 不代表内部是C语言代码，而是 **使用C程序** :command:`grt` 。后续的选项卡沿用此设定，不再解释。
 
 .. tabs:: 
 
diff --git a/docs/source/Tutorial/static/static_gfunc.rst b/docs/source/Tutorial/static/static_gfunc.rst
diff --git a/docs/source/Tutorial/static/static_strain_stress.rst b/docs/source/Tutorial/static/static_strain_stress.rst
diff --git a/docs/source/Tutorial/static/static_syn.rst b/docs/source/Tutorial/static/static_syn.rst
diff --git a/docs/source/index.rst b/docs/source/index.rst
diff --git a/docs/source/install.rst b/docs/source/install.rst

Original file line number	Diff line number	Diff line change
`@@ -1,5 +1,3 @@`
`1`		`-.. _static_uniform:`
`2`		`-`
`3`	`1`	`无限均匀介质中集中脉冲力产生的静态位移场`
`4`	`2`	`===========================================`
`5`	`3`