Plotting Sinyal Gempa Dengan GMT

Kali ini saya akan berbagi cara menampilkan sinyal gempa dengan GMT. Sebelumnya carilah event gempa yang anda inginkan, selanjutnya carilah sinyal dari event tersebut disini atau disini. Data yang anda download nantinya akan berformat SEED. Rubahlah menjadi format ASCII dan pisahkan menurut komponen-komponennya. Berikut adalah contoh data yang sudah saya rubah menjadi ASCII dan dipisah menurut komponenya. Data tersebut adalah data event gempa tanggal 1 Juli 2011 dari stasiun Lembang (LEM). Panjang data adalah 73 sekon. Episenter gempa terletak pada koordinat 108.82 BT - 8.19 LS. Selain itu dapat anda tambahkan pula data mekanisme fokal dari Global CMT. Berikut adalah data mekanisme fokal yang saya dapatkan :

108.82 -8.19 103 -3.72 5.55 -1.83 5.82 -1.70 0.29 24 104 -11 010177A          

Simpan dengan nama "jawa.gmt".

Episenter gempa terletak pada kolom 1 dan 2. Tetapi saya akan memplotnya pada koordinat di kolom 11 dan 12.

Sedangkan data episenter gempa yang saya dapatkan adalah :

2011    7    1    14    33    22.8  -8.19    108.82   103    4.6    4.6

Simpan dengan nama "maineq.dat".

Selanjutnya buatlah script berikut :

set F=signal.ps
set R=100/115/-15/0

#Membuat file warna
makecpt -Z -Cglobe > topo.cpt

#Plot peta dasar
psbasemap -JM12 -X4.5 -R%R% -Ba4g4WSne:."1 JULI 2011 EARTHQUAKE": --HEADER_FONT_SIZE=20 -P -K -Y15> %F%
pscoast -JM -R -Ggreen -Sblue -Dh -Wthin -K -P -O >> %F%
grdimage indo.nc -R -JM -Ctopo.cpt -K -O -Na -P >> %F%

#Plot kota Lembang
echo 107.6 -6.82 90 7.4 LEMBANG LT > lembang.d
psxy -R -J -O -K -Ss0.2 -Gred -Wthinnest lembang.d >>%F%
gawk "{print $1, $2, 9, 0, 1, $6, $5}" lembang.d | pstext -R%R% -J -O -K -Dj0.15/0 -N >> %F%

#Plot batas subduksi
psxy -R -JM -W10.0 -Sf0.8i/0.1ilt -Gblack -O -K trench.gmt -m>> %F%

#Plot garis
echo 108.82 -8.19 14 0 1 LT B> garis.dat
echo 104.0 -11 14 0 1 LT A>> garis.dat
psxy garis.dat -JM -R -W10 -P -O -K >> %F%

#Plot bola fokal
psmeca jawa.gmt -JM -R -Sm2.7/-1 -Gred -P -O -K -T0 -C>> %F%

#Plot gempa utama
gawk "{print $8, $7}" maineq.dat | psxy -JM -R -Sa0.9c -W1 -Gyellow -P -O -K >> %F%

#Plot sinyal komponen E-W
psbasemap -JX13.0c/3.0c -R0/80/-200000/200000 -Ba20f10:"":/a100000f100000:"":SW --LABEL_FONT_SIZE=16p -P -O -K -Y-5 >> %F%
gawk "{print $1, $2}" leme.dat | psxy -W2/0/0/255 -JX13.0c/3.0c -R0/80/-200000/200000 -P -O -K >> %F%
echo 0.2 4.0 16 0 0 LT E-W > arahEW.dat
pstext -JX13.0c/3.0c -R0/10/0/4 arahEW.dat -P -O -K >> %F%

#Plot sinyal komponen N-S
psbasemap -JX13.0c/3.0c -R0/80/-200000/200000 -Ba20f10:"":/a100000f100000:"Amplitude (count)":SW --LABEL_FONT_SIZE=16p -P -O -K -Y-4>> %F%
gawk "{print $1, $2}" lemn.dat | psxy -W2/0/255/0 -JX13.0c/3.0c -R0/80/-200000/200000 -P -O -K >> %F%
echo 0.2 4.0 16 0 0 LT N-S > arahNS.dat
pstext -JX13.0c/3.0c -R0/10/0/4 arahNS.dat -P -O -K >> %F%

#Plot sinyal komponen U-D
psbasemap -JX13.0c/3.0c -R0/80/-200000/200000 -Ba20f10:"Time (s)":/a100000f100000:"":SW --LABEL_FONT_SIZE=16p -P -O -K -Y-4 >> %F%
gawk "{print $1, $2}" lemv.dat | psxy -W2/255/0/0 -JX13.0c/3.0c -R0/80/-200000/200000 -P -O -K >> %F%
echo 0.2 4.0 16 0 0 LT U-D > arahUD.dat
pstext -JX13.0c/3.0c -R0/10/0/4 arahUD.dat -P -O >> %F%

#Print peta
ps2raster %F% -Gc:\programs\gs\gs8.53\bin\gswin32c -Tj -P -Fsignal

Simpan dengan ekstensi ".bat". Setelah dieksekusi, hasilnya akan seperti berikut :
 

Selamat mencoba dan semoga bermanfaat. Salam Orang Indonesia...