Skip to content

DOC: add grt.travt in Tutorial #21

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
merged 1 commit into from
Apr 20, 2025
Merged

DOC: add grt.travt in Tutorial #21

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
20 changes: 16 additions & 4 deletions docs/source/Tutorial/dynamic/gfunc.rst
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -43,16 +43,28 @@ Python中计算动态格林函数的主函数为 :func:`compute_grn() <pygrt.pym

一些基本信息（包括源点和场点的物性参数）保存在SAC头段变量中，其中 :c:var:`t0` 和 :c:var:`t1` 分别代表初至P波和初至S波的到时，由于程序中使用0作为参考时间，故其等价于走时。

如果你没有安装SAC软件，可以使用Python的ObsPy库读取数据，或者使用 :command:`grt.b2a` 工具临时将SAC格式文件转为如下的文本文件：
+ :command:`grt.travt` 命令可以显式地再计算初至波走时，

.. literalinclude:: run/run.sh
.. literalinclude:: run/run.sh
:language: bash
:start-after: BEGIN grt.travt
:end-before: END grt.travt

走时结果会输出到终端，

.. literalinclude:: run/travt
:language: text

+ 如果你没有安装SAC软件，可以使用Python的ObsPy库读取生成的SAC数据，或者使用 :command:`grt.b2a` 工具临时将SAC格式文件转为如下的文本文件：

.. literalinclude:: run/run.sh
:language: bash
:start-after: BEGIN grt.b2a
:end-before: END grt.b2a

输出的文本文件如下，两列分别为时间点和幅值。这种输出仅保留波形信息，缺失SAC文件中的头段变量。
输出的文本文件如下，两列分别为时间点和幅值。这种输出仅保留波形信息，缺失SAC文件中的头段变量。

.. literalinclude:: run/HFZ_head
.. literalinclude:: run/HFZ_head
:language: text

.. group-tab:: Python
Expand Down
8 changes: 8 additions & 0 deletions docs/source/Tutorial/dynamic/run/run.sh
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -8,6 +8,14 @@ grt -Mmilrow -D2/0 -N500/0.02 -OGRN -R5,8,10
# END GRN
# -----------------------------------------------------------------------------------

# -----------------------------------------------------------------------------------
# BEGIN grt.travt
grt.travt -Mmilrow -D2/0 -R5,8,10
# END grt.travt
# -----------------------------------------------------------------------------------
grt.travt -Mmilrow -D2/0 -R5,8,10 > travt


# -----------------------------------------------------------------------------------
# BEGIN grt.b2a
grt.b2a GRN/milrow_2_0_10/HFZ.sac > HFZ
Expand Down