# Save existing GMT defaults
echo Saving GMT defaults...
gmtdefaults -L > gmtdefaults$$
#
# Set new GMT defaults
echo Setting new GMT defaults...
gmtset PAPER_WIDTH 8.5
gmtset ANOT_FONT Helvetica
gmtset LABEL_FONT Helvetica
gmtset HEADER_FONT Helvetica
gmtset ANOT_FONT_SIZE 8
gmtset LABEL_FONT_SIZE 8
gmtset HEADER_FONT_SIZE 10
gmtset FRAME_WIDTH 0.075
gmtset TICK_LENGTH 0.075
gmtset PAGE_ORIENTATION LANDSCAPE
gmtset COLOR_BACKGROUND 0/0/0
gmtset COLOR_FOREGROUND 255/255/255
gmtset COLOR_NAN 255/255/255
#
# Make color pallette table file
echo Making color pallette table file...
echo  -37.5 255 255 255      0 255 186 133 > transit10m.grd.cpt
echo      0 255 186 133   37.5 255 161  68 >> transit10m.grd.cpt
echo   37.5 255 161  68     75 255 189  87 >> transit10m.grd.cpt
echo     75 255 189  87  112.5 240 236 121 >> transit10m.grd.cpt
echo  112.5 240 236 121    150 205 255 162 >> transit10m.grd.cpt
echo    150 205 255 162  187.5 138 236 174 >> transit10m.grd.cpt
echo  187.5 138 236 174    225 106 235 255 >> transit10m.grd.cpt
echo    225 106 235 255  262.5  50 190 255 >> transit10m.grd.cpt
echo  262.5  50 190 255    300  40 127 251 >> transit10m.grd.cpt
echo    300  40 127 251  337.5  37  57 175 >> transit10m.grd.cpt
#
# Make color image
echo Running grdimage...
grdimage transit10m.grd -Jm30.3030303030303 \
	-R-170.852/-170.555/-14.3785/-14.2261 -Ctransit10m.grd.cpt \
	-X1 -Y2.1320303030303 -K -V > transit10m.grd.ps
#
# Make contour plot
#echo Running grdcontour...
#grdcontour transit10m.grd -Jm30.3030303030303 \
#	-R-170.852/-170.555/-14.3785/-14.2261 \
#	-C25 \
#	-L-17.1338/302.693 -Wc1p \
#	-K -O -V >> transit10m.grd.ps
#
# Make color scale
echo Running psscale...
psscale -Ctransit10m.grd.cpt \
	-D4.5000/-0.5000/9.0000/0.1500h \
	-B":.Depth (m):" \
	-K -O -V >> transit10m.grd.ps
#
# Make basemap
echo Running psbasemap...
psbasemap -Jm30.3030303030303 \
	-R-170.852/-170.555/-14.3785/-14.2261 \
	-Bf2mg1m/f2mg1m:."Transit - 10-m Bathymetry Grid": \
	-O -V >> transit10m.grd.ps
#
# Delete surplus files
echo Deleting surplus files...
rm -f transit10m.grd.cpt
#
# Reset GMT default fonts
echo Resetting GMT fonts...
mv gmtdefaults$$ .gmtdefaults
#
# Run ghostview
echo Running ghostview in background...
#ghostview  transit10m.grd.ps &
#
# All done!
echo All done!