1、首先，撰写脚本，并将脚本命名为topo10.gs。备注：脚本中用到的topo10.ctl，topo10.bin文件，请到“忻稿痞烁气象家园”论坛搜索“绘制地形图并标出研究区域”，然后下载这些文件。脚本如下：'reinit'set vpage 0.0 11.0 0.0 8.5'set parea 1.0 10.5 1.5 8.0'set grads off'set grid off''open topo10.ctl'set dfile 1'set x 600 781''set y 630 781''set vrange 15 40''set vrange2 100 130''set ylint 5''set xlint 5''set ylopts 1 4 0.18''set xlopts 1 4 0.18''set gxout shaded''set rgb 20 24 116 205''set rgb 21 51 204 255''set rgb 22 0 255 255''set rgb 23 0 255 100''set rgb 24 0 255 0''set rgb 25 90 149 44''set rgb 26 155 207 52''set rgb 27 173 223 80''set rgb 28 197 239 116''set rgb 29 229 243 167''set rgb 30 249 242 191''set rgb 31 243 234 169''set rgb 32 231 216 133''set rgb 33 218 197 93''set rgb 34 210 174 56''set rgb 35 203 153 28''set rgb 36 167 116 37''set rgb 37 131 84 40''set rgb 38 100 64 30''set rgb 39 82 54 25''set rgb 40 122 100 79''set clevs 0 100 200 300 400 500 600 700 800 900 1000 1100 1200 1300 1400 1500''set ccols 20 21 22 23 24 25 26 27 28 29 30 33 35 36 37 38 39 40''d tp''cbarn 1 0 5.7 0.35'*添加经纬度标注'set string 1 c 5 90''set strsiz 0.24''draw string 1.5 4.9 `1Latitude''set string 1 c 5 0''set strsiz 0.24''draw string 6.0 0.9 `1Longitude''printim topo10.jpg x1000 y800 white''print -R topo10.eps''disable print';
2、其次，运行上述脚本。即runtopo10.gs。
3、第三，得到如下图形。
4、第四，用红框标出三个研究区域。三个研究区域分别为（108-119E 19-25N），（116-121E 26-29N）和（112-122E 29-33N）。只需要添加如下代码：'q w2xy 108 19'x1=subwrd(result,3)y1=subwrd(result,6)'q w2xy 119 25'x2=subwrd(result,3)y2=subwrd(result,6)'set line 2 1 6''draw rec 'x1' 'y1' 'x2' 'y2'''q w2xy 116 26'x3=subwrd(result,3)y3=subwrd(result,6)'q w2xy 121 29'x4=subwrd(result,3)y4=subwrd(result,6)'set line 2 1 6''draw rec 'x3' 'y3' 'x4' 'y4'''q w2xy 112 29'x5=subwrd(result,3)y5=subwrd(result,6)'q w2xy 122 33'x6=subwrd(result,3)y6=subwrd(result,6)'set line 2 1 6''draw rec 'x5' 'y5' 'x6' 'y6''
5、最后，给出完整的脚本和最终的图像。争犸禀淫'reinit''set vpage 0.0 11.0 0.0 8.5''set parea 1.0 10.5 1.5 8.0''set grads off''set grid off''open topo10.ctl''set dfile 1''set x 600 781''set y 630 781''set vrange 15 40''set vrange2 100 130''set ylint 5''set xlint 5''set ylopts 1 4 0.18''set xlopts 1 4 0.18''set gxout shaded''set rgb 20 24 116 205''set rgb 21 51 204 255''set rgb 22 0 255 255''set rgb 23 0 255 100''set rgb 24 0 255 0''set rgb 25 90 149 44''set rgb 26 155 207 52''set rgb 27 173 223 80''set rgb 28 197 239 116''set rgb 29 229 243 167''set rgb 30 249 242 191''set rgb 31 243 234 169''set rgb 32 231 216 133''set rgb 33 218 197 93''set rgb 34 210 174 56''set rgb 35 203 153 28''set rgb 36 167 116 37''set rgb 37 131 84 40''set rgb 38 100 64 30''set rgb 39 82 54 25''set rgb 40 122 100 79''set clevs 0 100 200 300 400 500 600 700 800 900 1000 1100 1200 1300 1400 1500''set ccols 20 21 22 23 24 25 26 27 28 29 30 33 35 36 37 38 39 40''d tp''cbarn 1 0 5.7 0.35'*添加经纬度标注'set string 1 c 5 90''set strsiz 0.24''draw string 1.5 4.9 `1Latitude''set string 1 c 5 0''set strsiz 0.24''draw string 6.0 0.9 `1Longitude'*用红框标出三个研究区域'q w2xy 108 19'x1=subwrd(result,3)y1=subwrd(result,6)'q w2xy 119 25'x2=subwrd(result,3)y2=subwrd(result,6)'set line 2 1 6''draw rec 'x1' 'y1' 'x2' 'y2'''q w2xy 116 26'x3=subwrd(result,3)y3=subwrd(result,6)'q w2xy 121 29'x4=subwrd(result,3)y4=subwrd(result,6)'set line 2 1 6''draw rec 'x3' 'y3' 'x4' 'y4'''q w2xy 112 29'x5=subwrd(result,3)y5=subwrd(result,6)'q w2xy 122 33'x6=subwrd(result,3)y6=subwrd(result,6)'set line 2 1 6''draw rec 'x5' 'y5' 'x6' 'y6'''printim topo10.jpg x1000 y800 white''print -R topo10.eps''disable print';
