20
1 FORMAT(//
*10X,’ Расчет коэффициента сжимаемости природного газа’//
*10Х,’ ----------------Метод расчета----------------- ’/
*10Х,’ ’/
*10Х,’ 1. Модифицированный метод NX 19 ’/
*10Х,’ ’/
*10Х,’ 2. Уравнение состояния GERG-91 ’/
*10Х,’ ’/
*10Х,’ 3. Уравнение состояния AGA8-92DC ’/
*10Х,’ ’/
*10Х,’ 4. Уравнение состояния ВНИЦ СМВ ’/
*10Х,’ ’/
*10Х,’---------------------------------------------------’/)
WRITE(*,5)
5 FORMAT(/,3X,
*’Введите порядковый номер метода расчета или 5 для выхода в ДОС’,
*\)
READ(*,*)NVAR
RETURN
END
SUBROUTINE RANGE(NRANGE)
IMPLICIT REAL*8(A-H,О-Z)
COMMON/Z/Z
WRITE(*,1)
1 FORMAT(//
*’ Выбранная Вами методика при заданных параметрах «не работает»’/
*’ Продолжить работу программы ? 0 - нет, 1 - да ’\)
READ(*,*)NRANGE
RETURN
END
SUBROUTINE TABL(YA,YY,PI,TI,ZP,NP,NT,NVAR,AR)
IMPLICIT REAL*8(A-H,О-Z)
CHARACTER*26 AR(25), FNAME
CHARACTER METH(4)*31,A*6,LIN1(5)*9,LIN2(5)*9,LIN3(6)*9,LIN4*9,
*AT(06)*28
CHARACTER*70 F,FZ(11,2)
DIMENSION PI(100),TI(100),ZP(100,100),ZPP(6)
COMMON/RON/RON/YI/YC(25)/NPR/NPR
DATA METH/
*’(модифицированный метод NX19)’,
*’(уравнение состояния GERG-91)’,
*’(уравнение состояния AGA8-92DC)’,
*’(уравнение состояния ВНИЦ СМВ)’/
DATA LIN1/5*’------’/,LIN2/5*’------’/,LIN3/6*’------’/,
*LIN4/’------’/,A/’ - ’/
DATA AT/
*’ T, K’,’ T, K’,’ T, K’,’ T,K’,
*’ T, K’,’ T, K’/
DATA FZ/
*’(3X,F5.2,2X,6(3X,F6.4))’,’(3X,F5.2,5X,A6,5(3X,F6.4))’,
*’(3X,F5.2,2X,2(3X,A6),4(3X,F6.4))’,’(3X,F5.2,2X,3(3X,A6),
*3(3X,F6.4))’,