具體代碼如下：其中changjaing.txt為流域邊界數(shù)據(jù)，為經(jīng)度、緯度兩列數(shù)據(jù)，negative.txt為GPS站點數(shù)據(jù)，包含經(jīng)度、緯度、東西向分量（為0）,南北向分量

REM 2021.3.29/ Z.Q.WEN ECJTU

set R=70/140/15/55

REM set J=j0.1i

set J = M8c

gmt set FORMAT_GEO_MAP=ddd:mm:ssF

gmt set MAP_FRAME_TYPE=plain

gmt set MAP_TICK_LENGTH=10p

gmt set FONT_LABEL 7p,Times-Roman

gmt set FONT_TITLE 7p,Times-Roman

gmt set MAP_TITLE_OFFSET=0.4c

gmt set PS_MEDIA=a3

gmt set FONT_ANNOT_PRIMARY 8p,Times-Roman

gmtset MAP_FRAME_TYPE plain

set PS="indus_ganges_.ps"

REM gmt grdcut ETOPO1_Bed_g_gdal.grd -R%R% -GcutTopo.grd

REM gmt grdgradient cutTopo.grd -Ne0.7 -A50 -GcutTopo_i.grd

REM gmt makecpt -Cglobe -T-8000/5000? -Di >t.cpt

gmt makecpt -Cgray -T1000/20000/100 -I >t.cpt

gmt grdimage cutTopo.grd -IcutTopo_i.grd? -R%R% -J%J% -Xc -Yc? -Ct.cpt? -K >%PS%

gmt pscoast -R%R% -J%J%? -Slightblue -A100? -W0.1p -O -K >>%PS%

gmt xyz2grd? speed.txt -R%R% -I1? -Gtmp.grd

gmt grdsample tmp.grd -Gtmp.grd -I0.01

gmt makecpt -Crainbow -T-106/7.7/1 -D >tem.cpt

REM gmt psclip load.txt? -J%J% -R%R% -Am? -O -K >>%PS%

REM gmt grdimage tmp.grd -R%R% -J%J% -Ctem.cpt -B5f1/5f1WenS -t80 -O -K>>%PS%

gmt psxy gps.txt? -J%J% -R%R%? -Sc0.01c -Gblue? -O -K>> %ps%

REM gmt psclip -C -O -K? >>%PS%

gmt psxy? changjiang.txt -R -J? -O -K -W2p,gray44 >>? %ps%

echo 108? 30? 1|gmt pstext? -R%R% -J%J% -O -K -F+f10p,Times-Roman >>? %ps%

gmt psxy hexi.txt -R -J -O -K -W2p,gray44>> %ps%

echo 100? 40? 2|gmt pstext? -R%R% -J%J% -O -K -F+f10p,Times-Roman >> %ps%

gmt psxy huanghe.txt -R -J? ?-O -K -W2p,gray44>>? %ps%

echo 109? 37? 3|gmt pstext? -R%R% -J%J% -O -K -F+f10p,Times-Roman >>? %ps%

gmt psxy tp.txt -R -J? -O -K -W2p,gray44>>? %ps%

echo 86? 33? 4|gmt pstext? -R%R% -J%J% -O -K -F+f10p,Times-Roman >> %ps%

gmt psxy neimenggu.txt -R -J -O -K -W2p,gray44>>? %ps%

echo 113? 43? 5|gmt pstext? -R%R% -J%J% -O -K -F+f10p,Times-Roman >>? %ps%

gmt psxy southwest.txt -R -J? -O -K -W2p,gray44 >>? %ps%

echo 95? 30? 6|gmt pstext? -R%R% -J%J% -O -K -F+f10p,Times-Roman >>? %ps%

gmt psxy southeast.txt -R -J? -O -K -W2p,gray44 >>? %ps%

echo? 119? 27? 7|gmt pstext? -R%R% -J%J% -O -K -F+f10p,Times-Roman >>? %ps%

gmt psxy tarim.txt -R -J? ? -O -K -W2p,gray44 >>? %ps%

echo? 85? 40? 8|gmt pstext? -R%R% -J%J% -O -K -F+f10p,Times-Roman >>? %ps%

gmt psxy northeast.txt -R -J? -O -K? -W2p,gray44>>? %ps%

echo? 125? 47? 9|gmt pstext? -R%R% -J%J% -O? -K -F+f10p,Times-Roman >>? %ps%

gmt psxy zj.txt -R -J? -O -K? -W2p,gray44>>? %ps%

