monochrome contour map - 等深線図をつくる　

---

  使うコマンド　　grdcontour coast   basemap



bash スクリプト例
#  plot bathymetry contours modern

# parameter setting
region=130/145/25/35                    # map region east/west/south/north
proj=M10                                # map projection and scale
ticks=a5f1g5                            # boundary tick info
frame=WSne+tContour                     # boundary frame info
cint=500                                # contour interval
aint=2000                               # annotation contour interval
limit=-10000/-500                       # contour min/max
grdfile=JTOPO1_30.grd         # input bathymetry grid fil                  
#
gmt begin contour        
       gmt grdcontour $grdfile -R$region -J$proj -C$cint -A$aint -L$limit
       gmt coast -Di -Ggray -Wthin,black
       gmt basemap -B$ticks -B$frame
gmt end show



Tips

• -Lオプションで等深線の範囲を明示的に決めている。普通は不要だが、エラー値のプロットや海岸線データと地形の０mコンターが厳密に合わずに交差してしまうことを避けるために指定。
• アノテーション数値の間隔はオプション指定できるが、重なりなどは制御できないので、最終的にはお絵かきソフトで動かしたり消したりするほうが楽。