#!/bin/tcsh
#
#  Example of using "pscontour" and "pscoast" to create a map rather
#  than "pscontour", "psbasemap", and "grdcontour" (see other example).
#  The "pscoast" program does the job of both "psbasemap" and
#  "grdcontour" (replacing the latter with its own internal
#  geographical information).
#
gmtset ANOT_FONT_SIZE 12 LABEL_FONT_SIZE 13
gmtset HEADER_FONT_SIZE 24
#
echo '-0.5  192   217   217    0.5  192   217   217' >! fig1.cpt
echo '0.5    50   153   204    1.5   50   153   204' >> fig1.cpt
echo '1.5    35   142    35    2.5   35   142    35' >> fig1.cpt
echo '2.5   153   204    50    3.5  153   204    50' >> fig1.cpt
echo '3.5    33    94    33    4.5   33    94    33' >> fig1.cpt
echo '4.5     0     0   156    5.5    0     0   156' >> fig1.cpt
echo '5.5   159    95   159    6.5  159    95   159' >> fig1.cpt
echo '6.5   255   255     0    7.5  255   255     0' >> fig1.cpt
echo '7.5   255     0     0    8.5  255     0     0' >> fig1.cpt
echo '8.5   166    42    42    9.5  166    42    42' >> fig1.cpt
echo '9.5   217   135    25   10.5  217   135    25' >> fig1.cpt
echo '10.5  255   255   255   11.5  255   255   255' >> fig1.cpt
echo 'B     255   255   255'                          >> fig1.cpt
echo 'F     255   255   255'                          >> fig1.cpt
#
echo '3.25 3.35 18 0 1 6 LGM VEGETATION BIOMES - Adams' >! temp1.title
pstext temp1.title -R0/6.5/0/9 -Jx1 -P -K >! ex10.ps
echo '3.25 8.8 18 0 1 6 LGM VEGETATION BIOMES - Crowley' >! temp1.title
pstext temp1.title -R0/6.5/0/9 -Jx1 -P -O -K >> ex10.ps
#
grdimage veg-LGM-GBC.grd -R-180/180/-90/90 -Jx0.016d -V -X0.5 -Y5.7 \
  -Cfig1.cpt -O -K >> ex10.ps
