TTT 計算メモ †
コマンド用法 †
ttt <inputfile> [-B] [-D<zmin>] [-E<lon/lat> or -e<file>] [-Ga|b|f|i] [-N<nodes>] [-O<mm/dd/yyyy/hh/mi>[/<ss>]] [-R<w/e/s/n>] [-S[<radius>][/<depth>]] [-T<ttt_prefix>] [-V[l]]
inputfile の部分は下記の ttt_topo_XXX.i2 のファイル（XXX は解像度）のうちのいずれかを選ぶ．
ttt_topo_10m.i2 ttt_topo_15m.i2 ttt_topo_1m.i2 ttt_topo_20m.i2 ttt_topo_2m.i2 ttt_topo_30m.i2 ttt_topo_30s.i2 ttt_topo_5m.i2 ttt_topo_60m.i2
デフォルトの通りインストールしていれば /usr/share/geoware/ttt/data にこれらのファイルが格納されているはず．
拡張子の .i2 を含めないように注意．
ttt ttt_topo_10m (...続き...)
のように書く．
その他のオプションは GMT と似たような使い方．R オプションは GMT と同じだと思って問題ない．
各オプションの詳細は以下の通り．
Usage: ttt <inputfile> [-B] [-D<zmin>] [-E<lon/lat> or -e<file>] [-Ga|b|f|i] [-N<nodes>]
[-O<mm/dd/yyyy/hh/mi>[/<ss>]] [-R<w/e/s/n>] [-S[<radius>][/<depth>]] [-T<ttt_prefix>] [-V[l]]
ttt attempts to decode <inputfile> using the following order:
1. If filename ends in ".b" it is read as bathymetry data [GMT binary float format].
2. If filename ends in ".i2" it is read as bathymetry data [GMT binary short format].
3. If $TTT_DIR/<inputfile>.i2 exists it will be used as bathymetry data [GMT binary short format].
OPTIONS:
-B Do NOT normalize the travel times [Default normalizes to avoid bias].
-D Specify a minimum water depth [0]. Depths shallower than 2*z0 will be adjusted
so that depth goes quadratically from z0 to 2*z0 instead of 0-2*z0.
-E Set the location of the epicenter.
-e Give filename with multiple "epicenters" to mimic a non-point source.
-G Store travel times in one of four formats. Choose from
-Ga Store travel times as ESRI Arc/Info ASCII Grid Interchange format (*.asc).
-Gb Store travel times as ESRI Gridfloat file (separate header *.hdr and binary *.flt files).
-Gf Store travel times as GMT native float grid (*.b) [Default].
-Gi Store travel times as GMT native 2-byte int grid with units of 10 sec (*.i2).
Do no supply a file extension with -T as -G will set it.
-N Number of Huygens nodes to use (8, 16, 32, 48, or 64) [64]
-O Set earthquake origin time (UTC). Use lower case -o if local time is used.
-R Specify a sub-region of the grid [Use entire grid].
-S Substitute nearest ocean node if epicenter is on land.
Optionally, append search radius in degrees [5].
Furthermore, you may append the minimum depth you want to place epicenter [0].
-T Set prefix for the output grid file with travel times in hrs [ttt].
Extension will be added automatically (see -G for file formats).
-V Give feedback while running [Default is quiet].
Append l to see relocated epicenter (if using -S).
計算例 †
TTT_DIR=/usr/share/geoware/ttt/data
epicenter="142.86E/38.10333N"
region="135/150/15/45"
topo="ttt_topo_10m"
output="test1"
#ttt64 ${topo} -E${epicenter} -R${region} -T${output} -V
ttt64 ${topo} -E${epicenter} -T${output} -V
## makecpt
#gmt makecpt -Chot -T-10/30/1 > tmp.cpt
#echo 'N gray' >> tmp.cpt
gmt makecpt -Cmagma -T0/30/1 > tmp.cpt
## makegrd
gmt grdcut ${TTT_DIR}/${topo}.i2 -R${output}.b=bf -Gtmp.i2
gmt grdgradient tmp.i2 -A0 -Nt1.2 -Gtmp.nc
## plot
gmt grdimage ${output}.b=bf -Itmp.nc -JA${epicenter}/12 -Rg -Ctmp.cpt -P -K -V > ${output}.ps
gmt psscale -DjMR+w10.0/0.25+o-1.0/-0.2 -Ctmp.cpt -Ba5f5 -J -R -P -K -V -O >> ${output}.ps
gmt grdcontour ${output}.b=bf -J -O -K -C1 -Wgray -GlZ-/Z+ >> ${output}.ps
gmt pscoast -R${output}.b=bf -J -O -Glightbrown -W0.1p -A500 -Baf >> ${output}.ps
rm tmp.*
ps2pdf ${output}.ps
pdfcrop -margins '5 5 5 5' ${output}.pdf
convert -density 300 ${output}-crop.pdf ${output}.png
rm ${output}*.pdf
