9
*’ Расчет физических свойств природного газа’/
*’ по уравнению состояния’/////)
WRITE(*,1)
1 FORMAT(’ Введите исходные данные для расчета.’/)
WRITE(*,35)
35 FORMAT(’ Введите 0, если состав задан в молярных долях’/
*’ или 1, если состав задан в объемных долях ’\)
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)