二、圖示「綠水觀測站」相關的基本資料

1.以GMT畫出綠水附近的地圖。（在此特別感謝嘉鈴學姊提供立霧溪、和平溪與美崙溪的位置資料點）
------------------------------------------------------------------------------------------------
psbasemap -JM6i -X1.5 -R121/122/24/25 -B0.2NWes:."LU-SHUI": -V -K > LU-SHUI.ps

pscoast -JM6i -R121/122/24/25 -Df -W2 -Na -Ia -G162/205/90 -T123/34/1 -L123/32.5/27/100k -V -K -O >> LU-SHUI.ps

psxy hp.txt -JM6i -R121/122/24/25 -MX -W -V -H -K -O >> LU-SHUI.ps
psxy lw.txt -JM6i -R121/122/24/25 -MX -W -V -H -K -O >>LU-SHUI.ps
psxy mk.txt -JM6i -R121/122/24/25 -MX -W -V -H -K -O >> LU-SHUI.ps

psxy temp1.txt -JM6i -R121/122/24/25 -Sa0.25 -G255/0/0 -V -H -K -O >> LU-SHUI.ps
------------------------------------------------------------------------------------------------

補充說明：紅色星號處為綠水觀測站的位置，以下就不再贅述了。

2.畫出綠水觀測站附近的等高線圖
------------------------------------------------------------------------------------------------
makecpt -Crelief -T300/3000/100 -V -Z > high.cpt

psbasemap -JM6i -X1.5 -R121.4/121.6/24.1/24.3 -B0.2NWes:."LU-SHUI": -V -K > contour_lu-shui.ps

pscoast -JM6i -R121.4/121.6/24.1/24.3 -Df -W2 -Na -Ia -T123/34/1 -L123/32.5/27/100k -V -K -O >> contour_lu-shui.ps

grdcontour twdtm_small.grd -JM6i -R121.4/121.6/24.1/24.3 -C500 -A500 -V -K -O >> contour_lu-shui.ps

psxy lw.txt -JM6i -R121.4/121.6/24.1/24.3 -Sc0.05 -W0/200/255 -MX -W -V -H -K -O >>contour_lu-shui.ps

psxy temp1.txt -JM6i -R121.4/121.6/24.1/24.3 -Sa0.4 -G255/0/0 -V -H -K -O >> contour_lu-shui.ps
------------------------------------------------------------------------------------------------

補充說明：藍色線為立霧溪，紅色星號為綠水觀測站。

3.畫出綠水觀測站附近的立體圖
------------------------------------------------------------------------------------------------
makecpt -Crainbow -T0/400/40 -V -Z > depth.cpt

grd2cpt twdtm_small.grd -Crelief -R121.4/121.7/24.0/24.3 -V -Z > lushui_test.cpt
grdgradient twdtm_small.grd -Gtwdtm_small_temp_shade.grd -A90 -Ne0.9 -V


psbasemap -JM6i -X1.5 -R121.4/121.7/24.0/24.3 -B0.2NWes:."LU-SHUI": -V -K > shade_LU-SHUI.ps

grdimage twdtm_small.grd -Itwdtm_small_temp_shade.grd -JM6i -R121.4/121.7/24.0/24.3 -Clushui_test.cpt -V -K -O >> shade_LU-SHUI.ps

psxy lw.txt -JM6i -R121.4/121.7/24.0/24.3 -Sc0.05 -W0/200/255 -MX -W -V -H -K -O >>shade_LU-SHUI.ps

psxy temp1.txt -JM6i -R121.4/121.7/24.0/24.3 -Sa0.4 -G255/0/0 -V -H -K -O >> shade_LU-SHUI.ps

pscoast -JM6i -R121.4/121.7/24.0/24.3 -Dh -W30 -Na -Ia -W10/0/0/255 -T123/34/1 -L123/32.5/27/100k -V -K -O >> shade_LU-SHUI.ps

psscale -D17.5c/11.6c/10c/0.5c -Clushui_test.cpt -B500 -V -O >> shade_LU-SHUI.ps

------------------------------------------------------------------------------------------------

以地形高低和立霧溪走向來看，原始立霧溪的資料點位置似乎有點偏差。

經過修正後，（經度增加0.018，緯度增加0.045），得到上圖，看起來比較能符合地形。但最終應該還是得靠現地測量或是有其他可靠來源，才能驗證立霧溪資料正確性。




4.以3D立體圖來呈現綠水觀測站附近的地形高低起伏
------------------------------------------------------------------------------------------------
makecpt -Crainbow -T0/3500/100 -V -Z > high.cpt

grdgradient twdtm_small.grd -Gtwdtm_small_temp_shade.grd -A0 -Ne0.9 -V

psbasemap -P -JM6i -JZ2i -R121.3/121.55/24.0/24.3/0/3000 -B0.1/0.08/1000:."LU-SHUI": -E180/60 -T121.45/24.32/0.5 -V -K > 3D_LU-SHUI.ps

grdview twdtm_small.grd -JM6i -JZ2i -R121.3/121.6/24.0/24.3/0/3000 -Chigh.cpt -Itwdtm_small_temp_shade.grd -Qs -E180/60 -N0 -V -K -O >> 3D_LU-SHUI.ps

psxyz temp1.txt -JM6i -JZ6i -R121.3/121.6/24.0/24.3/0/3000 -Sa0.4 -G255/0/0 -V -H -K -O >> 3D_LU-SHUI.ps

psscale -D3/-0.2/5/0.15h -Chigh.cpt -B500 -V -K -O >> 3D_LU-SHUI.ps

