Prosedur Perhitungan Magnitudo Mwp menggunakan SAC
SAC> bd x ! Mulai plot dengan suatu grafik
SAC> qdp off ! Mematikan plot yang salah
SAC> r file x ! Baca seismogram file x
SAC> p1 ! Plot pada layar
SAC> transfer from polezero subtype KMI.PZ to none ! Membagi data dengan
sensitifitas instrumen
SAC> cut cr1 cr2 ! Potong seismogram dalam kisaran waktu
SAC> r ! Plot kembali
SAC> xlim off ! Plot berdasarkan sumbu x
SAC> rmean ! Menghapus rata-rata
SAC> rtr ! Menghapus trend linier
SAC> taper ! Menerapkan simetris masing-masing data
SAC> int ! Integral untuk mengubah kecepatan ke pergeseran
SAC> int ! Integral pergeseran
SAC> abs ! Absolut masing-masing data
SAC> mul 4 ! Perkalian 4
SAC> mul 3.1415926 ! Perkalian π
SAC> mul 3.5E+03 ! Perkalian kerapatan (densitas) batuan
SAC> mul apv3 ! Perkalian data kecepatan gelombang P (rumus Kanjo,
2006)
SAC> mul distance ! Perkalian jarak sebenarnya (km)
SAC> log10 ! Logaritma 10
SAC> sub 9.1 ! Pengurangan 9.1
SAC> div 1.5 ! Pembagian 1.5
SAC> add 0.2 ! Penambahan 0.2
SAC> ylim lim1 lim2 ! Menetapkan sumbu y dari lim1 ke lim2
SAC> p1 ! Plot pada layar
Akan muncul display Mw untuk stasiun KAPI. Yang diambil adalah magnitudo maksimum
melalui perintah lh DEPMAX = 7.461437e+00
Untuk mendapatkan harga Mwp sebagaimana yang digunakan oleh PTWC untuk peringatan
dini tsunami dimana Mwp = Mw + 0.2 . Digunakan perintah allMWP.m menambahkan dengan
sintaks berikut :
SAC> r MWRAYN
SAC> add 0.2
SAC> ylim 7 8.5
SAC> w MWPRAYN
SAC> r MWPRAYN
SAC> p1
SAC> lh
Komentar
Posting Komentar