Mencari Rupture Time Function berdasarkan Radiasi Energi Seismik


  • Mengubah seismogram ke catatan frekwensi tinggi untuk mencari waktu radiasi energi

SAC > macro energy.m

r  *II*KAPI*00*LHZ*SAC
rmean
rtrend
int
write KAPI_INT.SAC
r *II*KAPI*00*LHZ*SAC
rmean
rtrend
hp bu co 0.001 0.005
write KAPI_1H.SAC
r *II*KAPI*00*LHZ*SAC
rmean
rtrend
hp bu co 0.001 0.005
envelope
write KAPI_ENV.SAC
r *II*KAPI*00*LHZ*SAC
rmean
rtrend
hp bu co 0.001 0.005
envelope
smooth mean halfwidth 8
write KAPI_E8.SAC
r *II*KAPI*00*LHZ*SAC
rmean
rtrend
hp bu co 0.001 0.005       !  Highpass filter Butterworth frek cor 0.001 – 0.005 Hz
envelope                          !  Fungsi Hilbert transform    e = sqrt (x^2 + y^2).
                                        !  Kadang envelope harus terpisah dari program macro

smooth mean halfwidth 128    !  Algoritma diperhalus

write KAPI_E128.SAC
bd x
qdp off
r  *II*KAPI*00*LHZ*SAC KAPI_INT.SAC KAPI_1H.SAC KAPI_ENV.SAC KAPI_E8.SAC KAPI_E128.SAC
p1



Komentar

Postingan populer dari blog ini

Install SAC

Metode Magnet Sebagai Prekursor Gempabumi

Download data IRIS Wilber orfeus II