ps=igs_cn.ps

gmtset HEADER_FONT_SIZE 10 X_ORIGIN 4c Y_ORIGIN 8c
   ANNOT_FONT_SIZE_PRIMARY 6 FRAME_WIDTH 0.05c

pscoast -R70/140/15/55 -JM5i -Glightyellow -Slightblue -U"Made by Feng Zhou" -A1000
    -Ba10f5/5:."IGS Stations in China": -K -Dh -N1/1p -Wthinnest -P > $ps

# Plot red circle at cities and plot names:
psxy -R -J -O -K -Sc0.15 -Gred -Wthinnest igs_sta.cn >> $ps
awk '{print $1, $2, 5, 0, 5, $4, $3}' igs_sta.cn | pstext -R -J -O -K -Dj0.15/0 -Gred -N >> $ps

# Clean up after ourselves:

rm -f igs_sta.cn .gmt*

转载本文请联系原作者获取授权，同时请注明本文来自周锋科学网博客。
链接地址：https://m.sciencenet.cn/blog-350278-495111.html

上一篇：fortran lapack的dgesvd和matlab SVD结果不一致，请大神指点！
下一篇：终于用GMT搞定全球部分IGS台站分布图了