夏の学校で作った絵の書き方 (GMT 3.2)
地形を3次元的に書く
必要なデータ
海のxyzファイル（ここでは名前をsea.xyzとする）
陸のxyzファイル（ここでは名前をjapan.xyzとする）
シェル
#!/bin/csh -f
surface japan.xyz -GTmp1.grd -R129/147/30/46 -I0.05/0.05 -V
surface sea.xyz -GTmp2.grd -R129/147/30/46 -I0.05/0.05 -V
grdmath Tmp1.grd 1.5 x = Tmp3.grd
grdmath Tmp3.grd Tmp2.grd ADD = all.grd
grdgradient all.grd -A290 -Gtmp2.grad -V
grdhisteq tmp2.grad -Gtmp3.hst -N -V
grdmath tmp3.hst 3. / = tmp4.hst
grdview all.grd -JX15.0/12.5 -Jz0.00015 -E160/15 -R129/147/30/46/-7000/4000 \
-Ccolor.cpt -Qs -Itmp4.hst -X5 -Y10 -V >! map1.ps
解説
1: 日本地図のグリットファイル（*.grd）を作る．
2：海のグリットファイルを作る．
3：陸の高度を高くして，山を強調する．
4：日本地図と海のグリットファイルを結合する．
5-7：光のグリットファイルを作る．
8：作ったグリットファイルをプロットする．
震源の3次元プロット
必要なデータ
震源データ(hypo.dat)
参考までに今回使った震源データファイルのフォーマットは，
年 月 日 時間（時分秒） 緯度 経度 深さ マグニチュード
1998 07 06 21 19 5.50 136.891 35.813 6.58 1.1
シェル
#!/bin/csh -f
awk '{print $7,$8,$9,$9}' hypo.dat | psxyz -JX15.0/12.5 -Jz-0.015 -E160/15 -R129
/147/30/46/0/400 -Su0.06 -Y2. -Ccolorh.cpt -K -B10/10/100WSneZ+ > map2.ps
pscoast -JX -E160/15 -R129/147/30/46 -B -W2.0 -O -Y5.8 -Dh >> map2.ps
説明無し，
協力者：中川茂樹