示例数据: 速度模型数据

本文展示了如何绘制地壳速度模型图。示例数据分为 6 列,分别为 地层厚度,S波速度,P波速度,密度,S波Q值和P波Q值, 绘图脚本只用到了前 3 列数据。

Bash 绘图脚本如下:

#!/bin/bash
gmt set MAP_TICK_LENGTH_PRIMARY -0.1c
model=model.txt
PS=model.ps
gmt psbasemap -R0/10/0/50 -JX8c/-10c -Bx1+l"Velocity (km/s)" -By10+l"Depth (km)" \
                                       -BWS+t'Crust model' -K -P > $PS
gmt psbasemap -R -J -Ben -B0 -K -O >> $PS
# 利用 awk 计算每一层上下边界的深度,提取该层的速度
# 每一层上边界的深度即为前几层的厚度累加,下边界的深度为上边界的深度加该层层厚
# Vs
awk 'BEGIN{d1=0;d2=0}{d2=d2+$1;print $2,d1;print $2,d2;d1=d2;}' $model | \
gmt psxy -R -J -W1.5p,blue -K -O >> $PS
# Vp
awk 'BEGIN{d1=0;d2=0}{d2=d2+$1;print $3,d1;print $3,d2;d1=d2;}' $model | \
gmt psxy -R -J -W1.5p,green,6_2:2 -K -O >> $PS
gmt pslegend -R -J -DjBL+w2c+o0.2c/0.2c -F -O >> $PS <<EOF
S 0.4c - 0.5c - 1p,blue 1c @;blue;Vs@;;
S 0.4c - 0.5c - 1p,green,6_2:0 1c @;green;Vp@;;
EOF
gmt psconvert $PS -A -Tg -Z
rm gmt*

绘图效果如下:

绘制速度模型