21
DO 1 J=l,8
1 BK(I)=BK(I)+B(I,J)/TR**(J-1)
7 ROR=RO*VCM
AO=0.D0
A1=0.D0
IF(NPRIZ.EQ.l) GO TO 5
A2=0.D0
A3=0.D0
5 DO 33 I=1,10
D=BK(I)*ROR**I
AO=AO+D
A1=A1+(I+1)*D
IF(NPRIZ.EQ.1) GO TO 33
DO 3 J=1,8
D1=B(I,J)*ROR**I/TR**(J-1)
A2=A2+(2-J)*D1
3 A3=A3+(J-1)*(2-J)*D1/I
33 CONTINUE
RETURN
END
С Подпрограмма расчета плотности, показателя адиабаты, скорости
С звука
SUBROUTINE TP(ROM)
IMPLICIT REAL*8(A-H,O-Z)
COMMON/BMM/BMM/AI/AO,A1,A2,A3/RM/RM/T/T/TS/RO,PA,W/Z/Z
CALL IDGFU(T,CVOS)
RO=BMM*ROM
R=RM/BMM
A11=1.D0+A1
A21=1.D0+A2
CV=R*(A3+CVOS)
CP=CV+R*A21**2/A11
W=DSQRT(DABS(1.DЗ*R*T*CP/CV))*DSQRT(DABS(A11))
PA=CP/CV*A11/Z
RETURN
END
С Подпрограмма расчета изохорной теплоемкости в идеально газовом
С состоянии
SUBROUTINE IDGFU(T,CVOS)
IMPLICIT REAL*8(A-H,O-Z)
DIMENSION CPO(8),CVO(8)
COMMON/IDGF/CPC(20,8),TOI(8),MCO(8),MCP(8)/Y/Y(8)/NC/NC
CVOS=0.D0
DO 21 I=1,NC
M=MCP(I)
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