GMT绘图笔记：用深度作为表面高度（Z 值），用其他物理量（泊松比）给表面着色

之前用GMT绘制莫霍面的三维示意图是用的莫霍面的深度作为表面着色。

GMT绘图笔记：绘制堆叠三维图_gmt画图-优快云博客

如果要换成其他的物理场，比如泊松比，则需要使用以下的代码

gmt begin BS_figures2 png,pdf E600
# 绘制底图
# Chistrong Wen
# University of Stuttgart
# 2025-3-12 11:14
gmt set FORMAT_GEO_MAP = ddd:mm:ssF
gmt set MAP_FRAME_TYPE = plain
gmt set MAP_FRAME_PEN = 1p,black
gmt set MAP_TICK_PEN = 1p,black
gmt set FONT_ANNOT_PRIMARY = 13p,Helvetica,black
gmt set FONT_ANNOT_SECONDARY = 13p,Times-Roman,black
gmt set MAP_TICK_LENGTH = 1.5p
gmt set LABEL_FONT = 12p
gmt set FONT_LABEL = 12p
gmt set LABEL_FONT_SIZE = 12p
gmt set ANNOT_FONT_SIZE_PRIMARY = 12p
gmt set MAP_FRAME_WIDTH = 2p 
gmt set FRAM_WIDTH = 1c
gmt set MAP_ANNOT_OFFSET_PRIMARY = 0.15c
gmt set MAP_TICK_LENGTH_PRIMARY = 0.2c
gmt set FORMAT_ANNOT_PRIMARY = 2p,Times-Roman,black

gmt grdsample @earth_relief_06m -Gdem.grd -I0.5
gmt xyz2grd depthtomoho1.txt -Rg -I1.1 -Gm1.grd
gmt grdsample m1.grd -Ga11.grd -I0.5
gmt xyz2grd depthtomoho.txt -Rg -I1.1 -Gm2.grd
gmt grdsample m2.grd -Gb11.grd -I0.5

gmt grd2cpt a11.grd -Cjet -T20/35/0.5
gmt basemap -R70/80/20/25 -JM13.3c -BWSne -Bxa2f2 -Bya2f2 -Xc -Yc
gmt grdimage a11.grd -R70/80/20/25 -C
gmt colorbar -Baf1
gmt basemap -R70/80/20/25 -JM13.3c -BWSne -Bxa2f2 -Bya2f2 -X20c -Yc
gmt grdimage b11.grd -R70/80/20/25 -Cjet
gmt colorbar -Baf1
gmt grd2cpt a11.grd -Cjet -T20/35/0.5
gmt grdview b11.grd -JX12c -JZ5c -R70/80/20/25 -Bxa2f2 -Bya2f2 -Bza500f500+l"Moho Depth (km)" \
                    -p135/40 -Qs -Ga11.grd -I+d -C -X-22c -Y-16c
gmt grdview b11.grd -JX12c -JZ5c -R70/80/20/25 -Bxa2f2 -Bya2f2 -Bza500f500+l"Moho Depth (km)" \
                    -p135/40 -Qs -I+d -Cjet -X20c -Y0c

gmt end 

绘制效果：图1是泊松比的2D图，图2是莫霍面的2D图，图3是采用泊松比作为表面着色，图4是使用莫霍面深度着色。

❤️欢迎点赞收藏❤️