#!/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 '-10000    0     0   156    -8000    0     0   156' >! fig1.cpt
echo '-8000    35   107   142    -6000   35   107   142' >> fig1.cpt
echo '-6000    50    50   205    -4000   50    50   205' >> fig1.cpt
echo '-4000    77    77   255    -2000   77    77   255' >> fig1.cpt
echo '-2000     0   127   255    -1000    0   127   255' >> fig1.cpt
echo '-1000    50   153   204    -500    50   153   204' >> fig1.cpt
echo '-500    143   143   189    -100   143   143   189' >> fig1.cpt
echo '-100    192   217   217    0      192   217   217' >> fig1.cpt
echo 'B     255   255   255'                          >> fig1.cpt
echo 'F     255   255   255'                          >> fig1.cpt
#
grd2xyz aust-med.grd > aust-med.xyz
pscontour aust-med.xyz -V -R124/136/-12/2 -JM6.5 -I -Cfig1.cpt -P \
  -X1.0 -Y2.0 -K >! fig.ps
psscale -Cfig1.cpt -D3.25/-.5/3.5/0.15h -O -K >> fig.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
pscoast -R -N1/1ta -N3 -JM6.5 -G0/230/0 \
  -Ba2f1SNEW \
  -Di -A40/3 -W -O -K >> fig.ps
#  -Ba2f1SNEW:."BALTIC SEA": \
#  -Di -A20/3 -K -W >! fig.ps
#
# -Di - data set choice: c(rude), l(ow), i(ntermediate), h(igh)
# -G200 - specify land fill attributes, see p. 13
# -N1 - draw national boundaries
# -N3 - draw marine boundaries
# -A20/3 - draw areas larger than 20 km2 and level 3 (see pscoast
#      manpage)
# -W - draw coastlines
#
#  Add geographical names.
#
pstext aust-med.txt -R124/136/-12/2 -JM6.5 -P -S -O >> fig.ps
#
mv fig.ps ex01.ps
rm aust-med.xyz