本文所说的 3D 数据是指 XYZD 数据。常见的情况是数据中包含四列：经度、纬度、深度、值。

GMT 几乎是不可以处理 3D 数据的，比如与网格相关的命令只支持 2D netCDF 文件。另外，GMT 绘制 3D 图的效果也很一般，因而通常情况下都需要从 3D 数据中截取剖面，绘制剖面图。

本文试着介绍如何从 3D 数据中截取自己想要的剖面。

在用 GMT 绘制底图时，通常 Y 轴的标注都是水平的，比如如下绘制底图的命令:

gmt psbasemap -R0/50/-10/10 -JM10c -Bx10 -By10 > test.ps

效果如下图

有些时候想要让 Y 轴的标注逆时针旋转 90 度，即 Y 轴的标注方向与 Y 轴平行而不是与 X 轴平行。想要实现的效果如下图。

fk 是 Prof. Lupei Zhu 写的一个用于计算 水平分层模型 下的理论格林函数并合成 理论地震图的代码包。代码是开源的，可以直接编译使用。

本文介绍几个专业名词：

• 慢度
• 水平慢度 / 垂直慢度
• 射线参数
• 视速度 / 视慢度

这些名词其实很简单，但是在大三一直到研一的几年间，这几个名词在我脑中的概念是 有些混乱的，我把混乱的原因归结于不同课本对同一个专业名词乱用了不同的词汇。 故而在这里总结一下。