. " -3 1" ( )

.. "122 1944 ." ( )

. " " ( )

.. " " ( )

.. " " ( )
- ..
.., .., .., .. . .. .: , 1982. 267 c.
( ): primenenevm1982.djvu
<< 1 .. 86 87 88 89 90 91 < 92 > 93 94 95 96 97 98 .. 103 >>

* 108,11)
105 FORMAT (2,1,4 1,7,
* 4,5,7(91,4),5 I ,
* 1)
106 FORMAT (2,1,4 I,
* 108X,4HI ,1)
225
1/2 15 .1554
107 K)RMAT (2X95H. 1,48X9
* 20195911 NO. ,
* 20A1 94X95HI ./691.9
* 108(1H-),1./ .,116X9
* 1./2,39(. )91.)
108 FORMAT (6X91HIt86X9
* 1192999
* IHI)
109 FORMAT (2X9AI94H I918X9
* 1 9HCOCTAB ,
* 2 [97A4,41X94HI ,Al)
110 FORMAT (2X9AI ,4H I9 4X9
* 12197(9194)95 I ,
* Al)
111 IORMAT (2X9AI94 1,18,
* 2OHOT ,
* 13 [,74,
* 29X94HI 91)
112 FORMAT (2XftA194H 1,18,
* 20HCOCTAb 1- ,
* 20ͳ [9
* 7A4922X,4HI 91)
113 FORMAT (29194 1,18X9
* 2OHCCJCTAB 2- ,
* 20? [,
* 7A4922X,4HI 91) END
SUBROUTINE SNDVEL
* (SV1CPP9CVP)
DOUBLE PRECISION AC9RZ9T9
* A9STR9X9EX1S9H9CP9SM9TV9
* EXPi9ANJ9EXMI 9EXM29TZ9
* SM19XE
COMMON /GEN/KDR( 1I)9
* NC(15)9AC(10)/REZ/
* RZ(100)9T(100)//A(1009
* 101)9STR(101)9X(350)9
* XE(350)9NJI(21)/ID/
* DNUO) 9K(7)9P1(1 O)9
* P2(10)9EL(20)9NO(20)9
* MX(80)/BIBVF(3509 9)9
* NM(350,4)/PR/S(350),
* H(350)fCP(350)/LR/
* G(2797)9TS(700)
.DIMENSION TZ(29100)9
* TV(200),EX(350)
*****#*****#********
.. ******************************
DO 99 1=1,350 99 EX(I)=XE(I)
DO 1 1=1,100
TV(I)=RZ(I)
1 TV(I+100)=T(I) L1=NC(4)+K(5) L3=NC(3) EXMI=O.
IM1=L1+NC(7)+7 IF (NC(5).NE.O)
* EXM1=DEXP(TV(IM1)) EXM2=0.
IF (NC(6).NE.O)
* EXM2=DEXP(TV(IM1+1)) IK=LI+7
IKS=O
DO 2 1=19L3
IF (NM(I93).EQ.5)
* EX(I)=EX(I)*EXM1 IF (NM(I,3).EQ.6)
* EX(I)=EX(I)EXM2
IF (NM(I,3).NE.1) GO TO 2
EX(I)=DEXP(TV(IK))
IK= IK+1
II (EX(I) .GT.1E-10)
* IKS=IKS+!
2 CONTINUE L2=L1+IKS M=L2+1
IF (EXM1.GT.1E-10) M=M+1 IF (EXM2.GT.1E-10) M=M+1 M1=M+1 L4=K(5)
226
IKK=LI+2
NR*NC( 9)
SMI=O.
DO 3 I=I,L3
II (NM(I,3).EQ.1.OR.
* NM(I,3).GE.5) GO TO 3 SMI=SMI+EX(I)
3 CONTINUE
DO 23 11=1 ,2 SM= (2-ID/SMI DO 4 I=I1M DO 4 J=I1MI
4 A(I1J)=O. A(I1I)=-1.
II (EXMI.GT.1E-10)
* A(L2+2,L2+2)=-1.0 II (EXM2.GT.1E-10)
* A(M,M)=-1.0 IK=IKK
DO 21 1=1,L3 EXPI=EX(I) JI=NM(I,3)
CALL STRIKS (IfTV(2)) II (J1.EQ.1.CR.J1.GE.5)
* GO TO 6 DO b J=I ,L1
5 A(1,J+1)=A(1,J+1)+
* STR(J+3)*EXPI A(1,1)=A(1,1)+EXPI*SM A(I,MI)=A(I,M1)+((II-1)/
* TV(2)-STR(1))*EXPI
6 DO 14 J2=1,L4
IF (NJI(J2).EQ.0)
* GO TO 14 JA=J2+1 ANJ=NJK J2)
GO TO (7,8,8,8,10,13), Jl
7 II (EXPI.LT.1E-10)
* GO TO 14 A(JA,IK)=A(JA,IK)+ANJ GO TO 14
8 DO 9 J=I ,Ll
9 A(JA,J+1)=A(JA,J+1) +
* STR(J+3)*EXPI*ANJ A(JA,I)=A(JA1I)+
* EXPI *SM ANJ
A(JA1MI )=A(JA,Ml )+((II-1)
* /TV(2)-STR(1))*EXPI*ANJ GO TO 14
10 II (EXM1.LT.IE-10)
* GO TO 14 A(JA,L2+2)=A(JA,L2+2) +
* EXPI*ANJ/EXM1
11 DO 12 J=1,L4
12 A(JA,J+1)=A(JA,J+1)+
* STR(J+3)*EXPI*ANJ A(JA, MI)=A(JA, Ml)-
* STR(1)*EXPI*ANJ GO TO 14
13 II (EXM2.LT.1E-10)
* GO TO 14
A( JA1M) =A(JA,M) +
* EXPI*ANJ/EXM2 GO TO 11
14 CONTINUE
II (NC(4).EQ.O) GO TO 15 II (NJl(LI) .EQ.0)
* GO TO 15 JAsL4+2 ANJ=NJI(LI) DO 141 J=I,Ll
141 A(JA,J+1)=A(JA,J+1) +
* STR(J+3)*EXPI*ANJ A(JA,I)=A(JA11)+
* EXPI *SM*ANJ
A(JA,M1 )=A( JA,M1 )+( (II-1 )
* AV(2)-STR(1) )*EXPI*ANJ
15 II (JI .NE 1 .OREX(I)
* .LE.1E-10) GO TO 17 DO 16 J=1,L4
16 A(IK,J+1)=A(IK,J+1)+
* NJI(J)
A( IK , Ml) =A( IK, Ml) -H( I)/
* TV(2)**2 IK= IK+1 GO TO 21
227
17 I* (J1.NE.5.OR.EXMI .LE,
* 1E-10) GO TO 20 JA= L2+2
A(JA1JA)=A(JA1JA) +
* EXPI/EXM1
18 A(JA1MI )=A( JA,Ml)-STR( 1 )*
* EXPI
DO 19 J=I,L4 1 9 A(JA,J+1 )=A( JA1J+1 ) +
* STR(J+3)*EXPI GO TO 21
20 II (J1.NE.6.GR.EXM2.LE.
* 1-) GO TO 21 JA=M
A(JA,JA)=A( JA1JA) +
* EXPI/E XM2 GO TO 18
21 CONTINUE
CALL SLEQ (M,M) DO 22 1=1,100
22 TZ(II,I)=T(I)
23 CONTINUE CPP=G(8,NR)*1000 II (KDR(2).EQ.0)
* CPP=CPP*0.239005 CVP=CPP/G(9,NR) IK=IKK
DO 30 1=1,L3 EXPI=EX(I)
CALL STRIKS (I,TV(2)) JI=NM(I,3)
II (J1.NE.1) GO TO 24 II (EXPI.LE.1E-10)
* GO TO 30 CPP=TZ(I ,IK)*4(I)+CPP CVP=TZ (2, IK )* H( I)+CVP IK=IK+1
GO TO 30
24 IF (J1.NE.5) GO TO 27 II (EXMI.LE.1E-10)
* GO TO 30 R1=(STR(1)+TZ(1,L2+2)/
* EXM1)*EXPI
R2=(STR(1)+TZ(2,L2+2)/
* EXM1)*EXPI 2 5 DO 26 J=I,L4
R1=R1+STR(J+3)*EXPI*
* TZ(1,J+1)
26 R2=R2+STR(J+3)*EXPI*
* TZ(2,J+1) CPP=R1*H(I)+CPP CVP=R2*H(I)+CVP GO TO 30
2 7 II (-J1.NE.6) GO TO 28 II (EXM2.LE.1E-10)
* GO TO 30 R1=(STR(1)+TZ(1,M)/
* EXM2)*EXPI R2=(STR(1)+TZ(2,M)/
* EXM2)*EXPI GO TO 25
28 R1=EXPI*(TZ(1,1)/SM1+
* STR(I))
R2=STR( 1 )*EXPI-EXPI/T V(2)
DO 29 J=1,L1
RI=RI +STR ( J+3) EXPI*
* TZ(1,J+1)
29 R2=R2+STR(J+3)*EXPI*
* TZ(2,J+1) CPP=RI*H(I)+CPP CVP=R2*(H(I)-1.98719*
* TV(2))+CVP
30 CONTINUE
DO 31 1=1,100 RZ(I)=TV(I)
31 T(I)=TV(I+100)
SV= (8.31 44* (CEP/CVP )*
* SM1*TV(2)*(SM1+TV(2)*
* TZ(2,1))/(SM1+TV(2)*
* TZ(1,1))) SV=SQRT(SV) RETURN
END
228
C PROGRAM MDCLR
DIMENSION H(75,10),
* K(75,2) ******************************


<< 1 .. 86 87 88 89 90 91 < 92 > 93 94 95 96 97 98 .. 103 >>
 
 
2010 PiroChem. .