ГОСТ Р 8.769—2011
СIF(RM.LT. 0.55.OR. RM.GT. 0.90)STOP ’REL. MASS OUT OF RANGE*
C IF(X3.LT. 0.0 .OR. X3.GT. 0.30) STOP C02 OUT OF RANGE’
C IF(HS.LT.20.0 OR. HS.GT.48.0 ) STOP ’CALOR. VALUE OUT OF RANGE’
IF((0.55«0.97*X3-0.45*X5).GT.RM)STOP ’CONFLICTING INPUT’
SM = RM’RL
X7 = X5*0.0964D0
X33 = X3*X3
X55 = X5*X5
X77 * X7*X7
BEFF= -0.06500
H = 1000.000
AMOL= 1.0D0/(FA*8EFF)
К = 0
KK = 0
1CALL SMBER(H.SMTI)
IF(ABS(SM-SMT1) .GT. 1.0-6) THEN
CALL SMBER(H*1.0D0.SMT2)
DH = (SM-SMT1 )/(SMT2-SMT1)
H=H+OH
KK = KK*1
IF(KK.GT.20)STOP 1NO CONVERGENCY »V
GO TO 1
END IF
X11 = X1*X1
X12 = X1*X2
X13 = X1*X3
X22 = X2*X2
X23 = X2*X3
X25 * X2*X5
X15 = X1*X5
X17»X1*X7
CALL B11BER(T0.H.B11)
CALL BBER(T0,B11.BEFF)
AMOL= 1.0D0/(FA-*BEFF)
HSBER = X1*H*AMOL+(X5*H5«-X7’ H7)*AMOL
IF(A8S(HS-HSBER).GT.1.0D-4)THEN
K*K+1
IF(K.GT.20}STOP ’NO CONVERGENCY #2’
GOTO 1
END IF
IF(X2.LT.-0.01 .OR. X2.GT.0.5) STOP ’CALC. N2 OUT OF RANGE*
IF{X2+X3.GT.0.5) STOP ’N2 + C02 OUT OF RANGE’
IF(0.55+.4,X2+0.97*X3-O.45*X5).GT.RM)
♦STOP’CONFLICTING RESULT FOR N2’
02 = X2
T = TC*T0
CALL B11BER(T,H.B11)
CALL BBER(T.BH.B)
CALL CBER(T.H.C)
CALL ITER(P.T.B.C.V.Z)
0 = 1.0D0/V
END
SUBROUTINE SMBER(H.SM)
IMPLICIT DOUBLE PRECISION (A-H.O-Z)
COMMON/RBLOK/AMOL.HS
COMMON /XBLOK/X1.X2.X3.X11.X12.X13.X22.X23.X33
>.
x s
.
x t
.
x is
.
x u
.
x z s
.
x s s
.
x t
?
COMMON /MBLOK/GM1R0.GM1R1.GM2.GM3.GM5.GM7.FA.FB.RL.T0.H5.H7.R
GM1" GM1R0*GM1R1*H
X I = (HS4X5*H5»X7*H7rAMOL);H/AMOL
X2 = 1.0D0-X1-X3-X5-X7
26