18
*’ Введите количество точек по температуре: ’
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)
IF(NVAR.EQ.2) CALL GERG2(ICALC,YA,YY)
IF(NVAR.EQ.3) CALL AGA8DC(ICALC)
IF(NVAR.EQ.4) CALL VNIC(ICALC)
IF(Z.NE.0D0) NTS=NTS+1
ZP(I,J)=Z/ZN
7 CONTINUE
IF(NTS.EQ.0) THEN
CALL RANGE(NRANGE)
IF (NRANGE) 134,134,200
ELSE
I=1
9 IС=0
DO 11 J=1,NT
IF(ZP(I,J).EQ.0D0)
IC=IC+1