2004年スマトラ島M9.0津波Sim備忘録Part1 27．Dec．2004−
（一度書いて消してしまったのでもう一度）
http://ofgs.ori.u-tokyo.ac.jp/~okino/gmtscripts/
を参考にした．
＜gmt4.0のインストール＞
本家よりinstall_gmtを落とし，これに実行属性をつ
けて実行するだけ．
自動的にインストールされる（常時ネット接続環境下で）．
あと
.gmtdefaultsの
PAPER_MEDIA = a1
にいつものように変更．
＜地形水深データの作成＞
http://www.gina.alaska.edu/page.xml?group=data&page=griddata
から
North
East GMT File bzipped
North
FarEast GMT File bzipped
South
East GMT File bzipped
South
FarEast GMT File bzipped
の４種類を落とす．（ちょうどまたがるため）
これをgrdcutで
grdcut North_East.grd -R70/90/0/25 -GNE.grd
grdcut North_farEast.grd -R90/110/0/25 -GNF.grd
grdcut South_East.grd -R70/90/-10/0 -GSE.grd
grdcut South_farEast.grd -R90/110/-10/0 -GSF.grd
そしてこれを統合
grdpaste NE.grd NF.grd -GNN.grd
grdpaste SE.grd SF.grd -GSS.grd
grdpaste NN.grd SS.grd -GINDIA.grd
これで東経70-110，南緯-10から25
の格子が出来上がる．
カラーcptは下記のものをエディタで調整INDIA_rel2.cptを作る
makecpt -Crelief >INDIA_rel.cpt
これは陸地の色が通常の地図と逆転していることに注意．しかし海と陸の区別がつきやすい．また水深もわかりやすい．
＜海底地形図script＞
#! /bin/csh
# make color contour map (gradation + contour)with scale bar
# http://ofgs.ori.u-tokyo.ac.jp/~okino/gmtscripts/colcnt_grad.htmlよ
り改変
# for INDIAN Ocean
region
29.Dec2004 Y.Okamoto
set path = /home/seagull/GMT/GMT4.0/bin
set region =
70/110/-10/25
# map region INDIAN Ocean
set proj =
M58
# map projection and scale
set boundary =
f1a10g10WSNE
# boundary tick info
set reso =
i
# coastline resolution
set grdfile =
INDIA.grd
# input grid bathymetry
set cptfile =
INDIA_rel2.cpt
# color table
set cont = 500
# contour interval
set anot =
1000t
# anotation interval (transparent)
set limit =
-8000/4000
# min an max value for contour
set psfile =
INDIA_rel3.ps
# output postscript file name
#set scaleloc =
10.5/3/6/0.2
# scale bar location
set scaleloc =
10/-1.3/16/0.5h
# scale bar location horizontal
set scaletick = 'f1000a2000g1000:Bathymetry-Topo:' # scale bar tick info
#
#makecpt -Chaxby -T-10000/0/250 -Z > $cptfile
#
grdimage $grdfile -Y5 -X7 -R$region -J$proj -C$cptfile -K -V >
$psfile
#-Y-X adjust the origin point
grdcontour $grdfile -R$region -J$proj -C$cont -A$anot -W1a2 -L$limit -K
-O -V >> $psfile
#psscale -D$scaleloc -Y-2 -B$scaletick -C$cptfile -K -O -V >>
$psfile
#pscoast -R$region -J$proj -D$reso -W2 -G220 -K -O -V >> $psfile
psbasemap -R$region -J$proj -B$boundary -K -O -V >> $psfile
psscale -D$scaleloc -B$scaletick -C$cptfile -O -V >> $psfile
#gv $psfile
※Nortilusおよびgvはps画像の右端が正方形エリアでクリップされて
出ないことに注意！gimpではちゃんと出る．
＜スマトラ島近海バージョン海底地形図＞は下記
grdcut INDIA.grd -R78/103/-2/23 -GSUMATORA.grd
これだと3000×3000の格子になる．
#! /bin/csh
# make color contour map (gradation + contour)with scale bar
# http://ofgs.ori.u-tokyo.ac.jp/~okino/gmtscripts/colcnt_grad.htmlよ
り改変
# for INDIAN Ocean
region
29.Dec2004 Y.Okamoto
set path = /home/seagull/GMT/GMT4.0/bin
set region =
78/103/-2/23
# map region INDIAN Ocean
set proj =
M52
# map projection and scale
set boundary =
f1a10g10WSNE
# boundary tick info
set reso =
i
# coastline resolution
set grdfile =
SUMATORA.grd
# input grid bathymetry
set cptfile =
INDIA_rel2.cpt
# color table
set cont = 500
# contour interval
set anot =
1000t
# anotation interval (transparent)
set limit =
-8000/4000
# min an max value for contour
set psfile =
SUMATORA.ps
# output postscript file name
#set scaleloc =
10.5/3/6/0.2
# scale bar location
set scaleloc =
10/-1.3/16/0.5h
# scale bar location horizontal
set scaletick = 'f1000a2000g1000:Bathymetry-Topo:' # scale bar tick info
#
#makecpt -Chaxby -T-10000/0/250 -Z > $cptfile
#
grdimage $grdfile -Y5 -X7 -R$region -J$proj -C$cptfile -K -V >
$psfile
#-Y-X adjust the origin point
grdcontour $grdfile -R$region -J$proj -C$cont -A$anot -W1a2 -L$limit -K
-O -V >> $psfile
#psscale -D$scaleloc -Y-2 -B$scaletick -C$cptfile -K -O -V >>
$psfile
#pscoast -R$region -J$proj -D$reso -W2 -G220 -K -O -V >> $psfile
psbasemap -R$region -J$proj -B$boundary -K -O -V >> $psfile
psscale -D$scaleloc -B$scaletick -C$cptfile -O -V >> $psfile
#gv $psfile