Plotting Centroid Moment Tensor

Data Global Centroid Moment Tensor bisa diunduh dari situs globalcmt.org untuk selanjutnya kita plot menggunakan Generic Mapping Tools. Screenshot hasil filtering data untuk Nusa Tenggara Timur (dengan pilihan format output sebagai input psmeca) bisa dilihat di bawah ini:

Hasil Pencarian Dari Global CMT Catalog
Hasil Pencarian Dari Global CMT Catalog

Setelah kita simpan data di atas pada file gcmt_ntt.gmt (atau nama pilihan Anda), maka selanjutnya perlu kita filter menggunakan perintah UNIX awk untuk mendapatkan data yang kita perlukan dan membuang data pada kolom lain yang tidak akan kita gunakan untuk pengeplotan ini.

$ awk '{print $1, $2, $3, $4, $5, $6, $7, $8, $9, $10}' gcmt_ntt.gmt > psm_ntt.gmt

Centroid Moment Tensor kita plot dengan perintah psmeca. Pastikan GMT Anda mempunyai perintah suplement ini. Distribusi GMT linux ubuntu versi 5.2.1 tidak menyertakan psmeca seperti ditulis pada forum ini.

Script untuk pengeplotan semua CMT  ada di bawah ini:

#!/usr/bin/bash

# Plot CMT
# GMT 5.2.1
# oleh: Hendro Nugroho

gmt set MAP_FRAME_TYPE plain

# Variabel
out='cmt_ntt.ps'
area='-R118.5/127.5/-11.5/-7.5'

# extraksi kolom yang kita perlukan
#awk '{print $1, $2, $3, $4, $5, $6, $7, $8, $9, $10}' gcmt_ntt.gmt > psm_ntt.gmt

# frame peta
gmt psbasemap $area -JM8.5i -Xc -Yc -Bxa2f1 -Bya1f1 -B+t"CMT Nusa Tenggara" -K > $out

# peta garis pantai
gmt pscoast -R -JM -Xc -Yc -B -Di -Ggray -O -K >> $out

# plot CMT
gmt psmeca psm_ntt.gmt -R -JM -Xc -Yc -Sm0.4 -h13 -O >> $out

Hasil pengeplotan CMT dengan zero trace (-Sm):

CMT Nusa Tenggara Timur
CMT Nusa Tenggara Timur (zero trace)

Dan hasil plot hanya bagian double couple dari moment tensor saja (-Sd):

cmt_ntt_sd
CMT Nusa Tenggara Timur (double couple)

Nah sekarang tugas Anda adalah membuat plot ini lebih informatif dan selektif.
Hal-hal yang bisa dikerjakan antara lain: 1) membedakan cmt untuk gempabumi, dangkal, sedang, dan dalam dengan menggunakan warna plot cmt yang berbeda, 2) melakukan pengeplotan menggunakan offset lokasi; gempa bumi diplot sebagai titik, lalu CMT tertentu diplot pada posisi off dengan referensi garis mengacu pada titik lokasi gempa.

Advertisements

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s