17
ПРИЛОЖЕНИЕ В
(рекомендуемое)
Листинг программы расчета коэффициента сжимаемости природного газа
C **********************************************************
C * *
С * Программа расчета коэффициента сжимаемости природного газа *
С * (основной модуль) *
С * *
C **********************************************************
IMPLICIT REAL*8(A-H,O-Z)
CHARACTER*26 AR(25)
DIMENSION PI(100),TI(100),ZP(100,100)
COMMON/P/P/T/T/RON/RON/YI/YC(25)/Z/Z/NPR/NPR
DATA AR/’ метана (СН4)’,’ этана (С2Н6)’,’ пропана (С3Н8)’,
*’ н-бутана (н-С4Н10)’,’ и-бутана (и-С4Н10)’,’ азота (N2)’,
*’ диоксида углерода (СO2)’,’ сероводорода (H2S)’,
*’ ацетилена (С2Н2)’,’ этилена (С2Н4)’,’ пропилена (С3Н6)’,
*’ н-пентана (н-С5Н12)’,’ и-пентана (и-C5H12)’,
*’ нео-пентана (нео-С5Н12)’,’ н-гексана (н-С6Н14)’,
*’ бензола (С6Н6)’,’ н-гептана (н-С7Н16)’,’ толуола (С7Н8)’,
*’ н-октана (н-С8Н18)’,’ н-нонана (н-С9Н20)’,
*’ н-декана (н-С10Н22)’,’ гелия (Не)’,’ водорода (Н2)’,
*’ моноксида углерода (СО)’,’ кислорода (О2)’/
200 WRITE(*,100)
CALL VAR(NVAR)
IF(NVAR.EQ.5) GO TO 134
WRITE(*,l00)
100 FORMAT(25(/))
WRITE(*,1)
1 FORMAT(’ Введите исходные данные для расчета.’/)
IF(NVAR.LE.2) THEN
WRITE(*,’(A\)’)
*’ Плотность при 293.15 К и 101.325 кПа, в кг/куб.м ’
READ(*,*)RON
WRITE(*,53)
53 FORMAT(’ Введите 0, если состав азота и диоксида углерода’,
*’ задан в молярных долях’/
*’ или 1, если состав этих компонентов задан’,
*’ в объемных долях ’\)
READ(*,*)NPR
IF(NPR.EQ.0) WRITE(*,3)
3 FORMAT (’ Значение молярной доли, в мол. %’)
IF(NPR.EQ.l) WRITE(*,33)
33 FORMAT(’ Значение объемной доли, в об. %’)
WRITE(*,’(A\)’) ’ азота (N2)
READ(*,*)YA
YA = YA/100.
WRITE(*,’(A\)’) ’ диоксида углерода (С02) ’
READ(*,*)YY
YY = YY/100.
ELSE