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
Posting Komentar