#! /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 LineD.grd -JX17.85/5 -G1 -V -O LineDProfile -Z0/10000 -B500:Trace:/0.1:Time::,sec::.LineD:WeSn -MGDANOT_FONT_SIZE/16 -MGDLABEL_FONT_SIZE/16 -MGDHEADER_FONT_SIZE/20
#
# Define shell variables used in this script:
set PS_FILE         = LineDProfile.ps
set CPT_FILE        = LineDProfile.cpt
set MAP_PROJECTION  = X
set MAP_SCALE       = 17.85/5
set MAP_REGION      = 20599.5/25900.5/-0.5/3.2e-05
set X_OFFSET        = 2.075
set Y_OFFSET        = 6.5
#
# 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
#
# Set user defined GMT defaults
echo Setting user defined GMT defaults...
gmtset ANOT_FONT_SIZE 16
gmtset LABEL_FONT_SIZE 16
gmtset HEADER_FONT_SIZE 20
#
# Make color pallette table file
echo Making color pallette table file...
echo      0  37  57 175   1250  40 127 251 > $CPT_FILE
echo   1250  40 127 251   2500  50 190 255 >> $CPT_FILE
echo   2500  50 190 255   3750 106 235 255 >> $CPT_FILE
echo   3750 106 235 255   5000 138 236 174 >> $CPT_FILE
echo   5000 138 236 174   6250 205 255 162 >> $CPT_FILE
echo   6250 205 255 162   7500 240 236 121 >> $CPT_FILE
echo   7500 240 236 121   8750 255 189  87 >> $CPT_FILE
echo   8750 255 189  87  10000 255 161  68 >> $CPT_FILE
echo  10000 255 161  68  11250 255 186 133 >> $CPT_FILE
echo  11250 255 186 133  12500 255 255 255 >> $CPT_FILE
#
# Define data files to be plotted:
set DATA_FILE        = LineD.grd
set INTENSITY_FILE   = 
#
# Make color image
echo Running grdimage...
grdimage $DATA_FILE -J$MAP_PROJECTION$MAP_SCALE \
	-R$MAP_REGION -C$CPT_FILE \
	-X$X_OFFSET -Y$Y_OFFSET -K -V >! $PS_FILE
#
# Make color scale
# echo Running psscale...
# psscale -C$CPT_FILE \
# 	-D8.9250/-1.0000/17.8500/0.3000h \
# 	-B":Data Values:" \
# 	-K -O -V >> $PS_FILE
#
# Make basemap
echo Running psbasemap...
psbasemap -J$MAP_PROJECTION$MAP_SCALE \
	-R$MAP_REGION \
	-B500:Trace:/0.1:Time::,sec::.LineD:WeSn \
	-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 -landscape -media BBox LineDProfile.ps &
#
# All done!
echo All done!