#! /bin/csh -f
#
# Shellscript to create Postscript plot of data in grd file
# Created by macro mbm_grdplot
#
# This shellscript created by following command line:
# mbm_grdplot -I206c.grd -G1 -C -D -V -LFile 206c.grd - Bathymetry Grid:Depth (m)
#
# Define shell variables used in this script:
set PS_FILE         = 206c.grd.ps
set CPT_FILE        = 206c.grd.cpt
set MAP_PROJECTION  = m
set MAP_SCALE       = 3.9456
set MAP_REGION      = -68.2221/-67.5285/-61.8031/-60.8903
set X_OFFSET        = 2.8817
set Y_OFFSET        = 2
#
# Save existing GMT defaults
echo Saving GMT defaults...
gmtdefaults -L > gmtdefaults$$
#
# Set new GMT defaults
echo Setting new GMT defaults...
gmtset MEASURE_UNIT inch
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   3375 255 255 255   3500 255 186 133 > $CPT_FILE
echo   3500 255 186 133   3625 255 161  68 >> $CPT_FILE
echo   3625 255 161  68   3750 255 189  87 >> $CPT_FILE
echo   3750 255 189  87   3875 240 236 121 >> $CPT_FILE
echo   3875 240 236 121   4000 205 255 162 >> $CPT_FILE
echo   4000 205 255 162   4125 138 236 174 >> $CPT_FILE
echo   4125 138 236 174   4250 106 235 255 >> $CPT_FILE
echo   4250 106 235 255   4375  50 190 255 >> $CPT_FILE
echo   4375  50 190 255   4500  40 127 251 >> $CPT_FILE
echo   4500  40 127 251   4625  37  57 175 >> $CPT_FILE
#
# Define data files to be plotted:
set DATA_FILE        = 206c.grd
set INTENSITY_FILE   = 
#
# Make color image
echo Running grdimage...
grdimage $DATA_FILE -J$MAP_PROJECTION$MAP_SCALE \
	-R$MAP_REGION -C$CPT_FILE \
	-P -X$X_OFFSET -Y$Y_OFFSET -K -V >! $PS_FILE
#
# Make contour plot
echo Running grdcontour...
grdcontour $DATA_FILE -J$MAP_PROJECTION$MAP_SCALE \
	-R$MAP_REGION \
	-C50 \
	-L3459.42/4477.87 -Wc1p \
	-P -K -O -V >> $PS_FILE
#
# Make color scale
echo Running psscale...
psscale -C$CPT_FILE \
	-D1.3683/-0.5000/2.7367/0.1500h \
	-B":.Depth (m):" \
	-P -K -O -V >> $PS_FILE
#
# Make basemap
echo Running psbasemap...
psbasemap -J$MAP_PROJECTION$MAP_SCALE \
	-R$MAP_REGION \
	-B10m/10m:."File 206c.grd - Bathymetry Grid": \
	-P -O -V >> $PS_FILE
#
# Delete surplus files
echo Deleting surplus files...
/bin/rm -f $CPT_FILE
#
# Reset GMT default fonts
echo Resetting GMT fonts...
/bin/mv gmtdefaults$$ .gmtdefaults
#
# Run /usr/local/packages/ImageMagick/display
echo Running /usr/local/packages/ImageMagick/display in background...
/usr/local/packages/ImageMagick/display  206c.grd.ps &
#
# All done!
echo All done!