GMT的操作有自己的邏輯和“命令”,但GMT是可以用Bash語言控制的,所以常常以.sh為后綴寫GMT程序。
GMT程序運行步驟如下:
- 采用cd ,定位到指定文件夾;
- 以sh ***.sh運行GMT,得到結果。
另外,遇到看不懂的GMT時,多學學Bash語言的命令,兩者結合可以滿足絕大多數需求。
#!/bin/bash# Output figure name
OUTFIG="world_eq_map"# Map region and projection
REGION="-180/180/-90/90"
PROJ="W15c"# Inline earthquake coordinates using cat (longitude latitude)
cat << EOF > eq_data.txt
138.5 35.6
142.9 38.3
95.9 3.3
-73.0 -36.1
EOF# Start GMT plotting
gmt begin $OUTFIG pngecho "Drawing base map..."gmt basemap -R$REGION -J$PROJ -Bafecho "Drawing coastlines..."gmt coast -R$REGION -J$PROJ -W0.25p -A1000 -Glightgray -Slightblueecho "Plotting earthquakes..."gmt plot eq_data.txt -R$REGION -J$PROJ -Sc0.3c -Gred -W0.25pgmt end show# Clean up (optional)
rm eq_data.txt