Membuat Peta Seismisitas

Untuk membuat peta seismisitas, diperlukan data lokasi gempa yang bisa kita peroleh dari USGS, GFZ, BMKG, atau situs-situs lainnya. Berikut adalah contoh data seismisitasnya :
96.46 4.49 
99.12 -0.71 
92.98 2.59 
93.18 2.45 
......................
kolom 1 = longitude
kolom 2 = latiude


disini saya simpan dengan nama "gempa.gmt". Selanjutnya buatlah baris-baris perintah berikut :

pscoast -R90/107/-6/7 -JM20c -Ba10g5 -Ggreen3 -S140/235/255 -Dh -N1-W1 -K > map.ps
psxy gempa.gmt -R -JM -Sc0.5 -G255/255/0 -W1 -O -K >> map.ps
echo 98.7 3.59 | psxy -R -JM -Sa0.5c -Gred -W1 -O -K >> map.ps
echo 98.7 3.59 14 0 1 LT TSI | pstext -R -JM -O -K >> map.ps



Setelah dieksekusi, hasilnya akan seperti ini :

Arti dari baris-baris perintah diatas adalah :
pscoast = file postscript untuk membuat garis pantai
-N1 = opsi pada pscoast untuk menampilkan batas negara
-W1 = opsi untuk menampilkan garis tepi dengan ketebalan garis 1 point
-K = opsi untuk menunjukkan bahwa masih ada perintah selanjutnya
-Sc0.5 = opsi pada psxy untuk menampilkan lingkaran dengan ukuran 0.5
-O = opsi untuk mengoverlay peta
psxy = file postscript untuk mem-plot titik
pstext = file postscript untuk menampilkan text
echo = disini dipakai untuk membuat titik yang akan diplot, setelah echo ada simbol pipeline (|) untuk memisahkannya dengan perintah pada GMT
-Sa0.5c = opsi untuk memplot simbol bintang dengan ukuran 0.5 cm
Perhatian, perintah pertama pada GMT selalu berakhiran ">", dan perintah selanjutnya berakhiran dengan ">>".

Semoga bermanfaat.