21
22 WRITE(*,44)
44 FORMAT(//’ Устройство вывода результатов расчета ?,’)
WRITE(*,’(A\)’)
*’ 0 - дисплей, 1 - принтер, 2 - файл на диске ’
READ(*,*)NYST
IF(NYST.EQ.0) OPEN(1,FILE=’CON’)
IF(NYST.EQ.l) OPEN(1,FILE=’PRN’)
IF(NYST.EQ.2) WRITE(*,’(A\)’) ’ Введите имя файла ’
IF(NYST.EQ.2) READ(*,’(A)’)FNAME
IF(NYST.EQ.2) OPEN(1,FILE=FNAME)
IF(NYST.EQ.0) WRITE(*,100)
100 FORMAT(25(/))
IF(NYST.EQ.l) PAUSE
*’ Включите принтер, вставьте бумагу и нажмите <ВВОД> ’
WRITE(1,88)METH(NVAR)
88 FORMAT(
*13X,’Коэффициент сжимаемости природного газа.’/
*18Х,А31/)
NW=3
IF(NVAR.LE.2) THEN
WRITE(1,1)RON
1 FORMAT(’ Плотность при 293.15 К и 101.325 кПа ’,F6.4,’ кг/куб.м’)
NW=NW+1
IF(YA.NE.0D0.OR.YY.NE.0D0) THEN
IF(NPR.EQ.0) WRITE(1,3)
3 FORMAT(’ Содержание в мол. %’)
IF(NPR.EQ.l) WRITE(1,33)
33 FORMAT(’ Содержание в об.%’)
NW=NW+1
IF(YA.NE.0D0) THEN
WRITE(1,5)AR(6),YA* 100.
5 FORMAT(2(A26,F7.4))
NW=NW+1
ENDIF
IF(YY.NE.0D0) THEN
WRITE(1,5)AR(7),YY*100.
NW=NW+1
ENDIF
ENDIF
ELSE
IF(NPR.EQ.0) WRITE(1,3)
IF(NPR.EQ.l) WRITE(1,33)
NW=NW+1
I=1
9 J=I+1
13 CONTINUE
IF(YC(J).NE.0D0) THEN
WRITE(1,5)AR(I),YC(I)*100.,AR(J),YC(J)*100.
NW=NW+1
DO 11 I=J+1,25
IF(YC(I).NE.0D0.AND.I.NE.25) GO TO 9
IF(YC(I).NE.0D0.AND.I.EQ.25) THEN