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の実行ができます.
●地下構造モデルの追加 (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