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 195.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...
Good....
ReplyDeleteMas, sy buat cross section qo ga keluar gambarnya ya? scriptnya sama, posisi long lat dan eventnya beda pakai data saya. utk peta dan fokalnya muncul tp cros sectionnya ga.. ini script CS sy:
ReplyDeleteset A=cross_sunda1.ps
echo 104 -8 A > teks.dat
echo 106 -8 B >> teks.dat
pstext teks.dat -R104/106/-7/-5 -JX-15 -K -N >> %A%
pscoupe focal.gmt -R104/106/-7/-5 -JX-15 -Aa104/-7/106/-5/90/10/0/200f -Ba200/a10 -Sm1 -G255/0/0 -N -T0 -O >> %A%
tolong dibantu mas, terima kasih
format datanya sama?
Deletesama semua mas utk format data, sy ikut spt yang mas buat.
Deleteini CS sj yg ga keluar, peta focalnya muncul
coba dirubah menjadi -Aa104/-7/106/-5/90/1/0/200f
DeleteHalo mas..saya coba membuat cross section, tetapi 0 nya mulai dari bawah, tidak dari atas seperti cross section yang mas buat, membaliknya pake apa ya? terimakasih
ReplyDeleteBisa diakali dengan -JY, jadi nilainya dibikin negatif, misalnya -JY-10, atau nilai kedalamannya yang dibuat negatif
ReplyDelete