set input1  = surfacedata
set input2  = highdata
set output  = surface.ps
set lonmin  = 99.0530701
set lonmax  = 101.4437801
set latmin  = 38.4404716
set latmax  = 40.2147616
set scale   = 6
set xshift  =   4.0
set yshift  =   4.0
set dlon1    =  0.00833
set dlat1    =  0.00833
set tmp1     =  land.grd
set fcolor1  =  eee.cpt
set dlon2    =  0.00833
set dlat2    =  0.00833
set tmp2     =  hight.grd

surface $input1 -G$tmp1 -R$lonmin/$lonmax/$latmin/$latmax -I$dlon1/$dlat1\
 -T0.0 -L0.0
grdimage $tmp1 -Jm$scale -R$lonmin/$lonmax/$latmin/$latmax -Ba0.5.f0.25\
 -C$fcolor1 -X$xshift -Y$yshift -K >! $output
\rm $tmp1

surface $input2 -G$tmp2 -R$lonmin/$lonmax/$latmin/$latmax -I$dlon2/$dlat2\
 -T0.0 -L0.0
grdcontour $tmp2 -Jm$scale -R$lonmin/$lonmax/$latmin/$latmax -Ba0.5.f0.25\
 -C500 -L0/5000 -A500 -W1a2   -O   >> $output
\rm $tmp2