#!/bin/csh

/bin/rm tmp.ps              tmp.grd tmp.int tmp.cpt

set lon1 =   -120.00
set lon2 =    -20.00
set lat1 =    -60.00
set lat2 =      0.00
set a =      20.0


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


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

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

grd2cpt tmp.grd -Ctopo       -Z >tmp.cpt

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

pscoast -R -JM -Df -W2 -Na/1 -O -V -K >>tmp.ps              

/bin/rm tmp.grd tmp.int tmp.cpt
echo "Your output postscript file is called tmp.ps              "

########### The above were created by topoimag
## All you need to do is adding a "-K" in the
#pscoast line above, and adding the following line:

awk '{print $4, $3}' /home/sgao/progs/evloc/evloc.neic | \
psxy -R -JM -O -Sc0.05 -G255/0/0 -W0.1 -L >>tmp.ps