#! /bin/csh -f
#
# Shellscript to create Postscript plot of swath sonar data
# Created by macro mbm_plot
#
# This shellscript created by following command line:
# mbm_plot -G5 -S -I 20040622_193147.s7k -V
#
# Define shell variables used in this script:
set PS_FILE         = 20040622_193147.s7k.ps
set CPT_FILE        = 20040622_193147.s7k.cpt
set MAP_PROJECTION  = m
set MAP_SCALE       = 282.62
set MAP_REGION      = -121.85519/-121.84452/36.79413/36.81547
set X_OFFSET        = 2.7422
set Y_OFFSET        = 2
set INPUT_FILE      = 20040622_193147.s7k
set INPUT_FORMAT    = 88
#
# 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_MEDIA archA+
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
gmtset DEGREE_FORMAT 3
#
# Make datalist file 
echo Making datalist file...
echo $INPUT_FILE $INPUT_FORMAT >! datalist$$
#
# Make color pallette table file
echo Making color pallette table file...
echo -15.8834 255 255 255 -7.63333 230 230 230 >! $CPT_FILE
echo -7.63333 230 230 230 0.61677 204 204 204 >> $CPT_FILE
echo 0.61677 204 204 204 8.86687 179 179 179 >> $CPT_FILE
echo 8.86687 179 179 179 18.4047 153 153 153 >> $CPT_FILE
echo 18.4047 153 153 153 36.5511 128 128 128 >> $CPT_FILE
echo 36.5511 128 128 128 57.8736 102 102 102 >> $CPT_FILE
echo 57.8736 102 102 102 88.3737  77  77  77 >> $CPT_FILE
echo 88.3737  77  77  77 157.733  51  51  51 >> $CPT_FILE
echo 157.733  51  51  51 428.421  26  26  26 >> $CPT_FILE
echo 428.421  26  26  26  31911   0   0   0 >> $CPT_FILE
#
# Run mbswath
echo Running mbswath...
mbswath -f-1 -Idatalist$$ \
	-J$MAP_PROJECTION$MAP_SCALE \
	-R$MAP_REGION \
	-C$CPT_FILE \
	-Z5 \
	-A1 \
	-p1 \
	-P -X$X_OFFSET -Y$Y_OFFSET -K -V >! $PS_FILE
#
# Make color scale
echo Running psscale...
psscale -C$CPT_FILE \
	-D1.5078/-0.5000/3.0155/0.1500h \
	-B":Sidescan Pixel Values:" \
	-L \
	-P -K -O -V >> $PS_FILE
#
# Make basemap
echo Running psbasemap...
psbasemap -J$MAP_PROJECTION$MAP_SCALE \
	-R$MAP_REGION \
	-B10c/10c:."Data File $INPUT_FILE": \
	-P -O -V >> $PS_FILE
#
# Delete surplus files
echo Deleting surplus files...
/bin/rm -f $CPT_FILE
/bin/rm -f datalist$$
#
# Reset GMT default fonts
echo Resetting GMT fonts...
/bin/mv gmtdefaults$$ .gmtdefaults
#
# Run gv
echo Running gv in background...
gv -portrait -media BBox $PS_FILE &
#
# All done!
echo All done!