#!/bin/csh
#
grdmath -R-2/2/-2/2 -I0.1 -V X Y R2 NEG EXP X MUL = z.grd
grdgradient z.grd -A270 -Gdzdx.grd
grdgradient z.grd -A180 -Gdzdy.grd
grdcontour dzdx.grd -JX3 -B1/1WSne -C0.1 -A0.5 -K -P -G2/10 -S4 -T0.1/0.03 -U"Example 13 in Cookbook" >! example_13.ps
grdcontour dzdy.grd -JX3 -B1/1WSne -C0.05 -A0.2 -O -K -G2/10 -S4 -T0.1/0.03 -X3.45 >> example_13.ps
grdcontour z.grd -JX3 -B1/1WSne -C0.05 -A0.1 -O -K -G2/10 -S4 -T0.1/0.03 -X-3.45 -Y3.45 >> example_13.ps
grdcontour z.grd -JX3 -B1/1WSne -C0.05 -O -K -G2/10 -S4 -X3.5 >> example_13.ps
grdvector dzdx.grd dzdy.grd -I0.2 -JX3 -O -K -Q0.03/0.10/0.09n0.25 -G0 -S5 >> example_13.ps
echo '3.2 3.6 40 0 6 2 z(x,y) = x * exp(-x@+2@+-y@+2@+)' | pstext -R0/8.5/0/5 -Jx1 -O -X-3.45 >> example_13.ps
\rm -f z.grd dzdx.grd dzdy.grd .gmtcommands