22
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,
*-1.15575D0,3.601316D0,-.7326041D0,-1.151685D0,.5403439D0,
*5*0.D0,.9060572D-1,-.5151915D0,.7622076D-1,7*0.D0,
*.4507142D-1,9*0.D0/
DATA BIJ/-.7187864D0,10.67179D0,-25.7687D0,17.13395D0,
*16.17303D0,-24.38953D0,7.156029D0,3.350294D0,-2.806204D0,
*.5728541D0,6.057018D0,-79.47685D0,216.7887D0,-244.732D0,
*78.04753D0,48.70601D0,-41.92715D0,10.00706D0,1.237872D0,