2次元量表示の最後のプログラム例 U2D7 は, 等高線とベクトル場の重 ね書きです. これまでに見てきた内容の描画を1つのフレームにまとめたもの です. 'UDLSET' ルーチンで等高線図に関する内部変数 'LABEL' と 'LMSG' を .FALSE. にして, ラベルとメッセージを描かない ように設定し, そのかわり, 右側にトーンバーをつけるようにしました. 61行 め以降で, もう一度正規化変換を行ない, UDCNTR と UETONE ルー チンでトーンバーを描いています.
PROGRAM U2D7
PARAMETER( NX=21, NY=21 )
PARAMETER( XMIN=-10, XMAX=10, YMIN=-10, YMAX=10 )
PARAMETER( DX1=1, DX2=5, DY1=1, DY2=5 )
PARAMETER( KMAX=5, PMIN=0, PMAX=1 )
REAL U(NX,NY), V(NX,NY), P(NX,NY), PI(2,KMAX+1)
DO 10 J=1,NY
DO 10 I=1,NX
X = XMIN + (XMAX-XMIN)*(I-1)/(NX-1)
Y = YMIN + (YMAX-YMIN)*(J-1)/(NY-1)
U(I,J) = X
V(I,J) = -Y
P(I,J) = EXP( -X**2/64 -Y**2/25 )
10 CONTINUE
WRITE(*,*) ' WORKSTATION ID (I) ? ;'
CALL SGPWSN
READ (*,*) IWS
CALL GROPN( IWS )
CALL SGLSET( 'LSOFTF', .TRUE. )
CALL GRFRM
CALL GRSWND( XMIN, XMAX, YMIN, YMAX )
CALL GRSVPT( 0.2, 0.8, 0.2, 0.8 )
CALL GRSTRN( 1 )
CALL GRSTRF
CALL UXAXDV( 'B', DX1, DX2 )
CALL UXAXDV( 'T', DX1, DX2 )
CALL UXSTTL( 'B', 'X-axis', 0. )
CALL UXMTTL( 'T', 'FIGURE TITLE', 0. )
CALL UYAXDV( 'L', DY1, DY2 )
CALL UYAXDV( 'R', DY1, DY2 )
CALL UYSTTL( 'L', 'Y-axis', 0. )
CALL UGLSET( 'LUNIT' , .TRUE. )
CALL UGRSET( 'VXUOFF', 0.04 )
CALL UGSUT( 'X', 'U' )
CALL UGSUT( 'Y', 'V' )
CALL UGVECT( U, NX, V, NX, NX, NY )
CALL UDLSET( 'LABEL', .FALSE. )
CALL UDLSET( 'LMSG' , .FALSE. )
CALL UDGCLB( P, NX, NX, NY, 0.1 )
CALL UDCNTR( P, NX, NX, NY )
DP = (PMAX-PMIN)/KMAX
DO 20 K=1,KMAX
TLEV1 = (K-1)*DP
TLEV2 = TLEV1 + DP
IPAT = 600 + K - 1
CALL UESTLV( TLEV1, TLEV2, IPAT )
20 CONTINUE
CALL UETONE( P, NX, NX, NY )
*-- トーンバー ----
CALL GRSWND( 0.0, 1.0, PMIN, PMAX )
CALL GRSVPT( 0.85, 0.9, 0.45, 0.75 )
CALL GRSTRN( 1 )
CALL GRSTRF
DO 30 K=1,KMAX+1
PI(1,K) = PMIN + (K-1)*DP
PI(2,K) = PMIN + (K-1)*DP
30 CONTINUE
CALL UDCNTR( PI, 2, 2, KMAX+1 )
CALL UETONE( PI, 2, 2, KMAX+1 )
CALL SLPVPR( 3 )
CALL UZLSET( 'LABELYR', .TRUE. )
CALL UZFACT( 0.8 )
CALL UYSFMT( '(F4.1)' )
CALL UYAXDV( 'R', DP, DP )
CALL GRCLS
END
