津波波形データベース構築の一例として，１月７日にプエルトリコで発生した地震による津波のデータ収集と計算をやってみた．
波源断層モデル
Moment Tensor (USGS)
(17.916, 293.187) 13.5 km, Dip: 43 deg., Strike: 268 deg., Rake: -58 deg.
M0 - S 関係 (cite) から，
(L, W) = (20.0073963517812 km, 10.00369817589060 km)
slip amount: 0.500816743451796 m
DEM作成
GEBCO 2019から，当該領域の切り出し (grdcut) を行った．GEBCO_2019.nc は， grdmath を用いて正負を反転してある．
% grdmath GEBCO/GEBCO_2019/GEBCO_2019.nc NEG = GEBCO_2019.nc
% grdcut -R290/300/15/20 $bathymetry/GEBCO_2019.nc -Ggebco.grd
GEBCO 2019 から切り出した地形データでJAGURSを走らせると，北西の領域端にかかっている島のところで発散が生じた．そのため，周囲３分角（15秒角の12メッシュ）は，水深100 mよりも浅い海と陸地を水深 100 mで固定する修正を grdmask と grdmath で行った．下記は， GMT 6.0.0 (Homebrew on OS X) を使用している（ -I15s は， -I15c に相当）．
% cat<<EOF | grdmask -A -Gmask.grd -I15s -Rgebco.grd -N100/100/NaN
290:3 15:3
290:3 19:57
299:57 19:57
299:57 15:3
EOF
（解説）
grdmask で，標準入力から縁取りの４点を入力して， -N 領域外/線上/内側 の値を指定する． -Rgebco.grd で元のDEMと範囲を合わせて，分解能も指定する（-I15s）．さらに，大円ではなく直線で四角にするため -A オプションを付ける．
-A Suppress connecting geographic points using great circle arcs, i.e., connect by straight lines,
unless m or p is appended to first follow meridian then parallel, or vice versa.
波源断層モデル
Moment Tensor (USGS)
(17.916, 293.187) 13.5 km, Dip: 43 deg., Strike: 268 deg., Rake: -58 deg.
M0 - S 関係 (cite) から，
(L, W) = (20.0073963517812 km, 10.00369817589060 km)
slip amount: 0.500816743451796 m
DEM作成
GEBCO 2019から，当該領域の切り出し (grdcut) を行った．GEBCO_2019.nc は， grdmath を用いて正負を反転してある．
% grdmath GEBCO/GEBCO_2019/GEBCO_2019.nc NEG = GEBCO_2019.nc
% grdcut -R290/300/15/20 $bathymetry/GEBCO_2019.nc -Ggebco.grd
GEBCO 2019 から切り出した地形データでJAGURSを走らせると，北西の領域端にかかっている島のところで発散が生じた．そのため，周囲３分角（15秒角の12メッシュ）は，水深100 mよりも浅い海と陸地を水深 100 mで固定する修正を grdmask と grdmath で行った．下記は， GMT 6.0.0 (Homebrew on OS X) を使用している（ -I15s は， -I15c に相当）．
% cat<<EOF | grdmask -A -Gmask.grd -I15s -Rgebco.grd -N100/100/NaN
290:3 15:3
290:3 19:57
299:57 19:57
299:57 15:3
EOF
% grdmath gebco.grd mask.grd LT mask.grd MUL 0 NAN gebco.grd AND = masked_gebco.grd=cf
（解説）
grdmask で，標準入力から縁取りの４点を入力して， -N 領域外/線上/内側 の値を指定する． -Rgebco.grd で元のDEMと範囲を合わせて，分解能も指定する（-I15s）．さらに，大円ではなく直線で四角にするため -A オプションを付ける．
-A Suppress connecting geographic points using great circle arcs, i.e., connect by straight lines,
unless m or p is appended to first follow meridian then parallel, or vice versa.
験潮所
ITIC (International Tsunami Information Center) のアラートで 0.02 mの津波が観測されたようなので，UNESCO IOC Sea Level Station Monitoring Facility の event からMagueyes Island PR staion のデータを "show data" から取得． 200〜2000秒のバンドパスフィルタをかけた．
(lat, lon) = (17.9701, -67.0464) だったが，DEMに重ね合わせると，ずれていたので， (17.9636, 292.9536) に移動した．
JAGURS計算
設定した断層パラメータで， dt=1.0 s で計算した．
