C
C
B.6. Source Programs
327
DAD,NPZON,RDC,DFU,NE,NZON,IHC)
CALL DATAW (WRAT,WCON,WLX,WLY,LEW,STEPW,NW,NTW,IHC)
CALL DATAS (SRAT,SCON,STEPS,NS,NTS,IHC)
CALL DATAO (OBWX,OBWY,LEOB,OBT,NOBW,NOBT)
CALL DATAC (DTI,DTF,DTL,ALF,EPS,PT,IHC)
WRITE (*, 102) NOBT
102 FORMAT (I, 5X, 'The total number of time periods is:', i5,/)
IF (IST .EQ. 1) THEN
WRITE (*, 103)
103 FORMAT (I, 5X, 'Finding the steady state .... ' ,I)
END IF
C Setting initial conditions
C
C
C
STEPW(NTW)=OBT(NOBT) +1.
STEPS(NTS) =OBT(NOBT) +1.
TSUM=O.
DTOO=DTI
DTO=DTI
DT=DTI
DO N=l, NP
HFOO(N)=O.
HFO(N)=O.
HO(N)=HOO(N)
IF (IHC .EQ. 2) THEN
CFOO(N)=O.
CFO (N) =0.
CO (N) =COO (N)
ENDIF
END DO
C Calculating the mass transport in transient flow field
C
C
C
L=l
IF (IST .EQ. 1) GOTO 45
10 TSUM=TSUM + DT
CALL FDB (TSUM,HBl,CBl,HO,CO,STEPB,NHBP,NCBP,
NP,NHB,NCB,NTB,IHC)
CALL FDW (TSUM,X,Y, IJM,WRAT,WCON,WLX,WLY, LEW, DWR,DWC,
STEPW,NP,NE,NW,NTW,IHC)
CALL FDS (TSUM,X,Y, IJM,SRAT,SCON,DSR,DSC, STEPS,
NSZON,NP,NE,NS,NTS,IHC)
C Solving the flow problem
C
C
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,NPl,NE,NB,ISO,IH,l)
CALL SNIM (MA,MCPN,IE,IEE,CFA,CFB,NHBP,HO,Hl,H2,HFOO,HFO,
1
DWR,DWC,DSR,DSC,TH,PS,PPR,PAD,EPS,ALF,RDC,
2
DTOO,DTO,DT,NP,NPl,NB,NHB,IH,l)
C Solving the mass transport problem
C
IF (IHC .EQ. 2) THEN
C
CALL FVD (X,Y, IJM,DTK,DPR,DAL,DAT,DAD,DFU, DXX, DXY, DYY,VX,VY,