#!/bin/csh
#
gmtset ANOT_FONT_SIZE 12 LABEL_FONT_SIZE 13
gmtset HEADER_FONT_SIZE 24
#
echo '-10000   30    30    30    -8000   30    30    30' >! fig1.cpt
echo '-8000    60    60    60    -6000   60    60    60' >> fig1.cpt
echo '-6000    90    90    90    -4000   90    90    90' >> fig1.cpt
echo '-4000   120   120   120    -2000  120   120   120' >> fig1.cpt
echo '-2000   150   150   150    -1000  150   150   150' >> fig1.cpt
echo '-1000   180   180   180    -500   180   180   180' >> fig1.cpt
echo '-500    210   210   210    -100   210   210   210' >> fig1.cpt
echo '-100    240   240   240    0      240   240   240' >> 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 -R116/122/-6/2 -JM6.0 -I -Cfig1.cpt -P \
  -X1.0 -Y1.5 -K >! fig.ps
psscale -Cfig1.cpt -L -D3.15/-.5/4.0/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.0 -G0/230/0 \
  -Ba2f1SNEW \
  -Di -A40/3 -W -O >> fig.ps
#  -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 -RR116/122/-6/2 -JM6.0 -P -S -O >> fig.ps
#
mv fig.ps makassar.ps
rm aust-med.xyz