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
test1.png

Front page   Edit Diff Attach Copy Rename Reload   New List of pages Search Recent changes   Help   RSS of recent changes
Last-modified: 2020-11-25 (Wed) 06:11:40 (4h)