Question

我在GMT为Southamerica的北部制作了一些焦点机制（FM），我需要将它们放在地图上。

对于FM我有一个85级地震的文件，这是代码：

psmeca eventos2.dat -R-80/-66/-5/14 -JM6i -Sm0.4/0/0.05u -P -V > coltopomap.ps
GMT psconvert coltopomap.ps -Tf

这是地图

GMT grdcut topo15.grd -R-80/-66/-5/14 -Gtopo15_cut_col.grd
GMT grdgradient topo15_cut_col.grd -Gtopo15_cut_shaded.grd -A310 -Ne0.6
GMT grdimage topo15_cut_col.grd -JM6i -R-80/-66/-5/14 -Cglobe -P -Ba --MAP_FRAME_TYPE=plain --MAP_GRID_CROSS_SIZE_PRIMARY=0p -Itopo15_cut_shaded.grd -K > coltopomap.ps 
GMT pscoast -R-80/-66/-5/14 -J -Df -N1/1,black -W1/0.5,black -Lf-68/-4/-68/-4/200k+l -Tdg-78/12+w0.3i+f3+lW,E,S,N --FONT_ANNOT=10p --FONT_LABEL=10p --FONT_TITLE=10p -O >> coltopomap.ps

要加入他们，我已尝试

GMT grdcut topo15.grd -R-80/-66/-5/14 -Gtopo15_cut_col.grd
GMT grdgradient topo15_cut_col.grd -Gtopo15_cut_shaded.grd -A310 -Ne0.6
GMT grdimage topo15_cut_col.grd -JM6i -R-80/-66/-5/14 -Cglobe -P -Ba --MAP_FRAME_TYPE=plain --MAP_GRID_CROSS_SIZE_PRIMARY=0p -Itopo15_cut_shaded.grd -K > coltopomap.ps 
GMT pscoast -R-80/-66/-5/14 -J -Df -N1/1,black -W1/0.5,black -Lf-68/-4/-68/-4/200k+l -Tdg-78/12+w0.3i+f3+lW,E,S,N --FONT_ANNOT=10p --FONT_LABEL=10p --FONT_TITLE=10p -O >> coltopomap.ps
GMT psmeca eventos2.dat -R-80/-66/-5/14 -JM6i -Sm0.4/0/0.05u -P -V >> coltopomap.ps
GMT psconvert coltopomap.ps -Tf

但我唯一得到的是一页中的地图和另一页中的FM，如下图所示：

https://i.stack.imgur.com/bAvx7.png

Answer 1

您似乎错过了-K命令中的GMT pscoast ...标记。

来自GMT documentation：

构建多层时，除了最后一个绘图命令之外的所有绘图命令都是必需的图

如何在通用映射工具（GMT）上进行FOCAL MECHANISMS？

1 个答案: