С * ля адиабаты, скорости звука и вязкости) природного газа по *
С * уравнению состояния ВНИЦ СМВ. *
С * *
С *******************************************************************
IMPLICIT REAL*8(A-H,O-Z)
CHARACTER*26 AR
DIMENSION PI(100),TI(100),ROP(100,100),PAP(100,100),
*WP(100,100),ETAP(100,100)
COMMON/P/P/T/T/RON/RON/YI/YC(25)/NPR/NPR/Z/Z/TS/RO,PA,W
*/ETA/ETA/AR/AR(25)
200 WRITE(*,300)
300 FORMAT(18(/))
WRITE(*,400)
400 FORMAT(
*' Расчет физических свойств природного газа'/
*' по уравнению состояния'/////)
WRITE(*,1)
1 FORMAT(' Введите исходные данные для расчета.'/)
WRITE(*,35)
35 FORMAT(' Введите 0, если состав задан в молярных долях'/
*' или 1, если состав задан в объемных долях '\)
READ(*,*)NPR
IF(NPR.EQ.l) 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\)')