问题：想要在地图上只绘制某条经线或 / 和某条纬线，该如何实现？

GMT 中的 -B 选项可以用于控制是否绘制网格线以及网格线的间隔。如果精心设计 -B 选项中的网格间隔以及第一条网格线的偏移量，是可以使得整张图中只有一条经线或 / 和纬线并保证该经线 / 纬线落在合适的位置，但是这种方法相对来说比较麻烦且很难用于批量绘图。

SOFI2D 是一个二维 PSV 交错网格有限差分代码。要使用该程序，至少需要阅读如下几篇参考文献：

• Virieux, J. 1986. P-SV wave propagation in heterogeneous media: velocity-stress finite- difference method. Geophysics, 51(4), 889–901.
• Levander, A.R. 1988. Fourth-order finite-difference P-SV seismograms. Geophysics, 53(11), 1425–1436.
• SOFI2D User Guide

具体的用法及每个参数的意义在官方文档已经介绍的很详细，本文是我使用 SOFI2D 过程中的经验总结。

注：本文使用的是 2015 年 10 月 28 日发布的 SOFI2D。

地震学中，经常需要绘制事件到台站的连线或台站与台站之间的连线。连线沿着大圆路径，理论上用 psxy 命令即可实现，实际绘图时要更复杂一些。

申请地震波形数据，按照具体的需求，大致可以分为两大类，即事件波形数据和连续波形数据。

这两者的主要区别在于如何确定要申请的数据时间窗，前者需要震相到时信息，后者则不需要。