GMT范例四-等值线

葡萄

gmtset PLOT_DEGREE_FORMAT ddd:mm:ssF
set filename=gmt2011_02d_contour.ps
grd2cpt tw.grd -Crainbow -Z > tmp.cpt
grdgradient tw.grd -A45 -Gtw.grd.int -Ne1 -M
grdimage tw.grd -Jm1.2i -R118.5/123.5/20/26.5 -Ctmp.cpt -Itw.grd.int -Y6c -P -K -V > %filename%
pscoast -R -Jm -Df -Ba1f0.5 -Ia/1p/100 -Lf122/26/24/200 -W3 -P -O -K -V >> %filename%
grdcontour tw.grd -R -Jm -L-8000/-500 -C500 -A1000 -Wa5/25 -Wc2/100 -S20 -V -O -P -K >> %filename%
grdcontour tw.grd -R -Jm -L500/4000 -C500 -Wc2/50 -S20 -V -O -P -K >> %filename%
psscale -D7.5/-1.5/17/0.4h -Ctmp.cpt -B1000:Bathymetry:/:meter: -I -E -S -O -V >> %filename%
del .gmt*
del *.int
del *.cpt