Membuat Focal Mechanism dari Global CMT

Kali ini saya akan berbagi cara mem-plot focal mechanism. Data seismisitas dapat kita peroleh dari Global CMT. Cara download datanya dengan mengisi parameter-parameter yang diminta, dan untuk output type pilihlah "GMT psmeca input". Jika sudah selesai, salinlah data yang muncul ke notepad dan disimpan dengan ekstensi ".gmt".

Berikut adalah contoh data yang didapat :

96.46 4.49 12 -0.05 -1.35 1.40 -0.87 -0.08 -0.38 24 X Y 201201011809A  
99.12 -0.71 25 0.95 -0.42 -0.52 4.46 -0.79 0.73 23 X Y 201201052314A  
92.98 2.59 24 -0.51 -2.96 3.47 -1.07 -0.90 -6.68 26 X Y 201201101837A  
.......................

disini saya simpan dengan nama "focal.gmt". Sebenarnya kalau kita hanya mem-plot titik saja, maka kolom yang diperlukan hanya kolom satu dan dua, yaitu longitude dan latitude. Tetapi karena saya juga mem-plot focal mechanism nya, maka diperlukan format lengkap seperti diatas.

Langkah selanjutnya untuk mem-plot focal mechanism, tulislah baris-barus perintah dibawah ini dan simpanlah dengan ekstensi ".bat" :

pscoast -R90/107/-6/7 -JM20c -Ba10g5 -Ggreen3 -S140/235/255 -Dh -N1 -W1 -K > map.ps
psmeca focal.gmt -R -JM -Sm0.5/-1 -G0/0/255 -T0 -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 berikut :

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 pantai dengan ketebalan garis 1 point
-K = opsi untuk menunjukkan bahwa masih ada perintah selanjutnya
psmeca = file postscript untuk membuat focal mechanism
-Sm0.5/-1 = opsi pada psmeca untuk menampilkan bola fokal dengan ukuran 0.5, arti dari -1 adalah label bola fokal tidak ditampilkan
-T0 = opsi pada psmeca untuk menampilkan kedua bidang nodal
-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 ">>".

Untuk pengertian dari perintah-perintah yang lain dapat dilihat di manual GMT.
Silakan dieksplorasi dan semoga bermanfaat.
Salam olahraga...