#!/bin/csh
#
grdmath -R-15/15/-15/15 -I0.3 -V X Y HYPOT DUP 2 MUL PI MUL 8 DIV COS EXCH NEG 10 DIV EXP MUL = sombrero.grd
echo '-5        128     128     128     5       128     128     128' >! gray.cpt
grdgradient sombrero.grd -A225 -G=1 -N -V | grdhisteq =1 -Gintensity.grd -N -V
grdview sombrero.grd -JX6 -JZ2 -B5/5/0.5SEwnZ -N-1/255/255/255 -Qs -Iintensity.grd -X1.5 -Cgray.cpt -R-15/15/-15/15/-1/1 -K -V -E120/30 -U/-1.25/-0.75/"Example 5 in Cookbook" >! example_05.ps
echo '4.1 5.5 50 0 33 2 z(r) = cos (2@~p@~r/8) * e@+-r/10@+' | pstext -R0/11/0/8.5 -Jx1 -O >> example_05.ps
\rm -f gray.cpt sombrero.grd intensity.grd .gmtcommands