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

  使うコマンド  grdcontour pscoast   psbasemap



bash スクリプト例

#  plot nathymetry contours

# 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
psfile=contour.ps                  # output postscript file name
#
gmt grdcontour $grdfile -R$region -J$proj -C$cint -A$aint -L$limit -K -V > $psfile
gmt pscoast -R$region -J$proj -Di -Ggray -Wthin,black -K -V -O >> $psfile
gmt psbasemap -R$region -J$proj -B$ticks -B$frame -O  -V >> $psfile
#


Tips