#!/bin/csh
gmt gmtset FORMAT_GEO_MAP DD
gmt gmtset MAP_FRAME_TYPE plain
gmt gmtset FONT_TITLE 24p
gmt gmtset MAP_TITLE_OFFSET -0.4c 
gmt gmtset FONT_LABEL 16p
gmt gmtset MAP_LABEL_OFFSET 0.3c
gmt gmtset FONT_ANNOT_PRIMARY 12p

set region =

set projection =

set open = "-K -V"

set add = "-O -K -V"

set close = "-O -V"

set psFile =

set pdfFile =

#*************************************Data examples*********************************************
#xyz2grd n35w115.bil -Gn35w115.grd -R-115/-114/35/36 -I3c -ZTLh -V

#To convert the Grid Float file downloaded from National Elevation Dataset (NED) Navigate to the directory that contains all of the files and use the following command:
#gdal_translate inputNED.flt -of GMT output.grd

#To convert .shp file to kml files
#ogr2ogr -f "KML" -where "NBRHOOD='Telegraph Hill'" realtor_neighborhoods.kml realtor_neighborhoods.shp
#To convert geo tif files:
#gdalwarp GMTtest.tiff GMTtest_northup.tiff
#gdalwarp -t_srs "+proj=latlong" GMTtest_northup.tiff GMTtest_latlon.tiff
#gdal_translate -b 1 -of GMT GMTtest_latlon.tiff GMTtest_latlon_red.grd
#gdal_translate -b 2 -of GMT GMTtest_latlon.tiff GMTtest_latlon_green.grd
#gdal_translate -b 3 -of GMT GMTtest_latlon.tiff GMTtest_latlon_blue.grd
#*************************************Data examples*********************************************


#gmt ps2raster $psFile -A -P -V -Tj -E300

ps2pdf $psFile $pdfFile

open -a Preview $pdfFile

echo "end of file"