#! /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 -I ../ReferenceGrids/MontereyA_bath.grd -R-122.25/-121.75/36.5/37 -O FigLocMARSFall2005 -LMBARI Mapping AUV Sea Tests - Monterey Bay:Topography (m) -G5 -D0/1 -A2 -PC -T -MXM -MXW2p/0/0/255 -MXI../MARSsurvey/MARSroute.prn -MXW1p -MXIMAUVSpring2005Nav.xy -MXW1p -MXI../Figures2004Spring/MAUV_Spring2004_Nav.xy -MXW1p -MXIMAUVMARS2005Nav.xy -MXW2p/255/0/0 -MXIC100mBounds.xy -MXIC300mBounds.xy -MXIC520mBounds.xy -MXIC1000mBounds.xy -MXIC1400mBounds.xy -MGLf-121.85/36.55/36.75/10 -V
#
# Define shell variables used in this script:
set PS_FILE         = FigLocMARSFall2005.ps
set CPT_FILE        = FigLocMARSFall2005.cpt
set MAP_PROJECTION  = m
set MAP_SCALE       = 24.141
set MAP_REGION      = -122.25/-121.75/36.5/37
set X_OFFSET        = 2.4646
set Y_OFFSET        = 4
#
# 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 archC+
gmtset ANOT_FONT Helvetica
gmtset LABEL_FONT Helvetica
gmtset HEADER_FONT Helvetica
gmtset ANOT_FONT_SIZE 16
gmtset LABEL_FONT_SIZE 16
gmtset HEADER_FONT_SIZE 20
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 color pallette table file
echo Making color pallette table file...
echo  -2700  37  57 175  -2400  40 127 251 > $CPT_FILE
echo  -2400  40 127 251  -2100  50 190 255 >> $CPT_FILE
echo  -2100  50 190 255  -1800 106 235 255 >> $CPT_FILE
echo  -1800 106 235 255  -1500 138 236 174 >> $CPT_FILE
echo  -1500 138 236 174  -1200 205 255 162 >> $CPT_FILE
echo  -1200 205 255 162   -900 240 236 121 >> $CPT_FILE
echo   -900 240 236 121   -600 255 189  87 >> $CPT_FILE
echo   -600 255 189  87   -300 255 161  68 >> $CPT_FILE
echo   -300 255 161  68      0 255 186 133 >> $CPT_FILE
echo      0 255 186 133    300 255 255 255 >> $CPT_FILE
#
# Define data files to be plotted:
set DATA_FILE        = ../ReferenceGrids/MontereyA_bath.grd
set INTENSITY_FILE   = 
#
# Get slope array
echo Getting slope array...
echo Running grdgradient to get x component of the gradient...
grdgradient $DATA_FILE -A90 -G$DATA_FILE.drvx -M
echo Running grdgradient to get y component of the gradient...
grdgradient $DATA_FILE -A0 -G$DATA_FILE.drvy -M
echo Running grdmath to get slope magnitude...
grdmath $DATA_FILE.drvx 2.0 POW \
	$DATA_FILE.drvy 2.0 POW ADD SQRT \
	-2 MUL \
	= $DATA_FILE.slope
/bin/rm -f $DATA_FILE.drvx $DATA_FILE.drvy 
#
# Make color image
echo Running grdimage...
grdimage $DATA_FILE -J$MAP_PROJECTION$MAP_SCALE \
	-R$MAP_REGION -C$CPT_FILE \
	-I$DATA_FILE.slope \
	-P -X$X_OFFSET -Y$Y_OFFSET -K -V >! $PS_FILE
#
# Make coastline data plot
echo Running pscoast...
pscoast \
	-J$MAP_PROJECTION$MAP_SCALE \
	-R$MAP_REGION \
	-Df \
	-G200 \
	-W1p \
	-P -K -O -V >> $PS_FILE
#
# Make xy data plot
echo Running psxy...
psxy ../MARSsurvey/MARSroute.prn \
	-J$MAP_PROJECTION$MAP_SCALE \
	-R$MAP_REGION \
	-M'>' \
	-W2p/0/0/255 \
	-P -K -O -V >> $PS_FILE
#
# Make xy data plot
echo Running psxy...
psxy MAUVSpring2005Nav.xy \
	-J$MAP_PROJECTION$MAP_SCALE \
	-R$MAP_REGION \
	-M'>' \
	-W1p \
	-P -K -O -V >> $PS_FILE
#
# Make xy data plot
echo Running psxy...
psxy ../Figures2004Spring/MAUV_Spring2004_Nav.xy \
	-J$MAP_PROJECTION$MAP_SCALE \
	-R$MAP_REGION \
	-M'>' \
	-W1p \
	-P -K -O -V >> $PS_FILE
#
# Make xy data plot
echo Running psxy...
psxy MAUVMARS2005Nav.xy \
	-J$MAP_PROJECTION$MAP_SCALE \
	-R$MAP_REGION \
	-M'>' \
	-W1p \
	-P -K -O -V >> $PS_FILE
#
# Make xy data plot
echo Running psxy...
psxy C100mBounds.xy \
	-J$MAP_PROJECTION$MAP_SCALE \
	-R$MAP_REGION \
	-M'>' \
	-W2p/255/0/0 \
	-P -K -O -V >> $PS_FILE
#
# Make xy data plot
echo Running psxy...
psxy C300mBounds.xy \
	-J$MAP_PROJECTION$MAP_SCALE \
	-R$MAP_REGION \
	-M'>' \
	-W2p/255/0/0 \
	-P -K -O -V >> $PS_FILE
#
# Make xy data plot
echo Running psxy...
psxy C520mBounds.xy \
	-J$MAP_PROJECTION$MAP_SCALE \
	-R$MAP_REGION \
	-M'>' \
	-W2p/255/0/0 \
	-P -K -O -V >> $PS_FILE
#
# Make xy data plot
echo Running psxy...
psxy C1000mBounds.xy \
	-J$MAP_PROJECTION$MAP_SCALE \
	-R$MAP_REGION \
	-M'>' \
	-W2p/255/0/0 \
	-P -K -O -V >> $PS_FILE
#
# Make xy data plot
echo Running psxy...
psxy C1400mBounds.xy \
	-J$MAP_PROJECTION$MAP_SCALE \
	-R$MAP_REGION \
	-M'>' \
	-W2p/255/0/0 \
	-P -K -O -V >> $PS_FILE
#
# Make color scale
echo Running psscale...
psscale -C$CPT_FILE \
	-D6.0354/-1.0000/12.0707/0.3000h \
	-B":Topography (m):" \
	-P -K -O -V >> $PS_FILE
#
# Make basemap
echo Running psbasemap...
psbasemap -J$MAP_PROJECTION$MAP_SCALE \
	-R$MAP_REGION \
	-B10m/10m:."MBARI Mapping AUV Sea Tests - Monterey Bay": \
	-Lf-121.85/36.55/36.75/10 \
	-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 gv
echo Running gv in background...
gv --orientation=portrait --media=BBox FigLocMARSFall2005.ps &
#
# All done!
echo All done!