Rabu, 15 Desember 2010
Sabtu, 04 Desember 2010
Cross of Sumatra
#!/bin/csh
# & rubah menjadi <
gmtset ANNOT_FONT_SIZE_PRIMARY 10p HEADER_FONT_SIZE 18p PLOT_DEGREE_FORMAT ddd:mm:ssF
set region = 90/107/-10/7
set size = M11c
set psfile = 4cross.ps
set grdfile=./sumatra.nc
set cptfile=color.cpt
makecpt -Cglobe -Z > $cptfile
set AZ = 45
#Cross box near subduction
set boxlon = 99.0
set boxlat = -4.0
set boxdep = 0
echo $boxlon $boxlat $boxdep >! box.d
#Draw bathymetry and basemap
grdgradient $grdfile -A30/270 -Gintens.grd -Nt0.30 -V
grdimage $grdfile -R$region -J$size -C$cptfile -Iintens.grd -K -P -Ba4f1NWse -X2 -Y17 >! $psfile
gmtset ANOT_FONT_SIZE 9
pscoast -R$region -J$size -C$cptfile -Dl -W1 -P -Lf92.0/-9/40/60 -Tf92.0/-8.65/10i/3 -O -K >> $psfile
#Draw cross centre
awk -F, '{ print $1, $2}' box.d | psxy -J$size -R$region -Sx.5/.5 -W2::0 -P -O -K >> $psfile
# Draw tectonical structures
awk -F, '{ print $1, $2}' ../focalproposal/trench.gmt | psxy -R$region -J$size -Gdarkblue -Sc0.05 -O -K >> $psfile
awk -F, '{ print $1, $2}' ../focalproposal/transform.gmt | psxy -R$region -J$size -Gdarkblue -Sc0.05 -O -K >> $psfile
awk -F, '{ print $1, $2}' ../focalproposal/ridge.gmt | psxy -R$region -J$size -Gdarkblue -Sc0.05 -O -K >> $psfile
#Draw usgs data
awk -F, '{ print $1, $2}' ../focalproposal/gempa_kecil.lst | psxy -J$size -R$region -Sc0.1c -W0.1 -H -P -O -K >> $psfile
awk -F, '{ print $1, $2}' ../focalproposal/gempa_sedang.lst | psxy -J$size -R$region -Sc0.15c -W0.15 -H -P -O -K >> $psfile
awk -F, '{ print $1, $2}' ../focalproposal/gempa_besar.lst | psxy -J$size -R$region -Sc0.2c -W0.2 -H -P -O -K >> $psfile
#Projection with the center
awk -F, '{ print $1, $2, $3}' ../focalproposal/trench.gmt | project -Q -C$boxlon/$boxlat -A$AZ -Fxyzpqrs -L-10/550 -W-50/50 >! xyzpqrs_trench
awk -F, '{ print $1, $2, $3}' ../focalproposal/transform.gmt | project -Q -C$boxlon/$boxlat -A$AZ -Fxyzpqrs -L-10/550 -W-50/50 >! xyzpqrs_transform
awk -F, '{ print $1, $2, $3}' ../focalproposal/ridge.gmt | project -Q -C$boxlon/$boxlat -A$AZ -Fxyzpqrs -L-10/550 -W-50/50 >! xyzpqrs_ridge
awk -F, '{ print $1, $2, $3}' ../focalproposal/gempa_kecil.lst | project -Q -C$boxlon/$boxlat -A$AZ -Fxyzpqrs -L-10/550 -W-50/50 >! xyzpqrs_kecil
awk -F, '{ print $1, $2, $3}' ../focalproposal/gempa_sedang.lst | project -Q -C$boxlon/$boxlat -A$AZ -Fxyzpqrs -L-10/550 -W-50/50 >! xyzpqrs_sedang
awk -F, '{ print $1, $2, $3}' ../focalproposal/gempa_besar.lst | project -Q -C$boxlon/$boxlat -A$AZ -Fxyzpqrs -L-10/550 -W-50/50 >! xyzpqrs_besar
#Plot projected earthquakes (red dots) and mainshock (yellow star) on map
awk '{print $1, $2}' xyzpqrs_kecil | psxy -J$size -R$region -Sc0.1c -W0.1/red -P -O -K >> $psfile
awk '{print $1, $2}' xyzpqrs_sedang | psxy -J$size -R$region -Sc0.15c -W0.15/red -P -O -K >> $psfile
awk '{print $1, $2}' xyzpqrs_besar | psxy -J$size -R$region -Sc0.2c -W0.2/red -P -O -K >> $psfile
#Mainshock will appear in front of all earthquakes
#BE CAREFULL !! MAINSHOCK FILE MUST CONSIST OF LON LAT DEPTH
awk '{print $1, $2}' ../focalproposal/cobi.d | psxy -J$size -R$region -Sa0.4c -Gyellow -Wthin -P -O -K >> $psfile
project -Q -C$boxlon/$boxlat -A$AZ ../focalproposal/cobi.d -Fxyzpqrs -L-10/550 -W-50/50 >! xyzpqrs_cobi
#Focalmec lines taken from true epicenter
sed -n '2,3p' ../focalproposal/test.txt | psxy -R$region -J$size -Ggray -W1 -O -K >> $psfile
sed -n '5,6p' ../focalproposal/test.txt | psxy -R$region -J$size -Ggray -W1 -O -K >> $psfile
sed -n '8,9p' ../focalproposal/test.txt | psxy -R$region -J$size -Ggray -W1 -O -K >> $psfile
....................................
#Modified location of Global CMT's focalmec
psmeca ../focalproposal/focaltop_cmt.gmt -R$region -J$size -Sm0.25c/-1 -Gred -W1 -O -K >> $psfile
pstext ../focalproposal/cibi.txt -R$region -J$size -Gblack -O -K >> $psfile
#Drawing box and line of cross section on map
psxy box.txt -J$size -R$region -W0.5 -P -O -K >> $psfile
awk '{print $1, $2}' line.txt | psxy -J$size -R$region -W0.5 -P -O -K >> $psfile
awk '{print $1, $2, 14, 0, 1, "LT", $3 }' line.txt | pstext -J$size -R$region -Gdarkblue -P -O -K >> $psfile
#DRAW CROSSEC DIAGRAM
psbasemap -JX9.5/-5 -R0/550/0/300 -B50/30NWes -P -O -K -Y-7 >> $psfile
pstext -JX -R -Gdarkblue -P -O -K <
0 0 14 0 1 LT C
520 0 14 0 1 LT C'
EOF
#Plot stars of structures, earthquakes, and mainshocks on crossec diagram
awk '{print $4, $3}' xyzpqrs_trench | psxy -JX -R -Gdarkblue -Sc0.2 -P -O -K >> $psfile
awk '{print $4, $3}' xyzpqrs_transform | psxy -JX -R -Gdarkblue -Sc0.2 -P -O -K >> $psfile
awk '{print $4, $3}' xyzpqrs_ridge | psxy -JX -R -Gdarkblue -Sc0.2 -P -O -K >> $psfile
awk '{print $4, $3}' xyzpqrs_kecil | psxy -JX -R -Sc0.1c -W0.1/red -P -O -K >> $psfile
awk '{print $4, $3}' xyzpqrs_sedang | psxy -JX -R -Sc0.15c -W0.15/red -P -O -K >> $psfile
awk '{print $4, $3}' xyzpqrs_besar | psxy -JX -R -Sc0.2c -W0.2/red -P -O -K >> $psfile
awk '{print $4, $3}' xyzpqrs_cobi | psxy -JX -R -Sa0.4c -Gyellow -Wthin -P -O -K >> $psfile
psscale -C$cptfile -Iintens.grd -D0.4c/1/10/0.15 -B1000 -P -O -K -X11.5c -Y12c >> $psfile
# Create legend of magnitude size
cat << EOF >! magnitude.legend
N 3
S 0.025i c 0.04i red 0.010p 0.1i mag 3.5-4.5
S 0.020i c 0.06i red 0.015p 0.1i mag 4.5-5.5
S 0.015i c 0.08i red 0.020p 0.1i mag 5.5-9.5
EOF
pslegend magnitude.legend -Dx0.1i/-0.2i/6.8i/1.5i/TC -Jx0.4i -R0/8/0/8 -O -X-4c -Y-5 >> $psfile
\rm -f *.d *.xy xyzpqrs_*
gs -sDEVICE=x11 $psfile
Label:
cross section,
focal mechanism,
gempa,
GMT,
sumatra
Jumat, 03 Desember 2010
KATALOG TSUNAMI INDONESIA 1990 - 2010
#!/bin/csh
gmtset ANNOT_FONT_SIZE_PRIMARY 10p HEADER_FONT_SIZE 18p PLOT_DEGREE_FORMAT ddd:mm:ssF
set area = 90/140/-15/9
set psfile = top.ps
set grdfile=./indo.nc
set cptfile=color.cpt
makecpt -Cglobe -Z > $cptfile
grdgradient $grdfile -A30/270 -Gintens.grd -Nt0.30 -V
grdimage $grdfile -R$area -JM18c -C$cptfile -Iintens.grd -P -K -X1.5 -Y8 >! $psfile
gmtset ANOT_FONT_SIZE 9
pscoast -R$area -JM18c -Df -W2,100 -P -Ba4f2/a4f2NWSE -O -K >> $psfile
awk -F, '{ print $1, $2}' ./trench.gmt | psxy -R$area -J -Gdarkblue -Sc0.05 -O -K >> $psfile
awk -F, '{ print $1, $2}' ./transform.gmt | psxy -R$area -J -Gdarkblue -Sc0.05 -O -K >> $psfile
awk -F, '{ print $1, $2}' ./ridge.gmt | psxy -R$area -J -Gdarkblue -Sc0.05 -O -K >> $psfile
awk -F, '{ print $1, $2}' ./cobi.d | psxy -R$area -J -Gyellow -W1 -Sa0.2 -O -K >> $psfile
sed -n '2,3p' test.txt | psxy -R$area -JM -Ggray -W1 -O -K >> $psfile
sed -n '5,6p' test.txt | psxy -R$area -JM -Ggray -W1 -O -K >> $psfile
sed -n '8,9p' test.txt | psxy -R$area -JM -Ggray -W1 -O -K >> $psfile
sed -n '11,12p' test.txt | psxy -R$area -JM -Ggray -W1 -O -K >> $psfile
sed -n '14,15p' test.txt | psxy -R$area -JM -Ggray -W1 -O -K >> $psfile
sed -n '17,18p' test.txt | psxy -R$area -JM -Ggray -W1 -O -K >> $psfile
sed -n '20,21p' test.txt | psxy -R$area -JM -Ggray -W1 -O -K >> $psfile
sed -n '23,24p' test.txt | psxy -R$area -JM -Ggray -W1 -O -K >> $psfile
sed -n '26,27p' test.txt | psxy -R$area -JM -Ggray -W1 -O -K >> $psfile
sed -n '29,30p' test.txt | psxy -R$area -JM -Ggray -W1 -O -K >> $psfile
sed -n '32,33p' test.txt | psxy -R$area -JM -Ggray -W1 -O -K >> $psfile
sed -n '35,36p' test.txt | psxy -R$area -JM -Ggray -W1 -O -K >> $psfile
sed -n '38,39p' test.txt | psxy -R$area -JM -Ggray -W1 -O -K >> $psfile
sed -n '41,42p' test.txt | psxy -R$area -JM -Ggray -W1 -O -K >> $psfile
sed -n '44,45p' test.txt | psxy -R$area -JM -Ggray -W1 -O -K >> $psfile
sed -n '47,48p' test.txt | psxy -R$area -JM -Ggray -W1 -O -K >> $psfile
sed -n '50,51p' test.txt | psxy -R$area -JM -Ggray -W1 -O -K >> $psfile
sed -n '53,54p' test.txt | psxy -R$area -JM -Ggray -W1 -O -K >> $psfile
sed -n '56,57p' test.txt | psxy -R$area -JM -Ggray -W1 -O -K >> $psfile
sed -n '59,60p' test.txt | psxy -R$area -JM -Ggray -W1 -O -K >> $psfile
sed -n '62,63p' test.txt | psxy -R$area -JM -Ggray -W1 -O -K >> $psfile
sed -n '65,66p' test.txt | psxy -R$area -JM -Ggray -W1 -O -K >> $psfile
sed -n '68,69p' test.txt | psxy -R$area -JM -Ggray -W1 -O -K >> $psfile
sed -n '71,72p' test.txt | psxy -R$area -JM -Ggray -W1 -O -K >> $psfile
sed -n '74,75p' test.txt | psxy -R$area -JM -Ggray -W1 -O -K >> $psfile
sed -n '77,78p' test.txt | psxy -R$area -JM -Ggray -W1 -O -K >> $psfile
sed -n '80,81p' test.txt | psxy -R$area -JM -Ggray -W1 -O -K >> $psfile
sed -n '83,84p' test.txt | psxy -R$area -JM -Ggray -W1 -O -K >> $psfile
sed -n '86,87p' test.txt | psxy -R$area -JM -Ggray -W1 -O -K >> $psfile
sed -n '89,90p' test.txt | psxy -R$area -JM -Ggray -W1 -O -K >> $psfile
psmeca focaltop_cmt.gmt -R$area -JM -Sm0.25c/-1 -Gred -W1 -O -K >> $psfile
pstext cibi.txt -R$area -JM -Gblack -O -K >> $psfile
psscale -C$cptfile -Iintens.grd -D0.4c/1/2/0.15 -B2000 -P -O -X20c -Y5c >> $psfile
Cut command linux
cut
cut options [files]Cut out selected columns or fields from one or more files.
In the following options, list is a sequence of integers. Use a comma between separate values, and a hyphen to specify a range (e.g., 1-10,15,20 or 50-). See also paste and join.
Options
-b list, --bytes listSpecify list of positions; only bytes in these positions will be printed.-c list, --characters list
Cut the column positions identified in list. Column numbers start with 1.-d c, --delimiter c
Use with -f to specify field delimiter as character c (default is tab); special characters (e.g., a space) must be quoted.-f list, --fields list
Cut the fields identified in list.-n
Don't split multibyte characters.-s, --only-delimited
Use with -f to suppress lines without delimiters.--output-delimiter=string
Use string as the output delimiter. By default, the output delimiter is the same as the input delimiter.--help
Print help message and exit.--version
Print version information and exit.
Examples
Extract usernames and real names from /etc/passwd:cut -d: -f1,5 /etc/passwd
Find out who is logged on, but list only login names:
who | cut -d" " -f1
Cut characters in the fourth column of file, and paste them back as the first column in the same file:
cut -c4 file | paste - file
-
cut -c1-3 text.txt
Output:
Thi
Cut the first three letters from the above line.
-
cut -d, -f1,2 text.txt
Output:This is, an example program
The above command is used to split the fields using delimiter and cut the first two fields.
Sed command linux
sed
(=stream editor) I use sed to filter text files. The pattern to match is typically included between a pair of slashes // and quoted.
sed G myfile.txt > newfile.txt
(=stream editor) I use sed to filter text files. The pattern to match is typically included between a pair of slashes // and quoted.
For example, to print lines containing the string "1024", I may use:
cat filename | sed -n '/1024/p'Here, sed filters the output from the cat command. The option "-n" tells sed to block all the incoming lines but those explicitly matching my expression. The sed action on a match is "p"= print.Another example, this time for deleting selected lines:
cat filename | sed '/.*o$/d' > new_file
In this example, lines ending the an "o" will be deleted. I used a regular expression for matching any string followed by an "o" and the end of the line. The output (i.e., all lines but those ending with "d") is directed to new_file.Another example. To search and replace, I use the sed 's' action, which comes in front of two expressions:
cat filename | sed 's/string_old/string_new/' > newfile
A shorter form for the last command is:To insert a text from a text file into an html file, I may use a script containing:
sed 's/string_old/string_new/' filename > newfile
sed '/text_which_is_a_placeholder_in__my_html_file/r text_file_to_insert.txt' index_master_file.html > index.htmll
sed G myfile.txt > newfile.txt
In the above example using the sed command with G would double space the file myfile.txt and output the results to the newfile.txt.
sed = myfile.txt | sed 'N;s/\n/\. /'The above example will use the sed command to output each of the lines in myfile.txt with the line number followed by a period and a space before each line. As done with the first example the output could be redirected to another file using > and the file name.
sed 's/test/example/g' myfile.txt > newfile.txtOpens the file myfile.txt and searches for the word "test" and replaces every occurrence with the word "example".
sed -n '$=' myfile.txtAbove this command count the number of lines in the myfile.txt and output the results.
Rabu, 01 Desember 2010
Cross section Sumbagsel
#!/bin/csh
gmtset ANNOT_FONT_SIZE_PRIMARY 10p HEADER_FONT_SIZE 18p PLOT_DEGREE_FORMAT ddd:mm:ssF
set cptfile = color.cpt
set psfile = south_sumatra.ps
set size = M11c
set region = 97/103/-4/3
makecpt -Cglobe -Z > $cptfile
set AZ = 45
#Padang_30/9/2010
set Alon = 98.0
set Alat = -2.5
set Adep = 0
echo $Alon $Alat $Adep >! A.d
set lon_X = 99.67
set lat_X = -0.79
set dep_X = 77.8
set lon_Xplus4 = 101.67
set lat_Xplus0 = -0.79
set lat_Xmin1 = -1.29
echo $lon_X $lat_X $dep_X >! mainshock_X.d
#Draw bathymetry and basemap
grdgradient indo.nc -A30/270 -Gintens.grd -Nt0.30 -V
grdimage indo.nc -R$region -J$size -C$cptfile -Iintens.grd -K -P -Ba1g0.5NWse -X2 -Y10 > $psfile
pscoast -R$region -J$size -P -C$cptfile -Dl -W1 -O -K -Lf108/-6/25/200 >> $psfile
#Draw cross centre
awk -F, '{ print $1, $2}' A.d | psxy -J$size -R$region -Sx.5/.5 -W2::0 -P -O -K >> $psfile
# Draw tectonical structures
awk -F, '{ print $1, $2}' trench.gmt | psxy -R$region -J$size -Gdarkblue -St0.2 -O -K >> $psfile
awk -F, '{ print $1, $2}' transform.gmt | psxy -R$region -J$size -Gdarkblue -St0.2 -O -K >> $psfile
awk -F, '{ print $1, $2}' ridge.gmt | psxy -R$region -J$size -Gdarkblue -St0.2 -O -K >> $psfile
#Draw usgs data
awk -F, '{ print $1, $2}' gempa_kecil.lst | psxy -J$size -R$region -Sc0.1c -W0.1 -H -P -O -K >> $psfile
awk -F, '{ print $1, $2}' gempa_sedang.lst | psxy -J$size -R$region -Sc0.15c -W0.15 -H -P -O -K >> $psfile
awk -F, '{ print $1, $2}' gempa_besar.lst | psxy -J$size -R$region -Sc0.2c -W0.2 -H -P -O -K >> $psfile
#Projection with the center
awk -F, '{ print $1, $2, $3}' trench.gmt | \
project -Q -C$Alon/$Alat -A$AZ -Fxyzpqrs -L-10/550 -W-50/50 >! xyzpqrs_trench
awk -F, '{ print $1, $2, $3}' transform.gmt | \
project -Q -C$Alon/$Alat -A$AZ -Fxyzpqrs -L-10/550 -W-50/50 >! xyzpqrs_transform
awk -F, '{ print $1, $2, $3}' ridge.gmt | \
project -Q -C$Alon/$Alat -A$AZ -Fxyzpqrs -L-10/550 -W-50/50 >! xyzpqrs_ridge
awk -F, '{ print $1, $2, $3}' gempa_kecil.lst | \
project -Q -C$Alon/$Alat -A$AZ -Fxyzpqrs -L-10/550 -W-50/50 >! xyzpqrs_kecil
awk -F, '{ print $1, $2, $3}' gempa_sedang.lst | \
project -Q -C$Alon/$Alat -A$AZ -Fxyzpqrs -L-10/550 -W-50/50 >! xyzpqrs_sedang
awk -F, '{ print $1, $2, $3}' gempa_besar.lst | \
project -Q -C$Alon/$Alat -A$AZ -Fxyzpqrs -L-10/550 -W-50/50 >! xyzpqrs_besar
project -Q -C$Alon/$Alat -A$AZ mainshock_X.d -Fxyzpqrs -L-10/550 -W-50/50 >! xyzpqrs_X
#Plot earthquakes projected (red dots) and mainshock (yellow star) on map
awk '{print $1, $2}' xyzpqrs_kecil | psxy -J$size -R$region -Sc0.1c -W0.1/red -P -O -K >> $psfile
awk '{print $1, $2}' xyzpqrs_sedang | psxy -J$size -R$region -Sc0.15c -W0.15/red -P -O -K >> $psfile
awk '{print $1, $2}' xyzpqrs_besar | psxy -J$size -R$region -Sc0.2c -W0.2/red -P -O -K >> $psfile
awk '{print $1, $2}' mainshock_X.d | psxy -J$size -R$region -Sa0.4c -Gyellow -Wthin -P -O -K >> $psfile
#Focalmec lines for true epicentres
psxy -R$region -J$size -P -O -K -W2 <> $psfile
$lon_X $lat_X
$lon_Xplus4 $lat_Xplus0
EOF
# Modified focalmec sphere from Global CMT
psmeca -R$region -J$size -Sm0.5c/-1 -Gred -Wthin -O -K <> $psfile
$lon_Xplus4 $lat_Xplus0 $dep_X 1.76 -0.76 -0.99 0.66 -0.99 -1.94 27 X Y 200909301016A
EOF
pstext -R$region -J$size -Gblack -Wwhite -O -K <> $psfile
$lon_Xplus4 $lat_Xmin1 9 0 1 LC 2009/09/30
EOF
#Create points of cross section line on map
psxy -J$size -R$region -W0.5 -P -O -K <> $psfile
98.0 -2.5
101.5 1.0
EOF
psxy -J$size -R$region -W0.5 -P -O -K <> $psfile
98.318 -2.818
101.813 0.677
101.177 1.313
97.682 -2.182
98.318 -2.818
EOF
pstext -J$size -R$region -Gdarkblue -P -O -K <> $psfile
98.0 -2.5 14 0 1 LT A
EOF
#DRAW CROSSEC DIAGRAM
psbasemap -JX9.5/-5 -R0/550/0/300 -B50/30NWes -P -O -K -Y-7 >> $psfile
pstext -JX -R -Gdarkblue -P -O -K <> $psfile
0 0 14 0 1 LT A
EOF
#Plot stars of structures, earthquakes, and mainshocks on crossec diagram
awk '{print $4, $3}' xyzpqrs_trench | \
psxy -JX -R -Gdarkblue -Sc0.2 -P -O -K >> $psfile
awk '{print $4, $3}' xyzpqrs_transform | \
psxy -JX -R -Gdarkblue -Sc0.2 -P -O -K >> $psfile
awk '{print $4, $3}' xyzpqrs_ridge | \
psxy -JX -R -Gdarkblue -Sc0.2 -P -O -K >> $psfile
awk '{print $4, $3}' xyzpqrs_kecil | \
psxy -JX -R -Sc0.1c -W0.1/red -P -O -K >> $psfile
awk '{print $4, $3}' xyzpqrs_sedang | \
psxy -JX -R -Sc0.15c -W0.15/red -P -O -K >> $psfile
awk '{print $4, $3}' xyzpqrs_besar | \
psxy -JX -R -Sc0.2c -W0.2/red -P -O -K >> $psfile
awk '{print $4, $3}' xyzpqrs_X | \
psxy -JX -R -Sa0.4c -Gyellow -Wthin -P -O -K >> $psfile
psscale -C$cptfile -Iintens.grd -D0.4c/1/10/0.15 -B1000 -P -O -K -X11.5c -Y12c >> $psfile
# Create legend of magnitude size
cat << EOF >! magnitude.legend
N 3
S 0.025i c 0.025i red 0.01p 0.050i mag 3-5
S 0.050i c 0.050i red 0.01p 0.075i mag 5-6
S 0.075i c 0.075i red 0.01p 0.100i mag 6-9
>
EOF
pslegend magnitude.legend -Dx0.1i/-0.2i/6.8i/1.5i/TC -Jx0.4i -R0/8/0/8 -O -X-4c -Y-5 >> $psfile
\rm -f *.d *.xy xyzpqrs_*
gs -sDEVICE=x11 $psfile
gmtset ANNOT_FONT_SIZE_PRIMARY 10p HEADER_FONT_SIZE 18p PLOT_DEGREE_FORMAT ddd:mm:ssF
set cptfile = color.cpt
set psfile = south_sumatra.ps
set size = M11c
set region = 97/103/-4/3
makecpt -Cglobe -Z > $cptfile
set AZ = 45
#Padang_30/9/2010
set Alon = 98.0
set Alat = -2.5
set Adep = 0
echo $Alon $Alat $Adep >! A.d
set lon_X = 99.67
set lat_X = -0.79
set dep_X = 77.8
set lon_Xplus4 = 101.67
set lat_Xplus0 = -0.79
set lat_Xmin1 = -1.29
echo $lon_X $lat_X $dep_X >! mainshock_X.d
#Draw bathymetry and basemap
grdgradient indo.nc -A30/270 -Gintens.grd -Nt0.30 -V
grdimage indo.nc -R$region -J$size -C$cptfile -Iintens.grd -K -P -Ba1g0.5NWse -X2 -Y10 > $psfile
pscoast -R$region -J$size -P -C$cptfile -Dl -W1 -O -K -Lf108/-6/25/200 >> $psfile
#Draw cross centre
awk -F, '{ print $1, $2}' A.d | psxy -J$size -R$region -Sx.5/.5 -W2::0 -P -O -K >> $psfile
# Draw tectonical structures
awk -F, '{ print $1, $2}' trench.gmt | psxy -R$region -J$size -Gdarkblue -St0.2 -O -K >> $psfile
awk -F, '{ print $1, $2}' transform.gmt | psxy -R$region -J$size -Gdarkblue -St0.2 -O -K >> $psfile
awk -F, '{ print $1, $2}' ridge.gmt | psxy -R$region -J$size -Gdarkblue -St0.2 -O -K >> $psfile
#Draw usgs data
awk -F, '{ print $1, $2}' gempa_kecil.lst | psxy -J$size -R$region -Sc0.1c -W0.1 -H -P -O -K >> $psfile
awk -F, '{ print $1, $2}' gempa_sedang.lst | psxy -J$size -R$region -Sc0.15c -W0.15 -H -P -O -K >> $psfile
awk -F, '{ print $1, $2}' gempa_besar.lst | psxy -J$size -R$region -Sc0.2c -W0.2 -H -P -O -K >> $psfile
#Projection with the center
awk -F, '{ print $1, $2, $3}' trench.gmt | \
project -Q -C$Alon/$Alat -A$AZ -Fxyzpqrs -L-10/550 -W-50/50 >! xyzpqrs_trench
awk -F, '{ print $1, $2, $3}' transform.gmt | \
project -Q -C$Alon/$Alat -A$AZ -Fxyzpqrs -L-10/550 -W-50/50 >! xyzpqrs_transform
awk -F, '{ print $1, $2, $3}' ridge.gmt | \
project -Q -C$Alon/$Alat -A$AZ -Fxyzpqrs -L-10/550 -W-50/50 >! xyzpqrs_ridge
awk -F, '{ print $1, $2, $3}' gempa_kecil.lst | \
project -Q -C$Alon/$Alat -A$AZ -Fxyzpqrs -L-10/550 -W-50/50 >! xyzpqrs_kecil
awk -F, '{ print $1, $2, $3}' gempa_sedang.lst | \
project -Q -C$Alon/$Alat -A$AZ -Fxyzpqrs -L-10/550 -W-50/50 >! xyzpqrs_sedang
awk -F, '{ print $1, $2, $3}' gempa_besar.lst | \
project -Q -C$Alon/$Alat -A$AZ -Fxyzpqrs -L-10/550 -W-50/50 >! xyzpqrs_besar
project -Q -C$Alon/$Alat -A$AZ mainshock_X.d -Fxyzpqrs -L-10/550 -W-50/50 >! xyzpqrs_X
#Plot earthquakes projected (red dots) and mainshock (yellow star) on map
awk '{print $1, $2}' xyzpqrs_kecil | psxy -J$size -R$region -Sc0.1c -W0.1/red -P -O -K >> $psfile
awk '{print $1, $2}' xyzpqrs_sedang | psxy -J$size -R$region -Sc0.15c -W0.15/red -P -O -K >> $psfile
awk '{print $1, $2}' xyzpqrs_besar | psxy -J$size -R$region -Sc0.2c -W0.2/red -P -O -K >> $psfile
awk '{print $1, $2}' mainshock_X.d | psxy -J$size -R$region -Sa0.4c -Gyellow -Wthin -P -O -K >> $psfile
#Focalmec lines for true epicentres
psxy -R$region -J$size -P -O -K -W2 <
$lon_X $lat_X
$lon_Xplus4 $lat_Xplus0
EOF
# Modified focalmec sphere from Global CMT
psmeca -R$region -J$size -Sm0.5c/-1 -Gred -Wthin -O -K <
$lon_Xplus4 $lat_Xplus0 $dep_X 1.76 -0.76 -0.99 0.66 -0.99 -1.94 27 X Y 200909301016A
EOF
pstext -R$region -J$size -Gblack -Wwhite -O -K <
$lon_Xplus4 $lat_Xmin1 9 0 1 LC 2009/09/30
EOF
#Create points of cross section line on map
psxy -J$size -R$region -W0.5 -P -O -K <
98.0 -2.5
101.5 1.0
EOF
psxy -J$size -R$region -W0.5 -P -O -K <
98.318 -2.818
101.813 0.677
101.177 1.313
97.682 -2.182
98.318 -2.818
EOF
pstext -J$size -R$region -Gdarkblue -P -O -K <
98.0 -2.5 14 0 1 LT A
EOF
#DRAW CROSSEC DIAGRAM
psbasemap -JX9.5/-5 -R0/550/0/300 -B50/30NWes -P -O -K -Y-7 >> $psfile
pstext -JX -R -Gdarkblue -P -O -K <
0 0 14 0 1 LT A
EOF
#Plot stars of structures, earthquakes, and mainshocks on crossec diagram
awk '{print $4, $3}' xyzpqrs_trench | \
psxy -JX -R -Gdarkblue -Sc0.2 -P -O -K >> $psfile
awk '{print $4, $3}' xyzpqrs_transform | \
psxy -JX -R -Gdarkblue -Sc0.2 -P -O -K >> $psfile
awk '{print $4, $3}' xyzpqrs_ridge | \
psxy -JX -R -Gdarkblue -Sc0.2 -P -O -K >> $psfile
awk '{print $4, $3}' xyzpqrs_kecil | \
psxy -JX -R -Sc0.1c -W0.1/red -P -O -K >> $psfile
awk '{print $4, $3}' xyzpqrs_sedang | \
psxy -JX -R -Sc0.15c -W0.15/red -P -O -K >> $psfile
awk '{print $4, $3}' xyzpqrs_besar | \
psxy -JX -R -Sc0.2c -W0.2/red -P -O -K >> $psfile
awk '{print $4, $3}' xyzpqrs_X | \
psxy -JX -R -Sa0.4c -Gyellow -Wthin -P -O -K >> $psfile
psscale -C$cptfile -Iintens.grd -D0.4c/1/10/0.15 -B1000 -P -O -K -X11.5c -Y12c >> $psfile
# Create legend of magnitude size
cat << EOF >! magnitude.legend
N 3
S 0.025i c 0.025i red 0.01p 0.050i mag 3-5
S 0.050i c 0.050i red 0.01p 0.075i mag 5-6
S 0.075i c 0.075i red 0.01p 0.100i mag 6-9
>
EOF
pslegend magnitude.legend -Dx0.1i/-0.2i/6.8i/1.5i/TC -Jx0.4i -R0/8/0/8 -O -X-4c -Y-5 >> $psfile
\rm -f *.d *.xy xyzpqrs_*
gs -sDEVICE=x11 $psfile
Senin, 22 November 2010
Ibu - Sulis
Ibu..ibu..ibu..ibu
ibu..ibu..ibu..ibu
ibu..ibu..ibu..ibu
Ibu senandung laguku
kunyanyikan slalu untukmu
engkau penyejuk hatiku
tanpamu tiadalah aku
Ku selalu kumohon doamukunyanyikan slalu untukmu
engkau penyejuk hatiku
tanpamu tiadalah aku
bahagia hidup matiku
ibu kaulah harapanku
tanpamu tiadalah aku
Sorga ditelapak kakimu
ridho Allah dengan ridhomu
begitulah sabda nabiku
tanpamu tiadalah aku
ridho Allah dengan ridhomu
begitulah sabda nabiku
tanpamu tiadalah aku
Ibu perisai hidupku
pertaruhkan nyawa bagiku
cintamu terangi jalanku
tanpamu tiadalah aku
pertaruhkan nyawa bagiku
cintamu terangi jalanku
tanpamu tiadalah aku
ibu…ibu…ibu…ibu…
ibu…ibu…ibu…ibu
ibu…ibu…ibu…ibu
gudanglagu.com Free Download Lagu Sulis Ibu MP3 Lirik 4shared Gratis Chord Video Album
Jumat, 19 November 2010
Mekanisme Sumber Gempa (II)
Polaritas Gelombang Seismik
Seismogram yang terekam bisa digunakan dalam studi geometri sesar atau mekanisme fokus karena tergantung radiasi gelombang gempabumi. Yang paling mudah ditentukan dengan menganalisa gelombang body dengan mengambil impuls awal atau polaritas awal gelombang seismik. Polaritas ini tergantung tipe gelombang dan posisi stasiun terhadap jarak gerakan awal pada sumber gempabumi (hypocentre).
Pada gambar sebelah kiri (a), sumber gaya tunggal S yang bisa dianalogikan dengan sumber palu yang dipukulkan secara horisontal. Apabila arah S menuju Δ1 maka di stasiun ini akan terobservasi gelombang P dengan polaritas kompresi (+) artinya menekan stasiun tersebut, di stasiun Δ4 terakam dilatasi yang artinya menarik stasiun tersebut, lain halnya dengan stasiun Δ2 tidak terekam gelombang P sama sekali.
Sebaliknya gelombang S akan terpolarisasi paralel dengan arah dispalcement statis S dan tegaklurus terhadap arah propagasinya dan terekam di Δ2, tapi tidak di Δ1 dan Δ4. Pada stasiun Δ3 gelombang P dan S akan terekam. Kalau letak sumber S di permukaan, maka sinyal gelombang P akan terekam pada komponen radial dan tidak ada pada komponen vertikal.
Gambar sebelah kanan (b) menceritakan pola polarisasi pada sumber shear dislocation sepanjang bidang sesar F yang berbeda dengan kasus gaya tunggal (Gambar 8.7.a). Di stasiun Δ1 dan Δ5 yang terletak pada arah strike, sinyal gelombang P tak terobservasi. Demikian juga di stasiun Δ3 yang tegak lurus sesar di titik F. Sebaliknya di Δ2 dan Δ4 yang bersudut 45° terhadap F, terekam amplitudo gerak gelombang P namun dengan arah berlawanan.
First motion gelombang P dari titik double couple ini, akan menggiring kita pada untuk mendefinisikan empat kuadran, yaitu dua kuadran kompresi dan dua kuadran dilatasi yang dipisahkan oleh dua bidang yang saling tegak lurus yang disebut bidang-bidang nodal, yaitu bidang sesar dan bidang bantu.
Radiasi gelombang P dari sumber sumber gempabumi mempunyai empat kurva daerah konsentrasi. Kita dapat memperkirakan bidang nodal 1 dan 2 (nodal bidang sesar dan bidang bantu) dengan menggunakan polaritas gelombang P dan atau amplitudonya.
Mekanisme Sumber Gempa (I)
Mekanisme terjadinya suatu gempabumi di dalam perut bumi sering dikaitkan dengan kombinasi gaya atau stress yang bekerja pada suatu batuan. Kombinasi stress, kompresi (tekanan kedalam) dan dilatasi (tarikan keluar), yang menyebabkan terjadinya suatu gempabumi dapat dimodelkan dengan mempelajari polarisasi gelombang gempabumi yang terekam pada komponen vertikal.
Model idealisasi dari mekanisme terjadinya suatu gempabumi dalam seismologi disebut dengan mekanisme fokus (focal mechanism). Melalui data seismogram bisa didapatkan banyak informasi gempabumi sehingga diketahui parameter gempabumi seperti : magnitude, kedalaman, lokasi, waktu asal gempabumi, termasuk juga mekanisme fokus. Dengan menganalisis mekanime fokus, kita bisa menganalisis sistem gaya-gaya tektonik yang bekerja pada suatu daerah (Puspito, 1997).
Sekarang kita tinjau bagaimana proses terjadinya sebuah gempabumi. Seorang ahli seismologi Amerika yang bernama Reid pada tahun 1906 mengadakan penelitian untuk membahas tentang proses pemecahan di sebuah sumber gempabumi pada gempabumi yang terjadi di San Andreas Fault. Pergeseran (displacement) dari Fault San Andres ini kebanyakan horizontal, dimana pada bagian timur (kanan) yang menghadap ke daratan Amerika bergerak ke selatan terhadap yang di sebelah barat (kiri) yang menghadap ke Pasifik.
Pada mulanya di dalam bumi terjadi gerakan yang terus-menerus, maka akan terdapat stress yang lama kelamaan akan terakumulasi dan mampu merubah bentuk geologi dari lapisan batuan. Berikutnya suatu lapisan batuan dikenai gaya stress sehingga terjadi perubahan bentuk geologi. Untuk daerah sebelah kiri mendapat stress ke atas, sedangkan daerah kanan mendapat stress ke bawah. Proses ini berjalan terus sampai stress yang terjadi (dikandung) di daerah ini cukup besar untuk merubahnya menjadi gesekan antara kedua daerah. Lama kelamaan karena lapisan batuan sudah tidak mampu lagi untuk menahan stress, maka akan terjadi suatu pergerakan atau perpindahan yang tiba-tiba sehingga terjadilah patahan. Peristiwa pergerakan secara tiba-tiba ini disebut gempabumi. Teori ini dikenal dengan nama Teori Reid atau juga Elastic Rebound Theory.
Mekanisme fokus memberikan tambahan informasi mengenai parameter gempa bumi seperti jenis sesar gempabumi. Parameter sesar terdiri dari ukuran sesar yang dinyatakan dalam km (kilometer) yaitu panjang dan lebar. Selain itu terdapat jarak pergeseran, momen seismik, stress drop, serta source process atau prose pecahnya batuan saat terjadi gempa atau rupture process.
Strike adalah jurus sesar diambil dari permukaan sesar yang arahnya ditarik searah jarum jam (utara mata angin). Kisaran derajatnya dari arah 0° – 360°. Dip adalah sudut kemiringan sesar dari blok yang tegak (foot-wall block) diukur dari bidang mendatar horisontal. Ukuran sudut nilainya dari 0° – 90°. Sedangkan rake atau slip adalah arah pergerakan sesar tersebut diukur dari penampang muka sesar dengan arah diukur dari arah strike ke arah mana slip bergerak (berlawanan arah strike dan dip). Ukuran sudutnya dari arah -180° – 180°. Jarak pergeseran slip atau dislocation dinyatakan dalam besaran m (meter).
Dalam keadaan yang sebenarnya permukaan sesar (patahan) atau fault dapat mempunyai keadaan yang berbeda dan demikian pula dengan gerakannya dapat mempunyai arah yang berlainan sepanjang permukaannya. Dapat dibedakan atas tiga bentuk gerakan dasar dari sesar, yaitu : sesar mendatar, turun, dan naik.
Gerakan sejajar jurus sesar, disebut sesar mendatar atau strike slip fault. Stress yang terbesar adalah stress horisontal dan stress vertikal kecil sekali. Sesar relatif ke bawah terhadap blok dasar, disebut sesar turun / sesar normal atau gravity fault. Gerakan relatif ke atas terhadap blok dasar, disebut sesar naik atau thrust fault / reverse fault.
Gerakan sejajar jurus sesar, disebut sesar mendatar atau strike slip fault. Stress yang terbesar adalah stress horisontal dan stress vertikal kecil sekali. Sesar relatif ke bawah terhadap blok dasar, disebut sesar turun / sesar normal atau gravity fault. Gerakan relatif ke atas terhadap blok dasar, disebut sesar naik atau thrust fault / reverse fault.
Kamis, 18 November 2010
Gelombang gempa
Gelombang seismik adalah gelombang elastik gempabumi yang menjalar ke seluruh bagian dalam bumi dan melalui permukaan bumi, akibat adanya lapisan batuan yang patah secara tiba – tiba atau adanya suatu ledakan. Dapat juga dianalogikan sebagai gelombang yang menjalar seperti pada suatu kolam air yang dijatuhkan di atasnya sebutir batu. Air mengalami gangguan dan gelombangnya terpancar keluar dari pusat awalnya mencapai jarak terjauh kolam. Akan tetapi partikel air yang terganggu tersebut tak bergeser dalam arah pergerakan gelombang.
Gelombang utama gempabumi bumi terdiri dari dua tipe yaitu gelombang bodi (body wave) dan gelombang permukaan (surface wave). Gelombang seismik merambat dalam lapisan bumi sesuai dengan prinsip yang berlaku pada perambatan gelombang cahaya: pembiasan dengan koefisien bias, pemantulan dengan koefisien pantul, hukum-hukum Fermat, Huygens, Snellius dan lain-lain.
Gelombang bodi menjalar melalui bagian dalam bumi dan biasa disebut free wave karena dapat menjalar ke segala arah di dalam bumi. Gelombang bodi terdiri atas gelombang primer dan gelombang sekunder. Gelombang primer P merupakan gelombang longitudinal atau gelombang kompresional, gerakan partikelnya sejajar dengan arah perambatannya. Sedang gelombang sekunder S merupakan gelombang transversal atau shear, gerakan partikelnya terletak pada suatu bidang yang tegak lurus dengan arah penjalarannya. Kecepatan gelombang P lebih tinggi dari gelombang S.
Gelombang S terdiri dari dua komponen, yaitu gelombang SH dengan gerakan partikel horizontal dan gelombang SV dengan gerakan partikel vertikal. Gelombang P mampu menembus lapisan inti bumi sedangkan gelombang S tidak bisa dikarenakan sifatnya yang tak bisa menembus media cair pada inti bumi.
Sifat penjalaran gelombang P yang langsung adalah bahwa gelombang ini akan menjadi hilang pada jarak lebih besar dari 130º, dan tidak terlihat sampai dengan jarak kurang dari 140º. Hal tersebut disebabkan karena adanya inti bumi. Gelombang langsung P akan menyinggung permukaan inti bumi pada jarak 103º dan pada jarak yang akan mengenai inti bumi pada jarak 144º. Gelombang P akan timbul kembali yaitu gelombang yang menembus inti bumi dengan dua kali mengalami refraksi. Menghilangnya gelombang P pada jarak 103º memungkinkan untuk menghitung kedalaman lapisan inti bumi.
Guttenberg (1913) mendapatkan kedalaman inti bumi 2900 km. Telah didapatkan pula bahwa batas mantel dengan inti bumi merupakan suatu diskontinuitas yang tajam. Daerah antara 103º - 144º disebut sebagai Shadow zone, walaupun sebenarnya fase yang lemah dapat pula terlihat di daerah ini.
Walaupun gelombang bodi dapat menjalar ke segala arah di permukaan bumi, namun tetap tidak dapat menembus inti bumi sebagai gelombang transversal. Keadaan ini membuktikan bahwa inti luar bumi berupa fluida. Untuk penelitian tetap diasumsikan keadaan homogen, yaitu bagian luar bumi dan inti bumi (dua media homogen yang berbeda).
Kadang – kadang juga ditemui suatu fase yang kuat di daerah Shadow zone sampai ke jarak kurang lebih 110º. Karena adanya fase inilah pada tahun 1930 ditemukan media lain yaitu inti dalam. Batas dari inti dalam ini terdapat pada kedalaman 5100 km . Diperkirakan kecepatan gelombang seismik di inti dalam lebih tinggi dari pada di inti luar. Untuk membedakan dan identifikasi, maka perlu pemberian nama untuk gelombang seismik yang melalui inti bumi luar dan dalam.
Gelombang permukaan (surface) merupakan gelombang elastik yang menjalar sepanjang permukaan bumi dan biasa disebut sebagai tide waves. Karena gelombang ini terikat harus menjalar melalui suatu lapisan atau permukaan. Gelombang permukaan terdiri dari gelombang Love (L) dan Rayleigh (R) yang menjalar melalui permukaan bebas dari bumi.
Gelombang L gerakan partikelnya sama dengan gelombang SH dan memerlukan media yang berlapis. Gelombang R lintasan gerak partikelnya merupakan suatu ellips. Bidang ellips ini vertikal dan berimpit dengan arah penjalarannya. Gerakan partikelnya ke belakang (bawah maju atas mundur). Gelombang R menjalar melalui permukaan media yang homogen.
Gelombang Stonely, arah penjalarannya seperti gelombang R tetapi menjalar melalui batas antara dua lapisan di dalam bumi. Gelombang Channel, yaitu gelombang yang menjalar melalui lapisan yang berkecepatan rendah (low velocity layer) di dalam bumi.
Gelombang Love dan Rayleigh ada juga yang memberi simbul LQ dan LR dimana L singkatan dari Long karena gelombang permukaan mempunyai sifat periode panjang dan Q adalah singkatan dari Querwellen yaitu nama lain dari Love seorang Jerman yang menemukan gelombang ini.
Gelombang LQ dan LR menjalar sepanjang permukaan bebas dari bumi atau lapisan batas diskontinuitas antara kerak dan mantel bumi. Amplitude gelombang LQ dan LR adalah yang terbesar pada permukaan dan mengecil secara eksponensial terhadap kedalaman. Dengan demikian pada gempa-gempabumi dangkal amplitude gelombang LQ dan LR akan mendominasi.
Dari hasil pengamatan gelombang permukaan ini diperoleh dua ketentuan utama baru yang menunjukkan bahwa bagian bumi berlapis-lapis dan tidak homogen. Ditemukan juga adanya perubahan dispersi kecepatan (velocity dispersion). Fakta menyebutkan bahwa gelombang L tidak dapat menjalar pada permukaan suatu media yang kecepatannya naik terhadap kedalaman. Oleh karena itu gelombang L dan R tidak datang bersama-sama pada suatu stasiun, tetapi gelombang yang mempunyai periode lebih panjang akan datang lebih dahulu. Dengan kata lain gelombang yang panjang periodenya mempunyai kecepatan yang tinggi.
Sabtu, 13 November 2010
Sejarah Seismograph Indonesia Era 2006
Setelah Tsunami Aceh 2004, Pemerintah Republik Indonesia bekerjasama dengan 14 negara donor serta institusi dalam dan luar negeri (diantaranya ; UNESCO, CTBTO, America, French, Japan, German, China) bersama-sama untuk membangun sistem baru peringatan dini tsunami atau Tsunami Early Warning System (TEWS). Tujuannya untuk mengurangi korban jiwa lebih besar diakibatkan oleh bahaya tsunami. Mulai tahun 2005 untuk mewujudkan program tersebut akan diinstalasi sekitar 160 seismograf, 500 akselerograf, dan 15 digital strong-motion akselerograf.
Pada tahun 2010 BMG berganti nama menjadi Badan Meteorologi Klimatologi dan Geofisika (BMKG). Operasional monitoring seismik seluruh wilayah Indonesia dikembangkan menjadi 10 PGR. Sedangkan untuk observasi muka laut BPPT dan RISTEK sebagai mitra kerja dalam negeri BMKG bertugas menangani operasional 60 tide gauge, dan 15 Dart-Buoy, secara berurutan disebar ke seluruh wilayah Indonesia. Telekomunikasi yang digunakan adalah 5 in 1 terdiri dari ; internet (web, mail), sms and mobile-phone, radio-internet, faximile, telepon. Oleh karenanya penentuan parameter gempabumi bumi disertai diseminasi peringatan gempabumi dan tsunami sekarang ini bisa dicapai dalam tempo 5 – 10 menit ke tangan pengguna.
Mulai tahun 2006, BMKG mengadopsi software analisa SeiscomP dari GFZ Jerman untuk menentukan parameter gempabumi. Hal ini adalah bentuk implementasi kerjasama bilateral Indonesia – Jerman. Institusi yang terbentuk adalah GITEWS (German Indonesia - Tsunami Early Warning System). Tahun 2007, China tidak mau ketinggalan untuk berkecimpung dalam Ina-TEWS, software analisa epicentre MSDP CEA di-instalasi untuk membandingkan hasil analisa SeiscomP pada saat penentuan lokasi pusat gempa.
Pembangunan Ina-TEWS secara masif diteruskan, sejak 2006 sampai 2008, BMKG terus mengusulkan penambahan sensor seismograf untuk melengkapi sebaran pusat gempabumi di daerah-daerah rawan tektonik. Pembangunan itu meliputi satu pusat nasional, 10 pusat regional, 160 seismometer broadband, dan 500 akselerometer
CTBT mempunyai beberapa jenis jaringan seismik di dunia, yaitu : sistem primary dan auxiliary. Sistem primary terdiri dari 50 stasiun, 30 stasiun array, 19 stasiun 3-komponen. Sistem auxiliary terdiri dari 120 stasiun, 7 stasiun array, 112 stasiun 3-komponen. Jaringan stasiun seismograf auxiliary CTBT, 6 stasiun berada di Indonesia yaitu : Kappang, Parapat, Lembang, Kupang, Sorong dan Jayapura.
GITEWS secara bertahap membangun sistem peringatan tsunami berbasis database pemodelan tsunami yang diverifikasi dengan observasi permukaan air laut. Sistem yang rencananya diluncurkan tahun 2010 dinamakan DSS (Decision Support System) bertujuan untuk membantu operator gempabumi untuk menentukan keputusan peringatan tsunami. DSS memilah-milah segmen pantai tingkat kecamatan menurut tingkatan peringatan (mayor, tsunami, saran)
berdasarkan nilai perkiraan ketinggian tsunami, kecepatan waktu tiba, dan proporsi populasi geografis di tiap segmen pantai rawan tsunami.
Berdasarkan informasi detil peringatan tsunami yang disampaikan tersebut, pemerintah daerah di daerah bencana tersebut akan mampu memutuskan tindakan mitigasi yang diperlukan, misalkan ; evakuasi total, sebagian atau hanya waspada. Informasi diteruskan oleh pemda ke masyarakat melalui sirine atau alat telekomunikasi setempat.
Sekalipun jumlah jaringan seismik BMKG telah mengalami peningkatan cukup siginifikan dari tahun-tahun sebelumnya, namun hal itu masih sangat dirasakan kurang dibandingkan dengan luas daerah Indonesia dengan aktifitas gempabumi yang tinggi, karenanya dilakukan kerjasama dengan jaringan seismograf luar negeri agar bisa menambah dan saling bertukar data gempa. Saat ini BMKG baru dapat menerima data seismik yang real time dari Australia, Malaysia dan beberapa jaringan seismik internasional seperti Geofon dan IRIS.
IRIS (Incorporated Reserach Institutions for Seismology) http://www.iris.edu/hq/ adalah suatu konsorsium nasional negara-negara Eropa dalam pengoperasian fasilitas ilmiah, manajemen dan distribusi data seismik global. IRIS telah berperan besar dalam rangka memajukan infrastruktur dan penelitian ilmiah tentang bencana gempabumi, eksplorasi sumber daya alam, dan monitoring percobaan ledakan nuklir, melalui jaringan nasional dan internasional seismik GSN (Global Seismografic Network), IRIS PASSCAL, dan IRIS DMS. IRIS telah melakukan hubungan kemitraan dan kolaborasi dengan hampir seluruh negara di dunia dengan membantu pengembangan infrastruktur teknis, dan kapasitas SDM.
Aktifitas USGS (United States of Geological Surveys) http://earthquake.usgs.gov/earthquakes/recenteqsww/ disamping memantau aktifitas getaran gempabumi yang terjadi di negara-negara bagian Amerika, namun juga di dunia lainnya. Perjanjian kerjasama dan kontrak kerja USGS dalam hal studi gempabumi bumi diterapkan melalui program hibah dengan pihak perguruan tinggi, negara, regional dan lokal instansi pemerintah, swasta dan industri yang bertujuan untuk mengembangkan informasi, pengetahuan, dan metode yang relevan dalam program bencana gempabumi. Data dan produk USGS beberapa diantaranya dapat diakses melalui internet seperti katalog gempa, waveform data, data bahaya gempa, getaran tanah, dan informasi kerak bumi.
BMKG bekerjasama juga dengan organisasi PBB yang membidangi pengawasan percobaan senjata nuklir yaitu CTBTO (Commision Nuclear Test-Ban Treaty Organization) dalam hal pertukaran data gempabumi dengan pengawasan IDC (International Data Centre) untuk keperluan sistem peringatan dini tsunami Indonesia. Saat ini CTBT sedang mengembangkan teknologi untuk memonitor ledakan nuklir dengan menggunakan metode seismik, infrasound, hydroakustic dan radiasi nuklir. Data Hydroakustik sangat berguna untuk tujuan penelitian seperti perambatan retakan. Data seismik tambahan dapat diminta untuk akses data real time.
CTBT mempunyai beberapa jenis jaringan seismik di dunia, yaitu : sistem primary dan auxiliary. Sistem primary terdiri dari 50 stasiun, 30 stasiun array, 19 stasiun 3-komponen. Sistem auxiliary terdiri dari 120 stasiun, 7 stasiun array, 112 stasiun 3-komponen. Jaringan stasiun seismograf auxiliary CTBT, 6 stasiun berada di Indonesia yaitu : Kappang, Parapat, Lembang, Kupang, Sorong dan Jayapura.
Label:
Broadband,
GITEWS,
Seismograp BMKG,
seismograph
Sejarah Seismograph Indonesia 1990 - 2003
Tidak mudah untuk menempatkan instrumen atau sensor gempabumi di setiap wilayah Indonesia. Ada beberapa persyaratan yang harus dilalui dalam hal pemilihan lokasi sensor diantaranya : berada pada batuan keras, jauh dari kebisingan jalan berkendaraan, keamanan terjamin, tersedia fasilitas listrik, solar panel, dan komunikasi.
Era Tahun 1900 – 1930. BMG (Badan Meteorologi Geofisika) sebagai lembaga pemerintah Indonesia bertugas untuk memonitor aktifitas gempabumi bumi di Indonesia sejak zaman kolonial Belanda pada tahun 1898 dengan mengoperasikan seismograf mekanik Ewing. Pada tahun 1908 telah mulai mengoperasikan stasiun pemantau gempabumi permanen, yakni dengan memasang seismograf Wichert komponen horisontal di Jakarta. Sedangkan komponen vertikal seismograph tersebut dipasang pada tahun 1928 di beberapa kota yaitu Jakarta, Medan, Bengkulu dan Ambon.
Era Tahun 1950 – 1980. Pada tahun 1953 dengan nama PMG (Pusat Meteorologi dan Geofisika) sebagai instansi yang terkait dengan pengamatan gempabumi bumi memasang seismograf Elektromagnetik Sprengnether di Lembang - Bandung yang disusul dengan pemasangan seismograf bertipe sama di Jakarta, Medan, Tangerang, Denpasar, Ujungpandang, Kupang, Jayapura, Manado dan Ambon sehingga terbentuk jaringan seismograf yang pertama kali di Indonesia. Seismograf 3 komponen ini beroperasi di sepuluh kota tersebut sampai dengan tahun 1980-an.
Pada tahun 1964 di stasiun Lembang dipasang Seismograf Teledyne Geotech yang termasuk dalam jaringan WWSSN (World Wide Standard Seismololgical Network). Seismograf ini memiliki 6 komponen dan mengalami modifikasi pada tahun 1978. UNDP-Unesco pada tahun 1974 mengadakan proyek pengembangan seismologi di Indonesia yang antara lain meliputi standarisasi seismograf dan proses pengolahan data gempabumi bumi serta pengembangan jaringan pemantau. Salah satu bentuknya adalah pemasangan seismograf periode pendek (Short Period Seismograf - Kinemetric) komponen Z di 27 stasiun seluruh Indonesia. Tiap-tiap stasiun dilengkapi dengan seismograf 1 komponen vertikal periode pendek, dan sinyal seismik direkam pada kertas seismogram.
Terdapat 30 stasiun geofisika konvensional dan 28 stasiun telemetri serta 10 stasiun diantaranya telah ditingkatkan menjadi 3 komponen periode pendek. BMG mempunyai sebuah PGN (Pusat Gempabumi Nasional) dan 5 PGR (Pusat Gempabumi Regional) atau wilayah besar pemantauan gempabumi sekaligus cuaca, yakni BMG Wilayah I di Medan, BMG Wilayah II di Ciputat, BMG Wilayah III di Denpasar, BMG Wilayah IV di Ujung Pandang dan BMG Wilayah V di Jayapura. Seluruh stasiun ini pada tahun 1998 dilengkapi dengan fasilitas GARNET.
Jaringan tersebut masih beroperasi hingga saat ini dan merupakan jaringan pemantau seismik utama BMG. Sejak tahun ini pula dapat dikatakan bahwa BMG memiliki dua tipe stasiun pemantau gempabumi bumi di Indonesia. Pertama adalah stasiun telemetri yang tidak berawak atau telemetry dan lainnya adalah stasiun geofisika konvensional. Di stasiun geofisika konvensional, data gempabumi bumi diobservasi dengan bantuan operator kemudian dilanjutkan dengan pengolahan data dan analisis parameter gempabumi bumi sementara.
Untuk pengolahan data gempabumi di Balai Wilayah, data gempabumi dari stasiun seismograf dikirim ke Balai Wilayah dan PGN setiap 3 jam melalui SSB, telex, internet, atau sarana telekomunikasi lain, bersama-sama dengan data meteorologi. Sekarang ini fasilitas komunikasi sudah dilengkapi dengan sarana VSAT untuk komunikasi stasiun dengan Balai wilayah dan Pusat.
BMG juga menginstalasi stasiun telemetri 1 komponen berasal dari Laboratoire De Geophysique (LDG) of France. Saat itu proses signal seismik di stasiun dan PGR masih dirasakan terbatas karena masih direkam secara analog pada recorder grafis. Hanya pada saat kejadian gempa, signal yang mengandung gempabumi diseleksi dan direkam secara digital dengan 50 sampel per detik dan 12 bit dikirimkan ke PGN di BMG Jakarta menggunakan PT. TELKOMSEL menggunakan medium kecepatan gelombang 4800 bps.
Pada tahun 1993 di Stasiun Geofisika Tretes Jawa Timur dipasang seismograf periode panjang (Long Period Seismograf) 3 komponen yang dilengkapi dengan TREMORS. Di tahun ini pula dipasang seismograf periode pendek 3 komponen SPS-3 (Kinemetrics) di 9 stasiun geofisika konvensional di seluruh Indonesia yaitu di Banda Aceh, Padang Panjang, Kepahyang, Kotabumi, Tanjungpandan, Kupang, Palu, Ambon dan Sorong.
Perkembangan lain dari sistem pemantau seismik BMG adalah dimulainya era broadband sejak tahun 1992 pada saat dioperasikannya seismograf 3 komponen tipe Broadband di stasiun Parapat dan Jayapura. Keduanya hingga saat ini masih beroperasi. Menyusul pada kurun waktu 1997-2001 dengan adanya proyek kerjasama Indonesia dan Jepang yaitu Joint Operation of Japan - Indonesia Seismik Network (JISNET) dipasang seismograf jenis broadband di 23 stasiun di seluruh Indonesia. Tahun 1996 sistem monitoring real-time telah ditingkatkan dengan menambah data akuisisi dengan fasilitas database komputer, display signal, dan perangkat pemetaan untuk prosesing digital di PGN dan waktu GPS di semua regional. Pada saat itu BMG telah mampu melalukan deteksi pusat gempabumi dalam waktu 15 menit sampai 1 jam dengan besaran gempabumi terendah mulai skala magnitude 4.
Sementara itu, pada tahun 1999 di Kappang (Sulawesi Selatan) dipasang seismograf 3 komponen jenis broadband yang merupakan kerjasama BMG-UCSD/USA. Pada tahun 2002 di stasiun yang sama kembali dipasang seismograf bertipe broadband yang merupakan salah satu dari 6 stasiun seismik CTBTO (Comprehensive Nuclear Test Ban Treaty Organization). Lima stasiun lainnya adalah Parapat, Lembang, Kupang, Sorong dan Jayapura. Proyek kerjasama ini dilanjutkan kembali antara NIED Jepang dan BMG untuk periode 2001-2006 dengan nama Operation and Data Exchange of Japan - Indonesia Seismik Network (JISNET continued). Pelaksanaan proyek ini meliputi pemasangan seismograf jenis Broadband di 22 stasiun seluruh Indonesia.
Sistem yang digunakan adalah PAC (Phasing and Analog Converter) yang telah sinkron dengan waktu GPS receiver. Software analisa dinamakan ARTDAS (Automatic Real Time Data Acquisition Software) dan XIDAS (X-windows Interactive Data Analysis Software) yang dioperasikan oleh Sunwork stations LDG France dengan melakukan akuisisi real time, rekaman kontinyu, dan proses semi otomatis menghasilkan parameter seismik seperti waktu tiba gelombang, amplitudo, periode, waktu asal gempa, lintang, bujur, dll yang semuanya disimpan di ORACLE database.
Era Tahun 2000. Pada tahun 2003 dibentuk Sistem Pemantauan Seismik Nasional (National Seismik Monitoring System) dengan penambahan seismograf broadband di 27 stasiun-stasiun seismik seluruh Indonesia. Seismograf ini terintegrasi dengan jaringan yang telah ada dan mempunyai sistem pengolahan data real time berlokasi di Jakarta dengan 3 Pusat Seismik Regional Mini (Mini Regional Seismik Center) yang berlokasi di Padangpanjang, Kepahyang, Palu.
Era Tahun 1900 – 1930. BMG (Badan Meteorologi Geofisika) sebagai lembaga pemerintah Indonesia bertugas untuk memonitor aktifitas gempabumi bumi di Indonesia sejak zaman kolonial Belanda pada tahun 1898 dengan mengoperasikan seismograf mekanik Ewing. Pada tahun 1908 telah mulai mengoperasikan stasiun pemantau gempabumi permanen, yakni dengan memasang seismograf Wichert komponen horisontal di Jakarta. Sedangkan komponen vertikal seismograph tersebut dipasang pada tahun 1928 di beberapa kota yaitu Jakarta, Medan, Bengkulu dan Ambon.
Era Tahun 1950 – 1980. Pada tahun 1953 dengan nama PMG (Pusat Meteorologi dan Geofisika) sebagai instansi yang terkait dengan pengamatan gempabumi bumi memasang seismograf Elektromagnetik Sprengnether di Lembang - Bandung yang disusul dengan pemasangan seismograf bertipe sama di Jakarta, Medan, Tangerang, Denpasar, Ujungpandang, Kupang, Jayapura, Manado dan Ambon sehingga terbentuk jaringan seismograf yang pertama kali di Indonesia. Seismograf 3 komponen ini beroperasi di sepuluh kota tersebut sampai dengan tahun 1980-an.
Pada tahun 1964 di stasiun Lembang dipasang Seismograf Teledyne Geotech yang termasuk dalam jaringan WWSSN (World Wide Standard Seismololgical Network). Seismograf ini memiliki 6 komponen dan mengalami modifikasi pada tahun 1978. UNDP-Unesco pada tahun 1974 mengadakan proyek pengembangan seismologi di Indonesia yang antara lain meliputi standarisasi seismograf dan proses pengolahan data gempabumi bumi serta pengembangan jaringan pemantau. Salah satu bentuknya adalah pemasangan seismograf periode pendek (Short Period Seismograf - Kinemetric) komponen Z di 27 stasiun seluruh Indonesia. Tiap-tiap stasiun dilengkapi dengan seismograf 1 komponen vertikal periode pendek, dan sinyal seismik direkam pada kertas seismogram.
Terdapat 30 stasiun geofisika konvensional dan 28 stasiun telemetri serta 10 stasiun diantaranya telah ditingkatkan menjadi 3 komponen periode pendek. BMG mempunyai sebuah PGN (Pusat Gempabumi Nasional) dan 5 PGR (Pusat Gempabumi Regional) atau wilayah besar pemantauan gempabumi sekaligus cuaca, yakni BMG Wilayah I di Medan, BMG Wilayah II di Ciputat, BMG Wilayah III di Denpasar, BMG Wilayah IV di Ujung Pandang dan BMG Wilayah V di Jayapura. Seluruh stasiun ini pada tahun 1998 dilengkapi dengan fasilitas GARNET.
Jaringan tersebut masih beroperasi hingga saat ini dan merupakan jaringan pemantau seismik utama BMG. Sejak tahun ini pula dapat dikatakan bahwa BMG memiliki dua tipe stasiun pemantau gempabumi bumi di Indonesia. Pertama adalah stasiun telemetri yang tidak berawak atau telemetry dan lainnya adalah stasiun geofisika konvensional. Di stasiun geofisika konvensional, data gempabumi bumi diobservasi dengan bantuan operator kemudian dilanjutkan dengan pengolahan data dan analisis parameter gempabumi bumi sementara.
Untuk pengolahan data gempabumi di Balai Wilayah, data gempabumi dari stasiun seismograf dikirim ke Balai Wilayah dan PGN setiap 3 jam melalui SSB, telex, internet, atau sarana telekomunikasi lain, bersama-sama dengan data meteorologi. Sekarang ini fasilitas komunikasi sudah dilengkapi dengan sarana VSAT untuk komunikasi stasiun dengan Balai wilayah dan Pusat.
BMG juga menginstalasi stasiun telemetri 1 komponen berasal dari Laboratoire De Geophysique (LDG) of France. Saat itu proses signal seismik di stasiun dan PGR masih dirasakan terbatas karena masih direkam secara analog pada recorder grafis. Hanya pada saat kejadian gempa, signal yang mengandung gempabumi diseleksi dan direkam secara digital dengan 50 sampel per detik dan 12 bit dikirimkan ke PGN di BMG Jakarta menggunakan PT. TELKOMSEL menggunakan medium kecepatan gelombang 4800 bps.
Pada tahun 1993 di Stasiun Geofisika Tretes Jawa Timur dipasang seismograf periode panjang (Long Period Seismograf) 3 komponen yang dilengkapi dengan TREMORS. Di tahun ini pula dipasang seismograf periode pendek 3 komponen SPS-3 (Kinemetrics) di 9 stasiun geofisika konvensional di seluruh Indonesia yaitu di Banda Aceh, Padang Panjang, Kepahyang, Kotabumi, Tanjungpandan, Kupang, Palu, Ambon dan Sorong.
Perkembangan lain dari sistem pemantau seismik BMG adalah dimulainya era broadband sejak tahun 1992 pada saat dioperasikannya seismograf 3 komponen tipe Broadband di stasiun Parapat dan Jayapura. Keduanya hingga saat ini masih beroperasi. Menyusul pada kurun waktu 1997-2001 dengan adanya proyek kerjasama Indonesia dan Jepang yaitu Joint Operation of Japan - Indonesia Seismik Network (JISNET) dipasang seismograf jenis broadband di 23 stasiun di seluruh Indonesia. Tahun 1996 sistem monitoring real-time telah ditingkatkan dengan menambah data akuisisi dengan fasilitas database komputer, display signal, dan perangkat pemetaan untuk prosesing digital di PGN dan waktu GPS di semua regional. Pada saat itu BMG telah mampu melalukan deteksi pusat gempabumi dalam waktu 15 menit sampai 1 jam dengan besaran gempabumi terendah mulai skala magnitude 4.
Sementara itu, pada tahun 1999 di Kappang (Sulawesi Selatan) dipasang seismograf 3 komponen jenis broadband yang merupakan kerjasama BMG-UCSD/USA. Pada tahun 2002 di stasiun yang sama kembali dipasang seismograf bertipe broadband yang merupakan salah satu dari 6 stasiun seismik CTBTO (Comprehensive Nuclear Test Ban Treaty Organization). Lima stasiun lainnya adalah Parapat, Lembang, Kupang, Sorong dan Jayapura. Proyek kerjasama ini dilanjutkan kembali antara NIED Jepang dan BMG untuk periode 2001-2006 dengan nama Operation and Data Exchange of Japan - Indonesia Seismik Network (JISNET continued). Pelaksanaan proyek ini meliputi pemasangan seismograf jenis Broadband di 22 stasiun seluruh Indonesia.
Sistem yang digunakan adalah PAC (Phasing and Analog Converter) yang telah sinkron dengan waktu GPS receiver. Software analisa dinamakan ARTDAS (Automatic Real Time Data Acquisition Software) dan XIDAS (X-windows Interactive Data Analysis Software) yang dioperasikan oleh Sunwork stations LDG France dengan melakukan akuisisi real time, rekaman kontinyu, dan proses semi otomatis menghasilkan parameter seismik seperti waktu tiba gelombang, amplitudo, periode, waktu asal gempa, lintang, bujur, dll yang semuanya disimpan di ORACLE database.
Era Tahun 2000. Pada tahun 2003 dibentuk Sistem Pemantauan Seismik Nasional (National Seismik Monitoring System) dengan penambahan seismograf broadband di 27 stasiun-stasiun seismik seluruh Indonesia. Seismograf ini terintegrasi dengan jaringan yang telah ada dan mempunyai sistem pengolahan data real time berlokasi di Jakarta dengan 3 Pusat Seismik Regional Mini (Mini Regional Seismik Center) yang berlokasi di Padangpanjang, Kepahyang, Palu.
Label:
Seismograp BMKG
Sabtu, 25 September 2010
Sabtu, 21 Agustus 2010
SWIFT FOCAL INVRSI DAN SOURCE TIME FUNCTION
CARA BUKA SWIFT
1) Buka Auto
cd seis/cmt/evt/2009...../09..../fs0.005.....auto/
less search.log ----------> untuk mencari ID
2) Buka Manual Script
cd seis/cmt/script.manual
source envset.sh
./manual.prep.sh ID ------> copy paste
3) Buka Manual Evt
cd seis/cmt/2009.../09..../fs......manual/
cd param
kghostview wplot........ps & -------------> untuk pilih sinyal baik. Yang buruk dikres ##
emacs sthls.prm.manual &
Save
4) Buka Manual Script
cd seis/cmt/script.manual
./manual 1st.search.sh -p ID
5) Buka Manual Evt
gs wffits.ps
1) Buka Auto
cd seis/cmt/evt/2009...../09..../fs0.005.....auto/
less search.log ----------> untuk mencari ID
2) Buka Manual Script
cd seis/cmt/script.manual
source envset.sh
./manual.prep.sh ID ------> copy paste
3) Buka Manual Evt
cd seis/cmt/2009.../09..../fs......manual/
cd param
kghostview wplot........ps & -------------> untuk pilih sinyal baik. Yang buruk dikres ##
emacs sthls.prm.manual &
Save
4) Buka Manual Script
cd seis/cmt/script.manual
./manual 1st.search.sh -p ID
5) Buka Manual Evt
gs wffits.ps
Selasa, 01 Juni 2010
Instalasi GMT
Pada dasarnya GMT dibuat berbasis bahasa pemrograman C-shell, Fortran yang cocok dioperasikan pada Unix/Linux. Tetapi kemudian si pencipta – Paul Wessel dan Walter Smith, membuat jembatan antara GMT C-shell ini dengan aplikasi Windows seperti yang popular, yaitu dengan software tambahan NetDF dan Cygwin, dan Ghost view serta Ghost script. Urutan instalasi NetCDF, Cygwin, Ghostview dan Ghostscript, dan terakhir GMT.
1.NetCDF (Network Common Data Form) adalah software untuk mengatur akses pemrograman data ilmiah.
http://www.unidata.ucar.edu/downloads/netcdf/index.jsp
2.Cygwin untuk tampilan linux pada windows dengan operasi C-shell. Untuk instalasi pakai perintah setup.exe
3.Kemudian mengubah tampilan Field Default pada All category menjadi Install.
4.Proses instalasi Cygwin dapat mencapai waktu 30 – 60 menit (sabar aja ya, jangan diinterupsi !)
http://cygwin.com/
5.Instal GMT
http://gmt.soest.hawaii.edu/
6.Ghostview Ghostscript untuk tampilan gambar.
http://pages.cs.wisc.edu/~ghost/
7.Sesuaikan HOME PATH pada Control panel sesuai directory cygwin (biasanya di C:).
1.NetCDF (Network Common Data Form) adalah software untuk mengatur akses pemrograman data ilmiah.
http://www.unidata.ucar.edu/downloads/netcdf/index.jsp
2.Cygwin untuk tampilan linux pada windows dengan operasi C-shell. Untuk instalasi pakai perintah setup.exe
3.Kemudian mengubah tampilan Field Default pada All category menjadi Install.
4.Proses instalasi Cygwin dapat mencapai waktu 30 – 60 menit (sabar aja ya, jangan diinterupsi !)
http://cygwin.com/
5.Instal GMT
http://gmt.soest.hawaii.edu/
6.Ghostview Ghostscript untuk tampilan gambar.
http://pages.cs.wisc.edu/~ghost/
7.Sesuaikan HOME PATH pada Control panel sesuai directory cygwin (biasanya di C:).
Langganan:
Postingan (Atom)