set input  = hight.date
set output = hight.ps
set lonmin = 96.5
set lonmax = 102.0166645
set latmin =  38.34
set latmax =  42.506665
set scale  = 3
set xshift = 4.0
set yshift = 4.0

set dlon   =  0.00833333
set dlat   =  0.00833333
set angle  = 90.0
set tmp1   = 3D.grd
set tmp2   = 3D.grad
set tmp3   = 3D.hst1
set tmp4   = 3D.hst2
set fcolor = color1.cpt

surface $input -G$tmp1 -R$lonmin/$lonmax/$latmin/$latmax -I$dlon/$dlat\
 -T0.0 -Ll-9.0

grdgradient $tmp1 -A$angle -G$tmp2 -V
grdhisteq $tmp2 -G$tmp3 -N -V
grdmath $tmp3 3. / = $tmp4

grdimage $tmp1 -Jm$scale -R$lonmin/$lonmax/$latmin/$latmax -Ba1/0.50.f0.25\
 -C$fcolor -I$tmp4 -X$xshift -Y$yshift >! $output

\rm $tmp1
\rm $tmp2
\rm $tmp3
\rm $tmp4