330
Appendix B
C
C
CALL SNIM
1
2
(MA,MCPN,IE,IEE,CFA,CFB,NHBP,HO,H1,H2,HFOO,HFO,
DWR,DWC,DSR,DSC,TH,PS,PPR,PAD,EPS,ALF,RDC,
DTOO,DTO,DT,NP,NP1,NB,NHB,IH,1)
C Testing the steady condition
C
C
C
C
DO N=l, NP
DD=HFO(N)
IF ( ABS(DD)/DTO .GT. EPS(l) ) THEN
IF (DT .LT. 10.) DT=DT*DTF(l)
GOTO 50
ENDIF
END DO
IF (IPR .LT. 4) THEN
WRITE (6, 52)
52 FORMAT (/,lX, 'THE STEADY STATE HEAD DISTRIBUTION' ,1,6X, 'NODE',
*
11X, ' HEAD' , /)
DO N=l, NP
WRITE (6, 54) N, HO(N)
54 FORMAT (I10, F15.2)
END DO
END IF
WRITE (11, 56)
56 FORMAT (/,lX, 'THE STEADY STATE HEAD DISTRIBUTION' ,I)
WRITE (11, 30)
DO N=l, NP
WRITE (11, 32) X(N), Y(N), HO(N)
END DO
C calculating mass transport in the steady flow field
C
C
C
C
IF (IHC .EQ. 2) THEN
CALL FVD (X,Y, IJM,DTK,DPR,DAL,DAT,DAD,DFU, DXX, DXY, DYY,VX,VY,
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)
L=l
TSUM=O.
DTOO=DTI
DTO=DTI
DT=DTI
60
TSUM=TSUM + DT
CALL SN IM (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. 0. ) THEN
NA=MA(N+1)
NC=MA(N)
ND=NA-NC
SS=O.
DD=ND
DO K=2, ND
MC=MCPN(K,N)