(上）光を北からあてた
（中）光を東からあてた
（下）光を南からあてた
シェルスクリプトはこちら

#!/bin/csh
#######################
set input=d1-ok.pmean
set output=sample3d.ps
set scale=m0.05
#set scale=A119/90/6.5
set range=50/180/40/70
set x=1
set y=8
set dlon=1
set dlat=1
set angle1=0
set angle2=90
set angle3=180
set grd1=data.grd1
set grd2=data.grd2
set grd3=data.hst1
set grd4=data.hst2
set fcolor=sample.cpt
set fontsize=25
set fonttype=7
###########################
surface $input -G$grd1 -R$range -I$dlon/$dlat -T0.0
grdgradient $grd1 -A$angle1 -G$grd2 -V
grdhisteq $grd2 -G$grd3 -N -V
grdmath $grd3 3. / = $grd4
grdimage $grd1 -J$scale -R$range -Ba10f10:."":neWS -C$fcolor -X$x -Y$y -P -I$grd4 -K >! $output
grdcontour $grd1 -J$scale -R$range -A4f9 -C4 -L900/1100 -W6 -P -O -K >> $output
######################
surface $input -G$grd1 -R$range -I$dlon/$dlat -T0.0
grdgradient $grd1 -A$angle2 -G$grd2 -V
grdhisteq $grd2 -G$grd3 -N -V
grdmath $grd3 3. / = $grd4
grdimage $grd1 -J$scale -R$range -Ba10f10:."":neWS -C$fcolor -X0 -Y-3.3 -P -I$grd4 -O -K >> $output
grdcontour $grd1 -J$scale -R$range -A4f9 -C4 -L900/1100 -W6 -P -O -K >> $output
######################
surface $input -G$grd1 -R$range -I$dlon/$dlat -T0.0
grdgradient $grd1 -A$angle3 -G$grd2 -V
grdhisteq $grd2 -G$grd3 -N -V
grdmath $grd3 3. / = $grd4
grdimage $grd1 -J$scale -R$range -Ba10f10:."":neWS -C$fcolor -X0 -Y-3.3 -P -I$grd4 -O -K >> $output
grdcontour $grd1 -J$scale -R$range -A4f9 -C4 -L900/1100 -W6 -P -O -K >> $output
####################
pstext -JX10/10 -R-5/10/-5/10 -X-1 -Y-4 -P -O << END >> $output
0 0 $fontsize 0 $fonttype 5 Sample-3A
END
\rm $grd1
\rm $grd2
\rm $grd3
\rm $grd4
gs $output
# もどります