#! /bin/csh -f
#
# Shellscript to create TIFF image of data in grd file
# Created by macro mbm_grdtiff
#
# This shellscript created by following command line:
# mbm_grdtiff -I MARSEM300SSLower.grd -G1 -W1/4 -D0/1 -OSSGeoTiff_geo -V
#
# Define shell variables used in this script:
set TIFF_FILE       = SSGeoTiff_geo.tif
set CPT_FILE        = SSGeoTiff_geo.cpt
set MAP_REGION      = -122.243/-122.115/36.6939/36.8573
#
# Save existing GMT defaults
echo Saving GMT defaults...
gmtdefaults -L > gmtdefaults$$
#
# Set new GMT defaults
echo Setting new GMT defaults...
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      0   0   0   0   6.25  26  26  26 > $CPT_FILE
echo   6.25  26  26  26   12.5  51  51  51 >> $CPT_FILE
echo   12.5  51  51  51  18.75  77  77  77 >> $CPT_FILE
echo  18.75  77  77  77     25 102 102 102 >> $CPT_FILE
echo     25 102 102 102  31.25 128 128 128 >> $CPT_FILE
echo  31.25 128 128 128   37.5 153 153 153 >> $CPT_FILE
echo   37.5 153 153 153  43.75 179 179 179 >> $CPT_FILE
echo  43.75 179 179 179     50 204 204 204 >> $CPT_FILE
echo     50 204 204 204  56.25 230 230 230 >> $CPT_FILE
echo  56.25 230 230 230   62.5 255 255 255 >> $CPT_FILE
#
# Define data files to be plotted:
set DATA_FILE        = MARSEM300SSLower.grd
set INTENSITY_FILE   = 
#
# Make tiff image
echo Running mbgrdtiff...
mbgrdtiff -I $DATA_FILE \
	-O $TIFF_FILE \
	-R$MAP_REGION \
	-C $CPT_FILE \
	-V 
#
# 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 xv
echo Running xv in background...
xv SSGeoTiff_geo.tif &
#
# All done!
echo All done!