GMT example - Gravity Map

pstext -R-119/-100/30/54 -Jt-110/1:11500000 -X1.0 -K <<END> dpgravmap.ps
-108 53.5 15 0 1 2 Bouguer Gravity for Deep Probe 1995
END
#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
#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
#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
#Bring in coord files
psxy -: /users/snelson/deep_probe/gravity/dpgrav.prof2 -R -Jt -P \
\-Sd0.06/255/0/0 -G200/0/0 -O -K >> dpgravmap.ps
psxy -: /users/snelson/deep_probe/gravity/106grav.prof -R -Jt -P \
\-Sa0.05/255/0/0 -G0/0/0 -O -K>> dpgravmap.ps
#Put state names in white
pstext -R -Jt -G255/255/255 -K -O <<END>> dpgravmap.ps
-108 50.5 15 0 1 2 Canada
-110 47 15 0 1 2 Montana
-114 44 15 0 1 2 Idaho
-106 43 15 0 1 2 Wyoming
-112 40 15 0 1 2 Utah
-106 39 15 0 1 2 Colorado
-111 36 15 0 1 2 Arizona
-106 35 15 0 1 2 New 
-106 34.5 15 0 1 2 Mexico
END
#Begin plotting legend
psxy -Jx1.0 -R0/3/0/7.5 -G255/255/255 -L -W10 -X-1 -Y0.1 -K -O <<END>> dpgravmap.ps
0 0
0 2.50
1.40 2.50
1.40 0
0 0
END
psxy -Jx -R -L -W3 -K -O <<END>> dpgravmap.ps
0.05 0.05
0.05 2.45
1.36 2.45
1.36 0.05
0.05 0.05
END
pstext -Jx -R -K -O <<END>> dpgravmap.ps
.7 .4 8 0 1 2 Scale 1:11,500,000
.6 .6 8 0 1 2 Contour Interval: 
1.0 .5 8 0 1 2 50 mGals
.7 2.3 10 0 1 2 Bouguer Gravity
.8 .2 8 0 7 2 Deep Probe Gravity
.8 .1 8 0 7 2 106 W Gravity
END
psxy -Jx -R -Sd0.06/255/0/0 -G200/0/0 -K -O <<END>> dpgravmap.ps
.2 .25
END
psxy -Jx -R -Sa0.05/0/0/0 -G0/0/0 -K -O <<END>> dpgravmap.ps
.2 .15
END
#Plot scale
psscale -Cdpgrav.cpt4 -D0.4/1.5/1.4/0.2 -I -O >> dpgravmap.ps
# Start ghostview
ghostview dpgravmap.ps
#
rm -i  dpgravmap.ps