#!/bin/csh

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

set lon1 =   -127.00
set lon2 =    -60.00
set lat1 =     25.00
set lat2 =     50.00


psbasemap -R$lon1/$lon2/$lat1/$lat2 -JM7 -Ba5f1g5 \
-V  -K -X0.5 -P -Y2 >tmp.ps              


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

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

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

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

pscoast -R -JM -B -S255 -Df -A100 -W12/255/0/0 -Na/12/255/0/0 -O -V >>tmp.ps