德克萨斯大学埃尔帕索地质研究所GMT实例之二

huolongguo

bat代码：

echo -108 53.5 15 0 1 2 Bouguer Gravity for Deep Probe 1995 | pstext -R-119/-100/30/54 -Jt-110/1:11500000 -X1.0 -K -P> dpgravmap.ps

rem Generate color image
grdimage dpgrav.grd -Cdpgrav.cpt4 -R-118/-102/31/53 -Jt-110/1:11500000 -Idpgravinten.grd -X1.5 -O -K >> dpgravmap.ps
rem Put in contours
grdcontour dpgrav.grd -C50 -A50f10 -Jt-110/1:11500000 -R-118/-102/31/53 -S25 -L-350/50 -Wc2/0/0/0 -K -O >> dpgravmap.ps
rem Put coastlines, national boundaries, rivers and scales and grid
pscoast -R -Jt -Lf-116/50/35/200 -Dh -B5 -Na/5/0/205/96 -S168/255/255 -W3 -K -O >> dpgravmap.ps
rem #Bring in coord files
psxy -: dpgrav.prof2 -R -Jt -P -Sd0.06/255/0/0 -G200/0/0 -O -K >> dpgravmap.ps
psxy -: 106grav.prof -R -Jt -P -Sa0.05/255/0/0 -G0/0/0 -O -K>> dpgravmap.ps
rem Put state names in white
echo -108 50.5 15 0 1 2 Canada >tmp
echo -110 47 15 0 1 2 Montana >>tmp
echo -114 44 15 0 1 2 Idaho >>tmp
echo -106 43 15 0 1 2 Wyoming >>tmp
echo -112 40 15 0 1 2 Utah >>tmp
echo -106 39 15 0 1 2 Colorado>>tmp
echo -111 36 15 0 1 2 Arizona>>tmp
echo -106 35 15 0 1 2 New >>tmp
echo -106 34.5 15 0 1 2 Mexico >>tmp
pstext tmp -R -Jt -G255/255/255 -K -O >> dpgravmap.ps
del tmp

rem #Plot scale
psscale -Cdpgrav.cpt4 -D14.5/3.5/6/0.3 -E -Ba100g50/:"(mgal)":  -I -O >> dpgravmap.ps
ps2raster dpgravmap.ps -A -P -Tg
del .gmt*