. " -3 1" ( )

.. "122 1944 ." ( )

. " " ( )

.. " " ( )

.. " " ( )
- ..
.., .., .., .. . .. .: , 1982. 267 c.
( ): primenenevm1982.djvu
<< 1 .. 91 92 93 94 95 96 < 97 > 98 99 100 101 102 .. 103 >>


DATA RO /1.98719/ =1.0
IF(ID(D.EQ.0) =0.239006
DO 99 1=1,100
IF (ABS(F(I)-.77777)
* .GT.1E-4) F(I)*F(I)*A IF (ABS(CP(I)-.77777)
* .GT.1E-4) CP(I)=CP(I)*A IF (ABS(H(I)-.77777)
* .GT.1E-4)H(I)=H(I)*A IF (ABS(S(I)-.77777)
* .GT.1E-4) S(I)=S(I)*A 99 CONTINUE
M1=0
DO 98 1=5,12
IF (ABS(R(I)-.77777)
* .GT.1E-4) MI=I 98 CONTINUE
IF (ABS(Ed )-.77777)
* .GT.1E-4) M1=1
IF (ABS(OMd )-.77777)
* .GT.1E-4) M1=1
IF (ABS(F(D-.77777)
* .GT.1E-4) M1=1
IF (ABS(CP(I)-.77777)
* .GT.1E-4) M1=1
IF (ABS(H(1)-.77777)
* .GT.1E-4) MI=I
IF (ABS(Sd )-.77777)
* .GT.1E-4) M1=1 NP=O
bo 97 1=1,5
IF (CO(I9I)-.77777 .GT.
* 1E-4) NP=NP+1 DO 97 J=3,9
IF (ABS(CO(I1J)-.77777)
* .GT.1E-4)
* CO(I,J)=CO(I,J)*A 97 CONTINUE
IF (ABS(RO)-.77777)
* .GT.1E-4) R(3)=R(3)*A IF (ABS(RU)-.77777)
* .GT.1E-4) R(4)=R(4)*A IF (ABS(R(13)-.77777)
* .GT.1E-4) R(13)=R(13)*A IF (ABS(15)-.77777)
* .GT.1E-4) R(15)-R(15)*A IF (ABS(R(16)-.77777)
* .GT.1E-4) R(16)=R(16)*A DO 96 1=19,25,2
IF (ABS(R(I)-.77777)
* .GT.1E-4) R(I)=R(I)*A
243
96 C(XITINUE
CALL MMASS (K1,K2,R(2)) IF (M1.EQ.1.0R.NP.GT.O)
* QO TO 1 PRINT 106 QO TO 999
C***
1 IF (T(1)-.77777.GT.1E-4)
* QO TO 3 T(1)=298.15 T(2)=250. DO 2 1=3,81
2 T(I)=T(I-I)+250. N=81
QO TO 13
C***
3 IF (T(2)-.77777.GT.1E-4)
* QO TO 4 PRINT 100 QO TO 999
0***
4 IF (ABS(F(I)-.77777)
* .GT.1E-4.0R.ABS(S(1)-
* .77777).GT.1E-4.GR.
ABS(H(1)-.77777).GT.1E-4
* .QR.CP(1)-.77777.GT.
* 1E-4) QO TO 10
IF (T(3)-.77777.GT.1E-4)
* GO TO 8
5 DT1=(T(2)-T(1))/98. DT=(T(2)-?(1))/6. N=8
IF (DT.LE.100.) GO TO 6 N=(DT1+100.)/100. OT=N*100
N=(T(2)-T(1))/DT+.5 N=N+2
6 DT1=T(2) T(2)=T<1) T(I)-298.15
, DO 7 I*3,N
7 T(I)=T(I-I )+W T(N)=DTI
GO TO 13
8 IF (T(4)-.77777.GT.1E-4)
* GO TO 9 N=(T(2)-tf(1))/T(3)+2. IF (N.LE.2) GO TO 9 IF (N.GT.100) GO TO 5 DT=T(3)
GO TO 6
9 N=I DT=298.15
DO 91 1*1,100
IF (T(I)-.77777.GT.1E-4)
* N=N+1 DTI=T(I) T(I)=DT
91 DT=DTI GO TO 13
C**
10 N=O
DO 11 1=1,100
IF (T(I)-.77777.GT.1E-4)
* N=N+1
11 CONTINUE
IF (ABS(F(N)-.77777)
* .GT. 1 E-4 .OR-ABS(S(N)-
* .77777).GT.1E-4.0R.
* ABS(H(N)-.77777).GT.1E-4
* .OR.ABS(CP(N)-.77777)
* .GT.1E-4) QO TO 12 PRINT 101
GO TO 999
12 IF (ABS(F(N+!)-.77777)
* .LE. 1E-4.AND.ABS(S(N+!)-
* .77777).LE. 1B-4.AND.
* ABS(H(N+1)-.77777).LE.
* 1E-4. AND . ABS(CP (N+D-
* .77777).LE.1 E-4)
* GO TO 600 PRINT 102 GO TO 999
?***
13 IF (Ml .EQ.O) QO TO 610
N1=0 N2=0
DO 14 1=1,500
IF (ABS(E(I)-.77777)
* .GT.1E-4) N1=N1+1
IF (G(I)-.77777.GT.1E-4)
* N2=N2+1
14 CONTINUE
IF (N1.EQ.N2) GO TO 15 PRINT 103 GO TO 999
15 IF (R(5)-.77777.GT.1E-4)
* GO TO 19 J1 = 1
W=1.5*AL0G(R(2))
16 DO 18 1=1,N DT*-11604.5/T(I) A=0.
B=0. C=0.
DO 17 J=I,N1
D=G(J)*EXP(DT*E(J))
A=A+D
D=D*E(J)
B=B+D
17 C=D*E(J)+C WI=B/A
IF (WI .LT. IE-20) W1=0* CPU>=( 2 . 5+DT*DT ( C/A-
Wf*2))*R0
H(IW 2 . 5-DTW1 )*T( I)*R0 S(I)=(?*2 . 5AL0G(T( I) ) +
ALOG(A) -DTW1)R*-
2.31482
18 f)sS(I)-H(I)/TCl*
m 500
1ST
DO ߻ 1=1,,100
?8(011(1)-.77777) * .GVfB-4) N2=N3*1? 20
IF (ABS(R( 10)-.77777)
* .GT.1E-4) GO TO 23 IF(ABS(R(1 1)-.77777)
* .GT.1E-4) GO TO 30 IF (ABS(R(9)-.77777)
* .LE.1E-4) GO TO 21 K3=K1/100000000 K4=K1-K3*100000000 M1=(K4+1)/1000000 K4=K4-M1*1000000 M2=(K4+1)/10000
IF (M2.EQ.0) M2=M1 R(8)=AM(M1)*AM(M2)/
* (AM(MI )+AM(M2))
* (R(9)*2)*1 .6596
21 IF (ABS(R(8)-.77777)
* .LE.1E-4) GO TO 22 R(7)=27.9932/R(8)
22 IF (ABS{R(7)-.77777)
* .LE.1B-4) GO TO 23 R(6)=1.43879*R(7)
23 J2=1
24 DO 29 I=,K
IF (ABS(R( 10)-.77777)
* .GT.1E-4) GO TO 25 DT=R(6)/T(I)
DT=(.2DT+1.)*DT/3.+1.0 CP(I)=(2.0-1.0/DT)R0/DT F( I)=ALOGK T (I) DT/
* (R(5)R(6)))*R0 S(I)F(I)+R0/DT GO TO 26
25 CP(I)=I.5*R0 DT=SQHT(T(I)**3*R( 10)/
* 20796.8)/R(5) F(I)=R0*ALOG(DT) S(I)=F(I)+CP(I)
26 DT=R0*( 1 .5*ALOG(R( 2) ) +
* 2.5ALOG(T(I))) F(I)=DT+F(I)-7.28279 S(I)=DT+S(I)-2.31482 CP(I)=CP(I)+2.5*R0 DO 27 J=1,N2
245
A=I.43879OM(J)A(I) B=EXP (-A) W=I.0-B WI=ALOG(W)
CP(I )=RO*B*(A/W)** 2+CP(I) F(I)=F(I)-R0*W1
27 S(I)=R0*(A*B/W-W1)+S(I) A=O.
B=O. C=O.
DT=-11604.5/T(I)
DO 28 J=I,N1
D=G(J)*EXP(DT*E(J))
A=A+D
D=D*E(J)
B=B+D
28 C=D*E(J)+C D=B/A
W1=R0*ALOG(A) F(I)=F(I)+W1 S(I)=S(I)+W1-R0*DT*D CP(I)=RO*DT*DT*(C/A-D**2) - +CP(I)
29 H(I)=(S(I)-F(I))*T(I) GO TO 500
C****
30 J3=1
IF (ABS(R(12)-.77777)
* .LT.1E-4) R(12)=R(7)**2*
* 5.0*(SQRT(R(11)/R(7))-
* 1.0)/OM(1) BB=R(7)-R(12)*0.5 BM=1.5*ALOG(R(2))
31 DO 33 1=1,N A=BB/T(I) C=A*A
CP(I)=O.5+.046C)*R0 H(I)=(3.5+.47958*A-
* .046*C)*T(I)*R0 S(I)-(BM+2.5ALOG(T(I))-
* ALOG(A*R(5))-0.023*C-
* 0.528733)*R0 =1 .43879()
B=A*(0M(1)-2.*R(11)) C=EXP(-B) DT=LO-C D=BC/DT
CP(I)=D*E/DT*RO +CP(I)
H(I)=RO^T(I)*D+H(I)
S(I)=(D-ALOG(DT))#RO+S(I)
C=EXP(B)
DT=R(I1)/0M(1)
D=R(12)/R(7)
DT1=R(7)/0M(1)
A1=8.0*DT1/B
A2=C-1 .0
A3=B*C
A4=A3*B/A2**2 A5=4.0*A4*DT/A2 CP(I) =( A1+A1-A4*D+
(2.0*A4*(D*C-2.0DT*(B+
2.0) )+3.0A5*A3)/A2)*R0+
CP(I) H(I)=(A1+B*(D*C-2.0*DT)
<< 1 .. 91 92 93 94 95 96 < 97 > 98 99 100 101 102 .. 103 >>
 
 
2010 PiroChem. .