#
# Draw continental outlines and box around figure.
#
#  -G240 for light gray shaded land
#  -S0/0/250 for blue shaded ocean
#  -G0/230/0 for green shaded land
psbasemap -R-180/180/-90/90 -Jx0.016d -Ba60f30/a30f15SWne -O -K >> ex10.ps
#pscoast -R-180/180/-90/90 -Jx0.024d -Ba60f30SNEW -Dc -W0 -O -K >> ex10.ps
#
echo '0.5  2.25 ' >! sym1.xy
psxy sym1.xy -R0/6.5/0/9 -Jx1 -Ss0.20 -L -G50/153/204 -O -K -Y-3.0 >> ex10.ps
echo '0.75  2.25  14  0  1  5  Trop. RF' >! sym1.xy
pstext sym1.xy -R0/6.5/0/9 -Jx1 -O -K >> ex10.ps 
#
echo '0.5  1.95 ' >! sym1.xy
psxy sym1.xy -R0/6.5/0/9 -Jx1 -Ss0.20 -L -G35/142/35 -O -K >> ex10.ps
echo '0.75  1.95  14  0  1  5  Deciduous' >! sym1.xy
pstext sym1.xy -R0/6.5/0/9 -Jx1 -O -K >> ex10.ps
#
echo '0.5  1.65 ' >! sym1.xy
psxy sym1.xy -R0/6.5/0/9 -Jx1 -Ss0.20 -L -G153/204/50 -O -K >> ex10.ps
echo '0.75  1.65  14  0  1  5  Mixed Decid.' >! sym1.xy
pstext sym1.xy -R0/6.5/0/9 -Jx1 -O -K >> ex10.ps
#
echo '0.5  1.35 ' >! sym1.xy
psxy sym1.xy -R0/6.5/0/9 -Jx1 -Ss0.20 -L -G33/94/33 -O -K >> ex10.ps
echo '0.75  1.35  14  0  1  5  Conifer' >! sym1.xy
pstext sym1.xy -R0/6.5/0/9 -Jx1 -O -K >> ex10.ps
#
#
echo '1.6  2.55 ' >! sym1.xy
psxy sym1.xy -R0/6.5/0/9 -Jx1 -Ss0.20 -L -G255/255/255 -O -K >> ex10.ps
echo '1.85  2.55  14  0  1  5  Trop. Grass' >! sym1.xy
pstext sym1.xy -R0/6.5/0/9 -Jx1 -O -K >> ex10.ps
#
echo '3.4  2.55 ' >! sym1.xy
psxy sym1.xy -R0/6.5/0/9 -Jx1 -Ss0.20 -L -G255/255/0 -O -K >> ex10.ps
echo '3.65  2.55  14  0  1  5  Temp. Grass' >! sym1.xy
pstext sym1.xy -R0/6.5/0/9 -Jx1 -O -K >> ex10.ps
#
echo '2.5  1.05 ' >! sym1.xy
psxy sym1.xy -R0/6.5/0/9 -Jx1 -Ss0.20 -L -G255/255/0 -O -K >> ex10.ps
echo '2.75  1.05  14  0  1  5  Grass' >! sym1.xy
pstext sym1.xy -R0/6.5/0/9 -Jx1 -O -K >> ex10.ps
#
#
echo '2.5  2.10 ' >! sym1.xy
psxy sym1.xy -R0/6.5/0/9 -Jx1 -Ss0.20 -L -G159/95/159 -O -K >> ex10.ps
echo '2.75  2.10  14  0  1  5  Savanna' >! sym1.xy
pstext sym1.xy -R0/6.5/0/9 -Jx1 -O -K >> ex10.ps
#
echo '2.5  1.80 ' >! sym1.xy
psxy sym1.xy -R0/6.5/0/9 -Jx1 -Ss0.20 -L -G255/0/0 -O -K >> ex10.ps
echo '2.75  1.80  14  0  1  5  Scrub/Grass' >! sym1.xy
pstext sym1.xy -R0/6.5/0/9 -Jx1 -O -K >> ex10.ps
#
echo '2.5  1.50 ' >! sym1.xy
psxy sym1.xy -R0/6.5/0/9 -Jx1 -Ss0.20 -L -G192/217/217 -O -K >> ex10.ps
echo '2.75  1.50  14  0  1  5  Ocean/Lakes' >! sym1.xy
pstext sym1.xy -R0/6.5/0/9 -Jx1 -O -K >> ex10.ps
#
#
echo '4.5  2.25 ' >! sym1.xy
psxy sym1.xy -R0/6.5/0/9 -Jx1 -Ss0.20 -L -G166/42/42 -O -K >> ex10.ps
echo '4.75  2.25  14  0  1  5  Semi-desert' >! sym1.xy
pstext sym1.xy -R0/6.5/0/9 -Jx1 -O -K >> ex10.ps
#
echo '4.5  1.95 ' >! sym1.xy
psxy sym1.xy -R0/6.5/0/9 -Jx1 -Ss0.20 -L -G0/0/156 -O -K >> ex10.ps
echo '4.75  1.95  14  0  1  5  Tundra' >! sym1.xy
pstext sym1.xy -R0/6.5/0/9 -Jx1 -O -K >> ex10.ps
#
echo '4.5  1.65 ' >! sym1.xy
psxy sym1.xy -R0/6.5/0/9 -Jx1 -Ss0.20 -L -G217/135/25 -O -K >> ex10.ps
echo '4.75  1.65  14  0  1  5  Desert' >! sym1.xy
pstext sym1.xy -R0/6.5/0/9 -Jx1 -O -K >> ex10.ps
#
echo '4.5  1.35 ' >! sym1.xy
psxy sym1.xy -R0/6.5/0/9 -Jx1 -Ss0.20 -L -G255/255/255 -O -K >> ex10.ps
echo '4.75  1.35  14  0  1  5  Land Ice' >! sym1.xy
pstext sym1.xy -R0/6.5/0/9 -Jx1 -O -K >> ex10.ps
#
#  Add second image of Global Biogeo. Cycles version.
#
grdimage veg-adams.grd -R-180/180/-90/90 -Jx0.016d -V -X0.0 -Y-2.5 \
  -Cfig1.cpt -O -K >> ex10.ps
psbasemap -R-180/180/-90/90 -Jx0.016d -Ba60f30/a30f15SWne -O >> ex10.ps

rm *.xy
rm *.cpt
rm *.title
mv ex10.ps veg-LGM-GBC-adams.ps