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 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