------------------------------------------------------------------------------------------------

可以看出綠水觀測站附近的地形高低起伏落差相當大。

5.檢視綠水觀測站附近相關的地震，圓點的顏色代表深度（單位：公里），圓點大小代表地震規模。可以看出地震的發生的立體空間分布與規模。

------------------------------------------------------------------------------------------------
5.1 3D圖
從（http://tecdc.earth.sinica.edu.tw/TEC/db/）設定好選取位置的經緯度，選csv輸出
貼到txt檔，再用excel轉檔（排好欄位）後儲存為hijmsus.txt，qawk才能使用
經緯度範圍121/122/24/25，時間從1900到2009.02.28


gawk "{print $6, $5, -$7, -$7, 0.0003*$8^4}" hijmsus.txt > eq_lu-shui.txt

makecpt -Crainbow -T-150/0/3 -V -Z > eq.cpt
psbasemap -JM8i -Jz0.06c -R121/122/24/25/-150/0 -E175/15 -B0.2:Longitude:/0.2:Latitude:/30:Depth::."LU-SHUI EARTHQUAKE 3-D PLOT":10 -V -K > eq.ps
pscoast -JM8i -Jz0.06c -R121/122/24/25/-150/0 -G162/205/90 -E175/15 -Df -W -V -K -O >> eq.ps
psxyz eq_lu-shui.txt -JM8i -Jz0.06c -R121/122/24/25/-150/0 -E175/15 -Sc -Ceq.cpt -V -K -O >> eq.ps
psscale -D6c/-0.2/4/0.2h -Ceq.cpt -B50 -V -O >> eq.ps

上圖左下角有五個圓形，是為了比對地震規模之用，從上到下的地震規模分別為4、5、6、7、8。
------------------------------------------------------------------------------------------------
5.2 2D俯視圖
gawk "{print $6, $5, -$7, -$7, 0.0003*$8^4}" hijmsus2.txt > eq_lu-shui.txt

makecpt -Crainbow -T-150/0/3 -V -Z > eq.cpt
psbasemap -JM5i -Jz0.06c -X3 -R121/122/24/25/-150/0 -E180/90 -B0.2:Longitude:/0.2:Latitude:/30:Depth::."LU-SHUI EARTHQUAKE 2-D PLOT":6 -V -K > eq2.ps
pscoast -JM5i -Jz0.06c -R121/122/24/25/-150/0 -G162/205/90 -E180/90 -Df -W -V -K -O >> eq2.ps

psxyz eq_lu-shui.txt -JM5i -Jz0.06c -R121/122/24/25/-150/0 -E180/90 -Sc -Ceq.cpt -V -K -O >> eq2.ps

rem 以下三行標註和平溪、立霧溪、美崙溪的位置
psxy hp.txt -JM5i -R121/122/24/25 -MX -W -V -H -K -O >> eq2.ps
psxy lw.txt -JM5i -R121/122/24/25 -MX -W -V -H -K -O >> eq2.ps
psxy mk.txt -JM5i -R121/122/24/25 -MX -W -V -H -K -O >> eq2.ps

rem 以下標註綠水觀測站的位置
psxy temp1.txt -JM5i -R121/122/24/25 -Sa0.4 -G0/0/0 -V -H -K -O >> eq2.ps

psscale -D15c/2.5/5/0.2 -Ceq.cpt -B50 -V -O >> eq2.ps

從以上兩圖可以看出：該地的地震多發生在海陸交界處，且多為淺層地震。其餘地震發生在中央山脈地區，但相較而言，為數不多。

------------------------------------------------------------------------------------------------




6.希望由「震源機制解」圖形來觀察綠水觀測站附近地區的地震型態。
------------------------------------------------------------------------------------------------
EQ_CMT_Harvard_2006_Taiwan.dat由 裡設定選取區域的經緯度範圍
output格式選GMT psmeca input，貼在檔案裡，不用再轉檔，就直接可以用在psmeca裡
有關震源機制解可以參考裡的解說

psbasemap -JM6i -R121/122/24/25 -X3 -Y0.7 -B0.5:."LU-SHUI FOCAL MECHANISM SOLUTION PLOT": -V -K> focal_mech_LU-SHUI.ps
pscoast -JM6i -R121/122/24/25 -G162/205/90 -Df -W1 -V -K -O >> focal_mech_LU-SHUI.ps

rem -Sd0.7c/7 前面的0.7表示海灘球的大小尺寸，後面的/7表示海灘球註解文字的大小
psmeca eq_CMT_LU-SHUI.txt -JM6i -R121/122/24/25 -Sd1c/9 -H -V -K -O >> focal_mech_LU-SHUI.ps

rem 以下三行標註和平溪、立霧溪、美崙溪的位置
psxy hp.txt -JM6i -R121/122/24/25 -MX -W -V -H -K -O >> focal_mech_LU-SHUI.ps
psxy lw.txt -JM6i -R121/122/24/25 -MX -W -V -H -K -O >> focal_mech_LU-SHUI.ps
psxy mk.txt -JM6i -R121/122/24/25 -MX -W -V -H -K -O >> focal_mech_LU-SHUI.ps

rem 以下標註綠水觀測站的位置
psxy temp1.txt -JM6i -R121/122/24/25 -Sa0.4 -G255/0/0 -V -H -O >> focal_mech_LU-SHUI.ps

從上圖可以看出，此一地區的地震型態相當複雜，包含了正斷層、逆斷層、平移斷層、平移與逆斷層...等等的複合地震型態。