C
C
C
C
C 20
C
C
C 25
C
C
C
C
B.6. Source Programs
343
DXX(N);(AL(N)*VX(N)*VX(N)+AT(N)*VY(N)*VY(N))/V + DFU
DYY(N);(AT(N) *VX(N) *VX(N)+AL(N) *VY(N) *VY(N) )/V + DFU
DXY(N);(AL(N)-AT(N))*(VX(N)*VY(N))/V
END DO
END1F
WR1TE (6, 20)
FORMAT (11, 3X, 'ELEMEN' ,8X, 'VX' ,8X, 'VY' ,7X, 'DXX' ,7X, 'DXY' ,7X, 'DYY' )
DO M;l, NE
WRITE (6, 25) M, VX(M) , VY(M) , DXX(M) , DXY(M), DYY(M)
FORMAT (110, 5El0.5)
END DO
RETURN
END
C********************************************************************
C
C
Subroutine COEFF forms coefficient matrices of MCB equations
C
C********************************************************************
C
C
C
C
C
C
C
C
SUBROUTINE COEFF (X,Y,IJM,MA,MCPN,VX,VY,DXX,DXY,DYY,TR,BR,TH,HO,
1
PS,TK,SE,PR,AD,PSE,PPR,PAD,CFA,CFB,PT,NP,NP1,NE,NB,ISO,IH,IHC)
IMPLICIT REAL*8 (A-H, O-Z)
DIMENSION X(NP),Y(NP) ,IJM(3,NE) ,MA(NP1) ,MCPN(NB,NP),
1
VX(NE) ,VY(NE) ,DXX(NE) ,DXY(NE) ,DYY(NE) ,TR(NP) ,BR(NP),
2
TH(NP) ,TK(NE) ,SE(NE) ,PR(NE) ,AD(NE) ,PSE(NP) ,PPR(NP),
3
PAD(NP) ,PS (NP) ,HO (NP) ,CFA(IH) ,CFB(IH) ,ME(3)
DO K;l, IH
CFA(K);O.
CFB(K);O.
END DO
DO 1;1, NP
PSE(I);O.
PPR(I);O.
PAD(I);O.
END DO
DO 1;1, NP
IF (HO(I) .GT. TR(I)) THEN
TH(I);TR(I)-BR(I)
ELSE
TH(I);HO(I)-BR(I)
ENDIF
END DO
DO N;l, NE
MI;IJM(l,N)
MJ;IJM(2,N)
MK;IJM(3,N)
BI;Y(MJ)-Y(MK)
BJ;Y(MK) -Y(MI)