#! /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 -f-41 -i 206a.list -L1 -C10/100/500/20/0.005/.011 -N -O 206a.list -V -Q -MGDBASEMAP_TYPE/PLAIN -MGDDEGREE_FORMAT/3 -MGDPAPER_WIDTH/36 -MGDANOT_FONT/Times-Roman -MGDANOT_FONT_SIZE/8 -MGDHEADER_FONT/Times-Roman -MGDHEADER_FONT_SIZE/10 -MGDLABEL_FONT/Times-Roman -MGDLABEL_FONT_SIZE/8
#
# Define shell variables used in this script:
set PS_FILE         = 206a.list.ps
set CPT_FILE        = 206a.list.cpt
set MAP_PROJECTION  = m
set MAP_SCALE       = 4.0183
set MAP_REGION      = -67.33398/-66.69862/-60.201635/-59.259265
set X_OFFSET        = 2.9735
set Y_OFFSET        = 2
set INPUT_FILE      = 206a.list
set INPUT_FORMAT    = -41
#
# 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
#
# Set user defined GMT defaults
echo Setting user defined GMT defaults...
gmtset BASEMAP_TYPE PLAIN
gmtset DEGREE_FORMAT 3
gmtset PAPER_WIDTH 36
gmtset ANOT_FONT Times-Roman
gmtset ANOT_FONT_SIZE 8
gmtset HEADER_FONT Times-Roman
gmtset HEADER_FONT_SIZE 10
gmtset LABEL_FONT Times-Roman
gmtset LABEL_FONT_SIZE 8
#
# Make datalist file 
echo Making datalist file...
echo $INPUT_FILE $INPUT_FORMAT >! datalist$$
#
# Run mbcontour
echo Running mbcontour...
mbcontour -f-1 -Idatalist$$ \
	-J$MAP_PROJECTION$MAP_SCALE \
	-R$MAP_REGION \
	-A10/100/500/20/0.005/0.1 \
	-D0.25/1/4/0.15 \
	-p1 \
	-P -X$X_OFFSET -Y$Y_OFFSET -K -V >! $PS_FILE
#
# Make basemap
echo Running psbasemap...
psbasemap -J$MAP_PROJECTION$MAP_SCALE \
	-R$MAP_REGION \
	-B10m/10m:."1": \
	-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 xpsview
#echo Running xpsview in background...
#xpsview -ps a -maxp 4m $PS_FILE &
#
# All done!
echo All done!