GMT范例三

葡萄

gmtset PLOT_DEGREE_FORMAT ddd:mm:ssF
set filename=gmt2011_01c.ps
grd2cpt tw.grd -Crelief -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%
psscale -D7.5/-1.5/15/0.4h -Ctmp.cpt -Bf200a1000:Topo_Bathymetry_relief:/:meter: -I -O -V >> %filename%
del .gmt*
del *.int
del *.cpt