Membuat Cross Section Mekanisme Fokal dengan GMT

Mulai postingan kali ini dan selanjutnya, saya akan menggunakan GMT versi tebaru yaitu GMT versi 5.1.0. Bagi yang masih menggunakan GMT versi 4.x.x kemungkinan akan mendapatkan hasil yang berbeda dan perlu sedikit modifikasi pada command anda agar mendapatkan hasil yang sama.

Kali ini saya akan berbagi cara untuk membuat cross section dari mekanisme fokal dengan GMT. Perintah yang digunakan untuk membt cross section ini adalah pscoupe. Langkah pertama downloadlah data mekanisme fokal, pada contoh ini saya mendapatkan data dari GlobalCMT.

Berikut adalah file yang berisi data mekanisme fokal :

95.48    3.57    31    9.56    -7.17    -2.39    -0.09    0.83    5.24    23    0    0    1
95.52    3.37    42    2.57    -2.05    -0.52    -0.02    -0.09    1.51    24    0    0    2
95.68    3.48    15    3.49    -2.03    -1.46    -3.05    0.3    -0.83    23    0    0    3
104.69    -6.38    15    4.87    -3.77    -1.1    5.74    -1.25    1.99    23    0    0    4
100.9    -4.44    15    3.14    -1.29    -1.85    2.45    -1.92    2.59    23    0    0    5
102.13    -5.13    44    -0.06    -9.59    9.66    -2.99    -6.54    -8.23    23    0    0    6
98.44    -0.71    35    1.97    -0.98    -1    2.33    -3.34    1.22    23    0    0    7
103.04    -4.19    15    2.92    -1.24    -1.68    2.56    -4.36    -1.19    23    0    0    8
101.93    -5.48    46    0.51    -0.28    -0.24    -0.25    -3.77    -7.86    23    0    0    9
107.87    -7.64    15    2.41    -3.3    0.89    7.07    -0.23    0.66    23    0    0    10

Simpan dengan nama foc.gmt

Selanjutnya untuk menampilkan mekanisme fokal pada peta dan garis cross sectionnya, ketikkan perintah berikut :

set F=psmeca.ps
pscoast -JM15 -R90/110/-10/10 -G100/255/0 -S200/255/255 -Dh -W0.1 -K -P > %F%
echo 93 5 A > garis.dat
echo 109 -9 B >> garis.dat
psxy garis.dat -JM -R -W10 -P -O -K -W2 >> %F%
pstext garis.dat -JM -D0.3 -R -P -O -K >> %F%
psmeca foc.gmt -JM15 -R90/110/-10/10 -Ba5f5/a5f5 -Sm1c -T0 -Gred -N -O >> %F%

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

 Selanjutnya untuk membuat cross sectionnya, ketikkan perintah berikut :

set F=pscoupe.ps
echo 93 -10 A > teks.dat
echo 109 -10 B >> teks.dat
pstext teks.dat -R93/109/-9/5 -JX-15 -K -N > %F%
pscoupe foc.gmt -R93/109/-9/5 -JX-15 -Aa93/5/109/-9/90/1/0/50f -Ba200/a10 -Sm1 -G255/0/0 -N -T0 -O >> %F%

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

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