From 5abd1e655dceda921f60a6eb0356f5bd8a0a44b5 Mon Sep 17 00:00:00 2001 From: Dongdong Tian Date: Fri, 1 Jan 2021 02:06:57 -0500 Subject: [PATCH 01/24] taup: Add tutorial for taup_time --- source/taup/taup_time.rst | 108 ++++++++++++++++++++++++++++++++++++++ 1 file changed, 108 insertions(+) create mode 100644 source/taup/taup_time.rst diff --git a/source/taup/taup_time.rst b/source/taup/taup_time.rst new file mode 100644 index 00000000..18df8195 --- /dev/null +++ b/source/taup/taup_time.rst @@ -0,0 +1,108 @@ +taup_time +========= + +TauP 提供的 ``taup_time`` 命令可用于计算震相走时、射线参数、出射角、 +入射角等信息。 + +基本用法 +-------- + +计算震源深度为 300 km，震中距为 60°时，P、S、PcP 等震相的走时信息:: + + taup_time -mod prem -ph P,S,PcP,ScS,PKiKP -deg 60 -h 300 + +其中： + +- ``-mod`` 指定速度模型为 TauP 内置的 prem 模型； +- ``-ph`` 后跟一系列震相名； +- ``-h`` 震源深度； +- ``-deg`` 后接震中距，单位为度； + +在震中距比较小的情况下，用千米表示震中距更方便，可以使用 ``-km`` +选项。比如计算一个浅源地震在 100 km 处常见 P 波震相的走时信息:: + + taup_time -mod prem -ph ttp+ -h 2 -km 100 + +当然，也可以通过给定地震位置以及台站位置:: + + taup_time -mod prem -ph P,PcP -h 100 -sta 20 130 -evt 40 150 + +输出信息 +-------- + +下面的命令计算了P和pP的走时信息，以它为例:: + + $ taup_time -mod prem -h 100 -deg 60 -ph pP,P + + Model: prem + Distance Depth Phase Travel Ray Param Takeoff Incident Purist Purist + (deg) (km) Name Time (s) p (s/deg) (deg) (deg) Distance Name + ----------------------------------------------------------------------------------- + 60.00 100.0 P 595.40 6.815 30.14 20.82 60.00 = P + 60.00 100.0 pP 618.89 6.890 149.49 21.06 60.00 = pP + +输出共 9 列，从左到右分别为： + +1. 震中距（度）； +2. 震源深度（km）； +3. 震相名； +4. 震相走时（秒）； +5. 射线参数（秒 / 度）； +6. 出射角（deg）：即射线与垂直向下方向的夹角，取值范围为0到180度。 + 夹角在 0 到 90 度内表示射线向下出射（例如P震相）； 夹角 90 度到 180 + 度内表示射线向上出射（例如pP震相)。 +7. 入射角（deg）：射线在台站（位于地表）处入射时与垂直向上方向的夹角。应该是恒正的。 +8. Purist Distance 和 Purist Name 不知为何物； + +关于射线参数 + +.. math:: + + p = r \frac{\sin \theta}{v_s} = 6371 km \frac{\sin 20.82}{5.8 km/s} + = 390.425 s/radian + +.. math:: = 390.425 * 2\pi/360 s/deg = 6.815 s/deg + +**老版本的TauP对出射角的定义稍有不同： 若夹角在 0 到 90 +度内，则取正值；若夹角在 90 度到 180 度内，则先减去 90，再取负值。 +因而，一般情况下，P 波出射者为正值，而 pP 出射角为负值。** + +信息提取 +-------- + +上面介绍的输出中包含了很多信息，想要在脚本中提取出想要的信息有些麻烦。所以 +``taup_time`` 提供了一些选项，使得在脚本中计算走时什么的更加方便。 + +time +~~~~ + +该选项使得输出中只包含走时信息:: + + $ taup_time -mod prem -ph P,S,PcP -deg 60 -h 100 --time + 595.3896 640.09875 1081.2472 + +需要注意的是，输出中走时不是按照 ``-ph`` +选项中震相的顺序进行排序的，而是按照走时递增的顺序排序。比如这里，S +震相放在震相列表的第二位，而走时 1081.2472 却放在第三位。 + +rayp +~~~~ + +该选项使得输出中只包含射线参数信息:: + + $ taup_time -mod prem -ph P,S,PcP -deg 60 -h 100 --rayp + 6.8185554 3.9990747 12.785722 + +这里，射线参数的输出是按照震相进行排序的。 + +time 和 rayp +~~~~~~~~~~~~ + +这两个选项是不能一起用的，也就是说不能同时输出震相走时和射线参数信息。当两个选项同时使用时，会以后出现的选项为准。 + +当同时需要走时和射线参数时，只能多次执行命令了。 + +rel +~~~ + +该选项用于计算震相的相对走时，似乎有 bug，无法正常工作。 From 3df3f2cda6d2a7a4324fa6d1cbe8b08112c83997 Mon Sep 17 00:00:00 2001 From: Dongdong Tian Date: Fri, 1 Jan 2021 17:31:43 -0500 Subject: [PATCH 02/24] Add to toc --- source/taup/index.rst | 1 + 1 file changed, 1 insertion(+) diff --git a/source/taup/index.rst b/source/taup/index.rst index aabdca58..dc58d421 100644 --- a/source/taup/index.rst +++ b/source/taup/index.rst @@ -27,3 +27,4 @@ TauP introduction install + taup_time From 9d1f7ecb5f54e6be31dbfeeeffdd1c3f40df2e94 Mon Sep 17 00:00:00 2001 From: Dongdong Tian Date: Fri, 1 Jan 2021 19:23:53 -0500 Subject: [PATCH 03/24] Updates --- source/taup/taup_time.rst | 67 +++++++++++++++++++++++++-------------- 1 file changed, 43 insertions(+), 24 deletions(-) diff --git a/source/taup/taup_time.rst b/source/taup/taup_time.rst index 18df8195..acbdb6b2 100644 --- a/source/taup/taup_time.rst +++ b/source/taup/taup_time.rst @@ -7,30 +7,39 @@ TauP 提供的 ``taup_time`` 命令可用于计算震相走时、射线参数、 基本用法 -------- -计算震源深度为 300 km，震中距为 60°时，P、S、PcP 等震相的走时信息:: +使用 ``taup_time --help`` 可查询该命令的完整用法和选项，此处仅介绍最常用用法。 - taup_time -mod prem -ph P,S,PcP,ScS,PKiKP -deg 60 -h 300 +计算震源深度为 300 km，震中距 60° 处，P、S、PcP 等震相的走时信息:: -其中： + $ taup_time -ph P,S,PcP,ScS,PKiKP -deg 60 -h 300 -- ``-mod`` 指定速度模型为 TauP 内置的 prem 模型； -- ``-ph`` 后跟一系列震相名； -- ``-h`` 震源深度； -- ``-deg`` 后接震中距，单位为度； +- ``-ph`` 指定了一系列以逗号分隔的震相名 +- ``-h`` 指定震源深度（单位 km） +- ``-deg`` 指定震中距（单位为°） -在震中距比较小的情况下，用千米表示震中距更方便，可以使用 ``-km`` -选项。比如计算一个浅源地震在 100 km 处常见 P 波震相的走时信息:: +TauP 默认使用 iasp91 模型，可以使用 ``-mod`` 选项设置要使用的速度模型:: - taup_time -mod prem -ph ttp+ -h 2 -km 100 + $ taup_time -ph P,S,PcP,ScS,PKiKP -deg 60 -h 300 -mod prem -当然，也可以通过给定地震位置以及台站位置:: +在震中距比较小的情况下，可以使用 ``-km`` 选项指定震中距（单位为千米）。 +比如计算一个浅源地震在 100 km 处常见 P 波震相的走时信息:: - taup_time -mod prem -ph P,PcP -h 100 -sta 20 130 -evt 40 150 + $ taup_time -mod prem -ph P,S -h 2 -km 100 -输出信息 +若已知震源和台站位置，也可以使用 ``-sta`` 和 ``-evt`` 分别指定台站和震源的 +纬度和经度。例如，台站位于 (20°N,130°E)，震源位于 (40°N,150°E)，则:: + + $ taup_time -mod prem -ph P,PcP -h 100 -sta 20 130 -evt 40 150 + +TauP 中默认台站位于地表处，也可使用 ``--stadepth`` 指定台站深度（单位为km）:: + + $ taup_time -ph P,S,PcP,ScS,PKiKP -deg 60 -h 300 --stadepth 0.5 + +输出格式 -------- -下面的命令计算了P和pP的走时信息，以它为例:: +下面的命令计算了 PREM 模型下，震源深度 100 km，震中距 60° 处，P 和 pP 震相的走时信息。 +其输出为:: $ taup_time -mod prem -h 100 -deg 60 -ph pP,P @@ -43,16 +52,26 @@ TauP 提供的 ``taup_time`` 命令可用于计算震相走时、射线参数、 输出共 9 列，从左到右分别为： -1. 震中距（度）； -2. 震源深度（km）； -3. 震相名； -4. 震相走时（秒）； -5. 射线参数（秒 / 度）； -6. 出射角（deg）：即射线与垂直向下方向的夹角，取值范围为0到180度。 - 夹角在 0 到 90 度内表示射线向下出射（例如P震相）； 夹角 90 度到 180 - 度内表示射线向上出射（例如pP震相)。 -7. 入射角（deg）：射线在台站（位于地表）处入射时与垂直向上方向的夹角。应该是恒正的。 -8. Purist Distance 和 Purist Name 不知为何物； +1. 震中距（度） +2. 震源深度（km） +3. 震相名 +4. 震相走时（秒） +5. 射线参数（秒/度） +6. 出射角（度）：即射线从震源出射时与\ **垂直向下**\ 方向的夹角，取值范围为0到180°。 + + - 出射角 0° 到 90° 内表示射线向下出射（例如P震相） + - 出射角 90° 到 180° 内表示射线向上出射（例如pP震相) + +7. 入射角（度）：射线在台站（位于地表）处入射时与垂直向上方向的夹角。应该是恒正的。 +8. Purist Distance +9. Purist Name + + +.. note:: + + 射线参数的单位 + + 关于射线参数 From d989102461819c2d440e60fc8db5e15bb1296e80 Mon Sep 17 00:00:00 2001 From: Dongdong Tian Date: Fri, 1 Jan 2021 22:47:25 -0500 Subject: [PATCH 04/24] Update --- source/taup/introduction.rst | 2 +- source/taup/taup_time.rst | 121 ++++++++++++++++++++++------------- 2 files changed, 76 insertions(+), 47 deletions(-) diff --git a/source/taup/introduction.rst b/source/taup/introduction.rst index e5b78542..25683464 100644 --- a/source/taup/introduction.rst +++ b/source/taup/introduction.rst @@ -29,7 +29,7 @@ TauP 提供了多个工具: - 图形界面： ``taup`` - 命令行工具： - - ``taup_time`` ：计算震相的走时、射线参数、出射角、入射角等 + - :doc:`taup_time` ：计算震相的走时、射线参数、出射角、入射角等 - ``taup_pierce`` ：计算震相在波速不连续面和指定深度的反射点和穿透点 - ``taup_path`` ：计算射线路径 - ``taup_wavefront`` ：计算波前面 diff --git a/source/taup/taup_time.rst b/source/taup/taup_time.rst index acbdb6b2..137be3a0 100644 --- a/source/taup/taup_time.rst +++ b/source/taup/taup_time.rst @@ -21,19 +21,19 @@ TauP 默认使用 iasp91 模型，可以使用 ``-mod`` 选项设置要使用的 $ taup_time -ph P,S,PcP,ScS,PKiKP -deg 60 -h 300 -mod prem -在震中距比较小的情况下，可以使用 ``-km`` 选项指定震中距（单位为千米）。 -比如计算一个浅源地震在 100 km 处常见 P 波震相的走时信息:: +TauP 中默认台站位于地表处，也可使用 ``--stadepth`` 指定台站深度（单位为km）:: - $ taup_time -mod prem -ph P,S -h 2 -km 100 + $ taup_time -ph P,S,PcP,ScS,PKiKP -deg 60 -h 300 --stadepth 0.5 若已知震源和台站位置，也可以使用 ``-sta`` 和 ``-evt`` 分别指定台站和震源的 纬度和经度。例如，台站位于 (20°N,130°E)，震源位于 (40°N,150°E)，则:: - $ taup_time -mod prem -ph P,PcP -h 100 -sta 20 130 -evt 40 150 + $ taup_time -mod prem -ph P,S,PcP,ScS,PKiKP -h 100 -sta 20 130 -evt 40 170 -TauP 中默认台站位于地表处，也可使用 ``--stadepth`` 指定台站深度（单位为km）:: +在震中距比较小的情况下，可以使用 ``-km`` 选项指定震中距（单位为千米）。 +比如计算一个浅源地震在 100 km 处常见 P 波震相的走时信息:: - $ taup_time -ph P,S,PcP,ScS,PKiKP -deg 60 -h 300 --stadepth 0.5 + $ taup_time -ph P,S,PcP,ScS,PKiKP -h 5 -km 100 输出格式 -------- @@ -52,39 +52,75 @@ TauP 中默认台站位于地表处，也可使用 ``--stadepth`` 指定台站 输出共 9 列，从左到右分别为： -1. 震中距（度） -2. 震源深度（km） -3. 震相名 -4. 震相走时（秒） -5. 射线参数（秒/度） -6. 出射角（度）：即射线从震源出射时与\ **垂直向下**\ 方向的夹角，取值范围为0到180°。 +1. 震中距（度）：地震震中与台站的大圆弧路径长度 +2. 震源深度（km） +3. 震相名 +4. 震相走时（秒） +5. 射线参数（秒/度） - - 出射角 0° 到 90° 内表示射线向下出射（例如P震相） - - 出射角 90° 到 180° 内表示射线向上出射（例如pP震相) + .. note:: -7. 入射角（度）：射线在台站（位于地表）处入射时与垂直向上方向的夹角。应该是恒正的。 -8. Purist Distance -9. Purist Name + 射线参数有三种常用单位：秒/度（s/deg ）、秒/弧度（s/radian）、秒/千米（s/km）。 + 三者之间的关系为（:math:`R_0` 为地球半径）： + .. math:: -.. note:: + 1 s/deg = \frac{1}{\pi*R_0/180} s/km = 1/111.2 s/km + + .. math:: + + 1 s/deg = \frac{1}{\pi / 180} s/radian = 57.3 s/radian + +6. 出射角（度）：即射线从震源出射时与\ **垂直向下**\ 方向的夹角，取值范围为0到180°。 - 射线参数的单位 + - 出射角 0° 到 90° 内表示射线向下出射（例如P震相） + - 出射角 90° 到 180° 内表示射线向上出射（例如pP震相) +7. 入射角（度）：射线在台站（位于地表）处入射时与垂直向上方向的夹角。应该是恒正的。 +8. Purist 距离：地震波实际走过的圆弧距离。该值可能与震中距不同，例如:: + $ taup_time -ph PKKKP -deg 50 -关于射线参数 + Model: iasp91 + Distance Depth Phase Travel Ray Param Takeoff Incident Purist Purist + (deg) (km) Name Time (s) p (s/deg) (deg) (deg) Distance Name + ----------------------------------------------------------------------------------- + 50.00 0.0 PKKKP 2519.56 2.242 6.72 6.72 410.00 = PKKKP -.. math:: +9. Purist 震相名：地震波的“真实”震相名。 - p = r \frac{\sin \theta}{v_s} = 6371 km \frac{\sin 20.82}{5.8 km/s} - = 390.425 s/radian + 例如，若想要计算 SS 震相在 660 km 不连续面的下表面的反射波走时，可指定其震相名为 + ``S^660S``。由于 PREM 模型中 660 km 不连续面的“真实”深度为670 km，此时 + TauP 实际返回的时从 670 km 处反射回的震相，此时 Purist 震相名会显示为 ``S^670S``:: -.. math:: = 390.425 * 2\pi/360 s/deg = 6.815 s/deg + $ taup_time -mod prem -ph S^660S -deg 120 -**老版本的TauP对出射角的定义稍有不同： 若夹角在 0 到 90 -度内，则取正值；若夹角在 90 度到 180 度内，则先减去 90，再取负值。 -因而，一般情况下，P 波出射者为正值，而 pP 出射角为负值。** + Model: prem + Distance Depth Phase Travel Ray Param Takeoff Incident Purist Purist + (deg) (km) Name Time (s) p (s/deg) (deg) (deg) Distance Name + ------------------------------------------------------------------------------------ + 120.00 0.0 S^660S 1975.75 12.336 20.79 20.79 120.00 * S^670S + +当使用 ``--rel`` 选项时，``taup_time`` 可输出相对于某个指定震相的走时差。例如， +下面的命令中，最后一列输出了所有震相相对于 PcP 震相的走时差:: + + $ taup_time -ph P,S,PcP,ScS,PKiKP -deg 60 -h 300 -mod prem --rel PcP + + Model: prem + Distance Depth Phase Travel Ray Param Takeoff Incident Purist Purist Relative to + (deg) (km) Name Time (s) p (s/deg) (deg) (deg) Distance Name PcP + ----------------------------------------------------------------------------------------------- + 60.00 300.0 P 574.72 6.739 33.65 20.58 60.00 = P -42.32 + 60.00 300.0 PcP 617.05 4.016 19.28 12.09 60.00 = PcP 0.00 + 60.00 300.0 PKiKP 993.45 1.245 5.88 3.72 60.00 = PKiKP 376.41 + 60.00 300.0 S 1043.98 12.659 34.16 21.37 60.00 = S 426.93 + 60.00 300.0 ScS 1134.89 7.480 19.38 12.43 60.00 = ScS 517.85 + +.. note:: + + 不同文章对出射角的定义不同，使用时应格外注意。 + + 老版本的TauP对出射角的定义不同，可视为BUG，应避免使用老版本TauP。 信息提取 -------- @@ -92,36 +128,29 @@ TauP 中默认台站位于地表处，也可使用 ``--stadepth`` 指定台站 上面介绍的输出中包含了很多信息，想要在脚本中提取出想要的信息有些麻烦。所以 ``taup_time`` 提供了一些选项，使得在脚本中计算走时什么的更加方便。 -time +走时 ~~~~ -该选项使得输出中只包含走时信息:: +``--time`` 选项使得输出中只包含走时信息:: $ taup_time -mod prem -ph P,S,PcP -deg 60 -h 100 --time 595.3896 640.09875 1081.2472 -需要注意的是，输出中走时不是按照 ``-ph`` -选项中震相的顺序进行排序的，而是按照走时递增的顺序排序。比如这里，S -震相放在震相列表的第二位，而走时 1081.2472 却放在第三位。 +需要注意的是，输出中走时不是按照 ``-ph`` 选项中震相的顺序进行排序的， +而是按照走时递增的顺序排序。比如这里，S 震相放在震相列表的第二位， +而 S 震相的走时 1081.2472 却放在第三位。 -rayp -~~~~ +射线参数 +~~~~~~~~ -该选项使得输出中只包含射线参数信息:: +``--rayp`` 选项使得输出中只包含射线参数信息:: $ taup_time -mod prem -ph P,S,PcP -deg 60 -h 100 --rayp 6.8185554 3.9990747 12.785722 -这里，射线参数的输出是按照震相进行排序的。 - -time 和 rayp -~~~~~~~~~~~~ +此时，射线参数的输出是按照震相进行排序的。 -这两个选项是不能一起用的，也就是说不能同时输出震相走时和射线参数信息。当两个选项同时使用时，会以后出现的选项为准。 - -当同时需要走时和射线参数时，只能多次执行命令了。 - -rel -~~~ +.. note:: -该选项用于计算震相的相对走时，似乎有 bug，无法正常工作。 + ``--rayp`` 和 ``--time`` 选项不能一起使用。 + 当两个选项同时使用时，会以后出现的选项为准。 From 5eb14550e0943b6544e8c24751561d84820b257b Mon Sep 17 00:00:00 2001 From: Dongdong Tian Date: Fri, 1 Jan 2021 23:17:58 -0500 Subject: [PATCH 05/24] fixes --- source/taup/taup_time.rst | 25 ++++++++++++++----------- 1 file changed, 14 insertions(+), 11 deletions(-) diff --git a/source/taup/taup_time.rst b/source/taup/taup_time.rst index 137be3a0..d0ea907e 100644 --- a/source/taup/taup_time.rst +++ b/source/taup/taup_time.rst @@ -7,7 +7,7 @@ TauP 提供的 ``taup_time`` 命令可用于计算震相走时、射线参数、 基本用法 -------- -使用 ``taup_time --help`` 可查询该命令的完整用法和选项，此处仅介绍最常用用法。 +使用 ``taup_time --help`` 可查询该命令的完整用法和选项。此处仅介绍最常用用法。 计算震源深度为 300 km，震中距 60° 处，P、S、PcP 等震相的走时信息:: @@ -21,7 +21,7 @@ TauP 默认使用 iasp91 模型，可以使用 ``-mod`` 选项设置要使用的 $ taup_time -ph P,S,PcP,ScS,PKiKP -deg 60 -h 300 -mod prem -TauP 中默认台站位于地表处，也可使用 ``--stadepth`` 指定台站深度（单位为km）:: +TauP 中默认台站位于地表处，也可使用 ``--stadepth`` 指定台站深度（单位为 km）:: $ taup_time -ph P,S,PcP,ScS,PKiKP -deg 60 -h 300 --stadepth 0.5 @@ -30,7 +30,7 @@ TauP 中默认台站位于地表处，也可使用 ``--stadepth`` 指定台站 $ taup_time -mod prem -ph P,S,PcP,ScS,PKiKP -h 100 -sta 20 130 -evt 40 170 -在震中距比较小的情况下，可以使用 ``-km`` 选项指定震中距（单位为千米）。 +在震中距比较小的情况下，可以使用 ``-km`` 选项指定震中距（单位为 km）。 比如计算一个浅源地震在 100 km 处常见 P 波震相的走时信息:: $ taup_time -ph P,S,PcP,ScS,PKiKP -h 5 -km 100 @@ -76,8 +76,17 @@ TauP 中默认台站位于地表处，也可使用 ``--stadepth`` 指定台站 - 出射角 0° 到 90° 内表示射线向下出射（例如P震相） - 出射角 90° 到 180° 内表示射线向上出射（例如pP震相) -7. 入射角（度）：射线在台站（位于地表）处入射时与垂直向上方向的夹角。应该是恒正的。 -8. Purist 距离：地震波实际走过的圆弧距离。该值可能与震中距不同，例如:: + .. note:: + + 不同文章对出射角的定义不同，使用时应格外注意。 + + 老版本的TauP对出射角的定义不同，可视为BUG，应避免使用老版本TauP。 + +7. 入射角（度）：射线在台站处入射时与\ **垂直向上**\ 方向的夹角。 +8. Purist 距离：地震波实际走过的圆弧距离。 + + 该值可能与震中距不同，例如 50° 震中距处记录到的 PKKKP 震相，实际走过的 + 圆弧距离为 410°:: $ taup_time -ph PKKKP -deg 50 @@ -116,12 +125,6 @@ TauP 中默认台站位于地表处，也可使用 ``--stadepth`` 指定台站 60.00 300.0 S 1043.98 12.659 34.16 21.37 60.00 = S 426.93 60.00 300.0 ScS 1134.89 7.480 19.38 12.43 60.00 = ScS 517.85 -.. note:: - - 不同文章对出射角的定义不同，使用时应格外注意。 - - 老版本的TauP对出射角的定义不同，可视为BUG，应避免使用老版本TauP。 - 信息提取 -------- From 302e4c204a50e0113dd59b96946160aded432c5c Mon Sep 17 00:00:00 2001 From: Dongdong Tian Date: Fri, 1 Jan 2021 23:26:10 -0500 Subject: [PATCH 06/24] Tiny fix --- source/taup/taup_time.rst | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/source/taup/taup_time.rst b/source/taup/taup_time.rst index d0ea907e..d178f9a2 100644 --- a/source/taup/taup_time.rst +++ b/source/taup/taup_time.rst @@ -7,7 +7,7 @@ TauP 提供的 ``taup_time`` 命令可用于计算震相走时、射线参数、 基本用法 -------- -使用 ``taup_time --help`` 可查询该命令的完整用法和选项。此处仅介绍最常用用法。 +使用 ``taup_time -help`` 可查询该命令的完整用法和选项。此处仅介绍最常用用法。 计算震源深度为 300 km，震中距 60° 处，P、S、PcP 等震相的走时信息:: From 6a7cbb7d49dce3408a34b409918b06228ab4fd4b Mon Sep 17 00:00:00 2001 From: Dongdong Tian Date: Sat, 2 Jan 2021 01:41:53 -0500 Subject: [PATCH 07/24] Apply suggestions from code review Co-authored-by: Jiayuan Yao --- source/taup/taup_time.rst | 20 ++++++++++---------- 1 file changed, 10 insertions(+), 10 deletions(-) diff --git a/source/taup/taup_time.rst b/source/taup/taup_time.rst index d178f9a2..29f5255e 100644 --- a/source/taup/taup_time.rst +++ b/source/taup/taup_time.rst @@ -7,23 +7,23 @@ TauP 提供的 ``taup_time`` 命令可用于计算震相走时、射线参数、 基本用法 -------- -使用 ``taup_time -help`` 可查询该命令的完整用法和选项。此处仅介绍最常用用法。 +使用 ``taup_time -help`` 可查询该命令的完整用法和选项，此处仅介绍常用用法。 -计算震源深度为 300 km，震中距 60° 处，P、S、PcP 等震相的走时信息:: +计算震源深度为 300 km，震中距 60° 处，P、S、PcP、ScS、PKiKP 震相的走时信息:: - $ taup_time -ph P,S,PcP,ScS,PKiKP -deg 60 -h 300 + $ taup_time -h 300 -deg 60 -ph P,S,PcP,ScS,PKiKP -- ``-ph`` 指定了一系列以逗号分隔的震相名 -- ``-h`` 指定震源深度（单位 km） -- ``-deg`` 指定震中距（单位为°） +- ``-h``：震源深度（单位 km） +- ``-deg``：震中距（单位 °） +- ``-ph``：以逗号分隔的震相名 -TauP 默认使用 iasp91 模型，可以使用 ``-mod`` 选项设置要使用的速度模型:: +TauP 默认使用 iasp91 模型，可以使用 ``-mod`` 选项设置其他 TauP 预定义速度模型:: - $ taup_time -ph P,S,PcP,ScS,PKiKP -deg 60 -h 300 -mod prem + $ taup_time -mod prem -h 300 -deg 60 -ph P,S,PcP,ScS,PKiKP -TauP 中默认台站位于地表处，也可使用 ``--stadepth`` 指定台站深度（单位为 km）:: +TauP 默认台站位于地表，可使用 ``--stadepth`` 选项设置台站深度（单位为 km）:: - $ taup_time -ph P,S,PcP,ScS,PKiKP -deg 60 -h 300 --stadepth 0.5 + $ taup_time -h 300 -deg 60 --stadepth 0.5 -ph P,S,PcP,ScS,PKiKP 若已知震源和台站位置，也可以使用 ``-sta`` 和 ``-evt`` 分别指定台站和震源的 纬度和经度。例如，台站位于 (20°N,130°E)，震源位于 (40°N,150°E)，则:: From 51366c9f02389806862bd58c9ed4da8142ae94a9 Mon Sep 17 00:00:00 2001 From: Dongdong Tian Date: Sat, 2 Jan 2021 02:01:25 -0500 Subject: [PATCH 08/24] Update source/taup/taup_time.rst Co-authored-by: Jiayuan Yao --- source/taup/taup_time.rst | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/source/taup/taup_time.rst b/source/taup/taup_time.rst index 29f5255e..5bda54b1 100644 --- a/source/taup/taup_time.rst +++ b/source/taup/taup_time.rst @@ -25,10 +25,10 @@ TauP 默认台站位于地表，可使用 ``--stadepth`` 选项设置台站深 $ taup_time -h 300 -deg 60 --stadepth 0.5 -ph P,S,PcP,ScS,PKiKP -若已知震源和台站位置，也可以使用 ``-sta`` 和 ``-evt`` 分别指定台站和震源的 -纬度和经度。例如，台站位于 (20°N,130°E)，震源位于 (40°N,150°E)，则:: +若已知震源和台站位置，也可以使用 ``-evt`` 和 ``-sta`` 选项分别指定震源和台站的 +纬度和经度。例如，震源位于 (80 °N，150 °E)，台站位于 (20 °N，130 °E):: - $ taup_time -mod prem -ph P,S,PcP,ScS,PKiKP -h 100 -sta 20 130 -evt 40 170 + $ taup_time -mod prem -h 300 -evt 80 150 -sta 20 130 -ph P,S,PcP,ScS,PKiKP 在震中距比较小的情况下，可以使用 ``-km`` 选项指定震中距（单位为 km）。 比如计算一个浅源地震在 100 km 处常见 P 波震相的走时信息:: From ceb9a7593f351420f13de961d39ecfb44005bd5f Mon Sep 17 00:00:00 2001 From: Dongdong Tian Date: Sat, 2 Jan 2021 02:01:43 -0500 Subject: [PATCH 09/24] Update source/taup/taup_time.rst Co-authored-by: Jiayuan Yao --- source/taup/taup_time.rst | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/source/taup/taup_time.rst b/source/taup/taup_time.rst index 5bda54b1..6634a9d4 100644 --- a/source/taup/taup_time.rst +++ b/source/taup/taup_time.rst @@ -31,9 +31,9 @@ TauP 默认台站位于地表，可使用 ``--stadepth`` 选项设置台站深 $ taup_time -mod prem -h 300 -evt 80 150 -sta 20 130 -ph P,S,PcP,ScS,PKiKP 在震中距比较小的情况下，可以使用 ``-km`` 选项指定震中距（单位为 km）。 -比如计算一个浅源地震在 100 km 处常见 P 波震相的走时信息:: +比如计算震源深度为 5 km，震中距 100 km 处，P,S,PcP,ScS,PKiKP 震相的走时信息:: - $ taup_time -ph P,S,PcP,ScS,PKiKP -h 5 -km 100 + $ taup_time -h 5 -km 100 -ph P,S,PcP,ScS,PKiKP 输出格式 -------- From 0d2b32e29a0f9cef8f23fc35048101a1b09652e6 Mon Sep 17 00:00:00 2001 From: Dongdong Tian Date: Sat, 2 Jan 2021 02:08:01 -0500 Subject: [PATCH 10/24] Update source/taup/taup_time.rst Co-authored-by: Jiayuan Yao --- source/taup/taup_time.rst | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/source/taup/taup_time.rst b/source/taup/taup_time.rst index 6634a9d4..a71e885e 100644 --- a/source/taup/taup_time.rst +++ b/source/taup/taup_time.rst @@ -61,11 +61,11 @@ TauP 默认台站位于地表，可使用 ``--stadepth`` 选项设置台站深 .. note:: 射线参数有三种常用单位：秒/度（s/deg ）、秒/弧度（s/radian）、秒/千米（s/km）。 - 三者之间的关系为（:math:`R_0` 为地球半径）： + 三者之间的关系为（:math:`R_0` 为地球半径，假设为 6371 km）： .. math:: - 1 s/deg = \frac{1}{\pi*R_0/180} s/km = 1/111.2 s/km + 1 s/deg = \frac{1}{\pi*R_0/180} s/km = 1/111.19492 s/km .. math:: From 61c5cbdd7189bd9094691b4f12720a14482e7666 Mon Sep 17 00:00:00 2001 From: Dongdong Tian Date: Sat, 2 Jan 2021 02:08:20 -0500 Subject: [PATCH 11/24] Update source/taup/taup_time.rst Co-authored-by: Jiayuan Yao --- source/taup/taup_time.rst | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/source/taup/taup_time.rst b/source/taup/taup_time.rst index a71e885e..57b6ad87 100644 --- a/source/taup/taup_time.rst +++ b/source/taup/taup_time.rst @@ -69,7 +69,7 @@ TauP 默认台站位于地表，可使用 ``--stadepth`` 选项设置台站深 .. math:: - 1 s/deg = \frac{1}{\pi / 180} s/radian = 57.3 s/radian + 1 s/deg = \frac{1}{\pi/180} s/radian = 57.29578 s/radian 6. 出射角（度）：即射线从震源出射时与\ **垂直向下**\ 方向的夹角，取值范围为0到180°。 From 2623e1ad1d20465aaff806f297aeaf74df96d8f5 Mon Sep 17 00:00:00 2001 From: Dongdong Tian Date: Sat, 2 Jan 2021 02:28:58 -0500 Subject: [PATCH 12/24] Apply suggestions from code review Co-authored-by: Jiayuan Yao --- source/taup/taup_time.rst | 16 ++++++++-------- 1 file changed, 8 insertions(+), 8 deletions(-) diff --git a/source/taup/taup_time.rst b/source/taup/taup_time.rst index 57b6ad87..4d5d00ae 100644 --- a/source/taup/taup_time.rst +++ b/source/taup/taup_time.rst @@ -13,7 +13,7 @@ TauP 提供的 ``taup_time`` 命令可用于计算震相走时、射线参数、 $ taup_time -h 300 -deg 60 -ph P,S,PcP,ScS,PKiKP -- ``-h``：震源深度（单位 km） +- ``-h``：震源深度（单位 km），默认深度为 0.0 km。 - ``-deg``：震中距（单位 °） - ``-ph``：以逗号分隔的震相名 @@ -71,24 +71,24 @@ TauP 默认台站位于地表，可使用 ``--stadepth`` 选项设置台站深 1 s/deg = \frac{1}{\pi/180} s/radian = 57.29578 s/radian -6. 出射角（度）：即射线从震源出射时与\ **垂直向下**\ 方向的夹角，取值范围为0到180°。 +6. 出射角（度）：即射线从震源出射时与\ **垂直向下**\ 方向的夹角，取值范围为 0 到 180°。 - - 出射角 0° 到 90° 内表示射线向下出射（例如P震相） - - 出射角 90° 到 180° 内表示射线向上出射（例如pP震相) + - 0° 到 90° 表示射线向下出射（例如 P 震相） + - 90° 到 180° 表示射线向上出射（例如 pP 震相) .. note:: 不同文章对出射角的定义不同，使用时应格外注意。 - 老版本的TauP对出射角的定义不同，可视为BUG，应避免使用老版本TauP。 + 老版本的 TauP 对出射角的定义不同，可视为 BUG，应避免使用老版本 TauP。 -7. 入射角（度）：射线在台站处入射时与\ **垂直向上**\ 方向的夹角。 +7. 入射角（度）：射线入射到台站时与\ **垂直向上**\ 方向的夹角。 8. Purist 距离：地震波实际走过的圆弧距离。 该值可能与震中距不同，例如 50° 震中距处记录到的 PKKKP 震相，实际走过的 圆弧距离为 410°:: - $ taup_time -ph PKKKP -deg 50 + $ taup_time -deg 50 -ph PKKKP Model: iasp91 Distance Depth Phase Travel Ray Param Takeoff Incident Purist Purist @@ -102,7 +102,7 @@ TauP 默认台站位于地表，可使用 ``--stadepth`` 选项设置台站深 ``S^660S``。由于 PREM 模型中 660 km 不连续面的“真实”深度为670 km，此时 TauP 实际返回的时从 670 km 处反射回的震相，此时 Purist 震相名会显示为 ``S^670S``:: - $ taup_time -mod prem -ph S^660S -deg 120 + $ taup_time -mod prem -deg 120 -ph S^660S Model: prem Distance Depth Phase Travel Ray Param Takeoff Incident Purist Purist From a4cee85995c9d493ca7de93473d4a29a9e36046a Mon Sep 17 00:00:00 2001 From: Dongdong Tian Date: Sat, 2 Jan 2021 02:33:59 -0500 Subject: [PATCH 13/24] Update source/taup/taup_time.rst Co-authored-by: Jiayuan Yao --- source/taup/taup_time.rst | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/source/taup/taup_time.rst b/source/taup/taup_time.rst index 4d5d00ae..f0932630 100644 --- a/source/taup/taup_time.rst +++ b/source/taup/taup_time.rst @@ -110,10 +110,10 @@ TauP 默认台站位于地表，可使用 ``--stadepth`` 选项设置台站深 ------------------------------------------------------------------------------------ 120.00 0.0 S^660S 1975.75 12.336 20.79 20.79 120.00 * S^670S -当使用 ``--rel`` 选项时，``taup_time`` 可输出相对于某个指定震相的走时差。例如， +可以使用 ``--rel`` 选项输出相对于某个指定震相的走时差。例如， 下面的命令中，最后一列输出了所有震相相对于 PcP 震相的走时差:: - $ taup_time -ph P,S,PcP,ScS,PKiKP -deg 60 -h 300 -mod prem --rel PcP + $ taup_time -mod prem -h 300 -deg 60 -ph P,S,PcP,ScS,PKiKP --rel PcP Model: prem Distance Depth Phase Travel Ray Param Takeoff Incident Purist Purist Relative to From c8a58dcfdb12fb25d49faaeb928d47348452d694 Mon Sep 17 00:00:00 2001 From: Dongdong Tian Date: Sat, 2 Jan 2021 02:34:05 -0500 Subject: [PATCH 14/24] Update source/taup/taup_time.rst Co-authored-by: Jiayuan Yao --- source/taup/taup_time.rst | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/source/taup/taup_time.rst b/source/taup/taup_time.rst index f0932630..a3dc54c5 100644 --- a/source/taup/taup_time.rst +++ b/source/taup/taup_time.rst @@ -129,7 +129,7 @@ TauP 默认台站位于地表，可使用 ``--stadepth`` 选项设置台站深 -------- 上面介绍的输出中包含了很多信息，想要在脚本中提取出想要的信息有些麻烦。所以 -``taup_time`` 提供了一些选项，使得在脚本中计算走时什么的更加方便。 +``taup_time`` 提供了一些选项，使得在脚本中计算走时、射线参数等更加方便。 走时 ~~~~ From 5ab811907fd593fc96eb1d070f589a3bd3fb91be Mon Sep 17 00:00:00 2001 From: Dongdong Tian Date: Sat, 2 Jan 2021 02:37:49 -0500 Subject: [PATCH 15/24] Update source/taup/taup_time.rst Co-authored-by: Jiayuan Yao --- source/taup/taup_time.rst | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/source/taup/taup_time.rst b/source/taup/taup_time.rst index a3dc54c5..004920d2 100644 --- a/source/taup/taup_time.rst +++ b/source/taup/taup_time.rst @@ -134,10 +134,10 @@ TauP 默认台站位于地表，可使用 ``--stadepth`` 选项设置台站深 走时 ~~~~ -``--time`` 选项使得输出中只包含走时信息:: +``--time`` 选项使得输出中只包含走时信息，并且小数点的位数更多:: - $ taup_time -mod prem -ph P,S,PcP -deg 60 -h 100 --time - 595.3896 640.09875 1081.2472 + $ taup_time -mod prem -h 300 -deg 60 -ph P,S,PcP,ScS,PKiKP --time + 574.7245 617.047 993.4524 1043.9805 1134.8923 需要注意的是，输出中走时不是按照 ``-ph`` 选项中震相的顺序进行排序的， 而是按照走时递增的顺序排序。比如这里，S 震相放在震相列表的第二位， From 5f9f21d4f63ac5de7498bc6bf2863aa44a7b8fab Mon Sep 17 00:00:00 2001 From: Dongdong Tian Date: Sat, 2 Jan 2021 02:39:48 -0500 Subject: [PATCH 16/24] Update source/taup/taup_time.rst Co-authored-by: Jiayuan Yao --- source/taup/taup_time.rst | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/source/taup/taup_time.rst b/source/taup/taup_time.rst index 004920d2..09a97ba6 100644 --- a/source/taup/taup_time.rst +++ b/source/taup/taup_time.rst @@ -141,7 +141,7 @@ TauP 默认台站位于地表，可使用 ``--stadepth`` 选项设置台站深 需要注意的是，输出中走时不是按照 ``-ph`` 选项中震相的顺序进行排序的， 而是按照走时递增的顺序排序。比如这里，S 震相放在震相列表的第二位， -而 S 震相的走时 1081.2472 却放在第三位。 +而 S 震相的走时 1043.9805 却放在第四位。 射线参数 ~~~~~~~~ From 77cb7d4cbbd58d194962c9dc3e7932012f9cfaa5 Mon Sep 17 00:00:00 2001 From: Dongdong Tian Date: Sat, 2 Jan 2021 02:48:19 -0500 Subject: [PATCH 17/24] Update source/taup/taup_time.rst Co-authored-by: Jiayuan Yao --- source/taup/taup_time.rst | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/source/taup/taup_time.rst b/source/taup/taup_time.rst index 09a97ba6..1ea0acf0 100644 --- a/source/taup/taup_time.rst +++ b/source/taup/taup_time.rst @@ -146,12 +146,12 @@ TauP 默认台站位于地表，可使用 ``--stadepth`` 选项设置台站深 射线参数 ~~~~~~~~ -``--rayp`` 选项使得输出中只包含射线参数信息:: +``--rayp`` 选项使得输出中只包含射线参数信息，并且小数点的位数更多:: - $ taup_time -mod prem -ph P,S,PcP -deg 60 -h 100 --rayp - 6.8185554 3.9990747 12.785722 + $ taup_time taup_time -mod prem -h 300 -deg 60 -ph P,S,PcP,ScS,PKiKP --rayp + 6.739042 4.0158706 1.2450149 12.65922 7.480467 -此时，射线参数的输出是按照震相进行排序的。 +此时，输出的射线参数是按照 ``-ph`` 选项中震相的顺序进行排序的。 .. note:: From 76c9ea3b8fc9dcf523e5a5f8f2ce90663e7637cf Mon Sep 17 00:00:00 2001 From: Dongdong Tian Date: Sat, 2 Jan 2021 02:48:27 -0500 Subject: [PATCH 18/24] Update source/taup/taup_time.rst Co-authored-by: Jiayuan Yao --- source/taup/taup_time.rst | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/source/taup/taup_time.rst b/source/taup/taup_time.rst index 1ea0acf0..19a1ac94 100644 --- a/source/taup/taup_time.rst +++ b/source/taup/taup_time.rst @@ -139,7 +139,7 @@ TauP 默认台站位于地表，可使用 ``--stadepth`` 选项设置台站深 $ taup_time -mod prem -h 300 -deg 60 -ph P,S,PcP,ScS,PKiKP --time 574.7245 617.047 993.4524 1043.9805 1134.8923 -需要注意的是，输出中走时不是按照 ``-ph`` 选项中震相的顺序进行排序的， +需要注意的是，输出的走时不是按照 ``-ph`` 选项中震相的顺序进行排序的， 而是按照走时递增的顺序排序。比如这里，S 震相放在震相列表的第二位， 而 S 震相的走时 1043.9805 却放在第四位。 From 0e7a255456daad1bf59580f395b13c4b7fbccdf6 Mon Sep 17 00:00:00 2001 From: Dongdong Tian Date: Sat, 2 Jan 2021 02:53:07 -0500 Subject: [PATCH 19/24] Update source/taup/taup_time.rst Co-authored-by: Jiayuan Yao --- source/taup/taup_time.rst | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/source/taup/taup_time.rst b/source/taup/taup_time.rst index 19a1ac94..296c3dba 100644 --- a/source/taup/taup_time.rst +++ b/source/taup/taup_time.rst @@ -132,7 +132,7 @@ TauP 默认台站位于地表，可使用 ``--stadepth`` 选项设置台站深 ``taup_time`` 提供了一些选项，使得在脚本中计算走时、射线参数等更加方便。 走时 -~~~~ +^^^^^ ``--time`` 选项使得输出中只包含走时信息，并且小数点的位数更多:: From 26f022f565ee779ef8323812bbe8d6bd68c2617c Mon Sep 17 00:00:00 2001 From: Dongdong Tian Date: Sat, 2 Jan 2021 02:53:17 -0500 Subject: [PATCH 20/24] Update source/taup/taup_time.rst Co-authored-by: Jiayuan Yao --- source/taup/taup_time.rst | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/source/taup/taup_time.rst b/source/taup/taup_time.rst index 296c3dba..0675045e 100644 --- a/source/taup/taup_time.rst +++ b/source/taup/taup_time.rst @@ -144,7 +144,7 @@ TauP 默认台站位于地表，可使用 ``--stadepth`` 选项设置台站深 而 S 震相的走时 1043.9805 却放在第四位。 射线参数 -~~~~~~~~ +^^^^^^^^^ ``--rayp`` 选项使得输出中只包含射线参数信息，并且小数点的位数更多:: From c58132c2434c1ce814d2b59103b27986f6407ca4 Mon Sep 17 00:00:00 2001 From: Dongdong Tian Date: Sat, 2 Jan 2021 02:53:32 -0500 Subject: [PATCH 21/24] Update source/taup/taup_time.rst Co-authored-by: Jiayuan Yao --- source/taup/taup_time.rst | 1 + 1 file changed, 1 insertion(+) diff --git a/source/taup/taup_time.rst b/source/taup/taup_time.rst index 0675045e..f68df0e0 100644 --- a/source/taup/taup_time.rst +++ b/source/taup/taup_time.rst @@ -20,6 +20,7 @@ TauP 提供的 ``taup_time`` 命令可用于计算震相走时、射线参数、 TauP 默认使用 iasp91 模型，可以使用 ``-mod`` 选项设置其他 TauP 预定义速度模型:: $ taup_time -mod prem -h 300 -deg 60 -ph P,S,PcP,ScS,PKiKP + $ taup_time -mod ak135 -h 300 -deg 60 -ph P,S,PcP,ScS,PKiKP TauP 默认台站位于地表，可使用 ``--stadepth`` 选项设置台站深度（单位为 km）:: From 302b982b6d2d993384fd03f4142bc8b49a7c95e6 Mon Sep 17 00:00:00 2001 From: Dongdong Tian Date: Sat, 2 Jan 2021 03:03:00 -0500 Subject: [PATCH 22/24] Update source/taup/taup_time.rst Co-authored-by: Jiayuan Yao --- source/taup/taup_time.rst | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/source/taup/taup_time.rst b/source/taup/taup_time.rst index f68df0e0..941b1b23 100644 --- a/source/taup/taup_time.rst +++ b/source/taup/taup_time.rst @@ -27,9 +27,9 @@ TauP 默认台站位于地表，可使用 ``--stadepth`` 选项设置台站深 $ taup_time -h 300 -deg 60 --stadepth 0.5 -ph P,S,PcP,ScS,PKiKP 若已知震源和台站位置，也可以使用 ``-evt`` 和 ``-sta`` 选项分别指定震源和台站的 -纬度和经度。例如，震源位于 (80 °N，150 °E)，台站位于 (20 °N，130 °E):: +纬度和经度。例如，震源位于 (40 °N，10 °E)，台站位于 (20 °S，10 °W):: - $ taup_time -mod prem -h 300 -evt 80 150 -sta 20 130 -ph P,S,PcP,ScS,PKiKP + $ taup_time -mod prem -h 300 -evt 40 10 -sta -20 -10 -ph P,S,PcP,ScS,PKiKP 在震中距比较小的情况下，可以使用 ``-km`` 选项指定震中距（单位为 km）。 比如计算震源深度为 5 km，震中距 100 km 处，P,S,PcP,ScS,PKiKP 震相的走时信息:: From 862e657d3c4a18fce51daa235b02b41d78c9ac26 Mon Sep 17 00:00:00 2001 From: Dongdong Tian Date: Sat, 2 Jan 2021 03:07:29 -0500 Subject: [PATCH 23/24] Apply suggestions from code review Co-authored-by: Jiayuan Yao --- source/taup/taup_time.rst | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/source/taup/taup_time.rst b/source/taup/taup_time.rst index 941b1b23..7ec84e68 100644 --- a/source/taup/taup_time.rst +++ b/source/taup/taup_time.rst @@ -29,7 +29,7 @@ TauP 默认台站位于地表，可使用 ``--stadepth`` 选项设置台站深 若已知震源和台站位置，也可以使用 ``-evt`` 和 ``-sta`` 选项分别指定震源和台站的 纬度和经度。例如，震源位于 (40 °N，10 °E)，台站位于 (20 °S，10 °W):: - $ taup_time -mod prem -h 300 -evt 40 10 -sta -20 -10 -ph P,S,PcP,ScS,PKiKP + $ taup_time -h 300 -evt 40 10 -sta -20 -10 -ph P,S,PcP,ScS,PKiKP 在震中距比较小的情况下，可以使用 ``-km`` 选项指定震中距（单位为 km）。 比如计算震源深度为 5 km，震中距 100 km 处，P,S,PcP,ScS,PKiKP 震相的走时信息:: @@ -42,7 +42,7 @@ TauP 默认台站位于地表，可使用 ``--stadepth`` 选项设置台站深 下面的命令计算了 PREM 模型下，震源深度 100 km，震中距 60° 处，P 和 pP 震相的走时信息。 其输出为:: - $ taup_time -mod prem -h 100 -deg 60 -ph pP,P + $ taup_time -mod prem -h 100 -deg 60 -ph P,pP Model: prem Distance Depth Phase Travel Ray Param Takeoff Incident Purist Purist From 325889c3e5fdb65b49034ab4f2fd6eec92ae1145 Mon Sep 17 00:00:00 2001 From: Dongdong Tian Date: Sat, 2 Jan 2021 03:16:14 -0500 Subject: [PATCH 24/24] Update source/taup/taup_time.rst Co-authored-by: Jiayuan Yao --- source/taup/taup_time.rst | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/source/taup/taup_time.rst b/source/taup/taup_time.rst index 7ec84e68..00fb7b12 100644 --- a/source/taup/taup_time.rst +++ b/source/taup/taup_time.rst @@ -130,7 +130,7 @@ TauP 默认台站位于地表，可使用 ``--stadepth`` 选项设置台站深 -------- 上面介绍的输出中包含了很多信息，想要在脚本中提取出想要的信息有些麻烦。所以 -``taup_time`` 提供了一些选项，使得在脚本中计算走时、射线参数等更加方便。 +``taup_time`` 提供了一些选项，使得在脚本中计算走时、射线参数更加方便。 走时 ^^^^^