#!/bin/csh
#
# 3-D mesh plot of Hawaiian topography and geoid and cross-section color
#
echo '-10     255     0       255' >! zero.cpt
echo '0       100     10      100' >> zero.cpt
grdcontour geoid.grd -Jm0.45 -E60/30 -R195/210/18/25 -C1 -A5 -V -K -P -X1.5 -Y1.5 -U/-1.25/-1.25/"Example 4 in Cookbook" >! example_04.ps
pscoast -Jm -E60/30 -R -V -B2/2NEsw -G0 -O -K >> example_04.ps
echo '205 26 0 0 1.1' | psxyz -Jm -E60/30 -R -SV0.2/0.5/0.4 -W4 -O -K -N >> example_04.ps
echo '205 29.2 36 -90 1 5 N' | pstext -Jm -E60/30 -R -V -O -K -N >> example_04.ps
grdview topo.grd -Jm -Jz0.34 -Czero.cpt -E60/30 -R195/210/18/25/-6/4 -N-6/200/200/200 -Qsm -O -K -V -B2/2/2:"Topo (km)":neswZ -Y2.2 >> example_04.ps
echo '3.25 5.75 60 0.0 33 2 H@#awaiian@# R@#idge' | pstext -R0/10/0/10 -Jx1 -O >> example_04.ps
\rm zero.cpt
job4c