GTOPO30，即global digital elevation model (DEM)，水平格网间隔为弧度30'（最大有1km）。数据块在全球的分布如下图：中国区域包含了四个DEM文件E060N90、E100N90、E060N40、E100N40。

    例如，对于地震多发的四川地区，经纬度包含在lat26-lat34,lon97-lon110，使用GTOPO30画四川地区的grd文件就需要将E60N40、E100N40两个文件制作的grd文件合并。GMT-4.5.12中可以使用grdpaste命令合并。grdpaste用法如下：

usage: grdpaste <file_a> <file_b> -G<outfile> [-V] [-f[i|o]<colinfo>]

    where file_a and file_b are to be combined into outfile.
    file_a and file_b must have same dx,dy and one edge in common.
    If in doubt, run grdinfo first and check your files.
    Use grdcut and/or grdsample to adjust files as necessary.

在运行paste时候，即使是两个相邻的有公共边的文件，有时会报错，例如：

#！ /bin/tcsh -f

set rangeEN="100/140/40/60"
set rangeES="100/140/0/40"
set rangeWS="60/100/0/40"
set rangeWN="60/100/40/60"
grdraster 6 -I0.5m -R$rangeEN -GChinaEN.grd
grdraster 7 -I0.5m -R$rangeES -GChinaES.grd
grdraster 8 -I0.5m -R$rangeWS -GChinaWS.grd
grdraster 9 -I0.5m -R$rangeWN -GChinaWN.grd

grdpaste ChinaES.grd ChinaEN.grd -GChinaE.grd              
grdpaste ChinaWN.grd ChinaWS.grd -GChinaW.grd          
grdpaste ChinaW.grd ChinaE.grd -GChina.grd                    

运行该gmt脚本，会报错：              

grdpaste:  Grids do not share a common edge!

此时，就很纳闷了，明明是有公共边的却报错。

尝试在grdpaste中 a.grd与b.grd的顺序调换，则程序不会报错。

此处改为：grdpaste ChinaEN.grd ChinaES.grd -GChinaE.grd              #ok

结论：这是gmt-4.5.12的一个bug

如果语句 grdpaste <file_a> <file_b> -G<outfile> 报错 Grids do not share a common edge!

可互换<file_a>与<file_b>的位置。

将文件E060N90、E100N90、E060N40、E100N40合并成一个大grd文件：

set rangeEN="100/140/40/60"
set rangeES="100/140/0/40"
set rangeWS="60/100/0/40"
set rangeWN="60/100/40/60"
grdraster 6 -I0.5m -R$rangeEN -GChinaEN.grd
grdraster 7 -I0.5m -R$rangeES -GChinaES.grd
grdraster 8 -I0.5m -R$rangeWS -GChinaWS.grd
grdraster 9 -I0.5m -R$rangeWN -GChinaWN.grd

grdpaste ChinaEN.grd ChinaES.grd -GChinaE.grd              #ok
grdpaste ChinaWN.grd ChinaWS.grd -GChinaW.grd           #ok
grdpaste ChinaW.grd ChinaE.grd -GChina.grd                  #ok