DO 3600 J=IP1,NMV
SU(J)=-G(J)
DO 3601 K=1,IP
SU(J)=SU(J)+A(K,J)*G(K)
3601
CONTINUE
IF(J.NE.IV) EX(J)=DEXP(SU(J)-ACT(J))
X(J)=EX(J)*X(IV)
3600
CONTINUE
END IF
END IF
IN=ICLOS+L
IN1=IN+1
5
IF(ICK.GT.400) RETURN
ICK=ICK+1
C..
EXAMINE FOR AQUEOUS SOLUTION PRESENCE.............
IF(ICLOS.LE.1) THEN
ISOLU=0
DO 3700 I=IP1,NL
X(I)=0.
3700
CONTINUE
GO TO 400
END IF
ISL=0
DO 3800 J=IP1,NL
SU(J)=G(J)
DO 3801 K=1,IP
SU(J)=SU(J)-A(K,J)*G(K)
3801
CONTINUE
EX(J)=X(J)/X(IV)
ACTT(J)=ACT(J)
IF(IS.EQ.1)ACTT(J)=0.
3800
CONTINUE
IS=IS+1
PHITT=PHI
CALL PITZER (T,2,EX)
PHI=(PHITT+PHI)/2
SUM=0.
DO 3900 J=IP1,NMV
ACT(J)=(ACT(J)+ACTT(J))/2.
JI=J-IP
P1(JI,JI)=1./X(J)
33
Previous Page
