#!/bin/csh
#
# This script will make two contour maps based on the data in the file osu91a1f_16.grd
#
gmtset GRID_CROSS_SIZE 0 ANOT_FONT_SIZE 10
psbasemap -R0/6.5/0/9 -Jx1 -B0 -P -K -U"Example 1 in Cookbook" >! example_01.ps
pscoast -R-180/180/-90/90 -JH0/6 -X0.25 -Y0.5 -O -K -Bg30 -Dc -G200 >> example_01.ps
grdcontour -R osu91a1f_16.grd -JH -C10 -A50f7 -L-1000/-1 -Wc1ta -Wa3t8_8:0 -O -K -T0.1/0.02 >> example_01.ps
grdcontour -R osu91a1f_16.grd -JH -C10 -A50f7 -L-1/1000 -O -K -T0.1/0.02 >> example_01.ps
pscoast -R0/360/-90/90 -JH180/6 -Y4 -O -K -Bg30:."Low Order Geoid": -Dc -G200 >> example_01.ps
grdcontour osu91a1f_16.grd -JH -C10 -A50f7 -L-1000/-1 -Wc1ta -Wa3t8_8:0 -O -K -T0.1/0.02:-+ >> example_01.ps
grdcontour osu91a1f_16.grd -JH -C10 -A50f7 -L-1/1000 -O -T0.1/0.02:-+ >> example_01.ps
\rm -f .gmtcommands