334
Appendix B
C
46 FORMAT {lX, 'THE MODIFIED MA IS: 'I, (lX,1015))
WRITE
(6, 48) IH
48 FORMAT (lX, 'THE LENTH OF COEFFICIENT MATRIX IS',I6)
ENDIF
C ca1culating the area of each exclusive subdomain
C
C
DO I=l,NP
PS{I)=O.
END DO
DO N=l,NE
DO J=l, 3
ME{J)=IJM{J,N)
END DO
MI=ME{l)
MJ=ME(2)
MK=ME(3)
BI=Y{MJ)-Y{MK)
BJ=Y{MK) -Y{MI)
CI=X{MK)-X{MJ)
CJ=X{MI)-X{MK)
TRS=ABS{CJ*BI-BJ*CI) 16.
PS{MI)=PS{MI)+TRS
PS{MJ) =PS{MJ) +TRS
PS (MK) =PS{MK) +TRS
END DO
IF (IPR .LT. 2) THEN
WRITE
(6, 56)
56 FORMAT (lX, 'THE AREA OF EXCLUSIVE SUBDOMAIN OF EACH NODE: ',I)
DO N=l, NP
WRITE (6, 58) N, PS{N)
58 FORMAT (I10, F15.2)
END DO
ENDIF
100 RETURN
END
C********************************************************************
C
C
C
Subroutine DATAI reads initial conditions
C********************************************************************
C
C
SUBROUTINE DATAI (HOO,COO,NP,IHC)
IMPLICIT REAL*8 (A-H, O-Z)
DIMENSION HOO{NP) ,COO{NP)
C Reading initial conditions
C
C
C
CALL JUMP (3, 3)
DO 1=1, NP
IF (IHC .EQ. 1) THEN
READ (3, 10) HOO{I)
ELSE
READ (3, 20) HOO{I), COO{I)
10 FORMAT (10X, F20.5)
20 FORMAT (10X, 2F20.5)
END IF
END DO