9
READ(*,*)NPR
IF(NPR.EQ.1)THEN
WRITE(*,’(A\)’)
*’ Плотность при 293.15 К и 101.325 кПа, в кг/куб.м’
READ(*,*)RON
WRITE(*,33)
33 FORMAT(’ Значение объемной доли, в об.%’)
ELSE
RON=0D0
WRITE(*,3)
3 FORMAT(’ Значение молярной доли, в мол.%’)
ENDIF
DO 5 I=1,25
WRITE(*,’(A\)’) AR(I)
READ(*,*)YC(I)
5 YC(I)=YC(I)/100.
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
CALL EOSVNIC(ICALC)
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)
DO 7 J=1,NT
T=TI(J)
CALL EOSVNIC(ICALC)
IF(Z.NE.0D0) THEN
NTS=NTS+1
ROP(I,J)=RO
PAP(I,J)=PA
WP(I,J)=W
ETAP(I,J)=ETA
ELSE