332
Appendix B
C
C
C
C
94
FORMAT (I, llX, 'TIME', 16X, 'HEAD', 7X, 'CONCENTRATION', I)
DO L=l, NOBT
WRITE (12, 96) OBT(L), HOB(N,L), COB(N,L)
96
FORMAT (F15.3, 2F20.3)
END DO
END DO
RETURN
END
C********************************************************************
C
C
C
C
Subroutine DATAG reads and processes geometrie and
geologie information
C*******************************************************************
C
SUBROUTINE DATAG (X,Y,TR,BR,TH,PS,IJM,MA,MCPN,NPZON,NSZON,
1
NP,NP1,NE,NB,IH,IPR)
IMPLICIT REAL*8 (A-H, O-Z)
DIMENSION X(NP) ,Y(NP) ,TR(NP) ,BR(NP),TH(NP),PS(NP),
1
IJM(3,NE) ,MA(NP1) ,MCPN(NB,NP) ,NPZON(NE) ,
2
NSZON(NE),ME(3),CHB(70)
C
C reading nodal coordinates
C
C
C
CALL JUMP (2, 2)
READ (2, 1) CHB
1 FORMAT (70A1)
DO
N=l, NP
READ (2, 2) X(N), Y(N), TR(N) , BR(N)
TH(N)=TR(N)-BR(N)
2 FORMAT (10X, 2F15.2, 2F10.2)
END DO
IF (IPR .LT. 1) THEN
WRITE (6, 4) CHB
4 FORMAT (li, 70A1, I)
DO N=l, NP
WRITE(6, 6) N, X(N), Y(N), TR(N) , BR(N)
6 FORMAT (110, 2F15.2, 2F10.2)
END DO
ENDIF
C reading node numbers of eaeh element and zonation patterns
C
C
CALL JUMP (2, 2)
READ (2, 1) CHB
DO M=l, NE
READ (2, 12) IJM(l,M) ,IJM(2,M) ,IJM(3,M) ,NPZON(M) ,NSZON(M)
12 FORMAT (10X, 5110)
END DO
IF (IPR .LT. 1) THEN
WRITE (6, 14) CHB
14 FORMAT (I, 70A1,/)
DO M=l, NE
WRITE (6, 16) M,IJM(1,M),IJM(2,M),IJM(3,M) ,NPZON(M) ,NSZON(M)
16 FORMAT (6110)