22
N=MCO(I)
TAU=T/TOI(I)
S1=0.D0
S2=0.D0
S3=0.D0
S1=CPC(1,I)
IF(M.EQ.0) GO TO 7
DO 9 J=1,M
9 S2=S2+СРС(J+1,I)*ТАU**J
7 IF(N.EQ.0) GO TO 11
DO 13 J=1,N
13 S3=S3+CPC(M+J+1,I)/TAU**J
11 CPO(I)=S1+S2+S3
CVO(I)=CPO(I)-1.D0
21 CVOS=CVOS+Y(I)*CVO(I)
RETURN
END
С Подпрограмма расчета вязкости
SUBROUTINE ETAS(ROM)
IMPLICIT REAL*8(A-H,O-Z)
COMMON/ETA/ETA/PARCM/TCM,VCM/BMM/BMM/T/T/PIM/PIM/PCM/PCM
DKSI=TCM*(1D0/6D0)/BMM**.5/PCM**(2D0/3D0)
ROR=VCM*ROM
TR=T/TCM
ETA=78.037D0+3.85612*PIM-29.0053*PIM**2-156.728/TR+145.519/TR**2
*-51.1082/TR**3+6.57895*ROR+(11.7452D0-95.7215*PIM**2/TR)*ROR**2+
*17.1027*ROR**3*PIM+.519623/TR**2*ROR**5
ETA=ETA/DKSI/10.
RETURN
END
BLOCK DATA BDVNIC
IMPLICIT REAL*8(A-H,O-Z)
CHARACTER*26 AR
COMMON/PARCD/VCD(8),TCD(8),PIID(8)/ABIJ/AIJ(10,8),BIJ(10,8)
COMMON/CPCI/CPC1(20,5),CPC2(20,3)/IDGFD/TOID(8),MCOD(8),MCPD(8)
*/AR/AR(25)
DATA TCD/190.67D0,305.57D0,369.96D0,425.4D0,407.96D0,
*125.65D0,304.11D0,373.18D0/
DATA VCD/163.03D0,205.53D0,218.54D0,226.69D0,225.64D0,
*315.36D0,466.74D0,349.37D0/
DATA PIID/0.0006467D0,0.1103D0,0.1764D0,0.2213D0,0.2162D0,
*0.04185D0,0.2203D0,0.042686D0/
DATA AIJ/.6087766D0,-.4596885D0,1.14934D0,-.607501D0,
*-.894094D0,1.144404D0,-.34579D0,-.1235682D0,.1098875D0,
*-.219306D-1,-1.832916D0,4.175759D0,-9.404549D0,10.62713D0,
*-3.080591D0,-2.122525D0,1.781466D0,-.4303578D0,-.4963321D-1,
*.347496D-1,1.317145D0,-10.73657D0,23.95808D0,-31.47929D0,
*18.42846D0,-4.092685D0,-.1906595D0,.4015072D0,-.1016264D0,
*-.9129047D-2,-2.837908D0,15.34274D0,-27.71885D0,35.11413D0,
*-23.485D0,7.767802D0,-1.677977D0,.3157961D0,.4008579D-2,0.D0,
*2.606878D0,-11.0б722D0,12.79987D0,-12.11554D0,7.580666D0,
*-1.894086D0,4*0.D0,