GMT范例注释07

REM GMT EXAMPLE 07

pscoast -R-50/0/-10/20 -JM9i -K -GP300/26 -Dl -W0.25p -B10 -U"Example 7 in Cookbook" > e_07.ps

-GP300/26：陆地填充用模板26，反相，深灰（300）底色。

psxy -R -JM -O -K -M fz.xy -W0.5pta >> e_07.ps

文件fz.xy内容如下：

> 629 OCEAN: 2 FEATURE: 6 NAME: KANE

284.74011 35.04543

284.82080 34.99857

284.90210 34.94939

…………………

> 445 OCEAN: 2 FEATURE: 6 NAME: FIFTEEN TWENTY

298.96030 18.06335

299.04568 18.01421

299.13681 17.96838

…………………

文件fz.xy为区内大洋转换断层的数据，为多段数据格式（-M）。每一段数据由一行“>”开头的文字和一系列x、y坐标组成，用来绘制一条断层，文字行包含了断层名字和x、y坐标组数等信息。

gawk "{print $1-360.0, $2, $3*0.01}" quakes.xym | psxy -R -JM -O -K -H1 -Sci -G255 -W0.25p >> e_07.ps

文件quakes.xym内容如下

lon lat mag1

338.799999 -8.4 4.8

346.68 -8.36 4.4

346.47 -8.3 4.5

346.57 -8.23 4.5

…………………

gawk命令由文件quakes.xym中读取大洋中脊地震的坐标和震级，绘制地震分布图。

psxy -R -JM -O -K -M isochron.xy -W0.75p >> e_07.ps

psxy -R -JM -O -K -M ridge.xy -W1.75p >> e_07.ps

文件isochron.xy和ridge.xy分别为区内“等时”断层和大洋中脊断层的数据，数据格式与fz.xy相同。

echo -14.5 15.2 > tmp

echo -2 15.2 >> tmp

echo -2 17.8 >> tmp

echo -14.5 17.8 >> tmp

psxy -R -JM -O -K -G255 -W1p -A tmp >> e_07.ps

4行echo语句给文件tmp赋值，作为矩形内框的4角坐标。

echo -14.35 15.35 > tmp

echo -2.15 15.35 >> tmp

echo -2.15 17.65 >> tmp

echo -14.35 17.65 >> tmp

psxy -R -JM -O -K -G255 -W0.5p -A tmp >> e_07.ps

4行echo语句给文件tmp重新赋值，作为矩形外框的4角坐标。

echo -13.5 16.5 | psxy -R -JM -O -K -Sc0.08i -G255 -W0.5p >> e_07.ps

在矩形框内绘制地震图例（圆圈），echo给出圆心位置。

echo -12.5 16.5 18 0 6 5 ISC Earthquakes | pstext -R -JM -O -K >> e_07.ps

在图例（圆圈）右面写“ISC Earthquakes”，echo后面参数依次为x坐标、y坐标、字符大小（字号）、角度（0）、字体（6）和对齐方式，5表示坐标值为字串的中心。

echo -43 -5 30 0 1 6 SOUTH > tmp

echo -43 -8 30 0 1 6 AMERICA >> tmp

echo -7 11 30 0 1 6 AFRICA >> tmp

pstext -R -JM -O -S0.75p -G255 tmp >> e_07.ps

3行echo语句给文件tmp重新赋值，分别写3个字串，对齐方式6表示坐标值为字串右边中点。

del .gmt*