18
WRITE(*,35)
35 FORMAT(’ Введите 0, если состав задан в молярных долях’/
*’ или 1, если состав задан в объемных долях ’\)
READ(*,*)NPR
IF(NPR.EQ.0) WRITE(*,3)
IF(NPR.EQ.l) WRITE(*,33)
DO 5 I=1,25
WRITE(*,’(A\)’) AR(I)
READ(*,*)YC(I)
5 YC(I) = YC(I)/100.
ENDIF
WRITE(*,’(A\)’)
*’ Введите количество точек по давлению: ’
READ(*,*)NP
WRITE(*,’(A\)’)
*’ Введите количество точек по температуре: ’
READ(*,*)NT
WRITE(*,’(A\)’)
*’ Введите значения давлений в МПа: ’
READ(*,*)(PI(I),I=1,NP)
WRITE(*,’(A\)’)
*’ Введите значения температур в К: ’
READ(*,*)(TI(I),I=1,NT)
WRITE(*,’(A\)’)
*’ Ввод исходных данных завершен. ’
P=.101325D0
T=293.15D0
ICALC=1
GO TO (10,20,30,40) NVAR
10 CALL NX19(YA,YY)
ZN=Z
GO TO 50
20 CALL GERG2(ICALC,YA,YY)
ZN=Z
GO TO 50
30 CALL AGA8DC(ICALC)
ZN=Z
GO TO 50
40 CALL VNIC(ICALC)
ZN=Z
50 CONTINUE
IF(Z.EQ.0D0) THEN
CALL RANGE(NRANGE)
IF(NRANGE) 134,134,200
ENDIF
ICALC=2
NTS=0
DO 7 I=1,NP
P=PI(I)
D07 J=1,NT
T=TI(J)
IF(NVAR.EQ.l) CALL NX19(YA,YY)