triangulate table_5.11 –Grawt.grd –R –I0.1 > NUL 执行Delauney三角型法构筑grd文件（尽可能等边）， grdfilter rawt.grd –Gfiltered.grd –D0 –Fc1
读入grd文件rawt.grd，按时间域，用多种方法进行筛选 –D0：x,y网格与宽度单位一致，采用笛卡
儿距离； –Fc1：采用余弦Arch进行筛选；输出文件filtered.grd
grdview 绘制并标注右下图框，由网格文件filtered.grd 及色谱文件ex16.cpt绘制等值线并填充。 echo 3.25 7 ?? 写字符串triangulate→grdfilter，底中对齐
echo 3.2125 7.5 ??写图名字符串Gridding of Data，底中对齐（疑为CT，顶中对齐） psscale 画图例，绘制系统时间标志。
–––––––––––––––––––––––––––––––––––––––– Figure 6.18 体积计算和特定地区的圈定
gmtset ELLIPSOID Sphere
REM Define location of Pratt seamount echo -142.65 56.25 > pratt.d
makecpt -Crainbow -T-60/60/10 -Z > grav.cpt
grdgradient ss_grav.grd -Nt1 -A45 -Gss_grav_i.grd
grdimage ss_grav.grd -Iss_grav_i.grd -JM5.5i -Cgrav.cpt -B2f1 -P -K -E100 -X1.5i -Y5.85i > example_18.ps pscoast1 -R-149/-135/52.5/58 -JM -O -K -Di -G160 -W0.25p >> example_18.ps psscale -D2.75i/-0.4i/4i/0.15ih -Cgrav.cpt -B20f10/:mGal: -O -K >> example_18.ps echo {print $1, $2, 12, 0, 1, \
gawk1 -f t pratt.d | pstext -R -JM -O -K -D0.1i/0.1i >> example_18.ps
gawk2 \
grdcontour1 ss_grav.grd -JM -C20 -B2f1WSEn -O -K -Y-4.85i -U/-1.25i/-0.75i/\example_18.ps
grdcontour2 ss_grav.grd -JM -C10 -L49/51 -O -K -Dsm -Wc0.75p/0/255/0 >> example_18.ps pscoast2 -R -JM -O -K -Di -G160 -W0.25p >> example_18.ps gawk3 \-O -K -SE -W0.25p >> example_18.ps REM Only consider closed contours
REM Make simple gawk script to calculate average position of locations
echo BEGIN { > center.awk echo x = 0 >> center.awk echo y = 0 >> center.awk echo n = 0 >> center.awk echo } >> center.awk echo { >> center.awk
echo x += $1 >> center.awk echo y += $2 >> center.awk echo n++ >> center.awk echo } >> center.awk
echo END { >> center.awk
echo print x/n, y/n >> center.awk echo } >> center.awk
36
REM Now determine centers of each enclosed seamount exceeding 50 mGal but REM only plot the ones within 200 km of Pratt seamount.
REM Now determine mean location of each closed contour and REM add it to the file centers.d using center.awk script
for %%f in (sm_*.xyz) do gawk -f center.awk %%f >> centers.d REM Only plot the ones within 200 km
gmtselect -R -JM -C200/pratt.d centers.d > tmp2
psxy1 tmp2 -R -JM -O -K -SC0.07i -G255/0/0 -W0.25p >> example_18.ps psxy2 -R -JM -O -K -ST0.1i -G255/255/0 -W0.25p pratt.d >> example_18.ps REM Then report the volume and area of these seamounts only REM by masking out data outside the 200 km-radius circle REM and then evaluate area/volume for the 50 mGal contour Grdmath1 -R -I2m -F -142.65 56.25 GDIST = mask.grd grdclip mask.grd -Sa200/NaN -Sb200/1 -Gmask.grd grdmath2 ss_grav.grd mask.grd MUL = tmp.grd echo -148.5 52.75 > tmp1 echo -140.5 52.75 >> tmp1 echo -140.5 53.75 >> tmp1 echo -148.5 53.75 >> tmp1
psxy3 -R -JM -A -O -K -L -W0.75p -G255 tmp1 >> example_18.ps
echo {printf \53.08 14 0 1 LM Areas: %%s km@+2@+\\n-148 53.42 14 0 1 LM Volumes: %%s km@+3@+\\n\
grdvolume tmp.grd -C50 -Sk | gawk -f t1 | pstext -R -JM -O >> example_18.ps del .gmt*
部分命令注释
gmtset ELLIPSOID Sphere 改缺省的ELLIPSOID WGS-84为Sphere
echo -142.65 56.25 > pratt.d 把Prett海山位置（经纬度）写入文件pratt.d
makecpt -Crainbow -T-60/60/10 -Z > grav.cpt 以色谱文件库的rainbow为标准建立连续色谱文件grav.cpt grdgradient ss_grav.grd -Nt1 -A45 -Gss_grav_i.grd
由ss_grav.grd计算45°方向的方向梯度文件ss_grav_i.grd。–Nt见 Figure 6.2注释。 ss_grav.xyz 文件如下
-148.983 57.9833 99.8354 -148.95 57.9833 98.0651 ???????
ss_grav_i.xyz文件如下
-148.983 57.9833 -0.449835 -148.95 57.9833 -0.520409 ???????
注意两个网格文件x的增量均为30分之1度，即2分。
grdimage用ss_grav.grd ss_grav_i.grd 和grav.cpt绘制影象图(上图)，设置grd的分辨率为100dpi。 pscoast1 绘制、标注上图图框。
psscale 绘制、标注图例条。以图左下角为原点，水平图例框（h）中心点X坐标（该图宽5.5i，图例
中心X坐标5.5/2=2.75i）；图框中心的Y坐标（图框下边界下方0.4i），图例的长、宽（4i、0.15i）。图例的标记刻度非标记刻度（20、10），图例的右边标“mGal”，。 echo {print $1, ??把“Pratt”写入文件t
37
gawk1 根据文件t和数据文件pratt.d写“Prett”海山名字。
gawk2 画椭圆命令，6个参数，中心经、纬度、角度、长轴、短轴长度（km）。例中，椭圆中心坐标由
pratt.d里读出，由于长轴、短轴长度一样，实际是以Prett海山为圆心，画半径200km 的圆。 grdcontour1 按照ss_grav.grd绘制全部等值线图，间距20。
grdcontour2 按照ss_grav.grd绘制等值线图，间距10，限定只画49-51范围内的等值线（因此只有50
一条），分别将每一条50等值线的x、y、z值建立一个xyz文件存盘，用sm_50_n_i来命名，其中n为序号，i表示为封闭的等值线。等值线为绿色（0/255/0）。 pscoast2 绘制、标注下图图框。
gawk3 以Prett海山位置为圆心，画半径200km 的大圆。
echo x = 0 ?? 用13行echo写一段gawk程序，存入center.awk。程序为 BEGIN { x = 0 y = 0 n = 0
} （变量赋初值） {
x += $1 （x坐标累加） y += $2 （y坐标累加） n++ （n记数器累加）
}
END { print x/n, y/n （求算并输出等值线中心经纬度坐标）
}
for %%f in (sm_*.xyz) do gawk -f center.awk %%f >> centers.d
顺序读取文件sm_*.xyz，执行gawk，即把所有封闭的超过50 mGal等值线中心坐标（由执行center.awk得）建立一个数据文件centers.d。本例中有33个等值线中心坐标满足条件。 gmtselect -R -JM -C200/pratt.d centers.d > tmp
筛选，从centers.d中选择距Pratt，即大圆圆心的距离在200公里以内的封闭50等值线。存入临时文件tmp。本例中有11个等值线中心坐标满足条件。
psxy1 以tem中的坐标为圆心绘制直径0.04i的红色（255/0/0）圆 psxy2 以Pratt海山为中心，画边长0.1i的黄色（255/255/0）三角形 grdmath -R -I2m -F -142.65 56.25 GDIST = mask.grd
以2分为网格间隔，计算由Pratt海山到每个节点的距离。生成一个新的网格文件mask.grd grdclip mask.grd -Sa200/NaN -Sb200/1 -Gmask.grd
对mask.grd进行处理，将>200的节点改为空白，将<200的节点改值为1。处理后的网格文件另存为mask.grd
grdmath ss_grav.grd mask.grd MUL = tmp.grd 两个网格文件相乘，生成临时网格文件tmp.grd echo -148.5 52.75 ??
psxy ?? 用上面4行echo给出的坐标画矩形框
echo {printf ?? 用echo写两个字串Areas: km2和Volumes: km3，存入ft，分别作为面积和体积的
标注。
grdvolume tmp.grd -C50 -Sk | gawk -f t | pstext -R -JM -O >> example_18.ps
计算给定条件，本例为距Pratt海山200km以内，等值线z=50以上的（山体）体积和表面积。单位度转换为km。计算结果通过管道命令作为gawk -f t的输入，即分别给体积和表面积赋值；最后通过pstext命令输出到.ps文件。
38
–––––––––––––––––––––––––––––––––––––––– Figure 6.19 彩色图案模版的使用
gmtset COLOR_MODEL rgb
grdmath1 -R-180/180/-90/90 -I1 -F Y COSD 2 POW = lat.grd grdmath2 -R-180/180/-90/90 -I1 -F X = lon.grd echo 0 255 255 255 1 0 0 255 > lat.cpt
makecpt -Crainbow -T-180/180/60 -Z > lon.cpt
grdimage1 lat.grd -JI0/6.5i -Clat.cpt -P -K -Y7.5i -B0 > example_19.ps pscoast1 -R -JI -O -K -Dc -A5000 -Gc >> example_19.ps grdimage2 lon.grd -JI -Clon.cpt -O -K >> example_19.ps pscoast2 -R -JI -O -K -Q >> example_19.ps
pscoast3 -R -JI -O -K -Dc -A5000 -W0.25p >> example_19.ps
echo 0 20 32 0 1 CM FIRST INTERNATIONAL | pstext -R -JI -O -K -G255/0/0 -S0.5p >> example_19.ps echo 0 -10 32 0 1 CM GMT CONFERENCE | pstext -R -JI -O -K -G255/0/0 -S0.5p >> example_19.ps
echo 0 -30 18 0 1 CM Honolulu, Hawaii, April 1, 2000 | pstext -R -JI -O -K -G0/255/50 -S0.25p >> example_19.ps
REM Then show example of color patterns
pscoast4 -R -JI -O -K -Dc -A5000 -Gp100/86:F255/0/0B255/255/0 -Sp100/7:F255/0/0B0/0/0 -B0 -Y-3.25i >> example_19.ps
echo 0 15 32 0 1 CM SILLY USES OF | pstext -R -JI -O -K -G50/255/50 -S0.5p >> example_19.ps
echo 0 -15 32 0 1 CM GMT COLOR PATTERNS | pstext -R -JI -O -K -G255/0/255 -S0.5p >> example_19.ps REM Finally repeat 1st plot but exchange the patterns
grdimage3 lon.grd -JI -Clon.cpt -O -K -Y-3.25i -B0 -U\pscoast5 -R -JI -O -K -Dc -A5000 -Gc >> example_19.ps grdimage4 lat.grd -JI -Clat.cpt -O -K >> example_19.ps pscoast6 -R -JI -O -K -Q >> example_19.ps
pscoast7 -R -JI -O -K -Dc -A5000 -W0.25p >> example_19.ps
echo 0 20 32 0 1 CM FIRST INTERNATIONAL | pstext -R -JI -O -K -G255/0/0 -S0.5p >> example_19.ps echo 0 -10 32 0 1 CM GMT CONFERENCE | pstext -R -JI -O -K -G255/0/0 -S0.5p >> example_19.ps
echo 0 -30 18 0 1 CM Honolulu, Hawaii, April 1, 2000 | pstext -R -JI -O -G0/255/50 -S0.25p >> example_19.ps del .gmt*
GMT EXAMPLE 19注释
Grdmath1 -R-180/180/-90/90 -I1 -F Y COSD 2 POW = lat.grd
按照表达式cos2 (y)计算每1度纬度（y）的换算值，作为z值，写入lat.grd -179.5 89.5 7.61524e-005 -178.5 89.5 7.61524e-005 ???????
grdmath2 -R-180/180/-90/90 -I1 -F X = lon.grd
把每1度经度（x）的值，作为z值，写入lon.grd
39
-179.5 89.5 -179.5 -178.5 89.5 -178.5 ???????
echo 0 255 255 255 1 0 0 255 > lat.cpt
makecpt -Crainbow -T-180/180/60 -Z > lon.cpt 建立2个.cpt文件
grdimage1 采用正弦曲线混合投影画上面影像图。由lat.grd和lat.cpt绘制彩色影象，z值0、1对应颜色白、蓝。
pscoast1 绘制海陆界线，-Gc 剪裁去陆地区，即只填充水体彩色影像。
grdimage2 由lon.grd和lon.cpt绘制陆地彩色影象，z值由-180度到180对应彩虹色紫到红。 pscoast2 -Q 恢复pscoast1剪裁去的陆地区 pscoast3 绘制海陆界线，填充陆地彩色影象。
echo 0 20 ??分别书写“FIRST INTERNATIONAL” 、“GMT CONFERENCE” 和“Honolulu, Hawaii,
April 1, 2000”3行文字
pscoast4 下移3.25i画中间图。陆地填充分辨率100dpi，采用彩色模板86，前景红色，背景黄色；水体
填充分辨率100dpi，采用模板7，前景红色，背景黑色。；
echo 0 15 ??分别书写“SILLY USES OF ”和|“GMT COLOR PATTERNS” 2行文字
grdimage3 再下移3.25i画下图。由lon.grd和lon.cpt绘制彩色影象，z值由-180度到180对应彩虹紫到红。
pscoast5 绘制海陆界线，-Gc 剪裁去陆地区，即只填充水体彩色影像。
grdimage4 由lat.grd和lat.cpt绘制陆地彩色影象，z值0、1对应颜色白、蓝。 pscoast6 -Q 恢复pscoast5剪裁去的陆地区
pscoast7 绘制海陆界线，填充陆地彩色影象。 echo 0 20?? 分别书写“FIRST INTERNATIONAL” 、“GMT CONFERENCE” 和“Honolulu, Hawaii, April 1, 2000”3行文字
40