Gempa Tsunami GMT SAC

$ ssh -X sugeng@192.168.15.102 password:123456 [sugeng@seiscomp ~]$ cd QUAKEDURATION/ [sugeng@seiscomp QUAKEDURATION]$ ls common.py   data.inp     duration.pyc     ttimes.lst common.pyc  duration.py  runduration.pyc [sugeng@seiscomp QUAKEDURATION]$ python2.7 Python 2.7 (r27:82500, May 20 2013, 20:31:27) [GCC 4.1.2 20080704 (Red Hat 4.1.2-52)] on linux2 Type "help", "copyright", "credits" or "license" for more information. >>> import duration as dr >>> hd = dr.Duration() >>> hd.setDirSAC('/home/sugeng/TITIPANSUGENG/1_bku2000/data') >>> hd.loadSAC() OK >>> hd.setDirSAC('/home/sugeng/TITIP/1_bku2000/data') >>> hd.loadSAC() . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . OK >>> hd.stage1() >>> hd.plot('all')  SEISMIC ANALYSIS CODE [01/04/2013 (Version 101.4)]  Copyright 1995 Reg

#!/bin/csh gmtset ANNOT_FONT_SIZE_PRIMARY 12p HEADER_FONT_SIZE 18p PLOT_DEGREE_FORMAT ddd:mm:ssF #set region = 92/106/-7.5/6 set region = 96/103/-6/1 set size = M12c set psfile = xxxx_all_3.ps set grdfile=./indonesia.nc set cptfile=color.cpt makecpt -Csealand -Z > $cptfile set AZ = 55 set pointboxlon_1 = 96.5 set pointboxlat_1 = -4.5 echo $pointboxlon_1 $pointboxlat_1 S1 >! S1.d awk '{print $1+5.5, $2+4, $3"`" }' S1.d >! S11.d cat S1.d S11.d >! garis_S11.d #Draw bathymetry and basemap grdgradient $grdfile -A30/270 -Gintens.grd -Nt0.30 -V gmtset ANOT_FONT_SIZE 10 grdimage $grdfile -R$region -J$size -C$cptfile -Iintens.grd -B2g1/2g1wnSE -K -P -X0.7 -Y3.5 >! $psfile pscoast -R$region -J$size -Dh -W1 -K -P -O >> $psfile # Draw tectonical structures  ;   trech3 contains jasum & nusampu awk -F, '{ print $1, $2}' trench.gmt    | psxy -R$region -J$size -Sf0.20i/0.04irt -W1.5/black -Gblack -M -N -O -K >> $psfile

#!/bin/csh -f pscoast -R122/128/-2/4 -JM15c -Ba2g1 -Ggreen3 -S140/235/255 -Dh -N1-W1 -P -K > focalbunaken19762010.ps psmeca bunaken.gmt -R -JM -Sm0.5c/-1 -G0/0/255 -P -O -K >> focalbunaken19762010.ps pstext -R -JM -P -O -K << EOF >> focalbunaken19762010.ps 124.8 1.53 16 0 1 LT Bunaken EOF psxy -R -JM -Sa0.7c -Gred -W1 -P -O << EOF >> focalbunaken19762010.ps 124.8 1.53 EOF

bd sgf # start sac graphics file (sgf) p1 sgftops f001.sgf COCO_enr.ps # convert to sgf to ps file quit from sac rm -f *.sgf # remove sgf files gs COCO_enr.ps Membuat label judul, absis, dan ordinat ================================ title on xlabel "Time @(s@)" ylabel "Velocity @(nm/s@)" * set the xlabel and ylabel ylim all * set the y-axis to be the same xlim o o +8 set the time axis to be from origin time to 8 sec after the origin time Gabungkan tiga komponen sekaligus ================================ rm -f *.sgf setbb sacfile "$sta$.*" #setbb sacfile "$sta$.[zne]" setbb psfile "$sta$_3c.eps" #set three component data qdp off r COCO.sac rmean bp p 2 n 4 c 1 10 #4th order bandpass butterworth filter with 1 - 10 Hz setbb ds '( concatenate ' &1,kstnm ' ' &1,kzdate ' ' &1,kztime ' ) ' evaluate to r &1,dist title '( concatenate ' %ds% ' ', dis

SAC > qdp off SAC > bd x SAC > synch # memperjelas sinyal to synchronice the time all components SAC > interpolate delta 0.5 SAC > p1 SAC > lh DELTA : 5E-1 sebelumnya DELTA : 1E+0 NyQIST frequency = 1/(2*DELTA) = 1/ (2*0.5) = 1.0 NyQIST frequency adalah maksimum frequency. DELTA adalah sampling per second SAC > decimate 2 Diperlukan untuk mengurangi harga filter untuk men-downsampling sehingga benar-benar mendapatkan harga frekwensi corner yang layak. Harga hanya diberikan 2 – 7. Prosedur untuk phase picking gelombang P dan S =============================================== SAC > r coco.sac SAC > lh DIST Ubah dalam bentuk degree DEG = DIST / 111.267 Lihat dalam Table/grafik Jeffrey Bullens untuk DEGREE berapa waktu kedatangan gel P dan S Misal : COCO DEG = 1.201442e+03 / 111.267 = 10.8 dalam tabel didapatkan waktu gel P = 3.70 dan S = 5.80 Kemudian observasi sinyal secara manual untuk mendapatkan nilai waktu sebenarnya SAC > bd

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 # Bandpass # BANDPASS BUTTER CORNER 0.1 0.4 NPOLES 2 PASSES 1 TRANBW 0.3 ATTEN 30.

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

1. Install SAC terlebih dahulu dengan mengetikkan : % echo $SACAUX /usr/local/geophysics/sac/sac/aux % echo $PATH 2. Tambahkan dalam lingkungan (environment) bashrc ${HOME}/.bashrc Atau ${HOME} vi bashrc # added for Seismic Analysis Code (SAC) export SACAUX=/usr/local/geophysics/sac/sac/aux export PATH=/usr/local/geophysics/bin:/usr/local/geophysics/sac/ local_bin:/usr/local/geophysics/sac/sac/bin:$PATH setenv SACHOME /usr/local/geophysics/sac/sac setenv SACAUX $SACHOME/aux set path=($SACHOME/bin $path) 3. Buka file package dari SAC % tar zxvf /usr/local/geophysics/sac/Sac_Tutorial_2009.tar.gz % cd SAC_Tutorial/SAC_Tutorial_WF % wget http://geophysics.eas.gatech.edu/people/zpeng/Teaching/ Sac_Tutorial_2009.tar.gz % tar zxvf Sac_Tutorial_2009.tar.gz % cd SAC_Tutorial/SAC_Tutorial_WF 4. Setelah sukses instalasi, ketikkan ls untuk menampilkan file SAC komponen LHZ jaringan II Perintah ls : sugeng@sugeng-laptop:~/DATE$ ls *II*00*LHZ*SAC 2010.298.14.43.17.0695_

Hasil-hasil W phase oleh Duputel dkk (2011) menggunakan jarak Δ<90° dari jaringan stasiun FDSN, GSN and STS1 global virtual networks (most stations are from II, IU and G). http://wphase.unistra.fr/MGE65/index.html

Lecture http://www.gm.univ-montp2.fr/CRYSTAL2PLATE/Crystal2Plate_Seismic_course.html Zhigang Peng http://geophysics.eas.gatech.edu/classes/SAC/ STA/LTA to minimize noise then easy define the phase http://psn.quake.net/sdrdocs/stalta.htm