1、自定义坐标轴

gmt basemap -R%R% -J%J% -Bycyannots.txt  -Bxa50f10+l"Distance" -BWES

-Bc使用自定义坐标轴，yannots.txt是标注的内容：

-4000 af  -20km

-2000 af  -10km

0     af    0m

2000  af    2km

4000  af    4km

2、awk筛选数据

awk "$1> 99 && $1< 101.9 &&  $2 > 37.3 &&  $2 < 38{print  $1, $2,$3}" earthquake_data.txt |
gmt plot -R%R% -J%J% -W0.3p,blue -Sc

3、剖线投影：-S 进行排序，-Fpz 只输出距离和深度

gmt project seismic_data.txt -C99.82/36.2 -E101.38/37.92 -S -Fpz -Q > seismic_profile.txt