328
Appendix B
C
C
C
C
1
HO,NP,NE,2)
CALL COEFF (X,Y,IJM,MA,MCPN,VX,VY,DXX,DXY,DYY,TR,BR,TH,HO,PS,
1
DTK,DSE,DPR,DAD,PSE,PPR,PAD,CFA,CFB,PT,NP,NP1,NE,NB,ISO,IH,2)
CALL SNIM (MA,MCPN,IE,IEE,CFA,CFB,NCBP,CO,C1,C2,CFOO,CFO,
1
DWR,DWC,DSR,DSC,TH,PS,PPR,PAD,EPS,ALF,RDC,
2
DTOO,DTO,DT,NP,NP1,NB,NCB,IH,2)
DO N=l, NP
IF ( CO(N) .LT. O. ) THEN
NA=MA(N+1)
NC=MA(N)
ND=NA-NC
SS=O.
DD=ND
DO K=2, ND
MC=MCPN(K,N)
SS=SS+CO(MC)
END DO
CO(N)=SS/DD
IF (CO(N) .LT. 0) CO(N)=O.
END IF
END DO
ENDIF
C writing head and concentration distributions
C
C
C
C
C
IF ( ABS(OBT(L)-TSUM) .LT. 1.E-5) THEN
IF (IPR .LT. 3) THEN
WRITE (6, 20) OBT(L)
20 FORMAT (/,lX, ' THE OBSERVATION TIME IS:', F10.2,/)
IF (IHC .EQ. 1) THEN
WRITE (6, 22)
22
FORMAT (I,6X,'NODE',11X, 'HEAD',/)
DO N=l,NP
WRITE (6,24) N, HO(N)
24
FORMAT (110, F15.2)
END DO
ELSE
WRITE (6, 26)
26
FORMAT (I,6X,'NODE',11X, 'HEAD',2X,'CONCENTRATION',/)
DO N=l, NP
WRITE (6, 28) N, HO(N), CO(N)
28
FORMAT (110, 2F15.2)
END DO
END IF
END IF
WRITE (11, 20) OBT(L)
IF (IHC .EQ. 1) THEN
WRITE (11, 30)
30
FORMAT (I,9X,'X-COOR', 9X, 'Y-COOR', 16X, 'HEAD',/)
DO N=l,NP
WRITE (11,32) X(N), Y(N), HO(N)
32
FORMAT (2F15.2, F20.3)
END DO
ELSE