征集令
维护者征集 下面列出了我发起的几个开源项目。挖了坑但是没有时间去填，所以需要征集志愿者参与维护。志愿者越多越好，项目开源的目的就是让所有人都可
判断 SAC 数据是否已去除仪器响应
有些时候，波形数据拿到手了，但是却不知道波形数据是否有被处理过，尤其是有没有去除仪器响应。
先说结论，判断一个 SAC 数据是否已经去除仪器响应的准则是：
若一个 SAC 数据中每个数据点的值都是整数，则这个数据的仪器响应没有被去除的概率为 99.9%
在极坐标下绘制深度剖面
在做大区域层析成像时，通常需要绘制深度剖面。一个简单的示例如下：
这样的深度剖面图实际上是在极坐标下绘制的。先直接给出绘制上图所使用的命令:
gmt psbasemap -R60/80/5571/6371 -JPa15c/70z -Bx5 -By200 -BWNse > test.ps
用 SAC 快速拾取震相
震相拾取，是地震学的一个基本工作。SAC 提供了 plotpk 命令，可以用于手动拾取
震相到时，用起来还算是比较简单方便。
绘制颜色渐变的线条
想要绘制一条颜色变化的线段，下面是演示代码：
#!/bin/bash
# GMT v5.2.1
gmt makecpt -Crainbow -T-2/2/1 > lines.cpt
gmt psxy -JX15c/4c -R0/6/0/4 -B1 -Clines.cpt -W2p > test.ps << EOF
> -Z-1.5
1 2
2 2
> -Z-0.5
2 2
3 2
> -Z0.5
3 2
4 2
> -Z1.5
4 2
5 2
EOF
绘图效果如下图：
