C
B.6. Source Programs
329
WRITE (11, 34)
34
FORMAT (/, 9X, 'X-COOR' ,9X, 'Y-COOR', 16X, 'HEAD',
* 7X, 'CONCENTRATION',/)
DO N=l, NP
WRITE (11, 36) X(N), Y(N), HO(N), CO(N)
36
FORMAT (2F15.2, 2F20.3)
END DO
END IF
C writing solutions corresponding to observation wells and times
C
C
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)
IF (IHC .EQ. 1) THEN
WRITE (6, 22)
DO N=l,NOBW
WRITE (6,24) N, HOB(N,L)
END DO
ELSE
WRITE (6, 38)
38
FORMAT (I,4X,'NUMBER',l1X, 'HEAD',2X,'CONCENTRATION',1)
DO N=l, NOBW
WRITE (6, 28) N, HOB(N,L), COB(N,L)
END DO
END IF
END IF
C turning to the next time step
C
WRITE (*, 110) L
110 FORMAT(/,3X, 'The simulation of time period',i5,2x, 'is completed')
C
L=L+1
C
END IF
C
IF (L .GT. NOBT) GOTO 90
IF (DT .LT. DTL(l» DT=DT*DTF(l)
IF (DT+TSUM .GT. OBT(L» DT=OBT(L)-TSUM
C
GOTO 10
C
C generating the steady flow field
C
C
C
45 CALL FDB (TSUM,HB1,CB1,HO,CO,STEPB,NHBP,NCBP,
NP,NHB,NCB,NTB,IHC)
(TSUM,X,Y,IJM,WRAT,WCON,WLX,WLY,LEW,DWR,DWC,
*
*
CALL FDW
CALL FDS
STEPW,NP,NE,NW,NTW,IHC)
(TSUM,X,Y,IJM,SRAT,SCON,DSR,DSC,STEPS,
NSZON,NP,NE,NS,NTS,IHC)
CALL FVD (X,Y, IJM,DTK,DPR,DAL,DAT,DAD,DFU,DXX,DXY, DYY,VX,VY,
1
HO,NP,NE,l)
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,l)
50 TSUM = TSUM + DT