#!/bin/tcsh

set outfile = shubao.ps

set file = shubao.txt

gmtset BASEMAP_TYPE plain

psbasemap -R93/111/32/42 -JM18 -Ba2f1/a2f1WeSn -P  -X2.0 -Y4.0 -K > $outfile

grdimage asia_topo30.grd -CGMT_topo.cpt -R -JM18 -P -K -O -V >> $outfile

pscoast -R -JM18 -Lf95/33/50/200+l -Df -N2 -S32/128/250 -W2 -K -O -V >> $outfile

psxy CN-border-La.dat -J -R -M -W0.2p -O -K >> $outfile

#psvelo velo_ver_xm.txt -R -JM18 -R -A0.05/0.30/0.095 -Se0.2/0.4/15 -G255/000/000 -W2/000/000/000 -L -O -K -V >> $outfile

#awk '{print $1,$2}' $file  | psxy -R -JM18 -St0.2 -G255/0/0 -W1 -O -K >> $outfile

makecpt -Crainbow -T900/3600/200  > shubao.cpt

psxy haiba.txt -R -JM18 -Sc0.2 -Cshubao.cpt -W1 -O -K >> $outfile

awk '{print $1,$2,6,0,0,"CM",substr($3,1,4)}' $file | pstext -R -D0/-0.25 -JM18 -G0/0/0 -O -K >> $outfile

echo " 94.15 37.9 10 0 0 CM m" | pstext -R -J -G0/0/0 -O -K >> $outfile

psscale -D1/5/4/0.3 -Cshubao.cpt -Ba500f500 -K -O >> $outfile

psbasemap -R72/135/16/55 -JM5 -B0 -Gwhite -X13 -Y8.6 -K -O >> $outfile

pscoast -R -JM5 -Df -N1 -W -A5000 -K -O >> $outfile

psxy -R -JM5 -W2p,red  -O -L <<END  >> $outfile

093 32

093 42

111 42

111 32

END

ps2raster -A -P $outfile

ps2raster -A -P -Tf $outfile