GMT绘制全球地震观测网络站点

huolongguo

bat代码：

pscoast -R-180/180/-90/90 -JR20 -Bg30 -Dc -A10000 -Gburlywood4 -Swheat1 -P -K -V> GMT_winkel.ps> hw02.ps
gawk "{print $3,$2}" gsn_stations.txt| psxy -R -J -Ss0.1 -Gred -O -V -K>>hw02.ps
gawk "{print $3-4,$2+1,6,0,1,1,$1}" gsn_stations.txt | pstext -R -J -O -V >>hw02.ps
del .gmt*

效果图：