数据下载: mt.dat

M-T图表示地震随时间发展的数量密度和震级大小。M-T图的横轴是时间,纵轴是震级。图中每一条竖线代表一次地震,竖线的位置表示了地震发生的时间,而竖线的长度代表了震级的大小。

绘图效果如下:

M-T图

绘图脚本如下:

#!/bin/bash
PS=mt.ps
J=X15c/5c
R=2017-02-10T00:00:00/2017-05-04T00:00:00/0/5

# 设置 SECONDARY TICK 长度为 0
gmt set MAP_TICK_LENGTH_SECONDARY 0p
gmt set FORMAT_DATE_MAP yyyy-mm

# -Sb0.005c 中 0.005c 代表每条线的宽度
# -BSW 只显示左边和下边坐标
# -Bsxa1O 设置X轴次级标注的间隔为1个月,注意最后一个字符是大写字母O不是数字0
# -Bpxa7Rf1d 设置X轴主级标注的间隔为7天,刻度间隔为1天
# -B+l"Date" 设置X轴标签
# -Bya1 设置Y轴标注间隔为1
# -B+l@%2%M@-L@-@%% 设置Y轴标签,使用2号字体并将L设置为下标
gmt psxy mt.dat -R$R -J$J -Sb0.005c -Gblack \
    -BSW -Bsxa1O -Bpxa7Rf1d -Bx+l"Date" -Bya1 -By+l@%2%M@-L@-@%% > $PS
rm gmt.*