#! /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 C1400m200411Topo.grd -O RM1400FigTopoNav -MGQ100 -Jm900 -G2 -D0/1 -A0.5/320 -V -MXW0p -MXIauvnavp.xy -MXW1p255/0/0 -MXI../WebArticle200504/FigLocSBPNav.xy -LMBARI Mapping AUV RM1400 Survey:Topography (m) -MGLfx3.0/0.5/36.83/0.2
#
# Define shell variables used in this script:
set PS_FILE         = RM1400FigTopoNav.ps
set CPT_FILE        = RM1400FigTopoNav.cpt
set MAP_PROJECTION  = m
set MAP_SCALE       = 900
set MAP_REGION      = 237.98/237.989/36.7225/36.7291
set X_OFFSET        = 6.95
set Y_OFFSET        = 5.3105
#
# 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 -1417.5  37  57 175  -1400  40 127 251 > $CPT_FILE
echo  -1400  40 127 251 -1382.5  50 190 255 >> $CPT_FILE
echo -1382.5  50 190 255  -1365 106 235 255 >> $CPT_FILE
echo  -1365 106 235 255 -1347.5 138 236 174 >> $CPT_FILE
echo -1347.5 138 236 174  -1330 205 255 162 >> $CPT_FILE
echo  -1330 205 255 162 -1312.5 240 236 121 >> $CPT_FILE
echo -1312.5 240 236 121  -1295 255 189  87 >> $CPT_FILE
echo  -1295 255 189  87 -1277.5 255 161  68 >> $CPT_FILE
echo -1277.5 255 161  68  -1260 255 186 133 >> $CPT_FILE
echo  -1260 255 186 133 -1242.5 255 255 255 >> $CPT_FILE
#
# Define data files to be plotted:
set DATA_FILE        = C1400m200411Topo.grd
set INTENSITY_FILE   = 
#
# Get shading array
echo Getting shading 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 adjusted x gradient...
grdmath $DATA_FILE.drvx -0.5 MUL 10 MUL = $DATA_FILE.magx
echo Running grdmath to get adjusted y gradient...
grdmath $DATA_FILE.drvy -0.5 MUL 10 MUL = $DATA_FILE.magy
echo Running grdmath to get normalization factor...
grdmath $DATA_FILE.magx 2.0 POW \
	$DATA_FILE.magy 2.0 POW ADD \
	1.0 ADD SQRT = $DATA_FILE.denom
echo Running grdmath to get normalized x gradient...
grdmath $DATA_FILE.magx $DATA_FILE.denom DIV = $DATA_FILE.normx
echo Running grdmath to get normalized y gradient...
grdmath $DATA_FILE.magy $DATA_FILE.denom DIV = $DATA_FILE.normy
echo Running grdmath to get normalized z gradient...
grdmath 1.0 $DATA_FILE.denom DIV = $DATA_FILE.normz
echo Running grdmath to apply lighting vector to normalized gradient...
grdmath $DATA_FILE.normx -0.556670465300929 MUL \
	$DATA_FILE.normy 0.663413898555546 MUL ADD \
	$DATA_FILE.normz 0.49999999226498 MUL ADD -0.5 ADD = $DATA_FILE.int
/bin/rm -f $DATA_FILE.drvx $DATA_FILE.drvy \
	$DATA_FILE.magx $DATA_FILE.magy \
	$DATA_FILE.denom $DATA_FILE.normx \
	$DATA_FILE.normy $DATA_FILE.normz
#
# Make color image
echo Running grdimage...
grdimage $DATA_FILE -J$MAP_PROJECTION$MAP_SCALE \
	-R$MAP_REGION -C$CPT_FILE \
	-I$DATA_FILE.int \
	-E100 \
	-X$X_OFFSET -Y$Y_OFFSET -K -V >! $PS_FILE
#
# Make xy data plot
echo Running psxy...
psxy auvnavp.xy \
	-J$MAP_PROJECTION$MAP_SCALE \
	-R$MAP_REGION \
	-W0p \
	-K -O -V >> $PS_FILE
#
# Make xy data plot
echo Running psxy...
psxy ../WebArticle200504/FigLocSBPNav.xy \
	-J$MAP_PROJECTION$MAP_SCALE \
	-R$MAP_REGION \
	-W1p255/0/0 \
	-K -O -V >> $PS_FILE
#
# Make color scale
echo Running psscale...
psscale -C$CPT_FILE \
	-D4.0500/-1.0000/8.1000/0.3000h \
	-B":Topography (m):" \
	-K -O -V >> $PS_FILE
#
# Make basemap
echo Running psbasemap...
psbasemap -J$MAP_PROJECTION$MAP_SCALE \
	-R$MAP_REGION \
	-B5c/5c:."MBARI Mapping AUV RM1400 Survey": \
	-Lfx3.0/0.5/36.83/0.2 \
	-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=landscape --media=BBox RM1400FigTopoNav.ps &
#
# All done!
echo All done!