#!/bin/csh

/bin/rm tmp.ps  tmp.cpt

set lon1 =      5.4 
set lon2 =      5.7
set lat1 =      22.65
set lat2 =      22.9

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

grdcut /home/sgao/Dbase/gtopo30/north.grd  \
 -Gtmp1.grd -R  -V

grdsample tmp1.grd -Gtmp.grd -I0.002

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

#grd2cpt tmp.grd -Csealand  -Z >tmp.cpt
makecpt -Cpanoply -T1200/1800/100 -Z >tmp.cpt

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

grdcontour tmp.grd -T+ -J -R -C50 -A50 -W2 -O -K >>tmp.ps
grdcontour tmp.grd -T- -J -R -C100 -A100 -W2 -O -K >>tmp.ps

psscale -D2.5/-0.4/5/0.2h -Ctmp.cpt -O -K >>tmp.ps

pstext -R -J -N -O -K <<END>>tmp.ps
5.55 22.59 14 0 5 6 Elevation (m)
END

### Draw circles
awk '{print $2, $1}' circ1.out | \
psxy  -R -J -W8/255 -O -K >>tmp.ps

awk '{print $2, $1}' circ2.out | \
psxy  -R -J -W8/255 -O -K >>tmp.ps

awk '{print $2, $1}' circ3.out | \
psxy  -R -J -W8/255 -O -K >>tmp.ps

##--------------

pscoast -R -J -Df -Lf5.675/22.675/22.7/4 -W2 -Na/1 -O -K >>tmp.ps

psxy pie.dat -R -J -Sx0.1 -W5/255/0/0 -O -K >>tmp.ps

psxy -R -J -B -O -St0.2 -G0/255/0 <<END>>tmp.ps
    5.5270    22.7910
END


ps2jpeg150 tmp.ps