#!/bin/csh
#get and project the files
gmtset HEADER_FONT_SIZE 30 OBLIQUE_ANOTATION 0 DEGREE_FORMAT 0
#grdraster 1 -R160/20/220/30r -JOc190/25.5/292/69/4.5 -Getopo5.grd=0/0.001/0
#grdraster 4 -R -JO -Ggeosat5.grd
# Make illustration
grdimage geosat5.grd -R160/20/220/30r -JOc190/25.5/292/69/4.5 -E50 -K -P -U/-1.25/-1/"Example 2 in Cookbook" -V -B10 -Cgeoid.cpt -X1.5 -Y1.25 >! example_02.ps
psscale -Cgeoid.cpt -D5.1/1.35/2.88/0.4 -O -K -L -B:m::.GEOID: >> example_02.ps
grdimage etopo5.grd -R -JO -E50 -B10:."H@#awaiian@# T@#opo and @#G@#eoid:" -O -K -V -Ctopo.cpt -Y4.5 >> example_02.ps
psscale -Ctopo.cpt -D5.1/1.35/2.88/0.4 -O -K -B:m::.TOPO: >> example_02.ps
echo "-0.4 7.5 36 0.0 1 2 a)" | pstext -R0/8.5/0/11 -Jx1 -O -K -Y-4.5 >> example_02.ps
echo "-0.4 3.0 36 0.0 1 2 b)" | pstext -R0/8.5/0/11 -Jx1 -O >>  example_02.ps
#rm -f *.grd .gmtcommands
# For Suns only
#
#ps2ras -w 9 -h 11.52 -d 24 -p 100 -b 255 255 255 example_02.ps example_02.32
#24to8 -v -large example_02.32 tmp
#compact tmp example_02.ras
#rm -f tmp example_02.32