set input   = hyoukou.data
set output  = hyoukou.ps
set lonmin  = 96.5
set lonmax  = 101.9978
set latmin  = 38.3499985
set latmax  = 42.4983385
set scale   = 2
set xshift  =   4.0
set yshift  =   4.0

set dlon    =  0.00833
set dlat    =  0.00833
set tmp     =  hight.grd
set fcolor  =  xxx.cpt　　　　　　　

surface $input -G$tmp -R$lonmin/$lonmax/$latmin/$latmax -I$dlon/$dlat\
 -T0.0 -L-1.0

grdimage $tmp -Jm$scale -R$lonmin/$lonmax/$latmin/$latmax -Ba1.0/0.5.f0.25\
 -C$fcolor -X$xshift -Y$yshift -K >! $output

grdcontour $tmp -Jm$scale -R$lonmin/$lonmax/$latmin/$latmax -Ba1.0/0.5.f0.25\
 -C500 -L0/5000 -A500 -W1a2   -O   >> $output

\rm $tmp