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://wap.sciencenet.cn/blog-350278-495111.html
上一篇:
fortran lapack的dgesvd和matlab SVD结果不一致,请大神指点!下一篇:
终于用GMT搞定全球部分IGS台站分布图了