NM(I1)=NMC(I)
NMM(I1)=NUMC(I)
NUMIN(I1)=I1
2500
CONTINUE
IP1=IP+1
L=NCAT+NANI+NNEI
N=NL+NSOL
N1=N+1
DO 2600 J=1,IP
AIN(J,N1)=BAL(J)
2600
CONTINUE
IJJ=0
DO 2700 I=1, N1
NUM(I)=NUMIN(I)
2700
CONTINUE
C...........................
READ TABLE 3
.......................
OPEN(5, FILE='TABLE3')
DO 2800 I=IP+1,N
2810
READ(5,*, END=2820) NT1
IF(NMM(I).EQ.NT1) THEN
BACKSPACE(5)
READ(5,*) NT1, (G0(I,K),K=1,6)
REWIND(5)
ELSE
GO TO 2810
END IF
2820
REWIND(5)
2800
CONTINUE
CLOSE (5)
CALL PITZER (T,0,EX)
READ(9,*) TINIT
READ(9,*) IPATH
IF (IPATH.EQ.1) THEN
READ(9,*) TFIN
READ(9,*) DT
ELSE
READ(9,*) WFIN
READ(9,*) DW
END IF
T=TINIT
BALWAT=1000.
1
CONTINUE
DO 2900 I=IP+1,N
GIN(I)= PF(T,G0(I,1),G0(I,2),G0(I,3),G0(I,4),G0(I,5),
$
G0(I,6))
2900 CONTINUE
14
Previous Page
