R=-R114/117/-4/-1
gmt begin coba png
	gmt basemap $R -B1 -JM20
	gmt grdimage topo30.grd -JM20 -Ctopo.cpt -I+a15+ne0.75
	gmt coast -W1 -Df -Na/thickest,black -Ir/0.5p,blue -Lf114.5/-3.8/17/100+l+jt -Tf114.5/-3/0.9i/2
	awk '{print $1,$2,$3,$4*0.03}' event.txt | gmt psxy -JM20 -Cgempa.cpt -Sci -Wthin
	echo 115.2 -3.2 A> line.dat
	echo 116.4 -1.1 B>> line.dat
	gmt psxy line.dat -JM20 $R -W2
	gmt text line.dat -JM20
	awk '{print $1,$2,$3}' event.txt | gmt project -C115.2/-3.2 -E116.4/-1.1 -Fpz -W-100/100 > cross.dat
	gmt basemap -JX15/-10 -X-1 -R0/250/0/70 -B10:'Distance (km)':/:'Depth (km)':50WSne -Y-12.5
	echo 0 0 14 0 1 BC A> teks.dat
	echo 250 0 14 0 1 BC B>> teks.dat
	gmt text teks.dat -JX -Gblue -N
	gmt psxy cross.dat -JX -Sc0.2 -W1 -Gred
gmt end show