#!/bin/csh

/bin/rm tmp.ps              
/bin/rm tmp.cpt

set lon1 =   -180.00
set lon2 =    180.00
set lat1 =    -80.00
set lat2 =     80.00
set a =      30.0


psbasemap -R$lon1/$lon2/$lat1/$lat2 -JM5 -Ba30f10 \
-V -P -K -X1.5 -Y2 >tmp.ps              


grdcut /home/sgao/Dbase/topo30min/data.grd       \
-Gtmp.grd -R  -V

grdgradient tmp.grd -A0 -Gtmp.int -Nt -V

makecpt -Cgray -T-10000/0/1000 -Z >tmp.cpt

grdimage tmp.grd -Itmp.int -Ctmp.cpt -R -JM -O -K >>tmp.ps              


pscoast -R -JM -Dh -W2 -A1000 -N1/1 -O -K >>tmp.ps              

awk '{print $1, $2}' routes.dat | \
psxy -R -M -W8/255/0/0 -JM -O -K >>tmp.ps

awk '{print $1, $2}' locs.dat | \
psxy -R -JM -O -K -Sc0.1 -G0/0/255 >>tmp.ps

pstext -R -N -JM -O <<END>>tmp.ps
0 85 19 0 0 6 Steve Gao's Work-related Trips Since 1999
END