#!/bin/sh
#  Based on GMT EXAMPLE 05
# $Id: wave,v 1.1 2001-08-27 13:20:16-06 braup Exp $

rm .gmtdefaults
outfile=`basename $0`.ps
grdfile=`basename $0`.grd
colorfile='color.cpt'
radius=7
grdmath -R-15/15/-15/15 -I0.2 X Y HYPOT DUP $radius GT EXCH DUP 2 MUL PI MUL 4 DIV COS EXCH NEG 6 DIV EXP MUL MUL \
   = $grdfile
#echo '-5	128	5	128' >| $colorfile
makecpt -Csealand -T-1/1/0.1 -Z -V >| $colorfile
grdgradient $grdfile -A225 -Gintensity.grd -Nt0.75
grdview $grdfile -JX13.0c -JZ5c -B5/5/0.5SEwnZ -N-1/255/255/255 -Qs -Iintensity.grd -X1.5c -K \
   -C$colorfile -R-15/15/-15/15/-1/1 -E120/30 -P >| $outfile
echo "3.6 5.5 30 0 33 CM Furuike ya, Kawazu Tobikomu, Mizu no Oto" | pstext -R0/11/0/8.5 -Jx2.5c -O \
   >> $outfile
rm -f $colorfile $grdfile intensity.grd