35
DATA GD/0D0,.0793D0,.141239D0,.281835D0,.256692D0,.027815D0,
*.189065D0,.0885D0,.366911D0,.332267D0,.432254D0,.512507D0,
*.576242D0,.648601D0,.716574D0,0D0,.034369D0,.038953D0,.021D0/
DATA QD/6*0D0,.69D0,12*0D0/,FD/16*0D0,1D0,2*0D0/
DATA AN/.1538326D0,1.341953D0,-2.998583D0,-.04831228D0,
*.3757965D0,-1.589575D0,-.05358847D0,2.29129D-9,1576724D0,
*-.4363864D0,-.04408159D0,-.003433888D0,.03205905D0,.02487355D0, -
*.07332279D0,-.001600573D0,.6424706D0,-.4162601D0,-.06689957D0,
*.2791795D0,-.6966051D0,-.002860589D0,-.008098836D0,3.150547D0,
*.007224479D0,-.7057529D0,.5349792D0,-.07931491D0-1.418465D0,
*-5.99905D-17,.1058402D0,.03431729D0,-.007022847D0,.02495587D0,
*.04296818D0,.7465453D0,-.2919613D0,7.294616D0,-9.936757D0,
*-.005399808D0,-.2432567D0,.04987016D0,.003733797D0,1.874951D0,
*.002168144D0,-.6587164D0,.000205518D0,.009776195D0,-.02048708D0,
*.01557322D0,.006862415D0,-.001226752D0,.002850906D0/
DATA BN/13*1,9*2,10*3,7*4,5*5,2*6,2*7,3*8,2*9/
DATA CN/7*0,6*1,2*0,7*1,0,9*1,2*0,5*1,0,4*1,0,1,0,6*1/
DATA KN/7*0,3,3*2,2*4,2*0,3*2,4*4,0,2*1,2*2,2*3,3*4,2*0,3*2,
*2*4,0,2*2,2*4,0,2,0,2,1,4*2/
DATA UN/0D0,.5D0,1D0,3.5D0,-.5D0,4.5D0,.5D0,-6D0,2D0,3D0,2*2D0,
*11D0,-.5D0,.5D0,0D0,4D0,6D0,21D0,23D0,22D0,-1D0,-.5D0,7D0,-1D0,
*6D0,4D0,1D0,9D0,-13D0,21D0,8D0,-.5D0,0D0,2D0,7D0,9D0,22D0,23D0,
*1D0,9D0,3D0,8D0,23D0,1.5D0,5D0,-.5D0,4D0,7D0,3D0,0D0,1D0,0D0/
DATA GN/4*0,2*1,13*0,1,3*0,1,2*0,3*1,16*0,1,2*0,1,0,1,2*0/
DATA QN/6*0,1,3*0,1,9*0,1,0,1,8*0,1,4*0,1,4*0,1,0,1,2*0,1,5*0,1/
DATA FN/7*0,1,13*0,1,2*0,1,4*0,1,23*0/
END
C ***********************************************************
С * *
С * Подпрограмма расчета коэффициента сжимаемости природного *
С * газа по уравнению состояния ВНИЦ СМВ. *
С * *
C ***********************************************************
SUBROUTINE VNIC(ICALC)
IMPLICIT REAL*8(A-H,O-Z)
REAL*8 LIJ(8,8)
DIMENSION VC(8),TC(8),PII(8),DIJ(8,8)
COMMON/PARCD/VCD(8),TCD(8),PIID(8)/ABIJ/AIJ(10,8),BIJ(10,8)
*/B/B(10,8)/RM/RM/Y/Y(8)/BM/BM(8)/NI/NI(8)/NC/NC/Z/Z
RM=8.31451DO
IF(ICALC.NE.1) GO TO 1
CALL COMPON
IF(Z.EQ.0D0) GO TO 133
CALL DDIJ(DIJ,LIJ)
DO 75 I=1,NC
TC(I)=TCD(NI(I))
VC(I)=BM(I)/VCD(NI(I))
PII(I)=PIID(NI(I))
DO 123 J=1,NC
IF(I.GE.J) GO TO 123
DIJ(I,J)=DIJ(NI(I),NI(J))
LIJ(I,J)=LIJ(NI(I),NI(J))