GMT演習 震源位置、マグニチュード、震源深さのプロット

【準備するもの】
震源位置(緯度・経度)、震源深さ、マグニチュードのデータ、4列のデータ(earth.dat)

pscoast -JM -R130/145/30/45 -Df -S128 -G255 -B1 -W.1/255 -P -X1 -Y3 -K > depth_Mw.eps

makecpt -Cno_green -T0/400/10 > rb1.cpt
awk '{print $1,$2,$3,$4*0.05}' earth.dat | psxy -JM -R -Sc -Crb1.cpt -O -K -L -W3 >> depth_Mw.eps

psscale -Crb1.cpt -D7.5/-1/15/0.3h -B10 -O >> depth_Mw.eps