348
Appendix B
C
C
IF (IHC .EQ. 1) THEN
DD=DD+DWR(N)
DD=DD+DSR(N)
ELSE
RET=l.+(l.-PPR(N))*PAD(N)/PPR(N)
IF (DWR(N) .GT. 0.) THEN
DD=DD+DWR(N)*(DWC(N)-HO(N))/TH(N)/PPR(N)/RET
END IF
DD=DD+DSR(N)*(DSC(N)-HO(N))/TH(N)/PPR(N)/RET
ENDIF
C The selecting node iteration procedure
C
C
C
C
C
C
C
H2(N)=(1.-ALF)*Hl(N)+ALF*DD/TRS
IF ( ABS(H2(N)-H1(N)) .LT. EPS(IHC)) IE(N)=l
Hl (N) =H2 (N)
ENDIF
END DO
DO N=l, NP
IEE(N)=l
END DO
DO N=l,NP
IF(IE(N) .EQ. 0) THEN
IEE(N)=O
NA=MA(N+1)
NC=MA(N)
ND=NA-NC
DO K=2,ND
MC=MCPN(K,N)
IF(IE(MC). LT. 2) IEE(MC)=O
END DO
ENDIF
END DO
DO N=l, NP
IF(IE(N) .LT. 2) IE(N)=IEE(N)
END DO
DO N=l ,NP
IF(IE(N) .EQ. 0) GO TO 40
END DO
C The solution is obtained for this time step.
C
C
C
DTOO=DTO
DTO=DT1
DO N=l,NP
IF(IE(N) .LT.2) THEN
DFOO(N)=DFO(N)
DFO(N)=HO(N)-H1(N)
HO (N) =H1 (N)
ENDIF
END DO
RETURN
END