C
C
C
C
C
C
SS=SS+CO(MC)
END DO
CO (N) =SS/DD
IF (CO(N) .LT. 0) CO(N)=O.
END IF
END DO
IF ( ABS(OBT(L)-TSUM) .LT. 1.E-5) THEN
IF (IPR .LT. 3) THEN
WRITE (6, 20) OBT(L)
WRITE (6, 62)
B.6. Source Programs
331
62
FORMAT (I,6X, 'NODE', 2X, 'CONCENTRATION',I)
DO N=l, NP
WRITE (6, 24) N, CO(N)
END DO
END IF
WRITE (11, 20) OBT(L)
WRITE (11, 64)
64
FORMAT (I,9X, 'X-COOR',9X,'Y-COOR', 7X, 'CONCENTRATION',I)
DO N=l, NP
WRITE (11, 32) X(N), Y(N), CO(N)
END DO
C writing solutions corresponding to observation wells and times
C
C
CALL FOB (X,Y,IJM,HO,CO,HOB,COB,OBWX,OBWY,LEOB,
L,NP,NE,NOBW,NOBT,IHC)
IF (IPR .LT. 4) THEN
WRITE (6, 20) OBT(L)
WRITE (6, 62)
DO N=l, NOBW
WRITE (6, 24) N, COB(N,L)
END DO
ENDIF
C turning to the next time per iod
C
C
C
C
C
C
WRITE (*, 110) L
L=L+1
ENDIF
IF (L .LE. NOBT) THEN
IF (DT .LT. DTL(2)) DT=DT*DTF(2)
IF (DT+TSUM .GT. OBT(L)) DT=OBT(L)-TSUM
GOTO 60
ENDIF
ENDIF
C Form data file DRAWT
C
90
DO N=l, NOBW
WRITE (12, 92) N
92
FORMAT (/, 2X, 'THE TIME-HEAD-CONCENTRATION FOR OBS. WELL', 13)
WRITE (12, 94)