#!/bin/csh # $Id: job13.csh,v 1.8 2010/07/09 00:21:22 guru Exp $ # GMT EXAMPLE 13 # # Purpose: Illustrate vectors and contouring # GMT progs: grdmath, grdcontour, grdvector, pstext # Unix progs: echo, rm grdmath -R-2/2/-2/2 -I0.1 X Y R2 NEG EXP X MUL = z.nc grdmath z.nc DDX = dzdx.nc grdmath z.nc DDY = dzdy.nc grdcontour dzdx.nc -JX3i -B1/1WSne -C0.1 -A0.5 -K -P -Gd2i -S4 -T0.1i/0.03i \ -U"Example 13 in Cookbook" >! example_13.ps grdcontour dzdy.nc -J -B1/1WSne -C0.05 -A0.2 -O -K -Gd2i -S4 -T0.1i/0.03i -X3.45i >> example_13.ps grdcontour z.nc -J -B1/1WSne -C0.05 -A0.1 -O -K -Gd2i -S4 -T0.1i/0.03i -X-3.45i -Y3.45i \ >> example_13.ps grdcontour z.nc -J -B1/1WSne -C0.05 -O -K -Gd2i -S4 -X3.45i >> example_13.ps grdvector dzdx.nc dzdy.nc -I0.2 -J -O -K -Q0.03i/0.1i/0.09in0.25i -G0 -S5i >> example_13.ps echo "3.2 3.6 40 0 6 BC z(x,y) = x * exp(-x@+2@+-y@+2@+)" | pstext -R0/6/0/4.5 -Jx1i -O -X-3.45i \ >> example_13.ps \rm -f z.nc dzdx.nc dzdy.nc .gmt*