GMT之Bash语言使用

GMT的操作有自己的逻辑和"命令"，但GMT是可以用Bash语言控制的，所以常常以.sh为后缀写GMT程序。

GMT程序运行步骤如下：

1. 采用cd ，定位到指定文件夹；
2. 以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 png

  echo "Drawing base map..."
  gmt basemap -R$REGION -J$PROJ -Baf

  echo "Drawing coastlines..."
  gmt coast -R$REGION -J$PROJ -W0.25p -A1000 -Glightgray -Slightblue

  echo "Plotting earthquakes..."
  gmt plot eq_data.txt -R$REGION -J$PROJ -Sc0.3c -Gred -W0.25p

gmt end show

# Clean up (optional)
rm eq_data.txt