CALL PITZER (T,2,EXTMP)
IF (AH2O.LE.SU(IV)) GO TO 10
2900
CONTINUE
10
X(IV)=X(IV)/XMN*1.2
DO 3000 I=1,ICLOS
J=JCLOS(I)
NJ=NUM(J)
EX(NJ)=EXTMP(NJ)/1.2
G(J)=G(NJ)+DLOG(EX(NJ))+ACT(NJ)
3000
CONTINUE
ELSE
DO 3200 I=1,IP
DO 3201 J=IP1, NL
IF(NUMIN(I).EQ.J) THEN
X(J)=X(J)-.95*A(IR,N1)*AIN(I,IK)
IF(X(J).LE.0.)X(J)=1.E-5
END IF
3201
CONTINUE
3200
CONTINUE
IF(ICE.EQ.1) THEN
SU(IV)=-G(IV)
DO 3300 K=1,IP
SU(IV)=SU(IV)+A(K,IV)*G(K)
3300
CONTINUE
III=0
END IF
11
DO 3400 J=IP1,NMV
EX(J)=X(J)/X(IV)
3400
CONTINUE
III=III+1
CALL PITZER (T,2,EX)
IF(ICE.EQ.1.AND.IPATH.EQ.1) THEN
IF(AH2O.GT.SU(IV))THEN
X(IV)=X(IV)/1.2
GO TO 11
ELSE
C
IF(III.GT.1)X(IV)=X(IV)*1.1
END IF
END IF
DO 3500 I=1,ICLOS
J=JCLOS(I)
NJ=NUM(J)
IF(NJ.LT.NL) G(J)=G(NJ)+DLOG(EX(NJ))+ACT(NJ)
3500
CONTINUE
32
Previous Page