echo? 110? 24? 10|gmt pstext? -R%R% -J%J% -O -K -F+f10p,Times-Roman >>? %ps%

gmt psxy zhungaer.txt -R -J -O -K? -W2p,gray44 >> %ps%

echo? 88? 45? 11|gmt pstext? -R%R% -J%J% -O? -K -F+f10p,Times-Roman>>? %ps%

gmt psxy haihe.txt -R -J? -O -K-W2p,gray44>>? %ps%

echo 116? 38? 12|gmt pstext? -R%R% -J%J% -O -K -F+f10p,Times-Roman >>? %ps%

gmt psxy huaihe1.txt -R -J? -O -K -W2p,gray44>> %ps%

echo 118? 34? 13|gmt pstext? -R%R% -J%J% -O -K -F+f10p,Times-Roman >> %ps%

REM gawk "{print $1,$2,$3}" kk.txt |

echo 70 15 > temp

echo 80 15 >> temp

echo 80 25 >> temp

echo 70 25 >> temp

gmt psxy temp -J -R -Gwhite -W0.1p,white -A -O -K >>%PS%?

gmt psvelo positive.txt -R%R% -J%J% -Se0.2c/0.95/0 -A0.4c+e+p0.1p -Gred -W0.2p,red -O -K -N >>%PS%

gmt psvelo negative.txt -R%R% -J%J% -Se0.2c/0.95/0 -A0.4c+e+p0.1p -Gblue -W0.2p,blue -O -K -N >>%PS%

gmt psvelo nc.txt -R%R% -J%J% -Se0.05c/0.95/0 -A0.4c+e+p0.1p -Ggreen -W0.2p,green -O -K -N >>%PS%

echo 75 23 4 \261 1 mm/yr | gmt pstext -J -R -F+f7p,4,red -O -K >>%PS%

echo 75 18 -4 \261 1 mm/yr | gmt pstext -J -R -F+f7p,4,blue -O -K >>%PS%

echo 75 16 -16 \261 1 mm/yr | gmt pstext -J -R -F+f7p,4,green -O -K >>%PS%

REM echo 73 19 0 5 0 0 0 | gmt psvelo -J -R -Se0.2c/0.95/0 -A0c+e+p0.1p -W0.1p,Gray45 -O -K >>%PS%

echo 73 19 0 4 0 0 0 | gmt psvelo -J -R -Se0.2c/0.95/0 -A0.2c+e+p0.1p -Gred -W0.2p,red -O -K >>%PS%

REM echo 75 22 0 -5 0 0 0 | gmt psvelo -J -R -Se0.2c/0.95/0 -A0c+e+p0.1p -W0.1p,Gray45 -O -K >>%PS%

echo 75 22 0 -4 0 0 0 | gmt psvelo -J -R -Se0.2c/0.95/0 -A0.2c+e+p0.1p -Gblue -W0.2p,blue -O -K >>%PS%

REM echo 77 22 0 -5 0 0 0 | gmt psvelo -J -R -Se0.2c/0.95/0 -A0c+e+p0.1p -W0.1p,Gray45 -O -K >>%PS%

echo 77 22 0 -4 0 0 0 | gmt psvelo -J -R -Se0.2c/0.95/0 -A0.2c+e+p0.1p -Ggreen -W0.2p,blue -O -K >>%PS%

gmt psbasemap -R%R% -J%J% -Bxa10 -Bya5 -BWesN? -Xc -Yc -O >>%PS%

REM gmt psscale -Dx0.5c/-01.2c+jTC+w6c/0.4c+v+e --FONT_ANNOT_PRIMARY=14p,Helvetica,black -X1c -Y12c? -O? -Ctem.cpt -Ba25 -By+lmm/a>> %ps%

REM 繪制圖中圖

REM set Rg=-180/180/-90/90

REM set Jg=G78/32/5c

REM gmt grdimage ETOPO1_Bed_g_gdal.grd? -R%Rg% -J%Jg%? -Ct.cpt? ?-X11.5c? -Y-5.5c? -O -K >>%PS%

REM gmt pscoast -R%Rg% -J%Jg% -N1 -A1000 -B0? -Wthinnest? -O -K -X0c -Y0c? --MAP_FRAME_TYPE=plain>> %PS%

REM gmt psbasemap -R%Rg% -J%Jg% -D%R%? -Bf30a30 -Bg60? -Fd+p2p,red+c0.1? -O>>%PS%

gmt psconvert? %PS% -A -Tg -P