DO 203 K=I,NP1
AT=P1(I1,K)
P1(I1,K)=P1(I,K)
P1(I,K)=AT/AM
203
CONTINUE
DO 241 J=INI,NP
AM=P1(J,I)
IF(J.EQ.I) GO TO 241
DO 204 K=I,NP1
P1(J,K)=P1(J,K)-P1(I,K)*AM
204
CONTINUE
241
CONTINUE
201
CONTINUE
RETURN
END
C------------------------------------------------------------------C
C------------------------------------------------------------------C
C..
THE ROUTINE FOR CALCULATION OF CHEMICAL EQUILIBRIA
C..
IN WATER-SALT SYSTEMS AT THE TEMPERATURE RANGE 25- -60 C.
DOUBLE PRECISION FUNCTION PF(T,A1,A2,A6,A9,A3,A4)
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
PF=A1+A2*T+A6*T**2+A9*T**3+A3/T+A4*LOG(T)
RETURN
END
C------------------------------------------------------------------C
C------------------------------------------------------------------C
SUBROUTINE SOL (ICK,EX,ISOLU,IPATH)
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
CHARACTER *40 NM
DIMENSION G(40),X(40),EX(20),EXTMP(20),A(10,40),ACT(20),NM(40)
DIMENSION JOPEN(10),JCLOS(10),SU(20),P1(20,21),AIN(10,40)
DIMENSION NUMIN(40),NUM(40),GIN(40),ACTT(20),IZ(20)
COMMON /NUMBER/NCAT,NANI,NNEI,NMM(40)
COMMON /MASS1/ A,G,NUM,X
COMMON /MASS2/ NM
COMMON /INTGER/ NL,N,N1
COMMON /MATRIX/ ACT,UM,PHI,AH2O,IZ,IP,NMV
COMMON /AIN/ AIN,GIN,T,NUMIN,IV
COMMON /COMGG/ P1
COMMON /ENTER/ ISIM,ISIMIN
COMMON /OUTPUT/JOPEN,IOPEN
DO 1000 I = 1, 20
ACT(I) = 0.D0
1000
CONTINUE
29
Previous Page
