END IF
2400
CONTINUE
CALL PITZER (T,2,EX)
G(IVV)=AH2O
DO 2500 I=1,IP
NI=NUM(I)
IF (I.NE.IVV) THEN
G(I)=G(NI)+DLOG(EX(NI))+ACT(NI)
END IF
2500
CONTINUE
DO 2600 J=IP1,NMV
SU(J)=-G(J)
DO 2601 K=1,IP
SU(J)=SU(J)+A(K,J)*G(K)
2601
CONTINUE
EX(J)=DEXP(SU(J)-ACT(J))
X(J)=X(IV)*EX(J)
2600
CONTINUE
CALL SIMPL (NL1,NL1,1,IP,N1,IK,IR)
DELL=1.
75
CONTINUE
76
ICLOS=0
IOPEN=0
DO 2700 I=1,IP
IF (NUM(I).GE.NL1.AND.NUM(I).LE.N) THEN
IOPEN=IOPEN+1
JOPEN(IOPEN)=I
ELSE
ICLOS=ICLOS+1
JCLOS(ICLOS)=I
END IF
2700
CONTINUE
IF (IOPEN.GT.0) THEN
IF (IOPEN.EQ.1.AND.NMM(NUM(JOPEN(1))).EQ.31) THEN
ICE=1
SU(IV)=-G(IV)
DO 2800 K=1,IP
SU(IV)=SU(IV)+A(K,IV)*G(K)
2800
CONTINUE
XMN=1.
DO 2900 K=1,10
XMN=XMN*1.5
DO 2901 J=IP1,NMV
EXTMP(J)=EX(J)*XMN
SUMS=SUM*XMN
2901
CONTINUE
31
Previous Page
