TauP Toolkitについて
Last update: 2020.8.7

●準備(Javaのインストール)

TauP Toolkitを利用するために Java が必要になります.
Java がインストールされていない場合のみ追加して下さい。

▼Cent OS 7.x の yum コマンドでのインストール例
# su -
# yum -y install java-1.8.0-*


●TauP Toolkitのインストール

▼取得
http://www.seis.sc.edu/taup/index.html
から
TauP-2.4.5.tgz
を取得.

▼Javaが利用できる環境かどうか確認
# java -version

インストールされている場合の表示の例

openjdk version "1.8.0_191"
OpenJDK Runtime Environment (build 1.8.0_191-b12)
OpenJDK 64-Bit Server VM (build 25.191-b12, mixed mode)
などと出力される.

▼unzip & untar
# tar xvfz TauP-2.4.5.tgz -C /usr/local

※※※※※※※※※※※※※※※※※※※※※
※ ここから下は,各ユーザーごとの設定 ※
※※※※※※※※※※※※※※※※※※※※※

▼環境変数の設定
bashの場合
下記の行を.bash_profileに追加する.

export TAUP_HOME=/usr/local/TauP-2.4.5
export PATH=${TAUP_HOME}/bin:$PATH


▼実行のテスト1
# taup_time -mod iasp91 -h 400 -deg 80 -ph P,S,ScS

上記の実行をすると,下記のような出力結果が出てきます.

Model: iasp91
Distance Depth Phase Travel Ray Param Purist Purist
(deg) (km) Name Time (s) p (s/deg) Distance Name
----------------------------------------------------------------
80.0 400.0 P 686.11 5.278 80.0 = P
80.0 400.0 S 1256.51 10.295 80.0 = S
80.0 400.0 ScS 1273.05 8.174 80.0 = ScS


▼実行のテスト2
# taup_table -ph P

上記の実行をすると,下記のような出力結果が出てきます.

iasp91 0.00 0.0 P 0.00 19.172 0.00 P
iasp91 0.10 0.0 P 1.92 19.172 0.10 P
iasp91 0.20 0.0 P 3.83 19.172 0.20 P
iasp91 0.30 0.0 P 5.75 19.171 0.30 P
iasp91 0.40 0.0 P 7.67 19.171 0.40 P
以下, 指定したモデル(iasp91など)での震源の深さ,震央距離別の走時表が表示される


▼実行のテスト3
# taup

上記の実行をすると,GUIで taup_times, taup-pierce, taup_pathの実行ができます.

taup screenshot 01


●地下構造モデルの追加 (2020.8.7)

自作の地下構造モデルを追加することができます。
.tvelフォーマットと.ndフォーマットの2種類のフォーマットがあります。
.tvelフォーマットの場合は、最初の2行がコメントで、3行目から各層のパラメータとなっています。
Depth(km) Vp(km/s) Vs(km/s) ρ(g/cm^3)
の4列で、層の数だけ行数を追加していきます。
.ndフォーマットはUser's Guideの5.1に説明があります。
TauP-2.4.5.tgzの展開後のStdModelsディレクトリにあるモデルのファイルも参考になります。

作成例
1..tvelフォーマットで地下構造モデルを作成
以下の7行のファイルを作成して、ファイル名は model01.tvel として保存する。
# model01  .tvel format for Taup Toolkit
# The first two lines of this file is ignored by TauP. taup_create -tvel model01.tvel
0.0       4.800  2.844  2.7
20.0     6.450  3.754  2.7
35.0     7.230  4.161  2.7
100.0   7.949  4.446  2.7
600.0 10.250  5.679  2.7
2.作成した地下構造モデルを使用する例
-mod オプションの後にファイル名を指定して利用できます。
taup_timeでの実行例
$ taup_time -mod ./model01.tvel -km 10 -h 20 -ph P,p
Model: model01
Distance   Depth   Phase   Travel    Ray Param  Takeoff  Incident  Purist    Purist
  (deg)     (km)   Name    Time (s)  p (s/deg)   (deg)    (deg)   Distance   Name
-----------------------------------------------------------------------------------
    0.09    20.0   p         3.585     0.085    176.91     2.22     0.09   = p
層の数は多くても問題ないので、JMA2001モデル用のファイルを作って走時計算することもできます。


●参考資料
Crotwell, H.P., T.J. Owens, and J. Ritsema (1999).
The TauP Toolkit: Flexible seismic travel-time and raypath utilities,
Seis. Res. Letters, 70, 154-170.

プログラム付属のドキュメント(TauP/doc/taup.pdf)
The TauP Tookkit: Flexible Seismic Travel-Time and Raypath Utilities Version 2.0 Documentation