> DE = , F10.5, 2x, DPR = , F10.5, 2x, E1 = , F10.5, 2x, 2E = , F10.5) WRITE (20, 119) GM1, GM2, W1, W2, NV, NK1, NK2 119 FORMAT (2x, GM1 = , F6.2, 2x, GM2 = , F6.2, 2x, W1 = , F6.3, 2x, W2 = , F6 .2, 2x, NV = , I3, 2x, NK1 = , I3, 2x, NK2 = , I3) IF (N. EQ. 0) CALL AER IF (N.EQ. 1) CALL N2 IF (N.EQ. 3) CALL O2 IF (N.EQ. 4) CALL HE1 IF (N.EQ. 5) CALL HE2 IF (N. EQ. 6) CALL CH4 WRITE (20, 71) 71 FORMAT (25x, RASTCHTET RTEPLO ) T = (T1 + T2 )/2. TT = (T3 + T4)/2. RO1 = FRO (P1, T) MU1 = FMU (RO1, T) LAMDA1 = FLA (RO1, T) CP1 = CRCP (RO1, T) * 1000. WRITE (6, 555) RO1, MU1, LAMDA1, CP1 555 FORMAT (2x, 4E15.7) TR = TRO * D2 TZ = TZO * D2 IF (NV.EQ.5) TZ = D2 +1.5 * DPR IF (NV . EQ.5) TR ​​= .866 * (D2 +1.5 * DPR) DELTA = TR-D2 IF (NV.EQ.6) TR = DR * TRO IF (NV.EQ.6) TZ = DR * TZO IF ( NV.EQ.6) DELTA = TR-TZ XN = .75 DO 65 I = 1, 25 65 NX (I) = 1 + ((2 * I-1) ** 2-1) * NX NT = GM1/ (.785 * W1 * RO1 * D1 ** 2) + .5 IF (NT.LE.1) NT = 1 W1 = GM1/(.785 * NT * RO1 * D1 ** 2) CALL ALFAT (W1, D1 , DC, MU1, RO1, CP1, LAMDA1, ALFA1, DZITA1, RE1) N = NK2 IF (N.EQ. 0) CALL AER IF (N.EQ. 1) CALL N2 IF (N.EQ. 3) CALL O2 IF (N.EQ. 4) CALL HE1 IF (N.EQ. 5) CALL HE2 IF (N.EQ. 6) CALL CH4 RO2 = FRO (P2, TT) MU2 = FMU (RO2, TT) LAMDA2 = FLA ( RO2, TT) CP2 = CRCP (RO2, TT) * 1000. 100 IF (NV.NE. 5) CALL ALFAM (NV, W2, TR, TZ, D2, MU2, LAMDA2, RO2, ALFA2, DZITA2, RE2) DTL = ((T1-T4) - (T2-T3))/ ALOG ((T1-T4) - (T2-T3)) DTL = DTCR KQ = 1./(D2/(ALFA1 * D1) +1 ./ALFA2) IF (NV.EQ.5) KQ = 1./( FI/ALFA1 +1 ./ALFA2) IF (NV.EQ.6) KQ = 1./(FIT/ALFA1 +1 ./ALFA2) F1 = 11.3 * Q * 1000./(KQ * DTL) FYD = 1. - .785 * (D2 ** 2.)/(TR * TZ) IF (NV.EQ.6) FYD = 1. - .785 * (D2 ** 2. + (DR ** 2.-D2 ** 2.) * DELTAR/TAUR)/(TR * TZ) FP = GM2/(W2 * RO2) IF (NV.EQ.5) FYD = E1 IF (NV.EQ.6) S = FP/E1 + .785 * DC ** 2 S = FP/АНВ + .785 * DC ** 2. DH = SQRT (S/0.785) EO = (DH-DC)/(2. * TR) NKK = E0 + .5 IF (NV.EQ.5) NKK = .5 + ((DH-DC)/2. - DPR)/TR DHAP = DC +2. * TR * NKK SF = .785 * DHAP ** 2.-DC ** 2.) FHS = 0.785 * DHAP ** 2. FPY = FHS * FYD-.785 * FYD * DC ** 2. IF (NV.EQ.5) FPY = .785 * (DHAP ** 2.-DC ** 2.) * E1 W2Y = GM2/(FPY * RO2) EPS = ABS (W2-W2Y) W2 = W2Y IF ( NV.EQ.5) H = F1/(FP * E2) LO = F1/(3.14 * D2 * NT) IF (NV.EQ.5) LO = F1/(3.14 * D1 * ...
