C----------------------------------------------------------------------C C C SUBROUTINE CKWYP (P, T, Y, ICKWRK, RCKWRK, WDOT) C C INPUT C P - Pressure. C cgs units - dynes/cm**2 C Data type - real scalar C T - Temperature. C cgs units - K C Data type - real scalar C Y - Mass fractions of the species. C cgs units - none C Data type - real array C Dimension Y(*) at least KK, the total number of C species. C ICKWRK - Array of integer workspace. C Data type - integer array C Dimension ICKWRK(*) at least LENIWK. C RCKWRK - Array of real work space. C Data type - real array C Dimension RCKWRK(*) at least LENRWK. C C OUTPUT C WDOT - Chemical molar production rates of the species. C cgs units - moles/(cm**3*sec) C Data type - real array C Dimension WDOT(*) at least KK, the total number of C species. C C END PROLOGUE C C*****double precision IMPLICIT DOUBLE PRECISION (A-H, O-Z), INTEGER (I-N) C*****END double precision C*****single precision C IMPLICIT REAL (A-H, O-Z), INTEGER (I-N) C*****END single precision C PARAMETER (IGLO=12,IREAC=175,KK=16,KSS=15) PARAMETER (NITER=200) DIMENSION ICKWRK(*), RCKWRK(*), Y(*), WDOT(*) DIMENSION XM(IREAC), RF(IREAC), RB(IREAC), W(IREAC) DIMENSION RKF(IGLO), XCON(KK), WT(KK) DIMENSION B(KSS), ABV(KSS), DEN(KSS) LOGICAL LITER, LBOUND(KSS) DATA SMALL/1.D-50/ DATA RU/8.314510D7/ DATA WT/ 2.016, 1.008, 31.999, 17.007, 18.015, 33.007, 34.015, & 15.035, 16.043, 28.011, 44.010, 30.026, 26.038, 28.054, & 30.070, 28.013/ C C CALL CKYTCP( P, T, Y, ICKWRK, RCKWRK, XCON ) C SUMYOW = 0.0 DO K = 1, KK SUMYOW = SUMYOW + Y(K)/WT(K) ENDDO SUMYOW = SUMYOW*T*RU DO K = 1, KK XCON(K) = P*Y(K) / ( SUMYOW*WT(K) ) ENDDO C C SET LOCAL VALUES FOR THE CONCENTRATIONS C XH2 = MAX( XCON(1), SMALL ) XH = MAX( XCON(2), SMALL ) XO2 = MAX( XCON(3), SMALL ) XOH = MAX( XCON(4), SMALL ) XH2O = MAX( XCON(5), SMALL ) XHO2 = MAX( XCON(6), SMALL ) XH2O2 = MAX( XCON(7), SMALL ) XCH3 = MAX( XCON(8), SMALL ) XCH4 = MAX( XCON(9), SMALL ) XCO = MAX( XCON(10), SMALL ) XCO2 = MAX( XCON(11), SMALL ) XCH2O = MAX( XCON(12), SMALL ) XC2H2 = MAX( XCON(13), SMALL ) XC2H4 = MAX( XCON(14), SMALL ) XC2H6 = MAX( XCON(15), SMALL ) XN2 = MAX( XCON(16), SMALL ) BIG = 0.0 DO 20 K = 1, KK IF( XCON(K) .GT. 0.0 ) THEN BIG = MAX( BIG, XCON(K) ) ENDIF WDOT(K) = 0.0 20 CONTINUE C C EFFECTIVE THIRD BODY FOR ALL REACTIONS C CALL THIRDB( XM, XH2, XH, XO2, XOH, XH2O, XHO2, XH2O2, XCH3, XCH4, & XCO, XCO2, XCH2O, XC2H2, XC2H4, XC2H6, XN2 ) C C SET THE ELEMENTARY RATES C CALL ELEMRATE( RF, RB, T, XM ) C C EXPRESSIONS FOR STEADY-STATE SPECIES C DO 40 K = 1, KSS B(K) = 0.0 LBOUND(K) = .TRUE. 40 CONTINUE CALL UPVALUE( 1, B, BIG, LBOUND, CONMAX, XC2H, XC2H5, XHCCOH, & XC2H3, XCH2S, XCH2OH, XCH3OH, XCH3O, XHCO, XCH2, & XCH, XCH2CO, XC, XO, XHCCO ) ADJ = 1.D0/BIG DO 30 N = 1, NITER LITER = .TRUE. CONMAX = 0.0 CALL STEADYE( ABV, DEN, RF, RB, XM, ADJ, CONMAX, SMALL, LBOUND, & LITER, XH2, XH, XO, XO2, XOH, XH2O, XHO2, XH2O2, & XC, XCH, XCH2, XCH2S, XCH3, XCH4, XCO, XCO2, XHCO, & XCH2O, XCH2OH, XCH3O, XCH3OH, XC2H, XC2H2, XC2H3, & XC2H4, XC2H5, XC2H6, XHCCO, XCH2CO, XHCCOH, XN2 ) IF( LITER .AND. CONMAX .LT. 1.0D-5 ) GO TO 35 30 CONTINUE 35 CONTINUE C C NET PRODUCTION RATES FOR SKELETAL MECHANISM C CALL NETRATE( W, RF, RB, XM, XH2, XH, XO, XO2, XOH, XH2O, XHO2, & XH2O2, XC, XCH, XCH2, XCH2S, XCH3, XCH4, XCO, XCO2, & XHCO, XCH2O, XCH2OH, XCH3O, XCH3OH, XC2H, XC2H2, & XC2H3, XC2H4, XC2H5, XC2H6, XHCCO, XCH2CO, XHCCOH, & XN2 ) C C GLOBAL RATES FOR REDUCED MECHANISM C RKF(1) = +W(1) +W(12) +W(14) +W(30) -W(37) -W(43) -W(45) -W(47) & +W(84) +W(98) -W(116) -W(118) -W(121) -W(124) -W(131) & -W(134) -W(142) -W(143) -W(151) -W(153) -W(154) -W(169) & -W(171) -W(174) RKF(2) = +W(2) +W(11) -W(14) +W(15) +W(18) +W(19) +W(22) +W(23) & +W(25) +W(27) +W(32) +W(33) +W(34) +W(35) +W(36) +W(37) & +W(38) +W(39) +W(40) +W(41) +W(45) +W(52) +W(55) +W(56) & +W(57) +W(67) +W(68) +W(70) +W(73) +W(74) +W(77) +W(78) & -W(83) -W(85) -W(86) -W(87) -W(88) -W(91) -W(92) -W(93) & -W(98) -W(99) -W(101) -W(102) -W(105) -W(110) -W(113) & -W(117) +W(118) +W(120) +W(121) +W(124) -W(127) +W(131) & +W(134) -W(135) -W(136) -W(137) -W(144) -W(147) +W(153) & +W(154) -W(155) -W(158) -W(164) -W(165) +W(169) -W(170) & +W(171) +W(174) RKF(3) = +W(3) +W(13) +W(14) -W(23) +W(24) -W(25) -W(30) -W(38) & -W(39) -W(40) -W(41) -W(44) -W(46) -W(52) -W(57) -W(59) & -W(64) -W(67) -W(68) -W(70) -W(73) -W(74) -W(77) -W(78) & -W(79) -W(80) +W(83) -W(89) -W(90) +W(92) +W(99) -W(109) & +W(110) -W(121) -W(124) -W(126) +W(127) -W(131) +W(132) & +W(135) +W(136) +W(137) +W(139) +W(142) +W(143) +W(144) & +W(147) +W(158) +W(164) +W(165) +W(166) -W(169) +W(170) & -W(171) +W(173) RKF(4) = +W(4) -W(32) -W(33) -W(34) -W(35) -W(36) +W(43) +W(44) & +W(45) +W(47) -W(84) +W(86) +W(114) +W(115) +W(116) & +W(117) +W(118) +W(119) -W(166) -W(167) -W(168) -W(173) RKF(5) = +W(5) +W(46) +W(47) -W(84) +W(87) +W(88) -W(114) -W(115) & -W(120) +W(155) RKF(6) = +W(6) +W(20) +W(21) + 2.00*W(23) -W(24) +W(25) +W(28) & +W(29) + 2.00*W(30) -W(48) -W(50) +W(70) -W(72) +W(73) & -W(76) +W(78) +W(79) + 2.00*W(80) + 2.00*W(89) +W(90) & -W(92) +W(105) + 2.00*W(109) -W(110) +W(113) & + 2.00*W(121) +W(124) +W(125) +W(126) -W(127) +W(131) & -W(132) - 2.00*W(136) -W(137) - 2.00*W(139) -W(147) & + 2.00*W(169) +W(171) -W(172) -W(173) +W(174) RKF(7) = +W(7) -W(23) +W(24) +W(25) -W(30) +W(48) +W(49) -W(70) & +W(72) +W(73) -W(76) -W(89) +W(91) +W(92) -W(95) +W(110) & +W(116) -W(121) -W(123) -W(125) +W(127) -W(128) -W(129) & +W(134) +W(135) + 2.00*W(136) +W(138) +W(139) -W(141) & -W(146) -W(147) -W(149) -W(150) +W(171) +W(172) -W(173) RKF(8) = +W(8) +W(13) +W(14) -W(23) +W(24) -W(25) -W(30) +W(48) & +W(50) +W(54) -W(70) +W(72) -W(73) +W(76) -W(78) -W(80) & -W(82) -W(89) -W(90) +W(92) +W(99) -W(109) +W(110) & -W(121) -W(124) -W(125) -W(126) +W(127) -W(131) +W(132) & + 2.00*W(136) +W(137) +W(139) +W(142) +W(143) +W(147) & +W(158) +W(164) +W(165) +W(166) -W(169) -W(171) +W(172) & +W(173) RKF(9) = +W(9) -W(11) -W(13) -W(14) -W(18) -W(19) +W(23) -W(24) & -W(27) +W(30) -W(48) +W(51) -W(52) -W(54) +W(58) -W(61) & +W(62) -W(66) -W(67) -W(68) +W(70) +W(71) -W(72) -W(74) & +W(75) -W(77) +W(80) +W(82) +W(89) +W(90) -W(92) +W(93) & +W(94) -W(96) -W(97) -W(99) -W(103) -W(104) +W(109) & -W(110) -W(111) -W(112) +W(117) +W(121) +W(124) +W(125) & +W(126) -W(127) +W(131) -W(132) - 2.00*W(136) -W(138) & -W(139) +W(141) +W(144) + 2.00*W(145) +W(146) +W(147) & +W(149) +W(150) +W(151) +W(155) +W(156) +W(159) -W(164) & -W(165) -W(166) +W(169) +W(171) - 2.00*W(172) RKF(10) = +W(10) +W(11) +W(18) +W(19) +W(27) -W(49) -W(51) +W(52) & -W(58) -W(60) -W(62) -W(65) +W(67) +W(68) -W(71) +W(74) & -W(75) +W(77) -W(80) +W(95) +W(96) +W(97) +W(103) & +W(104) -W(109) +W(111) +W(112) -W(117) +W(118) +W(123) & +W(128) +W(129) -W(135) -W(144) -W(145) +W(153) +W(154) & -W(155) +W(157) -W(158) -W(159) +W(172) RKF(11) = +W(26) -W(73) +W(76) -W(156) -W(157) +W(173) RKF(12) = +W(42) +W(43) +W(47) +W(61) +W(66) +W(83) +W(85) +W(86) & +W(87) +W(88) +W(92) +W(95) +W(96) +W(97) +W(99) & +W(100) +W(101) +W(102) +W(103) +W(104) +W(108) +W(110) & +W(111) +W(112) +W(113) -W(126) +W(143) -W(145) C C REDUCED MECHANISM C C O2 + 2CO = 2CO2 C H + O2 + CO = OH + CO2 C H2 + O2 + CO = H + OH + CO2 C HO2 + CO = OH + CO2 C O2 + H2O2 + CO = OH + HO2 + CO2 C O2 + 0.50 C2H2 = H + CO2 C O2 + CH3 + CO + C2H4 = CH4 + CO2 + CH2O + 0.50 C2H2 C O2 + 2CH3 = H2 + CH4 + CO2 C O2 + 2CH3 + CO = CH4 + CO2 + CH2O C O2 + CH3 + CO = H + CO2 + CH2O C O2 + CO + C2H6 = CH4 + CO2 + CH2O C H + OH = H2O C C SPECIES PRODUCTION RATES C C H2 WDOT(1) = - RKF(3) + RKF(8) C H WDOT(2) = - RKF(2) + RKF(3) + RKF(6) + RKF(10) - RKF(12) C O2 WDOT(3) = - RKF(1) - RKF(2) - RKF(3) - RKF(5) - RKF(6) - RKF(7) & - RKF(8) - RKF(9) - RKF(10) - RKF(11) C OH WDOT(4) = + RKF(2) + RKF(3) + RKF(4) + RKF(5) - RKF(12) C H2O WDOT(5) = + RKF(12) C HO2 WDOT(6) = - RKF(4) + RKF(5) C H2O2 WDOT(7) = - RKF(5) C CH3 WDOT(8) = - RKF(7) - 2.00*RKF(8) - 2.00*RKF(9) - RKF(10) C CH4 WDOT(9) = + RKF(7) + RKF(8) + RKF(9) + RKF(11) C CO WDOT(10) = - 2.00*RKF(1) - RKF(2) - RKF(3) - RKF(4) - RKF(5) & - RKF(7) - RKF(9) - RKF(10) - RKF(11) C CO2 WDOT(11) = + 2.00*RKF(1) + RKF(2) + RKF(3) + RKF(4) + RKF(5) & + RKF(6) + RKF(7) + RKF(8) + RKF(9) + RKF(10) & + RKF(11) C CH2O WDOT(12) = + RKF(7) + RKF(9) + RKF(10) + RKF(11) C C2H2 WDOT(13) = - 0.50*RKF(6) + 0.50*RKF(7) C C2H4 WDOT(14) = - RKF(7) C C2H6 WDOT(15) = - RKF(11) C N2 WDOT(16) = 0.0 DO K = 1, KK IF( XCON(K) .LE. 0.D0 .AND. WDOT(K) .LT. 0.D0 ) WDOT(K) = 0.D0 END DO RETURN END SUBROUTINE THIRDB( XM, XH2, XH, XO2, XOH, XH2O, XHO2, XH2O2, XCH3, & XCH4, XCO, XCO2, XCH2O, XC2H2, XC2H4, XC2H6, & XN2 ) C*****double precision IMPLICIT DOUBLE PRECISION (A-H, O-Z), INTEGER (I-N) C*****END double precision C*****single precision C IMPLICIT REAL (A-H, O-Z), INTEGER (I-N) C*****END single precision DIMENSION XM(*) C C EFFECTIVE THIRD BODY FOR ALL REACTIONS C XM(1) = + 2.400*XH2+XH+XO2+XOH+15.400*XH2O+XHO2+XH2O2+XCH3 & + 2.000*XCH4+ 1.750*XCO+ 3.600*XCO2+XCH2O+XC2H2+XC2H4 & + 3.000*XC2H6+XN2 XM(2) = + 2.000*XH2+XH+XO2+XOH+ 6.000*XH2O+XHO2+XH2O2+XCH3 & + 2.000*XCH4+ 1.500*XCO+ 2.000*XCO2+XCH2O+XC2H2+XC2H4 & + 3.000*XC2H6+XN2 XM(12) = + 2.000*XH2+XH+ 6.000*XO2+XOH+ 6.000*XH2O+XHO2+XH2O2+XCH3 & + 2.000*XCH4+ 1.500*XCO+ 3.500*XCO2+XCH2O+XC2H2+XC2H4 & + 3.000*XC2H6+XN2 XM(33) = +XH2+XH+XOH+XHO2+XH2O2+XCH3+XCH4+ 0.750*XCO+ 1.500*XCO2 & +XCH2O+XC2H2+XC2H4+ 1.500*XC2H6 XM(38) = +XH+XO2+XOH+XHO2+XH2O2+XCH3+ 2.000*XCH4+XCO+XCH2O+XC2H2 & +XC2H4+ 3.000*XC2H6+XN2 XM(42) = + 0.730*XH2+XH+XO2+XOH+ 3.650*XH2O+XHO2+XH2O2+XCH3 & + 2.000*XCH4+XCO+XCO2+XCH2O+XC2H2+XC2H4+ 3.000*XC2H6+XN2 XM(49) = + 2.000*XH2+XH+XO2+XOH+ 6.000*XH2O+XHO2+XH2O2+XCH3 & + 2.000*XCH4+ 1.500*XCO+ 2.000*XCO2+XCH2O+XC2H2+XC2H4 & + 3.000*XC2H6+XN2 XM(51) = + 2.000*XH2+XH+XO2+XOH+ 6.000*XH2O+XHO2+XH2O2+XCH3 & + 2.000*XCH4+ 1.500*XCO+ 2.000*XCO2+XCH2O+XC2H2+XC2H4 & + 3.000*XC2H6+XN2 XM(53) = + 2.000*XH2+XH+XO2+XOH+ 6.000*XH2O+XHO2+XH2O2+XCH3 & + 2.000*XCH4+ 1.500*XCO+ 2.000*XCO2+XCH2O+XC2H2+XC2H4 & + 3.000*XC2H6+XN2 XM(55) = + 2.000*XH2+XH+XO2+XOH+ 6.000*XH2O+XHO2+XH2O2+XCH3 & + 2.000*XCH4+ 1.500*XCO+ 2.000*XCO2+XCH2O+XC2H2+XC2H4 & + 3.000*XC2H6+XN2 XM(56) = + 2.000*XH2+XH+XO2+XOH+ 6.000*XH2O+XHO2+XH2O2+XCH3 & + 2.000*XCH4+ 1.500*XCO+ 2.000*XCO2+XCH2O+XC2H2+XC2H4 & + 3.000*XC2H6+XN2 XM(58) = + 2.000*XH2+XH+XO2+XOH+ 6.000*XH2O+XHO2+XH2O2+XCH3 & + 2.000*XCH4+ 1.500*XCO+ 2.000*XCO2+XCH2O+XC2H2+XC2H4 & + 3.000*XC2H6+XN2 XM(62) = + 2.000*XH2+XH+XO2+XOH+ 6.000*XH2O+XHO2+XH2O2+XCH3 & + 2.000*XCH4+ 1.500*XCO+ 2.000*XCO2+XCH2O+XC2H2+XC2H4 & + 3.000*XC2H6+XN2 XM(69) = + 2.000*XH2+XH+XO2+XOH+ 6.000*XH2O+XHO2+XH2O2+XCH3 & + 2.000*XCH4+ 1.500*XCO+ 2.000*XCO2+XCH2O+XC2H2+XC2H4 & + 3.000*XC2H6+XN2 XM(70) = + 2.000*XH2+XH+XO2+XOH+ 6.000*XH2O+XHO2+XH2O2+XCH3 & + 2.000*XCH4+ 1.500*XCO+ 2.000*XCO2+XCH2O+XC2H2+XC2H4 & + 3.000*XC2H6+XN2 XM(71) = + 2.000*XH2+XH+XO2+XOH+ 6.000*XH2O+XHO2+XH2O2+XCH3 & + 2.000*XCH4+ 1.500*XCO+ 2.000*XCO2+XCH2O+XC2H2+XC2H4 & + 3.000*XC2H6+XN2 XM(73) = + 2.000*XH2+XH+XO2+XOH+ 6.000*XH2O+XHO2+XH2O2+XCH3 & + 2.000*XCH4+ 1.500*XCO+ 2.000*XCO2+XCH2O+XC2H2+XC2H4 & + 3.000*XC2H6+XN2 XM(75) = + 2.000*XH2+XH+XO2+XOH+ 6.000*XH2O+XHO2+XH2O2+XCH3 & + 2.000*XCH4+ 1.500*XCO+ 2.000*XCO2+XCH2O+XC2H2+XC2H4 & + 3.000*XC2H6+XN2 XM(82) = + 2.000*XH2+XH+XO2+XOH+ 6.000*XH2O+XHO2+XH2O2+XCH3 & + 2.000*XCH4+ 1.500*XCO+ 2.000*XCO2+XCH2O+XC2H2+XC2H4 & + 3.000*XC2H6+XN2 XM(84) = + 2.000*XH2+XH+XO2+XOH+ 6.000*XH2O+XHO2+XH2O2+XCH3 & + 2.000*XCH4+ 1.500*XCO+ 2.000*XCO2+XCH2O+XC2H2+XC2H4 & + 3.000*XC2H6+XN2 XM(94) = + 2.000*XH2+XH+XO2+XOH+ 6.000*XH2O+XHO2+XH2O2+XCH3 & + 2.000*XCH4+ 1.500*XCO+ 2.000*XCO2+XCH2O+XC2H2+XC2H4 & + 3.000*XC2H6+XN2 XM(130) = + 2.000*XH2+XH+XO2+XOH+ 6.000*XH2O+XHO2+XH2O2+XCH3 & + 2.000*XCH4+ 1.500*XCO+ 2.000*XCO2+XCH2O+XC2H2+XC2H4 & + 3.000*XC2H6+XN2 XM(139) = + 2.000*XH2+XH+XO2+XOH+ 6.000*XH2O+XHO2+XH2O2+XCH3 & + 2.000*XCH4+ 1.500*XCO+ 2.000*XCO2+XCH2O+XC2H2+XC2H4 & + 3.000*XC2H6+XN2 XM(145) = + 2.000*XH2+XH+XO2+XOH+ 6.000*XH2O+XHO2+XH2O2+XCH3 & + 2.000*XCH4+ 1.500*XCO+ 2.000*XCO2+XCH2O+XC2H2+XC2H4 & + 3.000*XC2H6+XN2 XM(156) = + 2.000*XH2+XH+XO2+XOH+ 6.000*XH2O+XHO2+XH2O2+XCH3 & + 2.000*XCH4+ 1.500*XCO+ 2.000*XCO2+XCH2O+XC2H2+XC2H4 & + 3.000*XC2H6+XN2 XM(165) = + 2.000*XH2+XH+XO2+XOH+XHO2+XH2O2+XCH3+ 2.000*XCH4 & + 1.500*XCO+ 2.000*XCO2+XCH2O+XC2H2+XC2H4+ 3.000*XC2H6 & +XN2 XM(172) = + 2.000*XH2+XH+XO2+XOH+ 6.000*XH2O+XHO2+XH2O2+XCH3 & + 2.000*XCH4+ 1.500*XCO+ 2.000*XCO2+XCH2O+XC2H2+XC2H4 & + 3.000*XC2H6+XN2 RETURN END SUBROUTINE ELEMRATE( RF, RB, T, XM ) C*****double precision IMPLICIT DOUBLE PRECISION (A-H, O-Z), INTEGER (I-N) C*****END double precision C*****single precision C IMPLICIT REAL (A-H, O-Z), INTEGER (I-N) C*****END single precision DIMENSION RF(*), RB(*), XM(*) C RATE CONSTANTS FOR SKELETAL MECHANSIM RUC = 8.31451D0/4.184D0 ALOGT = DLOG(T) RTR = 1.0D3/(RUC*T) TM1 = 1.0/T TM2 = 1.0/T/T TM3 = 1.0/T/T/T TP1 = T TP2 = T*T TP3 = T*T*T RF(1) = 1.200D+17 * EXP( -1.000*ALOGT ) RB(1) = EXP(+(-2.51660D+06)*TM3+( 3.69147D+04)*TM2 & +(-5.95882D+04)*TM1+( 3.63993D+01)+(-1.05609D-03)*TP1 & +( 1.85706D-07)*TP2+(-1.85434D-11)*TP3) RF(2) = 5.000D+17 * EXP( -1.000*ALOGT ) RB(2) = EXP(+( 3.64734D+06)*TM3+(-4.88930D+04)*TM2 & +(-5.06755D+04)*TM1+( 3.38868D+01)+(-2.21767D-04)*TP1 & +(-4.13655D-08)*TP2+( 7.53188D-12)*TP3) RF(3) = 5.000D+04 * EXP( 2.670*ALOGT - 6.29000*RTR ) RB(3) = EXP(+(-2.23489D+07)*TM3+( 2.80492D+05)*TM2 & +(-3.94618D+03)*TM1+( 2.82285D+01)+( 2.13355D-03)*TP1 & +(-4.33965D-07)*TP2+( 4.54247D-11)*TP3) RF(4) = 2.000D+13 RB(4) = EXP(+( 5.92451D+05)*TM3+(-3.41134D+04)*TM2 & +(-2.65087D+04)*TM1+( 3.02310D+01)+( 4.08842D-04)*TP1 & +(-8.62910D-08)*TP2+( 9.70706D-12)*TP3) RF(5) = 9.630D+06 * EXP( 2.000*ALOGT - 4.00000*RTR ) RB(5) = EXP(+(-1.30325D+07)*TM3+( 1.53786D+05)*TM2 & +(-1.04088D+04)*TM1+( 2.69031D+01)+( 2.32010D-03)*TP1 & +(-4.64925D-07)*TP2+( 4.48793D-11)*TP3) RF(6) = 5.700D+13 RB(6) = EXP(+(-5.23695D+06)*TM3+( 6.29753D+04)*TM2 & +(-8.91655D+04)*TM1+( 3.60343D+01)+(-4.42073D-04)*TP1 & +( 1.49384D-07)*TP2+(-1.62519D-11)*TP3) RF(7) = 8.000D+13 RB(7) = EXP(+(-1.03571D+07)*TM3+( 1.21403D+05)*TM2 & +(-4.64114D+04)*TM1+( 3.49933D+01)+(-7.79879D-04)*TP1 & +( 2.10360D-07)*TP2+(-2.33865D-11)*TP3) RF(8) = 1.500D+13 RB(8) = EXP(+( 8.60710D+06)*TM3+(-9.91271D+04)*TM2 & +(-9.45738D+04)*TM1+( 3.24001D+01)+( 1.72519D-04)*TP1 & +( 1.39133D-08)*TP2+(-6.20876D-12)*TP3) RF(9) = 1.500D+13 RB(9) = EXP(+(-1.97482D+06)*TM3+( 3.81025D+04)*TM2 & +(-5.06441D+04)*TM1+( 3.23466D+01)+(-7.28854D-04)*TP1 & +( 2.06068D-07)*TP2+(-2.38317D-11)*TP3) RF(10) = 8.430D+13 RB(10) = EXP(+(-1.45276D+07)*TM3+( 1.77812D+05)*TM2 & +(-3.52415D+04)*TM1+( 3.59934D+01)+(-7.92692D-04)*TP1 & +( 2.06574D-07)*TP2+(-2.21450D-11)*TP3) RF(11) = 1.020D+09 * EXP( 1.500*ALOGT - 8.60000*RTR ) RB(11) = EXP(+( 1.95738D+07)*TM3+(-2.19284D+05)*TM2 & +(-2.40285D+03)*TM1+( 2.50874D+01)+( 1.87962D-03)*TP1 & +(-2.71623D-07)*TP2+( 2.25030D-11)*TP3) RF(12) = 6.020D+14 * EXP( - 3.00000*RTR ) RB(12) = EXP(+(-2.32390D+06)*TM3+( 8.09642D+04)*TM2 & +(-6.57572D+04)*TM1+( 4.16383D+01)+(-1.00275D-03)*TP1 & +( 1.69711D-07)*TP2+(-1.55578D-11)*TP3) RF(13) = 3.000D+13 RB(13) = EXP(+( 9.14681D+06)*TM3+(-1.21514D+05)*TM2 & +(-4.30054D+04)*TM1+( 3.03552D+01)+( 8.90916D-04)*TP1 & +(-2.05596D-07)*TP2+( 2.07550D-11)*TP3) RF(14) = 3.000D+13 RB(14) = EXP(+( 1.10085D+07)*TM3+(-9.08238D+04)*TM2 & +(-5.59387D+04)*TM1+( 3.80351D+01)+(-6.93068D-04)*TP1 & +( 1.62980D-07)*TP2+(-1.81746D-11)*TP3) RF(15) = 3.900D+13 * EXP( - 3.54000*RTR ) RB(15) = EXP(+( 1.73749D+07)*TM3+(-2.19790D+05)*TM2 & +(-7.96059D+03)*TM1+( 2.63453D+01)+( 6.78941D-04)*TP1 & +(-1.02404D-07)*TP2+( 6.90075D-12)*TP3) RF(16) = 1.000D+13 RB(16) = EXP(+(-1.52049D+07)*TM3+( 1.53537D+05)*TM2 & +(-3.70153D+04)*TM1+( 3.05236D+01)+( 5.82746D-04)*TP1 & +(-1.50232D-07)*TP2+( 1.73226D-11)*TP3) RF(17) = 1.000D+13 RB(17) = EXP(+( 1.19988D+07)*TM3+(-1.65341D+05)*TM2 & +(-3.94627D+04)*TM1+( 2.68081D+01)+( 1.12910D-03)*TP1 & +(-2.83082D-07)*TP2+( 3.38972D-11)*TP3) RF(18) = 3.880D+05 * EXP( 2.500*ALOGT - 3.10000*RTR ) RB(18) = EXP(+( 8.23274D+06)*TM3+(-9.96787D+04)*TM2 & +(-4.53584D+03)*TM1+( 2.44522D+01)+( 2.63451D-03)*TP1 & +(-4.77102D-07)*TP2+( 4.40625D-11)*TP3) RF(19) = 1.300D+05 * EXP( 2.500*ALOGT - 5.00000*RTR ) RB(19) = EXP(+(-1.89709D+07)*TM3+( 2.19199D+05)*TM2 & +(-3.04453D+03)*TM1+( 2.70742D+01)+( 2.08815D-03)*TP1 & +(-3.44251D-07)*TP2+( 2.74880D-11)*TP3) RF(20) = 5.000D+13 RB(20) = EXP(+(-8.19898D+06)*TM3+( 9.00509D+04)*TM2 & +(-3.97874D+04)*TM1+( 3.07746D+01)+( 7.39886D-04)*TP1 & +(-1.11002D-07)*TP2+( 6.76476D-12)*TP3) RF(21) = 1.020D+07 * EXP( 2.000*ALOGT - 1.90000*RTR ) RB(21) = EXP(+(-1.46891D+07)*TM3+( 1.84717D+05)*TM2 & +(-1.20731D+04)*TM1+( 3.00062D+01)+( 1.00615D-03)*TP1 & +(-8.90544D-08)*TP2+( 3.91193D-12)*TP3) RF(22) = 4.600D+19 * EXP( -1.410*ALOGT - 28.95000*RTR ) RB(22) = EXP(+( 1.01289D+07)*TM3+(-1.92905D+05)*TM2 & +( 2.19999D+03)*TM1+( 3.03345D+01)+(-1.97314D-04)*TP1 & +(-2.89392D-08)*TP2+( 7.84947D-12)*TP3) RF(23) = 1.020D+07 * EXP( 2.000*ALOGT - 1.90000*RTR ) RB(23) = EXP(+(-2.88069D+07)*TM3+( 2.98392D+05)*TM2 & +(-2.59020D+04)*TM1+( 2.56523D+01)+( 2.72769D-03)*TP1 & +(-5.32393D-07)*TP2+( 5.50080D-11)*TP3) RF(24) = 3.000D+13 RB(24) = EXP(+( 1.06166D+07)*TM3+(-1.35563D+05)*TM2 & +(-4.50278D+04)*TM1+( 3.52174D+01)+(-5.71421D-04)*TP1 & +( 1.66026D-07)*TP2+(-1.94609D-11)*TP3) RF(25) = 1.920D+07 * EXP( 1.830*ALOGT - 0.22000*RTR ) RB(25) = EXP(+(-1.15317D+07)*TM3+( 6.96782D+04)*TM2 & +(-1.41140D+04)*TM1+( 2.30488D+01)+( 2.51186D-03)*TP1 & +(-4.79217D-07)*TP2+( 4.75268D-11)*TP3) RF(26) = 1.320D+14 RB(26) = EXP(+(-5.58411D+06)*TM3+( 1.20929D+04)*TM2 & +(-3.94571D+04)*TM1+( 3.11836D+01)+( 9.58585D-04)*TP1 & +(-2.04178D-07)*TP2+( 2.15354D-11)*TP3) RF(27) = 8.980D+07 * EXP( 1.920*ALOGT - 5.69000*RTR ) RB(27) = EXP(+( 1.19559D+06)*TM3+(-1.20562D+04)*TM2 & +(-4.09241D+03)*TM1+( 2.49785D+01)+( 2.32012D-03)*TP1 & +(-4.30151D-07)*TP2+( 4.01341D-11)*TP3) RF(28) = 1.000D+14 RB(28) = EXP(+(-1.11426D+07)*TM3+( 6.32901D+04)*TM2 & +(-5.19314D+04)*TM1+( 3.02585D+01)+( 1.25135D-03)*TP1 & +(-2.39710D-07)*TP2+( 2.50928D-11)*TP3) RF(29) = 1.000D+13 * EXP( - 8.00000*RTR ) RB(29) = EXP(+( 1.13827D+07)*TM3+(-1.48793D+05)*TM2 & +(-1.51461D+03)*TM1+( 2.54177D+01)+( 6.88118D-04)*TP1 & +(-9.85504D-08)*TP2+( 5.79617D-12)*TP3) RF(30) = 1.750D+12 * EXP( - 1.35000*RTR ) RB(30) = EXP(+(-8.73406D+05)*TM3+(-4.42754D+03)*TM2 & +(-2.49305D+04)*TM1+( 2.70007D+01)+( 8.25676D-04)*TP1 & +(-1.73313D-07)*TP2+( 1.79626D-11)*TP3) RF(31) = 2.500D+12 * EXP( - 47.80000*RTR ) RB(31) = EXP(+( 8.02559D+06)*TM3+(-5.51176D+04)*TM2 & +(-2.80744D+04)*TM1+( 3.22875D+01)+(-7.49661D-04)*TP1 & +( 1.41505D-07)*TP2+(-1.28544D-11)*TP3) RF(32) = 1.000D+14 * EXP( - 40.00000*RTR ) RB(32) = EXP(+( 1.67825D+07)*TM3+(-1.85677D+05)*TM2 & +( 2.00793D+02)*TM1+( 2.76827D+01)+( 2.70099D-04)*TP1 & +(-1.61133D-08)*TP2+(-2.80631D-12)*TP3) RF(33) = 2.800D+18 * EXP( -0.860*ALOGT ) RB(33) = EXP(+( 1.95829D+06)*TM3+(-8.96287D+02)*TM2 & +(-2.42563D+04)*TM1+( 3.69564D+01)+(-5.18190D-04)*TP1 & +( 2.28756D-08)*TP2+( 4.24094D-14)*TP3) RF(34) = 3.000D+20 * EXP( -1.720*ALOGT ) RB(34) = EXP(+( 8.69457D+06)*TM3+(-8.61799D+04)*TM2 & +(-2.37070D+04)*TM1+( 3.57880D+01)+(-1.20877D-03)*TP1 & +( 1.58325D-07)*TP2+(-1.35799D-11)*TP3) RF(35) = 9.380D+18 * EXP( -0.760*ALOGT ) RB(35) = EXP(+( 1.17500D+06)*TM3+( 9.02041D+03)*TM2 & +(-2.43202D+04)*TM1+( 3.88448D+01)+(-4.37890D-04)*TP1 & +( 7.12568D-09)*TP2+( 1.62640D-12)*TP3) RF(36) = 3.750D+20 * EXP( -1.720*ALOGT ) RB(36) = EXP(+( 8.69457D+06)*TM3+(-8.61799D+04)*TM2 & +(-2.37070D+04)*TM1+( 3.60112D+01)+(-1.20877D-03)*TP1 & +( 1.58325D-07)*TP2+(-1.35799D-11)*TP3) RF(37) = 8.300D+13 * EXP( - 14.41300*RTR ) RB(37) = EXP(+( 6.16394D+06)*TM3+(-8.58077D+04)*TM2 & +( 1.65979D+03)*TM1+( 2.81102D+01)+( 8.34323D-04)*TP1 & +(-2.27071D-07)*TP2+( 2.60752D-11)*TP3) RF(38) = 1.000D+18 * EXP( -1.000*ALOGT ) RB(38) = EXP(+( 5.08245D+06)*TM3+(-6.46086D+04)*TM2 & +(-5.15999D+04)*TM1+( 3.53104D+01)+(-2.11311D-04)*TP1 & +(-2.79241D-08)*TP2+( 4.39979D-12)*TP3) RF(39) = 9.000D+16 * EXP( -0.600*ALOGT ) RB(39) = EXP(+( 1.94930D+06)*TM3+(-2.49418D+04)*TM2 & +(-5.18553D+04)*TM1+( 3.56199D+01)+( 1.09889D-04)*TP1 & +(-9.09238D-08)*TP2+( 1.07358D-11)*TP3) RF(40) = 6.000D+19 * EXP( -1.250*ALOGT ) RB(40) = EXP(+( 7.04068D+06)*TM3+(-8.94004D+04)*TM2 & +(-5.14402D+04)*TM1+( 3.77063D+01)+(-4.12061D-04)*TP1 & +( 1.14507D-08)*TP2+( 4.39807D-13)*TP3) RF(41) = 5.500D+20 * EXP( -2.000*ALOGT ) RB(41) = EXP(+( 1.29153D+07)*TM3+(-1.63776D+05)*TM2 & +(-5.09612D+04)*TM1+( 3.48266D+01)+(-1.01431D-03)*TP1 & +( 1.29575D-07)*TP2+(-1.14401D-11)*TP3) RF(42) = 2.200D+22 * EXP( -2.000*ALOGT ) RB(42) = EXP(+( 2.90282D+06)*TM3+(-4.02691D+04)*TM2 & +(-5.90676D+04)*TM1+( 4.05688D+01)+(-1.19039D-03)*TP1 & +( 1.28742D-07)*TP2+(-8.88854D-12)*TP3) RF(43) = 3.970D+12 * EXP( - 0.67100*RTR ) RB(43) = EXP(+(-1.82101D+06)*TM3+(-1.21302D+04)*TM2 & +(-2.69644D+04)*TM1+( 2.74581D+01)+( 1.07755D-03)*TP1 & +(-3.00754D-07)*TP2+( 3.52018D-11)*TP3) RF(44) = 2.800D+13 * EXP( - 1.06800*RTR ) RB(44) = EXP(+( 2.02756D+06)*TM3+(-4.98290D+04)*TM2 & +(-2.79704D+04)*TM1+( 3.12979D+01)+( 4.19299D-04)*TP1 & +(-7.28496D-08)*TP2+( 6.57497D-12)*TP3) RF(45) = 1.340D+14 * EXP( - 0.63500*RTR ) RB(45) = EXP(+( 6.75639D+06)*TM3+(-1.19921D+05)*TM2 & +(-1.79156D+04)*TM1+( 2.81934D+01)+( 1.24317D-03)*TP1 & +(-3.13362D-07)*TP2+( 3.57823D-11)*TP3) RF(46) = 1.210D+07 * EXP( 2.000*ALOGT - 5.20000*RTR ) RB(46) = EXP(+(-1.15974D+07)*TM3+( 1.38071D+05)*TM2 & +(-1.19370D+04)*TM1+( 2.78618D+01)+( 2.33056D-03)*TP1 & +(-4.51484D-07)*TP2+( 4.17473D-11)*TP3) RF(47) = 1.000D+13 * EXP( - 3.60000*RTR ) RB(47) = EXP(+( 8.12238D+05)*TM3+(-5.66781D+04)*TM2 & +(-3.55569D+04)*TM1+( 2.56172D+01)+( 1.79165D-03)*TP1 & +(-4.50681D-07)*TP2+( 4.84013D-11)*TP3) RF(48) = 1.100D+14 RB(48) = EXP(+( 4.09971D+06)*TM3+(-3.92623D+04)*TM2 & +(-1.17301D+04)*TM1+( 3.32169D+01)+( 7.79596D-05)*TP1 & +( 3.94975D-08)*TP2+(-7.78065D-12)*TP3) RFH49 = 2.500D+16 * EXP( -0.800*ALOGT ) RFL49 = 3.200D+27 * EXP( -3.140*ALOGT - 1.23000*RTR ) PR = RFL49*XM(49)/RFH49 PRLOG = DLOG10(PR) F4 = 0.6800 F5 = 78.00 F6 = 1995. F7 = 5590. FC = (1.0-F4)*EXP(-T/F5)+F4*EXP(-T/F6)+EXP(-F7/T) FCLOG = DLOG10(FC) CPRLOG = PRLOG - (0.4 + 0.67*FCLOG) X = CPRLOG/(0.75-1.27*FCLOG-0.14*CPRLOG) F = 10.0**( FCLOG/(1.0+X*X) ) PCOR = PR*F/(1.0+PR) RF(49) = RFH49*PCOR RB(49) = EXP(+(-1.11237D+07)*TM3+( 1.34322D+05)*TM2 & +(-5.57939D+04)*TM1+( 3.62512D+01)+(-7.27296D-04)*TP1 & +( 3.33251D-08)*TP2+( 2.55764D-12)*TP3) RB(49) = RB(49)*PCOR RF(50) = 3.000D+13 RB(50) = EXP(+( 1.38441D+07)*TM3+(-1.62102D+05)*TM2 & +(-5.40836D+03)*TM1+( 2.87330D+01)+( 6.14592D-04)*TP1 & +(-1.35471D-07)*TP2+( 1.00431D-11)*TP3) RFH51 = 1.270D+16 * EXP( -0.630*ALOGT - 0.38300*RTR ) RFL51 = 2.477D+33 * EXP( -4.760*ALOGT - 2.44000*RTR ) PR = RFL51*XM(51)/RFH51 PRLOG = DLOG10(PR) F4 = 0.7830 F5 = 74.00 F6 = 2941. F7 = 6964. FC = (1.0-F4)*EXP(-T/F5)+F4*EXP(-T/F6)+EXP(-F7/T) FCLOG = DLOG10(FC) CPRLOG = PRLOG - (0.4 + 0.67*FCLOG) X = CPRLOG/(0.75-1.27*FCLOG-0.14*CPRLOG) F = 10.0**( FCLOG/(1.0+X*X) ) PCOR = PR*F/(1.0+PR) RF(51) = RFH51*PCOR RB(51) = EXP(+(-3.05740D+07)*TM3+( 3.55834D+05)*TM2 & +(-5.39874D+04)*TM1+( 3.85737D+01)+(-5.99776D-04)*TP1 & +(-6.42667D-08)*TP2+( 1.46495D-11)*TP3) RB(51) = RB(51)*PCOR RF(52) = 6.600D+08 * EXP( 1.620*ALOGT - 10.84000*RTR ) RB(52) = EXP(+( 2.00690D+07)*TM3+(-2.23100D+05)*TM2 & +(-4.53102D+03)*TM1+( 2.61978D+01)+( 1.98643D-03)*TP1 & +(-2.77081D-07)*TP2+( 2.12717D-11)*TP3) RFH53 = 1.090D+12 * EXP( 0.480*ALOGT + 0.26000*RTR ) RFL53 = 1.350D+24 * EXP( -2.570*ALOGT - 1.42500*RTR ) PR = RFL53*XM(53)/RFH53 PRLOG = DLOG10(PR) F4 = 0.7824 F5 = 271.0 F6 = 2755. F7 = 6570. FC = (1.0-F4)*EXP(-T/F5)+F4*EXP(-T/F6)+EXP(-F7/T) FCLOG = DLOG10(FC) CPRLOG = PRLOG - (0.4 + 0.67*FCLOG) X = CPRLOG/(0.75-1.27*FCLOG-0.14*CPRLOG) F = 10.0**( FCLOG/(1.0+X*X) ) PCOR = PR*F/(1.0+PR) RF(53) = RFH53*PCOR RB(53) = EXP(+(-2.53202D+07)*TM3+( 3.17664D+05)*TM2 & +(-4.53107D+04)*TM1+( 3.58545D+01)+( 2.87730D-04)*TP1 & +(-1.72060D-07)*TP2+( 2.40742D-11)*TP3) RB(53) = RB(53)*PCOR RF(54) = 7.340D+13 RB(54) = EXP(+( 1.05819D+07)*TM3+(-1.37230D+05)*TM2 & +(-4.39297D+04)*TM1+( 3.19804D+01)+( 9.01373D-04)*TP1 & +(-1.92155D-07)*TP2+( 1.76229D-11)*TP3) RFH55 = 5.400D+11 * EXP( 0.454*ALOGT - 3.60000*RTR ) RFL55 = 1.270D+32 * EXP( -4.820*ALOGT - 6.53000*RTR ) PR = RFL55*XM(55)/RFH55 PRLOG = DLOG10(PR) F4 = 0.7187 F5 = 103.0 F6 = 1291. F7 = 4160. FC = (1.0-F4)*EXP(-T/F5)+F4*EXP(-T/F6)+EXP(-F7/T) FCLOG = DLOG10(FC) CPRLOG = PRLOG - (0.4 + 0.67*FCLOG) X = CPRLOG/(0.75-1.27*FCLOG-0.14*CPRLOG) F = 10.0**( FCLOG/(1.0+X*X) ) PCOR = PR*F/(1.0+PR) RF(55) = RFH55*PCOR RB(55) = EXP(+( 7.46324D+06)*TM3+(-5.82414D+04)*TM2 & +(-1.64005D+04)*TM1+( 2.94363D+01)+( 3.63048D-04)*TP1 & +(-1.20138D-07)*TP2+( 1.32405D-11)*TP3) RB(55) = RB(55)*PCOR RFH56 = 5.400D+11 * EXP( 0.454*ALOGT - 2.60000*RTR ) RFL56 = 2.200D+30 * EXP( -4.800*ALOGT - 5.56000*RTR ) PR = RFL56*XM(56)/RFH56 PRLOG = DLOG10(PR) F4 = 0.7580 F5 = 94.00 F6 = 1555. F7 = 4200. FC = (1.0-F4)*EXP(-T/F5)+F4*EXP(-T/F6)+EXP(-F7/T) FCLOG = DLOG10(FC) CPRLOG = PRLOG - (0.4 + 0.67*FCLOG) X = CPRLOG/(0.75-1.27*FCLOG-0.14*CPRLOG) F = 10.0**( FCLOG/(1.0+X*X) ) PCOR = PR*F/(1.0+PR) RF(56) = RFH56*PCOR RB(56) = EXP(+(-1.97404D+07)*TM3+( 2.60636D+05)*TM2 & +(-1.34499D+04)*TM1+( 3.31518D+01)+(-1.83310D-04)*TP1 & +( 1.27130D-08)*TP2+(-3.33405D-12)*TP3) RB(56) = RB(56)*PCOR RF(57) = 2.300D+10 * EXP( 1.050*ALOGT - 3.27500*RTR ) RB(57) = EXP(+( 1.05855D+07)*TM3+(-1.31380D+05)*TM2 & +(-9.42217D+03)*TM1+( 2.67734D+01)+( 1.53255D-03)*TP1 & +(-2.54337D-07)*TP2+( 2.04006D-11)*TP3) RFH58 = 1.800D+13 RFL58 = 3.000D+31 * EXP( -4.800*ALOGT - 3.30000*RTR ) PR = RFL58*XM(58)/RFH58 PRLOG = DLOG10(PR) F4 = 0.7679 F5 = 338.0 F6 = 1812. F7 = 5081. FC = (1.0-F4)*EXP(-T/F5)+F4*EXP(-T/F6)+EXP(-F7/T) FCLOG = DLOG10(FC) CPRLOG = PRLOG - (0.4 + 0.67*FCLOG) X = CPRLOG/(0.75-1.27*FCLOG-0.14*CPRLOG) F = 10.0**( FCLOG/(1.0+X*X) ) PCOR = PR*F/(1.0+PR) RF(58) = RFH58*PCOR RB(58) = EXP(+(-3.20005D+07)*TM3+( 3.97870D+05)*TM2 & +(-4.99350D+04)*TM1+( 3.58494D+01)+(-4.57780D-05)*TP1 & +(-1.15511D-07)*TP2+( 1.89091D-11)*TP3) RB(58) = RB(58)*PCOR RF(59) = 2.000D+13 RB(59) = EXP(+(-1.37698D+07)*TM3+( 1.37822D+05)*TM2 & +(-3.79396D+04)*TM1+( 3.19472D+01)+( 5.93202D-04)*TP1 & +(-1.36790D-07)*TP2+( 1.41905D-11)*TP3) RF(60) = 1.200D+13 RB(60) = EXP(+(-6.77346D+05)*TM3+(-2.42745D+04)*TM2 & +(-1.77373D+03)*TM1+( 2.67779D+01)+( 1.37544D-03)*TP1 & +(-3.56806D-07)*TP2+( 3.94676D-11)*TP3) RF(61) = 6.000D+12 RB(61) = EXP(+(-4.43259D+06)*TM3+( 3.43570D+03)*TM2 & +(-1.58107D+03)*TM1+( 2.58398D+01)+( 1.82492D-03)*TP1 & +(-4.46096D-07)*TP2+( 4.74746D-11)*TP3) RFH62 = 5.000D+13 RFL62 = 8.600D+28 * EXP( -4.000*ALOGT - 3.02500*RTR ) PR = RFL62*XM(62)/RFH62 PRLOG = DLOG10(PR) F4 = 0.8902 F5 = 144.0 F6 = 2838. F7 = 0.4557E+05 FC = (1.0-F4)*EXP(-T/F5)+F4*EXP(-T/F6)+EXP(-F7/T) FCLOG = DLOG10(FC) CPRLOG = PRLOG - (0.4 + 0.67*FCLOG) X = CPRLOG/(0.75-1.27*FCLOG-0.14*CPRLOG) F = 10.0**( FCLOG/(1.0+X*X) ) PCOR = PR*F/(1.0+PR) RF(62) = RFH62*PCOR RB(62) = EXP(+(-4.79682D+06)*TM3+( 7.89925D+04)*TM2 & +(-5.23825D+04)*TM1+( 3.31555D+01)+( 5.00579D-04)*TP1 & +(-2.48362D-07)*TP2+( 3.54837D-11)*TP3) RB(62) = RB(62)*PCOR RF(63) = 3.400D+06 * EXP( 1.600*ALOGT ) RB(63) = EXP(+( 1.46711D+07)*TM3+(-1.60211D+05)*TM2 & +(-3.46931D+03)*TM1+( 2.21937D+01)+( 1.83116D-03)*TP1 & +(-3.84850D-07)*TP2+( 4.19185D-11)*TP3) RF(64) = 2.000D+13 RB(64) = EXP(+( 1.34339D+07)*TM3+(-1.81056D+05)*TM2 & +(-4.03870D+04)*TM1+( 2.82317D+01)+( 1.13956D-03)*TP1 & +(-2.69641D-07)*TP2+( 3.07651D-11)*TP3) RF(65) = 3.200D+13 RB(65) = EXP(+( 2.65263D+07)*TM3+(-3.43152D+05)*TM2 & +(-4.22115D+03)*TM1+( 2.40432D+01)+( 1.92180D-03)*TP1 & +(-4.89656D-07)*TP2+( 5.60422D-11)*TP3) RF(66) = 1.600D+13 RB(66) = EXP(+( 2.27711D+07)*TM3+(-3.15442D+05)*TM2 & +(-4.02850D+03)*TM1+( 2.31051D+01)+( 2.37128D-03)*TP1 & +(-5.78947D-07)*TP2+( 6.40491D-11)*TP3) RF(67) = 1.700D+07 * EXP( 2.100*ALOGT - 4.87000*RTR ) RB(67) = EXP(+( 1.28010D+07)*TM3+(-1.55061D+05)*TM2 & +(-6.09539D+03)*TM1+( 2.62451D+01)+( 2.32377D-03)*TP1 & +(-4.00661D-07)*TP2+( 3.45945D-11)*TP3) RF(68) = 4.200D+06 * EXP( 2.100*ALOGT - 4.87000*RTR ) RB(68) = EXP(+(-1.44027D+07)*TM3+( 1.63817D+05)*TM2 & +(-3.64797D+03)*TM1+( 2.85624D+01)+( 1.77741D-03)*TP1 & +(-2.67810D-07)*TP2+( 1.80199D-11)*TP3) RFH69 = 1.000D+17 * EXP( -1.000*ALOGT ) RFL69 = 3.750D+33 * EXP( -4.800*ALOGT - 1.90000*RTR ) PR = RFL69*XM(69)/RFH69 PRLOG = DLOG10(PR) F4 = 0.6464 F5 = 132.0 F6 = 1315. F7 = 5566. FC = (1.0-F4)*EXP(-T/F5)+F4*EXP(-T/F6)+EXP(-F7/T) FCLOG = DLOG10(FC) CPRLOG = PRLOG - (0.4 + 0.67*FCLOG) X = CPRLOG/(0.75-1.27*FCLOG-0.14*CPRLOG) F = 10.0**( FCLOG/(1.0+X*X) ) PCOR = PR*F/(1.0+PR) RF(69) = RFH69*PCOR RB(69) = EXP(+( 4.56282D+06)*TM3+( 4.18635D+03)*TM2 & +(-6.65431D+04)*TM1+( 3.76389D+01)+(-1.15668D-03)*TP1 & +( 2.09648D-07)*TP2+(-2.26519D-11)*TP3) RB(69) = RB(69)*PCOR RFH70 = 5.600D+12 * EXP( - 2.40000*RTR ) RFL70 = 3.800D+40 * EXP( -7.270*ALOGT - 7.22000*RTR ) PR = RFL70*XM(70)/RFH70 PRLOG = DLOG10(PR) F4 = 0.7507 F5 = 98.50 F6 = 1302. F7 = 4167. FC = (1.0-F4)*EXP(-T/F5)+F4*EXP(-T/F6)+EXP(-F7/T) FCLOG = DLOG10(FC) CPRLOG = PRLOG - (0.4 + 0.67*FCLOG) X = CPRLOG/(0.75-1.27*FCLOG-0.14*CPRLOG) F = 10.0**( FCLOG/(1.0+X*X) ) PCOR = PR*F/(1.0+PR) RF(70) = RFH70*PCOR RB(70) = EXP(+(-2.52082D+07)*TM3+( 3.21012D+05)*TM2 & +(-1.98449D+04)*TM1+( 2.98925D+01)+(-1.35317D-04)*TP1 & +(-4.03962D-08)*TP2+( 9.26855D-12)*TP3) RB(70) = RB(70)*PCOR RFH71 = 6.080D+12 * EXP( 0.270*ALOGT - 0.28000*RTR ) RFL71 = 1.400D+30 * EXP( -3.860*ALOGT - 3.32000*RTR ) PR = RFL71*XM(71)/RFH71 PRLOG = DLOG10(PR) F4 = 0.7820 F5 = 207.5 F6 = 2663. F7 = 6095. FC = (1.0-F4)*EXP(-T/F5)+F4*EXP(-T/F6)+EXP(-F7/T) FCLOG = DLOG10(FC) CPRLOG = PRLOG - (0.4 + 0.67*FCLOG) X = CPRLOG/(0.75-1.27*FCLOG-0.14*CPRLOG) F = 10.0**( FCLOG/(1.0+X*X) ) PCOR = PR*F/(1.0+PR) RF(71) = RFH71*PCOR RB(71) = EXP(+(-2.35726D+07)*TM3+( 3.03061D+05)*TM2 & +(-5.71239D+04)*TM1+( 3.73454D+01)+( 3.68631D-05)*TP1 & +(-1.14473D-07)*TP2+( 1.76431D-11)*TP3) RB(71) = RB(71)*PCOR RF(72) = 3.000D+13 RB(72) = EXP(+( 2.24578D+07)*TM3+(-2.86454D+05)*TM2 & +(-3.36014D+04)*TM1+( 3.11510D+01)+( 7.27005D-04)*TP1 & +(-1.45027D-07)*TP2+( 1.09712D-11)*TP3) RFH73 = 1.080D+12 * EXP( 0.454*ALOGT - 1.82000*RTR ) RFL73 = 1.200D+42 * EXP( -7.620*ALOGT - 6.97000*RTR ) PR = RFL73*XM(73)/RFH73 PRLOG = DLOG10(PR) F4 = 0.9753 F5 = 210.0 F6 = 984.0 F7 = 4374. FC = (1.0-F4)*EXP(-T/F5)+F4*EXP(-T/F6)+EXP(-F7/T) FCLOG = DLOG10(FC) CPRLOG = PRLOG - (0.4 + 0.67*FCLOG) X = CPRLOG/(0.75-1.27*FCLOG-0.14*CPRLOG) F = 10.0**( FCLOG/(1.0+X*X) ) PCOR = PR*F/(1.0+PR) RF(73) = RFH73*PCOR RB(73) = EXP(+(-1.67300D+07)*TM3+( 1.91195D+05)*TM2 & +(-1.97183D+04)*TM1+( 3.08448D+01)+( 3.50640D-04)*TP1 & +(-1.54781D-07)*TP2+( 2.06667D-11)*TP3) RB(73) = RB(73)*PCOR RF(74) = 1.325D+06 * EXP( 2.530*ALOGT - 12.24000*RTR ) RB(74) = EXP(+(-1.10989D+06)*TM3+( 9.16469D+03)*TM2 & +(-3.20326D+03)*TM1+( 2.58675D+01)+( 2.80322D-03)*TP1 & +(-5.11949D-07)*TP2+( 4.69485D-11)*TP3) RFH75 = 5.210D+17 * EXP( -0.990*ALOGT - 1.58000*RTR ) RFL75 = 1.990D+41 * EXP( -7.080*ALOGT - 6.68500*RTR ) PR = RFL75*XM(75)/RFH75 PRLOG = DLOG10(PR) F4 = 0.8422 F5 = 125.0 F6 = 2219. F7 = 6882. FC = (1.0-F4)*EXP(-T/F5)+F4*EXP(-T/F6)+EXP(-F7/T) FCLOG = DLOG10(FC) CPRLOG = PRLOG - (0.4 + 0.67*FCLOG) X = CPRLOG/(0.75-1.27*FCLOG-0.14*CPRLOG) F = 10.0**( FCLOG/(1.0+X*X) ) PCOR = PR*F/(1.0+PR) RF(75) = RFH75*PCOR RB(75) = EXP(+(-1.26657D+07)*TM3+( 1.54555D+05)*TM2 & +(-5.14742D+04)*TM1+( 4.03744D+01)+(-9.92101D-04)*TP1 & +( 8.48113D-08)*TP2+(-2.03116D-12)*TP3) RB(75) = RB(75)*PCOR RF(76) = 2.000D+12 RB(76) = EXP(+( 1.04234D+07)*TM3+(-1.11615D+05)*TM2 & +(-3.37261D+04)*TM1+( 2.89292D+01)+( 6.05610D-04)*TP1 & +(-1.02147D-07)*TP2+( 6.76438D-12)*TP3) RF(77) = 1.150D+08 * EXP( 1.900*ALOGT - 7.53000*RTR ) RB(77) = EXP(+( 2.78736D+06)*TM3+(-2.97552D+04)*TM2 & +(-5.92988D+03)*TM1+( 2.58204D+01)+( 2.31452D-03)*TP1 & +(-4.13559D-07)*TP2+( 3.66852D-11)*TP3) RF(78) = 1.000D+14 RB(78) = EXP(+(-2.25001D+07)*TM3+( 1.96976D+05)*TM2 & +(-9.59617D+03)*TM1+( 2.88550D+01)+( 1.67052D-03)*TP1 & +(-4.39047D-07)*TP2+( 5.15413D-11)*TP3) RF(79) = 5.000D+13 * EXP( - 8.00000*RTR ) RB(79) = EXP(+( 1.28178D+07)*TM3+(-1.64508D+05)*TM2 & +(-2.43893D+03)*TM1+( 2.77576D+01)+( 6.98575D-04)*TP1 & +(-8.51090D-08)*TP2+( 2.66409D-12)*TP3) RF(80) = 1.130D+13 * EXP( - 3.42800*RTR ) RB(80) = EXP(+(-1.59395D+07)*TM3+( 1.28264D+05)*TM2 & +(-1.80334D+04)*TM1+( 2.51875D+01)+( 1.74353D-03)*TP1 & +(-4.35699D-07)*TP2+( 4.87500D-11)*TP3) RF(81) = 1.000D+13 RB(81) = EXP(+(-1.11016D+07)*TM3+( 1.42257D+05)*TM2 & +(-1.57232D+04)*TM1+( 2.98741D+01)+( 6.61450D-05)*TP1 & +( 5.98073D-09)*TP2+(-6.75568D-13)*TP3) RFH82 = 4.300D+07 * EXP( 1.500*ALOGT - 79.60000*RTR ) RFL82 = 5.070D+27 * EXP( -3.420*ALOGT - 84.35000*RTR ) PR = RFL82*XM(82)/RFH82 PRLOG = DLOG10(PR) F4 = 0.9320 F5 = 197.0 F6 = 1540. F7 = 0.1030E+05 FC = (1.0-F4)*EXP(-T/F5)+F4*EXP(-T/F6)+EXP(-F7/T) FCLOG = DLOG10(FC) CPRLOG = PRLOG - (0.4 + 0.67*FCLOG) X = CPRLOG/(0.75-1.27*FCLOG-0.14*CPRLOG) F = 10.0**( FCLOG/(1.0+X*X) ) PCOR = PR*F/(1.0+PR) RF(82) = RFH82*PCOR RB(82) = EXP(+(-4.38917D+07)*TM3+( 5.56044D+05)*TM2 & +(-4.22194D+04)*TM1+( 3.25901D+01)+( 2.05416D-04)*TP1 & +(-1.40555D-07)*TP2+( 2.26081D-11)*TP3) RB(82) = RB(82)*PCOR RF(83) = 2.160D+08 * EXP( 1.510*ALOGT - 3.43000*RTR ) RB(83) = EXP(+(-2.18402D+07)*TM3+( 2.73249D+05)*TM2 & +(-1.07968D+04)*TM1+( 3.15026D+01)+( 1.03645D-03)*TP1 & +(-2.38657D-07)*TP2+( 2.64699D-11)*TP3) RFH84 = 7.400D+13 * EXP( -0.370*ALOGT ) RFL84 = 2.300D+18 * EXP( -0.900*ALOGT + 1.70000*RTR ) PR = RFL84*XM(84)/RFH84 PRLOG = DLOG10(PR) F4 = 0.7346 F5 = 94.00 F6 = 1756. F7 = 5182. FC = (1.0-F4)*EXP(-T/F5)+F4*EXP(-T/F6)+EXP(-F7/T) FCLOG = DLOG10(FC) CPRLOG = PRLOG - (0.4 + 0.67*FCLOG) X = CPRLOG/(0.75-1.27*FCLOG-0.14*CPRLOG) F = 10.0**( FCLOG/(1.0+X*X) ) PCOR = PR*F/(1.0+PR) RF(84) = RFH84*PCOR RB(84) = EXP(+(-1.06770D+07)*TM3+( 1.78051D+05)*TM2 & +(-2.63633D+04)*TM1+( 3.64487D+01)+(-1.67315D-03)*TP1 & +( 3.22699D-07)*TP2+(-3.14707D-11)*TP3) RB(84) = RB(84)*PCOR RF(85) = 3.570D+04 * EXP( 2.400*ALOGT + 2.11000*RTR ) RB(85) = EXP(+(-2.73763D+07)*TM3+( 3.45792D+05)*TM2 & +(-9.50175D+03)*TM1+( 2.95716D+01)+( 1.76158D-03)*TP1 & +(-3.65390D-07)*TP2+( 3.74354D-11)*TP3) RF(86) = 2.900D+13 * EXP( + 0.50000*RTR ) RB(86) = EXP(+(-7.98495D+06)*TM3+( 7.36775D+04)*TM2 & +(-3.52878D+04)*TM1+( 3.33863D+01)+( 2.43223D-04)*TP1 & +(-7.36830D-08)*TP2+( 9.12658D-12)*TP3) RF(87) = 1.750D+12 * EXP( - 0.32000*RTR ) RB(87) = EXP(+(-5.94416D+06)*TM3+( 6.32430D+04)*TM2 & +(-1.63103D+04)*TM1+( 2.82097D+01)+( 5.48486D-04)*TP1 & +(-1.37318D-07)*TP2+( 1.26190D-11)*TP3) RF(88) = 5.800D+14 * EXP( - 9.56000*RTR ) RB(88) = EXP(+(-5.94416D+06)*TM3+( 6.32430D+04)*TM2 & +(-2.09600D+04)*TM1+( 3.40131D+01)+( 5.48486D-04)*TP1 & +(-1.37318D-07)*TP2+( 1.26190D-11)*TP3) RF(89) = 5.000D+13 RB(89) = EXP(+(-7.90154D+06)*TM3+( 8.65219D+04)*TM2 & +(-7.83597D+04)*TM1+( 3.57483D+01)+(-5.09576D-04)*TP1 & +( 1.23328D-07)*TP2+(-1.16033D-11)*TP3) RF(90) = 3.000D+13 RB(90) = EXP(+(-1.43838D+07)*TM3+( 1.84489D+05)*TM2 & +(-4.61601D+04)*TM1+( 3.60694D+01)+(-1.33299D-03)*TP1 & +( 3.54980D-07)*TP2+(-3.70069D-11)*TP3) RF(91) = 2.000D+13 RB(91) = EXP(+(-2.77320D+07)*TM3+( 3.41193D+05)*TM2 & +(-4.02322D+04)*TM1+( 3.85562D+01)+(-1.45882D-03)*TP1 & +( 3.12764D-07)*TP2+(-3.02872D-11)*TP3) RF(92) = 1.130D+07 * EXP( 2.000*ALOGT - 3.00000*RTR ) RB(92) = EXP(+(-2.02166D+07)*TM3+( 2.43039D+05)*TM2 & +(-1.20690D+04)*TM1+( 3.05545D+01)+( 1.99349D-03)*TP1 & +(-4.47011D-07)*TP2+( 4.47198D-11)*TP3) RF(93) = 3.000D+13 RB(93) = EXP(+(-1.93497D+07)*TM3+( 2.57892D+05)*TM2 & +(-4.44649D+04)*TM1+( 3.79890D+01)+(-1.40780D-03)*TP1 & +( 3.08472D-07)*TP2+(-3.07324D-11)*TP3) RFH94 = 6.300D+13 RFL94 = 2.700D+38 * EXP( -6.300*ALOGT - 3.10000*RTR ) PR = RFL94*XM(94)/RFH94 PRLOG = DLOG10(PR) F4 = 0.2105 F5 = 83.50 F6 = 5398. F7 = 8370. FC = (1.0-F4)*EXP(-T/F5)+F4*EXP(-T/F6)+EXP(-F7/T) FCLOG = DLOG10(FC) CPRLOG = PRLOG - (0.4 + 0.67*FCLOG) X = CPRLOG/(0.75-1.27*FCLOG-0.14*CPRLOG) F = 10.0**( FCLOG/(1.0+X*X) ) PCOR = PR*F/(1.0+PR) RF(94) = RFH94*PCOR RB(94) = EXP(+(-3.13232D+07)*TM3+( 4.22145D+05)*TM2 & +(-4.81613D+04)*TM1+( 4.04402D+01)+(-1.42122D-03)*TP1 & +( 2.41294D-07)*TP2+(-2.05585D-11)*TP3) RB(94) = RB(94)*PCOR RF(95) = 5.600D+07 * EXP( 1.600*ALOGT - 5.42000*RTR ) RB(95) = EXP(+(-7.90555D+06)*TM3+( 1.03077D+05)*TM2 & +(-7.78939D+03)*TM1+( 2.74931D+01)+( 1.78531D-03)*TP1 & +(-3.45581D-07)*TP2+( 3.29056D-11)*TP3) RF(96) = 2.501D+13 RB(96) = EXP(+(-3.75524D+06)*TM3+( 2.77102D+04)*TM2 & +( 1.92658D+02)*TM1+( 3.06053D+01)+( 4.49484D-04)*TP1 & +(-8.92907D-08)*TP2+( 8.00698D-12)*TP3) RF(97) = 1.000D+08 * EXP( 1.600*ALOGT - 3.12000*RTR ) RB(97) = EXP(+( 1.02131D+07)*TM3+(-1.01577D+05)*TM2 & +(-8.73980D+03)*TM1+( 2.62281D+01)+( 1.79430D-03)*TP1 & +(-2.74764D-07)*TP2+( 2.35065D-11)*TP3) RF(98) = 4.760D+07 * EXP( 1.228*ALOGT - 0.07000*RTR ) RB(98) = EXP(+(-7.75714D+06)*TM3+( 1.52467D+05)*TM2 & +(-1.37529D+04)*TM1+( 3.37009D+01)+(-5.97901D-04)*TP1 & +( 1.75167D-07)*TP2+(-1.94782D-11)*TP3) RF(99) = 5.000D+13 RB(99) = EXP(+( 5.69407D+05)*TM3+(-1.37231D+04)*TM2 & +(-5.20361D+04)*TM1+( 3.36498D+01)+( 7.25298D-04)*TP1 & +(-1.92988D-07)*TP2+( 2.01745D-11)*TP3) RF(100) = 3.430D+09 * EXP( 1.180*ALOGT + 0.44700*RTR ) RB(100) = EXP(+(-4.45301D+05)*TM3+( 5.01807D+03)*TM2 & +(-1.57386D+04)*TM1+( 2.78069D+01)+( 1.46086D-03)*TP1 & +(-2.75646D-07)*TP2+( 2.50114D-11)*TP3) RF(101) = 5.000D+12 RB(101) = EXP(+(-2.37823D+07)*TM3+( 2.61328D+05)*TM2 & +(-4.60460D+04)*TM1+( 3.26142D+01)+( 4.17127D-04)*TP1 & +(-1.37624D-07)*TP2+( 1.67421D-11)*TP3) RF(102) = 5.000D+12 RB(102) = EXP(+( 3.42137D+06)*TM3+(-5.75496D+04)*TM2 & +(-4.84934D+04)*TM1+( 2.88987D+01)+( 9.63484D-04)*TP1 & +(-2.70475D-07)*TP2+( 3.33167D-11)*TP3) RF(103) = 1.440D+06 * EXP( 2.000*ALOGT + 0.84000*RTR ) RB(103) = EXP(+( 3.57178D+06)*TM3+(-4.14713D+04)*TM2 & +(-1.12645D+04)*TM1+( 2.51505D+01)+( 2.06739D-03)*TP1 & +(-3.85744D-07)*TP2+( 3.55621D-11)*TP3) RF(104) = 6.300D+06 * EXP( 2.000*ALOGT - 1.50000*RTR ) RB(104) = EXP(+(-2.36319D+07)*TM3+( 2.77406D+05)*TM2 & +(-9.99464D+03)*TM1+( 3.03418D+01)+( 1.52103D-03)*TP1 & +(-2.52894D-07)*TP2+( 1.89875D-11)*TP3) RF(105) = 2.000D+13 RB(105) = EXP(+( 1.89217D+06)*TM3+( 3.94622D+04)*TM2 & +(-2.57072D+04)*TM1+( 3.62692D+01)+(-1.53477D-03)*TP1 & +( 4.76958D-07)*TP2+(-5.79517D-11)*TP3) RF(106) = 2.180D-04 * EXP( 4.500*ALOGT + 1.00000*RTR ) RB(106) = EXP(+(-4.56540D+07)*TM3+( 5.81427D+05)*TM2 & +(-1.47216D+04)*TM1+( 2.69374D+01)+( 2.32553D-03)*TP1 & +(-3.84252D-07)*TP2+( 3.77156D-11)*TP3) RF(107) = 5.040D+05 * EXP( 2.300*ALOGT - 13.50000*RTR ) RB(107) = EXP(+(-1.73200D+07)*TM3+( 2.21003D+05)*TM2 & +(-4.88993D+03)*TM1+( 3.36121D+01)+( 4.92784D-04)*TP1 & +(-4.37346D-08)*TP2+( 3.54330D-12)*TP3) RF(108) = 3.370D+07 * EXP( 2.000*ALOGT - 14.00000*RTR ) RB(108) = EXP(+(-2.51587D+07)*TM3+( 2.53046D+05)*TM2 & +(-1.48552D+03)*TM1+( 2.83426D+01)+( 2.37529D-03)*TP1 & +(-5.53404D-07)*TP2+( 6.12832D-11)*TP3) RF(109) = 4.830D-04 * EXP( 4.000*ALOGT + 2.00000*RTR ) RB(109) = EXP(+(-5.76771D+07)*TM3+( 6.60107D+05)*TM2 & +(-3.02074D+04)*TM1+( 1.94677D+01)+( 3.66756D-03)*TP1 & +(-7.41202D-07)*TP2+( 7.85456D-11)*TP3) RF(110) = 5.000D+12 RB(110) = EXP(+( 1.24452D+07)*TM3+(-1.62947D+05)*TM2 & +(-4.17078D+04)*TM1+( 3.14126D+01)+( 5.50930D-04)*TP1 & +(-1.45861D-07)*TP2+( 1.35228D-11)*TP3) RF(111) = 3.600D+06 * EXP( 2.000*ALOGT - 2.50000*RTR ) RB(111) = EXP(+(-6.97098D+06)*TM3+( 8.01127D+04)*TM2 & +(-6.06982D+03)*TM1+( 2.53197D+01)+( 2.20156D-03)*TP1 & +(-4.29308D-07)*TP2+( 4.11049D-11)*TP3) RF(112) = 3.540D+06 * EXP( 2.120*ALOGT - 0.87000*RTR ) RB(112) = EXP(+(-8.94839D+06)*TM3+( 1.15568D+05)*TM2 & +(-1.08253D+04)*TM1+( 2.58875D+01)+( 2.31510D-03)*TP1 & +(-4.49042D-07)*TP2+( 4.27216D-11)*TP3) RF(113) = 7.500D+12 * EXP( - 2.00000*RTR ) RB(113) = EXP(+( 2.80532D+06)*TM3+(-4.10017D+04)*TM2 & +(-7.52602D+03)*TM1+( 2.79138D+01)+( 5.22500D-04)*TP1 & +(-8.59425D-08)*TP2+( 5.21569D-12)*TP3) RF(114) = 1.300D+11 * EXP( + 1.63000*RTR ) RB(114) = EXP(+(-2.04080D+06)*TM3+( 1.04345D+04)*TM2 & +(-1.85699D+04)*TM1+( 2.79598D+01)+(-3.05263D-04)*TP1 & +( 6.36354D-08)*TP2+(-3.49240D-12)*TP3) RF(115) = 4.200D+14 * EXP( - 12.00000*RTR ) RB(115) = EXP(+(-2.04080D+06)*TM3+( 1.04345D+04)*TM2 & +(-2.54287D+04)*TM1+( 3.60403D+01)+(-3.05263D-04)*TP1 & +( 6.36354D-08)*TP2+(-3.49240D-12)*TP3) RF(116) = 2.000D+13 RB(116) = EXP(+(-2.09756D+07)*TM3+( 2.21272D+05)*TM2 & +(-5.78282D+04)*TM1+( 3.42208D+01)+(-2.15655D-04)*TP1 & +(-5.98129D-10)*TP2+( 5.49510D-12)*TP3) RF(117) = 1.000D+12 RB(117) = EXP(+(-3.07307D+07)*TM3+( 3.33922D+05)*TM2 & +(-2.93915D+04)*TM1+( 3.30815D+01)+(-2.66276D-04)*TP1 & +(-5.09175D-08)*TP2+( 1.09639D-11)*TP3) RF(118) = 2.000D+13 RB(118) = EXP(+(-1.97699D+07)*TM3+( 2.23231D+05)*TM2 & +(-1.33749D+04)*TM1+( 3.33448D+01)+(-6.78630D-04)*TP1 & +( 1.76294D-07)*TP2+(-2.02599D-11)*TP3) RF(119) = 1.500D+14 * EXP( - 23.60000*RTR ) RB(119) = EXP(+( 8.61804D+06)*TM3+(-8.92309D+04)*TM2 & +(-4.24053D+04)*TM1+( 3.59861D+01)+(-3.40819D-04)*TP1 & +( 5.52138D-08)*TP2+(-3.14730D-12)*TP3) RF(120) = 1.000D+12 * EXP( - 8.00000*RTR ) RB(120) = EXP(+( 1.47417D+07)*TM3+(-1.75242D+05)*TM2 & +(-3.08637D+03)*TM1+( 2.54465D+01)+(-3.51639D-05)*TP1 & +( 4.75221D-08)*TP2+(-6.29872D-12)*TP3) RF(121) = 5.800D+13 * EXP( - 0.57600*RTR ) RB(121) = EXP(+(-1.73760D+06)*TM3+( 7.14168D+02)*TM2 & +(-6.97369D+04)*TM1+( 3.19571D+01)+( 3.24747D-04)*TP1 & +(-1.03744D-07)*TP2+( 1.44719D-11)*TP3) RF(122) = 5.000D+13 RB(122) = EXP(+( 4.32406D+06)*TM3+(-6.66151D+04)*TM2 & +(-3.88236D+04)*TM1+( 3.44598D+01)+(-6.96352D-04)*TP1 & +( 8.97091D-08)*TP2+(-4.74764D-12)*TP3) RF(123) = 5.000D+13 RB(123) = EXP(+( 1.84440D+07)*TM3+(-1.76917D+05)*TM2 & +(-4.97005D+04)*TM1+( 3.58199D+01)+(-9.65138D-04)*TP1 & +( 2.34532D-07)*TP2+(-2.67892D-11)*TP3) RF(124) = 3.300D+13 RB(124) = EXP(+(-8.21982D+06)*TM3+( 9.86815D+04)*TM2 & +(-3.72474D+04)*TM1+( 3.22251D+01)+(-4.98666D-04)*TP1 & +( 1.27909D-07)*TP2+(-1.09317D-11)*TP3) RF(125) = 1.107D+08 * EXP( 1.790*ALOGT - 1.67000*RTR ) RB(125) = EXP(+(-1.94826D+07)*TM3+( 2.56311D+05)*TM2 & +(-8.07975D+02)*TM1+( 3.20096D+01)+( 8.73802D-04)*TP1 & +(-1.50745D-07)*TP2+( 1.78651D-11)*TP3) RF(126) = 5.710D+12 * EXP( + 0.75500*RTR ) RB(126) = EXP(+(-2.31813D+07)*TM3+( 2.96488D+05)*TM2 & +(-3.05702D+04)*TM1+( 3.65759D+01)+(-1.84631D-03)*TP1 & +( 4.44776D-07)*TP2+(-4.33272D-11)*TP3) RF(127) = 4.000D+13 RB(127) = EXP(+( 7.90413D+06)*TM3+(-3.70824D+04)*TM2 & +(-6.54969D+04)*TM1+( 3.97531D+01)+(-1.56376D-03)*TP1 & +( 3.66779D-07)*TP2+(-3.95800D-11)*TP3) RF(128) = 3.000D+13 RB(128) = EXP(+( 8.59572D+04)*TM3+( 7.02744D+04)*TM2 & +(-2.78293D+04)*TM1+( 3.60756D+01)+(-1.61418D-03)*TP1 & +( 4.19057D-07)*TP2+(-4.55410D-11)*TP3) RF(129) = 6.000D+13 RB(129) = EXP(+( 1.41369D+07)*TM3+(-7.17483D+04)*TM2 & +(-3.04427D+04)*TM1+( 3.70705D+01)+(-1.70024D-03)*TP1 & +( 5.10600D-07)*TP2+(-5.68035D-11)*TP3) RFH130 = 5.000D+13 RFL130 = 2.690D+28 * EXP( -3.740*ALOGT - 1.93600*RTR ) PR = RFL130*XM(130)/RFH130 PRLOG = DLOG10(PR) F4 = 0.5757 F5 = 237.0 F6 = 1652. F7 = 5069. FC = (1.0-F4)*EXP(-T/F5)+F4*EXP(-T/F6)+EXP(-F7/T) FCLOG = DLOG10(FC) CPRLOG = PRLOG - (0.4 + 0.67*FCLOG) X = CPRLOG/(0.75-1.27*FCLOG-0.14*CPRLOG) F = 10.0**( FCLOG/(1.0+X*X) ) PCOR = PR*F/(1.0+PR) RF(130) = RFH130*PCOR RB(130) = EXP(+( 5.90560D+06)*TM3+(-3.14793D+02)*TM2 & +(-3.72340D+04)*TM1+( 3.78810D+01)+(-1.69342D-03)*TP1 & +( 3.89094D-07)*TP2+(-4.13447D-11)*TP3) RB(130) = RB(130)*PCOR RF(131) = 3.400D+12 * EXP( - 0.69000*RTR ) RB(131) = EXP(+(-1.62454D+07)*TM3+( 1.53799D+05)*TM2 & +(-3.35739D+04)*TM1+( 2.62121D+01)+( 2.50995D-04)*TP1 & +(-1.35963D-08)*TP2+( 1.92271D-12)*TP3) RF(132) = 9.460D+13 * EXP( + 0.51500*RTR ) RB(132) = EXP(+( 2.52301D+07)*TM3+(-2.43100D+05)*TM2 & +(-3.73564D+04)*TM1+( 3.74812D+01)+(-1.39291D-03)*TP1 & +( 3.78509D-07)*TP2+(-4.28569D-11)*TP3) RF(133) = 5.000D+13 RB(133) = EXP(+(-6.21365D+06)*TM3+( 7.65925D+04)*TM2 & +(-7.93258D+04)*TM1+( 3.56224D+01)+( 1.57778D-04)*TP1 & +(-7.65605D-08)*TP2+( 1.15161D-11)*TP3) RF(134) = 1.320D+13 * EXP( - 1.50000*RTR ) RB(134) = EXP(+(-4.19319D+06)*TM3+( 3.55954D+04)*TM2 & +(-3.82535D+04)*TM1+( 2.92518D+01)+( 5.44439D-05)*TP1 & +(-1.67115D-08)*TP2+( 2.68878D-12)*TP3) RF(135) = 5.000D+05 * EXP( 2.000*ALOGT - 7.23000*RTR ) RB(135) = EXP(+(-3.03054D+07)*TM3+( 3.77431D+05)*TM2 & +(-8.98192D+03)*TM1+( 2.99807D+01)+( 9.29414D-04)*TP1 & +(-2.22250D-07)*TP2+( 2.66697D-11)*TP3) RF(136) = 3.200D+13 RB(136) = EXP(+( 1.33659D+07)*TM3+(-1.15884D+05)*TM2 & +(-6.66726D+04)*TM1+( 3.82034D+01)+(-1.00020D-03)*TP1 & +( 2.35600D-07)*TP2+(-2.90916D-11)*TP3) RF(137) = 4.000D+13 RB(137) = EXP(+(-1.31596D+07)*TM3+( 2.33201D+05)*TM2 & +(-3.35769D+04)*TM1+( 4.04543D+01)+(-1.82225D-03)*TP1 & +( 4.01354D-07)*TP2+(-4.19261D-11)*TP3) RF(138) = 2.460D+06 * EXP( 2.000*ALOGT - 8.27000*RTR ) RB(138) = EXP(+( 2.45288D+06)*TM3+(-6.31961D+03)*TM2 & +(-7.54677D+03)*TM1+( 2.64583D+01)+( 1.61499D-03)*TP1 & +(-2.44182D-07)*TP2+( 2.22808D-11)*TP3) RFH139 = 8.100D+11 * EXP( 0.500*ALOGT - 4.51000*RTR ) RFL139 = 2.690D+33 * EXP( -5.110*ALOGT - 7.09500*RTR ) PR = RFL139*XM(139)/RFH139 PRLOG = DLOG10(PR) F4 = 0.5907 F5 = 275.0 F6 = 1226. F7 = 5185. FC = (1.0-F4)*EXP(-T/F5)+F4*EXP(-T/F6)+EXP(-F7/T) FCLOG = DLOG10(FC) CPRLOG = PRLOG - (0.4 + 0.67*FCLOG) X = CPRLOG/(0.75-1.27*FCLOG-0.14*CPRLOG) F = 10.0**( FCLOG/(1.0+X*X) ) PCOR = PR*F/(1.0+PR) RF(139) = RFH139*PCOR RB(139) = EXP(+(-5.36694D+06)*TM3+( 1.34975D+05)*TM2 & +(-4.25853D+04)*TM1+( 3.96140D+01)+(-1.42693D-03)*TP1 & +( 2.64275D-07)*TP2+(-2.56005D-11)*TP3) RB(139) = RB(139)*PCOR RF(140) = 3.000D+13 RB(140) = EXP(+(-2.32097D+07)*TM3+( 2.84245D+05)*TM2 & +(-4.69001D+04)*TM1+( 3.36662D+01)+(-5.66037D-06)*TP1 & +(-6.27118D-08)*TP2+( 1.10332D-11)*TP3) RF(141) = 1.500D+13 * EXP( - 0.60000*RTR ) RB(141) = EXP(+( 8.38231D+06)*TM3+(-8.33007D+04)*TM2 & +(-4.53466D+03)*TM1+( 2.93664D+01)+( 5.10252D-05)*TP1 & +(-4.29182D-09)*TP2+(-4.45233D-13)*TP3) RF(142) = 2.800D+13 RB(142) = EXP(+( 1.75215D+07)*TM3+(-2.19493D+05)*TM2 & +(-3.34226D+04)*TM1+( 2.84270D+01)+( 4.15154D-04)*TP1 & +(-2.77347D-08)*TP2+(-3.73257D-13)*TP3) RF(143) = 1.200D+13 RB(143) = EXP(+( 4.75853D+06)*TM3+(-6.14284D+04)*TM2 & +(-9.37675D+04)*TM1+( 3.02906D+01)+( 8.30767D-04)*TP1 & +(-2.13992D-07)*TP2+( 2.24181D-11)*TP3) RF(144) = 7.000D+13 RB(144) = EXP(+(-6.25728D+06)*TM3+( 9.57963D+04)*TM2 & +(-8.29904D+03)*TM1+( 3.41778D+01)+(-6.25560D-04)*TP1 & +( 8.84572D-08)*TP2+(-5.45537D-12)*TP3) RFH145 = 2.000D+13 RFL145 = 2.700D+38 * EXP( -6.300*ALOGT - 3.10000*RTR ) PR = RFL145*XM(145)/RFH145 PRLOG = DLOG10(PR) F4 = 0.1507 F5 = 134.0 F6 = 2383. F7 = 7265. FC = (1.0-F4)*EXP(-T/F5)+F4*EXP(-T/F6)+EXP(-F7/T) FCLOG = DLOG10(FC) CPRLOG = PRLOG - (0.4 + 0.67*FCLOG) X = CPRLOG/(0.75-1.27*FCLOG-0.14*CPRLOG) F = 10.0**( FCLOG/(1.0+X*X) ) PCOR = PR*F/(1.0+PR) RF(145) = RFH145*PCOR RB(145) = EXP(+(-2.75679D+07)*TM3+( 3.94434D+05)*TM2 & +(-4.83540D+04)*TM1+( 3.95377D+01)+(-1.87070D-03)*TP1 & +( 3.30585D-07)*TP2+(-2.85654D-11)*TP3) RB(145) = RB(145)*PCOR RF(146) = 3.000D+13 RB(146) = EXP(+( 8.38231D+06)*TM3+(-8.33007D+04)*TM2 & +(-4.23273D+03)*TM1+( 3.00595D+01)+( 5.10252D-05)*TP1 & +(-4.29183D-09)*TP2+(-4.45232D-13)*TP3) RF(147) = 1.200D+13 * EXP( + 0.57000*RTR ) RB(147) = EXP(+(-4.77731D+06)*TM3+( 1.49900D+05)*TM2 & +(-3.75227D+04)*TM1+( 3.82776D+01)+(-1.77123D-03)*TP1 & +( 3.97062D-07)*TP2+(-4.23714D-11)*TP3) RF(148) = 1.600D+13 * EXP( + 0.57000*RTR ) RB(148) = EXP(+( 2.65010D+07)*TM3+(-2.87954D+05)*TM2 & +(-6.05370D+03)*TM1+( 2.75861D+01)+( 6.00154D-05)*TP1 & +( 6.65251D-08)*TP2+(-9.84434D-12)*TP3) RF(149) = 9.000D+12 RB(149) = EXP(+( 8.38231D+06)*TM3+(-8.33007D+04)*TM2 & +(-4.23273D+03)*TM1+( 2.88556D+01)+( 5.10252D-05)*TP1 & +(-4.29183D-09)*TP2+(-4.45232D-13)*TP3) RF(150) = 7.000D+12 RB(150) = EXP(+( 8.38231D+06)*TM3+(-8.33007D+04)*TM2 & +(-4.23273D+03)*TM1+( 2.86042D+01)+( 5.10252D-05)*TP1 & +(-4.29183D-09)*TP2+(-4.45233D-13)*TP3) RF(151) = 1.400D+13 RB(151) = EXP(+(-2.12114D+07)*TM3+( 2.27202D+05)*TM2 & +(-3.15315D+04)*TM1+( 2.95469D+01)+( 1.76189D-04)*TP1 & +(-6.01038D-08)*TP2+( 8.19717D-12)*TP3) RF(152) = 4.000D+13 * EXP( + 0.55000*RTR ) RB(152) = EXP(+( 1.14126D+07)*TM3+(-1.22376D+05)*TM2 & +(-8.94944D+03)*TM1+( 2.79701D+01)+( 1.63261D-04)*TP1 & +(-2.58532D-08)*TP2+( 1.13398D-12)*TP3) RF(153) = 2.675D+13 * EXP( - 28.80000*RTR ) RB(153) = EXP(+(-2.03624D+07)*TM3+( 2.57344D+05)*TM2 & +(-1.35883D+03)*TM1+( 3.40314D+01)+(-1.08747D-03)*TP1 & +( 2.62585D-07)*TP2+(-2.99669D-11)*TP3) RF(154) = 3.600D+10 * EXP( - 8.94000*RTR ) RB(154) = EXP(+(-8.36362D+06)*TM3+( 9.20040D+04)*TM2 & +(-3.08276D+04)*TM1+( 2.42951D+01)+( 4.16309D-05)*TP1 & +(-2.04976D-08)*TP2+( 3.93026D-12)*TP3) RF(155) = 2.450D+04 * EXP( 2.470*ALOGT - 5.18000*RTR ) RB(155) = EXP(+(-4.80371D+07)*TM3+( 5.68430D+05)*TM2 & +(-1.41856D+04)*TM1+( 2.99684D+01)+( 2.02240D-03)*TP1 & +(-5.03576D-07)*TP2+( 5.35810D-11)*TP3) RFH156 = 2.120D+16 * EXP( -0.970*ALOGT - 0.62000*RTR ) RFL156 = 1.770D+50 * EXP( -9.670*ALOGT - 6.22000*RTR ) PR = RFL156*XM(156)/RFH156 PRLOG = DLOG10(PR) F4 = 0.5325 F5 = 151.0 F6 = 1038. F7 = 4970. FC = (1.0-F4)*EXP(-T/F5)+F4*EXP(-T/F6)+EXP(-F7/T) FCLOG = DLOG10(FC) CPRLOG = PRLOG - (0.4 + 0.67*FCLOG) X = CPRLOG/(0.75-1.27*FCLOG-0.14*CPRLOG) F = 10.0**( FCLOG/(1.0+X*X) ) PCOR = PR*F/(1.0+PR) RF(156) = RFH156*PCOR RB(156) = EXP(+(-2.17658D+07)*TM3+( 3.22258D+05)*TM2 & +(-4.67883D+04)*TM1+( 4.25668D+01)+(-2.72732D-03)*TP1 & +( 4.92413D-07)*TP2+(-4.53947D-11)*TP3) RB(156) = RB(156)*PCOR RF(157) = 4.990D+12 * EXP( 0.100*ALOGT - 10.60000*RTR ) RB(157) = EXP(+(-9.72674D+06)*TM3+( 1.75635D+05)*TM2 & +(-1.18240D+03)*TM1+( 3.51761D+01)+(-1.67098D-03)*TP1 & +( 3.95001D-07)*TP2+(-4.20964D-11)*TP3) RF(158) = 2.648D+13 RB(158) = EXP(+(-2.21763D+07)*TM3+( 2.46521D+05)*TM2 & +(-4.58882D+04)*TM1+( 3.60767D+01)+( 2.15798D-04)*TP1 & +(-1.70223D-07)*TP2+( 2.20119D-11)*TP3) RF(159) = 3.320D+03 * EXP( 2.810*ALOGT - 5.86000*RTR ) RB(159) = EXP(+(-3.59586D+07)*TM3+( 4.26904D+05)*TM2 & +(-1.38056D+04)*TM1+( 2.80950D+01)+( 2.26025D-03)*TP1 & +(-5.09604D-07)*TP2+( 5.26679D-11)*TP3) RF(160) = 3.000D+07 * EXP( 1.500*ALOGT - 9.94000*RTR ) RB(160) = EXP(+(-1.52575D+07)*TM3+( 1.69189D+05)*TM2 & +(-1.02220D+04)*TM1+( 2.78526D+01)+( 1.15639D-03)*TP1 & +(-2.84229D-07)*TP2+( 2.94795D-11)*TP3) RF(161) = 1.000D+07 * EXP( 1.500*ALOGT - 9.94000*RTR ) RB(161) = EXP(+(-4.24612D+07)*TM3+( 4.88067D+05)*TM2 & +(-7.77457D+03)*TM1+( 3.04695D+01)+( 6.10033D-04)*TP1 & +(-1.51378D-07)*TP2+( 1.29049D-11)*TP3) RF(162) = 2.270D+05 * EXP( 2.000*ALOGT - 9.20000*RTR ) RB(162) = EXP(+(-2.97167D+07)*TM3+( 3.40357D+05)*TM2 & +(-3.29349D+03)*TM1+( 2.56184D+01)+( 1.69206D-03)*TP1 & +(-4.06542D-07)*TP2+( 4.29423D-11)*TP3) RF(163) = 6.140D+06 * EXP( 1.740*ALOGT - 10.45000*RTR ) RB(163) = EXP(+(-2.87176D+07)*TM3+( 3.38129D+05)*TM2 & +(-9.25559D+03)*TM1+( 2.69191D+01)+( 1.50046D-03)*TP1 & +(-3.66427D-07)*TP2+( 3.85398D-11)*TP3) RF(164) = 2.244D+18 * EXP( -1.000*ALOGT - 17.00000*RTR ) RB(164) = EXP(+( 2.11652D+07)*TM3+(-2.70955D+05)*TM2 & +( 3.92816D+02)*TM1+( 3.48570D+01)+(-4.93315D-04)*TP1 & +( 1.50768D-07)*TP2+(-1.84567D-11)*TP3) RF(165) = 1.870D+17 * EXP( -1.000*ALOGT - 17.00000*RTR ) RB(165) = EXP(+( 2.11652D+07)*TM3+(-2.70955D+05)*TM2 & +( 3.92816D+02)*TM1+( 3.23721D+01)+(-4.93315D-04)*TP1 & +( 1.50768D-07)*TP2+(-1.84567D-11)*TP3) RF(166) = 7.600D+12 * EXP( - 0.40000*RTR ) RB(166) = EXP(+( 8.55436D+06)*TM3+(-8.74006D+04)*TM2 & +(-1.66980D+04)*TM1+( 2.93780D+01)+( 4.82074D-04)*TP1 & +(-1.19305D-07)*TP2+( 1.10480D-11)*TP3) RF(167) = 1.800D+13 * EXP( - 0.90000*RTR ) RB(167) = EXP(+(-1.57974D+07)*TM3+( 1.87651D+05)*TM2 & +(-1.09595D+04)*TM1+( 3.15071D+01)+( 1.73904D-04)*TP1 & +(-6.39407D-08)*TP2+( 7.61556D-12)*TP3) RF(168) = 4.280D-13 * EXP( 7.600*ALOGT + 3.53000*RTR ) RB(168) = EXP(+(-4.81236D+07)*TM3+( 6.22442D+05)*TM2 & +(-1.60316D+04)*TM1+( 2.04228D+01)+( 6.82306D-03)*TP1 & +(-1.39379D-06)*TP2+( 1.44574D-10)*TP3) RF(169) = 5.000D+13 * EXP( - 1.50000*RTR ) RB(169) = EXP(+(-1.64188D+07)*TM3+( 1.88732D+05)*TM2 & +(-7.77896D+04)*TM1+( 3.18721D+01)+( 2.41220D-04)*TP1 & +( 1.69071D-08)*TP2+(-4.16689D-12)*TP3) RF(170) = 4.070D+05 * EXP( 2.400*ALOGT - 0.20000*RTR ) RB(170) = EXP(+(-1.93186D+07)*TM3+( 3.06796D+05)*TM2 & +(-1.65767D+04)*TM1+( 3.38526D+01)+( 9.81827D-04)*TP1 & +(-1.40427D-07)*TP2+( 1.09642D-11)*TP3) RF(171) = 3.980D+12 * EXP( + 0.24000*RTR ) RB(171) = EXP(+(-2.28333D+07)*TM3+( 2.06219D+05)*TM2 & +(-4.45389D+04)*TM1+( 2.89944D+01)+( 3.22826D-04)*TP1 & +(-8.45743D-08)*TP2+( 1.24644D-11)*TP3) RFH172 = 8.000D+12 * EXP( 0.440*ALOGT - 88.77000*RTR ) RFL172 = 7.000D+50 * EXP( -9.310*ALOGT - 99.86000*RTR ) PR = RFL172*XM(172)/RFH172 PRLOG = DLOG10(PR) F4 = 0.7345 F5 = 180.0 F6 = 1035. F7 = 5417. FC = (1.0-F4)*EXP(-T/F5)+F4*EXP(-T/F6)+EXP(-F7/T) FCLOG = DLOG10(FC) CPRLOG = PRLOG - (0.4 + 0.67*FCLOG) X = CPRLOG/(0.75-1.27*FCLOG-0.14*CPRLOG) F = 10.0**( FCLOG/(1.0+X*X) ) PCOR = PR*F/(1.0+PR) RF(172) = RFH172*PCOR RB(172) = EXP(+( 4.04690D+07)*TM3+(-5.19107D+05)*TM2 & +(-2.17424D+04)*TM1+( 2.67435D+01)+( 1.26027D-03)*TP1 & +(-1.42379D-07)*TP2+( 4.57447D-12)*TP3) RB(172) = RB(172)*PCOR RF(173) = 8.400D+11 * EXP( - 3.87500*RTR ) RB(173) = EXP(+( 8.39586D+06)*TM3+(-6.17862D+04)*TM2 & +(-8.24308D+03)*TM1+( 2.77271D+01)+( 1.86311D-04)*TP1 & +(-2.92971D-08)*TP2+( 1.89402D-13)*TP3) RF(174) = 1.600D+12 * EXP( - 0.85400*RTR ) RB(174) = EXP(+(-4.97861D+06)*TM3+(-2.25177D+04)*TM2 & +(-4.34485D+04)*TM1+( 2.21836D+01)+( 2.08567D-03)*TP1 & +(-4.66782D-07)*TP2+( 5.11680D-11)*TP3) RF(175) = 1.000D+13 RB(175) = EXP(+(-1.21193D+07)*TM3+( 7.69072D+04)*TM2 & +(-4.20918D+04)*TM1+( 2.76750D+01)+( 1.85120D-03)*TP1 & +(-4.65655D-07)*TP2+( 5.28608D-11)*TP3) RETURN END SUBROUTINE STEADYE( ABV, DEN, RF, RB, XM, ADJ, CONMX, SMALL, & LBOUND, LITER, XH2, XH, XO, XO2, XOH, XH2O, & XHO2, XH2O2, XC, XCH, XCH2, XCH2S, XCH3, XCH4, & XCO, XCO2, XHCO, XCH2O, XCH2OH, XCH3O, XCH3OH, & XC2H, XC2H2, XC2H3, XC2H4, XC2H5, XC2H6, & XHCCO, XCH2CO, XHCCOH, XN2 ) C*****double precision IMPLICIT DOUBLE PRECISION (A-H, O-Z), INTEGER (I-N) C*****END double precision C*****single precision C IMPLICIT REAL (A-H, O-Z), INTEGER (I-N) C*****END single precision DIMENSION ABV(*), DEN(*), RF(*), RB(*), XM(*) LOGICAL LBOUND(*), LITER C STEADY-STATE EXPRESSION FOR C2H ABV(1) = RB(69)*XC2H2 +RF(22)*XO*XC2H2 +RB(170)*XH*XC2H2 +RB(105) & *XH*XHCCO +RB(20)*XCH*XCO +RF(122)*XC*XCH2 +RF(108)*XOH & *XC2H2 +RB(169)*XCO*XHCO DEN(1) = RF(69)*XH +RB(22)*XOH +RF(170)*XH2 +RF(105)*XOH +RF(20) & *XO +RB(122)*XH +RB(108)*XH2O +RF(169)*XO2 IF( LITER .OR. .NOT. LBOUND(1) ) THEN IF(DEN(1) .LT. 1.0) DEN(1) = MAX( ADJ*ABV(1), DEN(1), SMALL ) VOLD = XC2H XC2H = ABV(1)/DEN(1) DIFF = ABS( (XC2H-VOLD)/MAX(XC2H,VOLD,SMALL) ) CONMX = MAX( CONMX, DIFF ) ENDIF C STEADY-STATE EXPRESSION FOR C2H5 ABV(2) = RF(163)*XCH3*XC2H6 +RF(152)*XCH2S*XC2H6 +RB(173)*XHO2 & *XC2H4 +RB(75)*XC2H6 +RB(76)*XH2*XC2H4 +RF(157)*XCH3*XCH3 & +RF(27)*XO*XC2H6 +RB(26)*XCH3*XCH2O +RF(112)*XOH*XC2H6 + & RF(77)*XH*XC2H6 +RF(73)*XH*XC2H4 DEN(2) = RB(163)*XCH4 +RB(152)*XCH3 +RF(173)*XO2 +RF(75)*XH + & RF(76)*XH +RB(157)*XH +RB(27)*XOH +RF(26)*XO +RB(112) & *XH2O +RB(77)*XH2 +RB(73) IF( LITER .OR. .NOT. LBOUND(2) ) THEN IF(DEN(2) .LT. 1.0) DEN(2) = MAX( ADJ*ABV(2), DEN(2), SMALL ) VOLD = XC2H5 XC2H5 = ABV(2)/DEN(2) DIFF = ABS( (XC2H5-VOLD)/MAX(XC2H5,VOLD,SMALL) ) CONMX = MAX( CONMX, DIFF ) ENDIF C STEADY-STATE EXPRESSION FOR HCCOH ABV(3) = RF(107)*XOH*XC2H2 +RB(81)*XCH2CO*XH DEN(3) = RB(107)*XH +RF(81)*XH IF( LITER .OR. .NOT. LBOUND(3) ) THEN IF(DEN(3) .LT. 1.0) DEN(3) = MAX( ADJ*ABV(3), DEN(3), SMALL ) VOLD = XHCCOH XHCCOH = ABV(3)/DEN(3) DIFF = ABS( (XHCCOH-VOLD)/MAX(XHCCOH,VOLD,SMALL) ) CONMX = MAX( CONMX, DIFF ) ENDIF C STEADY-STATE EXPRESSION FOR C2H3 ABV(4) = RF(162)*XCH3*XC2H4 +RF(140)*XCH2*XHCCO +RB(71)*XC2H4 + & RB(110)*XH2O*XC2H2 +RB(171)*XHCO*XCH2O +RB(24)*XH*XCH2CO & +RF(128)*XCH*XCH3 +RB(72)*XH2*XC2H2 +RF(111)*XOH*XC2H4 + & RF(74)*XH*XC2H4 +RF(70)*XH*XC2H2 DEN(4) = RB(162)*XCH4 +RB(140)*XCO +RF(71)*XH +RF(110)*XOH + & RF(171)*XO2 +RF(24)*XO +RB(128)*XH +RF(72)*XH +RB(111) & *XH2O +RB(74)*XH2 +RB(70) IF( LITER .OR. .NOT. LBOUND(4) ) THEN IF(DEN(4) .LT. 1.0) DEN(4) = MAX( ADJ*ABV(4), DEN(4), SMALL ) VOLD = XC2H3 XC2H3 = ABV(4)/DEN(4) DIFF = ABS( (XC2H3-VOLD)/MAX(XC2H3,VOLD,SMALL) ) CONMX = MAX( CONMX, DIFF ) ENDIF C STEADY-STATE EXPRESSION FOR CH2S ABV(5) = RB(152)*XCH3*XC2H5 +RB(147)*XH*XC2H4 +RF(66)*XH*XCH3O + & RB(148)*XCH3*XCH3 +RF(78)*XH*XHCCO +RB(145)*XCH3OH +RB(9) & *XH*XHCO +RB(8)*XH2*XCO +RF(61)*XH*XCH2OH +RB(149)*XCH2 & *XCO +RB(50)*XH2*XCH +RB(93)*XH*XCH2O +RB(143)*XH2O*XCO + & RB(150)*XCH2*XCO2 +RB(142)*XH*XOH*XCO +RB(151)*XCO*XCH2O & +RB(144)*XH*XCH3 +RB(146)*XCH2*XH2O +RB(141)*XCH2*XN2 + & RF(96)*XOH*XCH3 DEN(5) = RF(152)*XC2H6 +RF(147)*XCH3 +RB(66)*XH2O +RF(148)*XCH4 + & RB(78)*XCO +RF(145)*XH2O +RF(9)*XO +RF(8)*XO +RB(61)*XH2O & +RF(149)*XCO +RF(50)*XH +RF(93)*XOH +RF(143)*XO2 +RF(150) & *XCO2 +RF(142)*XO2 +RF(151)*XCO2 +RF(144)*XH2 +RF(146) & *XH2O +RF(141)*XN2 +RB(96)*XH2O IF( LITER .OR. .NOT. LBOUND(5) ) THEN IF(DEN(5) .LT. 1.0) DEN(5) = MAX( ADJ*ABV(5), DEN(5), SMALL ) VOLD = XCH2S XCH2S = ABV(5)/DEN(5) DIFF = ABS( (XCH2S-VOLD)/MAX(XCH2S,VOLD,SMALL) ) CONMX = MAX( CONMX, DIFF ) ENDIF C STEADY-STATE EXPRESSION FOR CH2OH ABV(6) = RF(160)*XCH3*XCH3OH +RB(58)*XCH3OH +RF(63)*XCH3O*XH + & RB(16)*XOH*XCH2O +RB(101)*XH2O*XCH2O +RF(103)*XOH*XCH3OH & +RB(61)*XH2O*XCH2S +RF(18)*XO*XCH3OH +RB(59)*XH2*XCH2O + & RB(167)*XHO2*XCH2O +RF(67)*XH*XCH3OH +RB(60)*XOH*XCH3 + & RF(55)*XH*XCH2O DEN(6) = RB(160)*XCH4 +RF(58)*XH +RB(63)*XH +RF(16)*XO +RF(101) & *XOH +RB(103)*XH2O +RF(61)*XH +RB(18)*XOH +RF(59)*XH + & RF(167)*XO2 +RB(67)*XH2 +RF(60)*XH +RB(55) IF( LITER .OR. .NOT. LBOUND(6) ) THEN IF(DEN(6) .LT. 1.0) DEN(6) = MAX( ADJ*ABV(6), DEN(6), SMALL ) VOLD = XCH2OH XCH2OH = ABV(6)/DEN(6) DIFF = ABS( (XCH2OH-VOLD)/MAX(XCH2OH,VOLD,SMALL) ) CONMX = MAX( CONMX, DIFF ) ENDIF C STEADY-STATE EXPRESSION FOR CH3OH ABV(7) = RB(161)*XCH4*XCH3O +RB(160)*XCH4*XCH2OH +RF(62)*XH*XCH3O & +RF(58)*XH*XCH2OH +RB(19)*XOH*XCH3O +RF(145)*XH2O*XCH2S + & RB(103)*XH2O*XCH2OH +RB(68)*XH2*XCH3O +RB(18)*XOH*XCH2OH & +RB(104)*XH2O*XCH3O +RB(67)*XH2*XCH2OH +RF(94)*XOH*XCH3 DEN(7) = RF(161)*XCH3 +RF(160)*XCH3 +RB(62) +RB(58) +RF(19)*XO + & RB(145) +RF(103)*XOH +RF(68)*XH +RF(18)*XO +RF(104)*XOH + & RF(67)*XH +RB(94) IF( LITER .OR. .NOT. LBOUND(7) ) THEN IF(DEN(7) .LT. 1.0) DEN(7) = MAX( ADJ*ABV(7), DEN(7), SMALL ) VOLD = XCH3OH XCH3OH = ABV(7)/DEN(7) DIFF = ABS( (XCH3OH-VOLD)/MAX(XCH3OH,VOLD,SMALL) ) CONMX = MAX( CONMX, DIFF ) ENDIF C STEADY-STATE EXPRESSION FOR CH3O ABV(8) = RF(161)*XCH3*XCH3OH +RB(62)*XCH3OH +RB(63)*XCH2OH*XH + & RB(17)*XOH*XCH2O +RB(102)*XH2O*XCH2O +RB(66)*XH2O*XCH2S + & RB(65)*XOH*XCH3 +RB(64)*XH2*XCH2O +RF(118)*XHO2*XCH3 + & RF(153)*XO2*XCH3 +RB(168)*XHO2*XCH2O +RF(19)*XO*XCH3OH + & RF(68)*XH*XCH3OH +RF(104)*XOH*XCH3OH +RF(56)*XH*XCH2O DEN(8) = RB(161)*XCH4 +RF(62)*XH +RF(63)*XH +RF(17)*XO +RF(102) & *XOH +RF(66)*XH +RF(65)*XH +RF(64)*XH +RB(118)*XOH + & RB(153)*XO +RF(168)*XO2 +RB(19)*XOH +RB(68)*XH2 +RB(104) & *XH2O +RB(56) IF( LITER .OR. .NOT. LBOUND(8) ) THEN IF(DEN(8) .LT. 1.0) DEN(8) = MAX( ADJ*ABV(8), DEN(8), SMALL ) VOLD = XCH3O XCH3O = ABV(8)/DEN(8) DIFF = ABS( (XCH3O-VOLD)/MAX(XCH3O,VOLD,SMALL) ) CONMX = MAX( CONMX, DIFF ) ENDIF C STEADY-STATE EXPRESSION FOR HCO ABV(9) = RF(120)*XHO2*XCH2O +RB(53)*XCH2O +RF(169)*XO2*XC2H + & RF(171)*XO2*XC2H3 +RF(32)*XO2*XCH2O +RF(159)*XCH3*XCH2O + & RB(158)*XCH4*XCO +RF(25)*XO*XC2H4 +RF(9)*XO*XCH2S + & RF(131)*XCH*XCO2 +RF(90)*XOH*XCH +RB(13)*XOH*XCO +RB(14) & *XH*XCO2 +RB(166)*XHO2*XCO +RF(134)*XO2*XCH2 +RF(7)*XO & *XCH2 +RF(124)*XO2*XCH +RB(99)*XH2O*XCO +RF(15)*XO*XCH2O & +RB(54)*XH2*XCO +RB(165)*XH*XCO*XM(165) +RF(57)*XH*XCH2O & +RF(100)*XOH*XCH2O +RB(164)*XH*XCO*XH2O DEN(9) = RB(120)*XH2O2 +RF(53)*XH +RB(169)*XCO +RB(171)*XCH2O + & RB(32)*XHO2 +RB(159)*XCH4 +RF(158)*XCH3 +RB(25)*XCH3 + & RB(9)*XH +RB(131)*XCO +RB(90)*XH +RF(13)*XO +RF(14)*XO + & RF(166)*XO2 +RB(134)*XOH +RB(7)*XH +RB(124)*XO +RF(99) & *XOH +RB(15)*XOH +RF(54)*XH +RF(165)*XM(165) +RB(57)*XH2 & +RB(100)*XH2O +RF(164)*XH2O IF( LITER .OR. .NOT. LBOUND(9) ) THEN IF(DEN(9) .LT. 1.0) DEN(9) = MAX( ADJ*ABV(9), DEN(9), SMALL ) VOLD = XHCO XHCO = ABV(9)/DEN(9) DIFF = ABS( (XHCO-VOLD)/MAX(XHCO,VOLD,SMALL) ) CONMX = MAX( CONMX, DIFF ) ENDIF C STEADY-STATE EXPRESSION FOR CH2 ABV(10) = RB(140)*XCO*XC2H3 +RB(122)*XH*XC2H +RB(127)*XH*XC2H2 + & RB(116)*XOH*XCH2O +RF(30)*XO*XCH2CO +RB(138)*XCH3*XCH3 + & RB(136)*XH2*XC2H2 +RB(139)*XCH2CO +RB(137)*XH*XC2H4 + & RF(23)*XO*XC2H2 +RB(49)*XCH3 +RF(149)*XCH2S*XCO +RB(135) & *XH*XCH3 +RF(150)*XCH2S*XCO2 +RB(92)*XH2O*XCH +RB(91)*XH & *XCH2O +RB(134)*XOH*XHCO +RF(95)*XOH*XCH3 +RB(7)*XH*XHCO & +RF(146)*XCH2S*XH2O +RF(141)*XCH2S*XN2 +RF(125)*XH2*XCH DEN(10) = RF(140)*XHCCO +RF(122)*XC +RF(127)*XCH +RF(116)*XHO2 + & RB(30)*XCO2 +RF(138)*XCH4 +RF(136)*XCH2 +RF(139)*XCO + & RF(137)*XCH3 +RB(23)*XCO +RF(49)*XH +RB(149)*XCO + & RF(135)*XH2 +RB(150)*XCO2 +RF(92)*XOH +RF(91)*XOH + & RF(134)*XO2 +RB(95)*XH2O +RF(7)*XO +RB(146)*XH2O + & RB(141)*XN2 +RB(125)*XH IF( LITER .OR. .NOT. LBOUND(10) ) THEN IF(DEN(10) .LT. 1.0) DEN(10) = MAX( ADJ*ABV(10), DEN(10), SMALL ) VOLD = XCH2 ACH2 = RF(136) BCH2 = DEN(10) - ACH2*XCH2 ZCH2 = ABV(10)/BCH2 IF( ACH2 .GT. SMALL ) THEN B2P4AC = BCH2*BCH2 + 4.0*ACH2*ABV(10) XCH2 = 0.5*( SQRT(B2P4AC) - BCH2 )/ACH2 IF( XCH2 .LT. SMALL ) XCH2 = ZCH2 ELSE XCH2 = ZCH2 ENDIF DIFF = ABS( (XCH2-VOLD)/MAX(XCH2,VOLD,SMALL) ) CONMX = MAX( CONMX, DIFF ) ENDIF C STEADY-STATE EXPRESSION FOR CH ABV(11) = RB(133)*XCO*XC2H2 +RF(20)*XO*XC2H +RB(127)*XH*XC2H2 + & RB(130)*XHCCO +RB(128)*XH*XC2H3 +RB(132)*XH*XCH2CO + & RB(129)*XH*XC2H4 +RF(50)*XH*XCH2S +RB(6)*XH*XCO +RB(131) & *XCO*XHCO +RB(90)*XH*XHCO +RF(92)*XOH*XCH2 +RB(124)*XO & *XHCO +RB(48)*XH2*XC +RB(126)*XH*XCH2O +RB(125)*XH*XCH2 DEN(11) = RF(133)*XHCCO +RB(20)*XCO +RF(127)*XCH2 +RF(130)*XCO + & RF(128)*XCH3 +RF(132)*XCH2O +RF(129)*XCH4 +RB(50)*XH2 + & RF(6)*XO +RF(131)*XCO2 +RF(90)*XOH +RB(92)*XH2O +RF(124) & *XO2 +RF(48)*XH +RF(126)*XH2O +RF(125)*XH2 IF( LITER .OR. .NOT. LBOUND(11) ) THEN IF(DEN(11) .LT. 1.0) DEN(11) = MAX( ADJ*ABV(11), DEN(11), SMALL ) VOLD = XCH XCH = ABV(11)/DEN(11) DIFF = ABS( (XCH-VOLD)/MAX(XCH,VOLD,SMALL) ) CONMX = MAX( CONMX, DIFF ) ENDIF C STEADY-STATE EXPRESSION FOR CH2CO ABV(12) = RF(81)*XHCCOH*XH +RF(106)*XOH*XC2H2 +RB(29)*XOH*XHCCO + & RB(30)*XCH2*XCO2 +RF(24)*XO*XC2H3 +RF(139)*XCH2*XCO + & RB(113)*XH2O*XHCCO +RB(80)*XCH3*XCO +RB(79)*XH2*XHCCO + & RF(132)*XCH*XCH2O DEN(12) = RB(81)*XH +RB(106)*XH +RF(29)*XO +RF(30)*XO +RB(24)*XH + & RB(139) +RF(113)*XOH +RF(80)*XH +RF(79)*XH +RB(132)*XH IF( LITER .OR. .NOT. LBOUND(12) ) THEN IF(DEN(12) .LT. 1.0) DEN(12) = MAX( ADJ*ABV(12), DEN(12), SMALL ) VOLD = XCH2CO XCH2CO = ABV(12)/DEN(12) DIFF = ABS( (XCH2CO-VOLD)/MAX(XCH2CO,VOLD,SMALL) ) CONMX = MAX( CONMX, DIFF ) ENDIF C STEADY-STATE EXPRESSION FOR C ABV(13) = RB(122)*XH*XC2H +RB(123)*XH*XC2H2 +RB(89)*XH*XCO + & RB(121)*XO*XCO +RF(48)*XH*XCH DEN(13) = RF(122)*XCH2 +RF(123)*XCH3 +RF(89)*XOH +RF(121)*XO2 + & RB(48)*XH2 IF( LITER .OR. .NOT. LBOUND(13) ) THEN IF(DEN(13) .LT. 1.0) DEN(13) = MAX( ADJ*ABV(13), DEN(13), SMALL ) VOLD = XC XC = ABV(13)/DEN(13) DIFF = ABS( (XC-VOLD)/MAX(XC,VOLD,SMALL) ) CONMX = MAX( CONMX, DIFF ) ENDIF C STEADY-STATE EXPRESSION FOR O ABV(14) = RB(22)*XOH*XC2H +RB(20)*XCH*XCO +RB(29)*XOH*XHCCO + & RB(30)*XCH2*XCO2 +RB(27)*XOH*XC2H5 +RB(26)*XCH3*XCH2O + & RB(17)*XOH*XCH2O +RF(31)*XO2*XCO +RB(24)*XH*XCH2CO + & RF(153)*XO2*XCH3 +RB(28)*XH*XCO*XCO +RB(19)*XOH*XCH3O + & RB(25)*XCH3*XHCO +RB(16)*XOH*XCH2O +RB(21)*XH*XHCCO + & RB(23)*XCH2*XCO +RB(1)*XO2*XM(1) +RB(9)*XH*XHCO +RB(8) & *XH2*XCO +RB(5)*XOH*XHO2 +RB(18)*XOH*XCH2OH +RB(2)*XOH & *XM(2) +RF(43)*XH*XHO2 +RB(12)*XCO2*XM(12) +RB(4)*XO2 & *XOH +RB(6)*XH*XCO +RB(13)*XOH*XCO +RB(14)*XH*XCO2 + & RF(121)*XO2*XC +RB(7)*XH*XHCO +RF(124)*XO2*XCH +RB(15) & *XOH*XHCO +RB(11)*XOH*XCH3 +RF(85)*XOH*XOH +RB(3)*XH*XOH & +RB(10)*XH*XCH2O +RF(37)*XH*XO2 DEN(14) = RF(22)*XC2H2 +RF(20)*XC2H +RF(29)*XCH2CO +RF(30)*XCH2CO & +RF(27)*XC2H6 +RF(26)*XC2H5 +RF(17)*XCH3O +RB(31)*XCO2 + & RF(24)*XC2H3 +RB(153)*XCH3O +RF(28)*XHCCO +RF(19)*XCH3OH & +RF(25)*XC2H4 +RF(16)*XCH2OH +RF(21)*XC2H2 +RF(23)*XC2H2 & +RF(1)*XO*XM(1) +RF(9)*XCH2S +RF(8)*XCH2S +RF(5)*XH2O2 + & RF(18)*XCH3OH +RF(2)*XH*XM(2) +RB(43)*XH2O +RF(12)*XCO & *XM(12) +RF(4)*XHO2 +RF(6)*XCH +RF(13)*XHCO +RF(14)*XHCO & +RB(121)*XCO +RF(7)*XCH2 +RB(124)*XHCO +RF(15)*XCH2O + & RF(11)*XCH4 +RB(85)*XH2O +RF(3)*XH2 +RF(10)*XCH3 +RB(37) & *XOH IF( LITER .OR. .NOT. LBOUND(14) ) THEN IF(DEN(14) .LT. 1.0) DEN(14) = MAX( ADJ*ABV(14), DEN(14), SMALL ) VOLD = XO AO = RF(1)*XM(1) BO = DEN(14) - AO*XO ZO = ABV(14)/BO IF( AO .GT. SMALL ) THEN B2P4AC = BO*BO + 4.0*AO*ABV(14) XO = 0.5*( SQRT(B2P4AC) - BO )/AO IF( XO .LT. SMALL ) XO = ZO ELSE XO = ZO ENDIF DIFF = ABS( (XO-VOLD)/MAX(XO,VOLD,SMALL) ) CONMX = MAX( CONMX, DIFF ) ENDIF C STEADY-STATE EXPRESSION FOR HCCO ABV(15) = RB(175)*XCO*XCO*XC2H2 +RB(133)*XCO*XC2H2 +RB(140)*XCO & *XC2H3 +RF(105)*XOH*XC2H +RF(29)*XO*XCH2CO +RB(174)*XOH & *XCO*XCO +RF(130)*XCH*XCO +RF(113)*XOH*XCH2CO +RF(79)*XH & *XCH2CO +RB(28)*XH*XCO*XCO +RF(21)*XO*XC2H2 +RB(78) & *XCH2S*XCO DEN(15) = RF(175)*XHCCO +RF(133)*XCH +RF(140)*XCH2 +RB(105)*XH + & RB(29)*XOH +RF(174)*XO2 +RB(130) +RB(113)*XH2O +RB(79) & *XH2 +RF(28)*XO +RB(21)*XH +RF(78)*XH IF( LITER .OR. .NOT. LBOUND(15) ) THEN IF(DEN(15) .LT. 1.0) DEN(15) = MAX( ADJ*ABV(15), DEN(15), SMALL ) VOLD = XHCCO AHCCO = RF(175) BHCCO = DEN(15) - AHCCO*XHCCO ZHCCO = ABV(15)/BHCCO IF( AHCCO .GT. SMALL ) THEN B2P4AC = BHCCO*BHCCO + 4.0*AHCCO*ABV(15) XHCCO = 0.5*( SQRT(B2P4AC) - BHCCO )/AHCCO IF( XHCCO .LT. SMALL ) XHCCO = ZHCCO ELSE XHCCO = ZHCCO ENDIF DIFF = ABS( (XHCCO-VOLD)/MAX(XHCCO,VOLD,SMALL) ) CONMX = MAX( CONMX, DIFF ) ENDIF RETURN END SUBROUTINE UPVALUE( IOPT, VNEW, BIG, LBOUND, CONMAX, XC2H, XC2H5, & XHCCOH, XC2H3, XCH2S, XCH2OH, XCH3OH, XCH3O, & XHCO, XCH2, XCH, XCH2CO, XC, XO, XHCCO ) C*****double precision IMPLICIT DOUBLE PRECISION (A-H, O-Z), INTEGER (I-N) C*****END double precision C*****single precision C IMPLICIT REAL (A-H, O-Z), INTEGER (I-N) C*****END single precision DIMENSION VNEW(*) LOGICAL LBOUND(*) IF( IOPT .EQ. 1 ) THEN XC2H = VNEW(1) XC2H5 = VNEW(2) XHCCOH = VNEW(3) XC2H3 = VNEW(4) XCH2S = VNEW(5) XCH2OH = VNEW(6) XCH3OH = VNEW(7) XCH3O = VNEW(8) XHCO = VNEW(9) XCH2 = VNEW(10) XCH = VNEW(11) XCH2CO = VNEW(12) XC = VNEW(13) XO = VNEW(14) XHCCO = VNEW(15) ELSEIF( IOPT .EQ. 2 ) THEN CALL EXBOUND( XC2H, VNEW(1), LBOUND(1), CONMAX, BIG ) CALL EXBOUND( XC2H5, VNEW(2), LBOUND(2), CONMAX, BIG ) CALL EXBOUND( XHCCOH, VNEW(3), LBOUND(3), CONMAX, BIG ) CALL EXBOUND( XC2H3, VNEW(4), LBOUND(4), CONMAX, BIG ) CALL EXBOUND( XCH2S, VNEW(5), LBOUND(5), CONMAX, BIG ) CALL EXBOUND( XCH2OH, VNEW(6), LBOUND(6), CONMAX, BIG ) CALL EXBOUND( XCH3OH, VNEW(7), LBOUND(7), CONMAX, BIG ) CALL EXBOUND( XCH3O, VNEW(8), LBOUND(8), CONMAX, BIG ) CALL EXBOUND( XHCO, VNEW(9), LBOUND(9), CONMAX, BIG ) CALL EXBOUND( XCH2, VNEW(10), LBOUND(10), CONMAX, BIG ) CALL EXBOUND( XCH, VNEW(11), LBOUND(11), CONMAX, BIG ) CALL EXBOUND( XCH2CO, VNEW(12), LBOUND(12), CONMAX, BIG ) CALL EXBOUND( XC, VNEW(13), LBOUND(13), CONMAX, BIG ) CALL EXBOUND( XO, VNEW(14), LBOUND(14), CONMAX, BIG ) CALL EXBOUND( XHCCO, VNEW(15), LBOUND(15), CONMAX, BIG ) ELSEIF( IOPT .EQ. 3 ) THEN VNEW(1) = XC2H VNEW(2) = XC2H5 VNEW(3) = XHCCOH VNEW(4) = XC2H3 VNEW(5) = XCH2S VNEW(6) = XCH2OH VNEW(7) = XCH3OH VNEW(8) = XCH3O VNEW(9) = XHCO VNEW(10) = XCH2 VNEW(11) = XCH VNEW(12) = XCH2CO VNEW(13) = XC VNEW(14) = XO VNEW(15) = XHCCO ENDIF RETURN END SUBROUTINE EXBOUND( VOLD, B, LBOUND, CONMAX, BIG ) C*****double precision IMPLICIT DOUBLE PRECISION (A-H, O-Z), INTEGER (I-N) C*****END double precision C*****single precision C IMPLICIT REAL (A-H, O-Z), INTEGER (I-N) C*****END single precision PARAMETER (VBOUND = 100.0) LOGICAL LBOUND AA = LOG(VOLD) - B IF( ABS(AA) .LT. VBOUND ) THEN IF( EXP(AA) .LT. BIG ) THEN DIFF = ABS( (EXP(AA)-VOLD)/MAX(EXP(AA),VOLD) ) CONMAX = MAX( CONMAX, DIFF ) VOLD = EXP(AA) LBOUND = .TRUE. ELSE LBOUND = .FALSE. ENDIF ELSE LBOUND = .FALSE. ENDIF RETURN END SUBROUTINE NETRATE( W, RF, RB, XM, XH2, XH, XO, XO2, XOH, XH2O, & XHO2, XH2O2, XC, XCH, XCH2, XCH2S, XCH3, XCH4, & XCO, XCO2, XHCO, XCH2O, XCH2OH, XCH3O, XCH3OH, & XC2H, XC2H2, XC2H3, XC2H4, XC2H5, XC2H6, & XHCCO, XCH2CO, XHCCOH, XN2 ) C*****double precision IMPLICIT DOUBLE PRECISION (A-H, O-Z), INTEGER (I-N) C*****END double precision C*****single precision C IMPLICIT REAL (A-H, O-Z), INTEGER (I-N) C*****END single precision DIMENSION W(*), RF(*), RB(*), XM(*) C NET PRODUCTION RATES FOR SKELETAL MECHANSIM W(1) = RF(1)*XO*XO*XM(1) - RB(1)*XO2*XM(1) W(2) = RF(2)*XH*XO*XM(2) - RB(2)*XOH*XM(2) W(3) = RF(3)*XH2*XO - RB(3)*XH*XOH W(4) = RF(4)*XO*XHO2 - RB(4)*XO2*XOH W(5) = RF(5)*XO*XH2O2 - RB(5)*XOH*XHO2 W(6) = RF(6)*XO*XCH - RB(6)*XH*XCO W(7) = RF(7)*XO*XCH2 - RB(7)*XH*XHCO W(8) = RF(8)*XO*XCH2S - RB(8)*XH2*XCO W(9) = RF(9)*XO*XCH2S - RB(9)*XH*XHCO W(10) = RF(10)*XO*XCH3 - RB(10)*XH*XCH2O W(11) = RF(11)*XO*XCH4 - RB(11)*XOH*XCH3 W(12) = RF(12)*XO*XCO*XM(12) - RB(12)*XCO2*XM(12) W(13) = RF(13)*XO*XHCO - RB(13)*XOH*XCO W(14) = RF(14)*XO*XHCO - RB(14)*XH*XCO2 W(15) = RF(15)*XO*XCH2O - RB(15)*XOH*XHCO W(16) = RF(16)*XO*XCH2OH - RB(16)*XOH*XCH2O W(17) = RF(17)*XO*XCH3O - RB(17)*XOH*XCH2O W(18) = RF(18)*XO*XCH3OH - RB(18)*XOH*XCH2OH W(19) = RF(19)*XO*XCH3OH - RB(19)*XOH*XCH3O W(20) = RF(20)*XO*XC2H - RB(20)*XCH*XCO W(21) = RF(21)*XO*XC2H2 - RB(21)*XH*XHCCO W(22) = RF(22)*XO*XC2H2 - RB(22)*XOH*XC2H W(23) = RF(23)*XO*XC2H2 - RB(23)*XCH2*XCO W(24) = RF(24)*XO*XC2H3 - RB(24)*XH*XCH2CO W(25) = RF(25)*XO*XC2H4 - RB(25)*XCH3*XHCO W(26) = RF(26)*XO*XC2H5 - RB(26)*XCH3*XCH2O W(27) = RF(27)*XO*XC2H6 - RB(27)*XOH*XC2H5 W(28) = RF(28)*XO*XHCCO - RB(28)*XH*XCO*XCO W(29) = RF(29)*XO*XCH2CO - RB(29)*XOH*XHCCO W(30) = RF(30)*XO*XCH2CO - RB(30)*XCH2*XCO2 W(31) = RF(31)*XO2*XCO - RB(31)*XO*XCO2 W(32) = RF(32)*XO2*XCH2O - RB(32)*XHO2*XHCO W(33) = RF(33)*XH*XO2*XM(33) - RB(33)*XHO2*XM(33) W(34) = RF(34)*XH*XO2*XO2 - RB(34)*XHO2*XO2 W(35) = RF(35)*XH*XO2*XH2O - RB(35)*XHO2*XH2O W(36) = RF(36)*XH*XO2*XN2 - RB(36)*XHO2*XN2 W(37) = RF(37)*XH*XO2 - RB(37)*XO*XOH W(38) = RF(38)*XH*XH*XM(38) - RB(38)*XH2*XM(38) W(39) = RF(39)*XH*XH*XH2 - RB(39)*XH2*XH2 W(40) = RF(40)*XH*XH*XH2O - RB(40)*XH2*XH2O W(41) = RF(41)*XH*XH*XCO2 - RB(41)*XH2*XCO2 W(42) = RF(42)*XH*XOH*XM(42) - RB(42)*XH2O*XM(42) W(43) = RF(43)*XH*XHO2 - RB(43)*XO*XH2O W(44) = RF(44)*XH*XHO2 - RB(44)*XH2*XO2 W(45) = RF(45)*XH*XHO2 - RB(45)*XOH*XOH W(46) = RF(46)*XH*XH2O2 - RB(46)*XH2*XHO2 W(47) = RF(47)*XH*XH2O2 - RB(47)*XOH*XH2O W(48) = RF(48)*XH*XCH - RB(48)*XH2*XC W(49) = RF(49)*XH*XCH2 - RB(49)*XCH3 W(50) = RF(50)*XH*XCH2S - RB(50)*XH2*XCH W(51) = RF(51)*XH*XCH3 - RB(51)*XCH4 W(52) = RF(52)*XH*XCH4 - RB(52)*XH2*XCH3 W(53) = RF(53)*XH*XHCO - RB(53)*XCH2O W(54) = RF(54)*XH*XHCO - RB(54)*XH2*XCO W(55) = RF(55)*XH*XCH2O - RB(55)*XCH2OH W(56) = RF(56)*XH*XCH2O - RB(56)*XCH3O W(57) = RF(57)*XH*XCH2O - RB(57)*XH2*XHCO W(58) = RF(58)*XH*XCH2OH - RB(58)*XCH3OH W(59) = RF(59)*XH*XCH2OH - RB(59)*XH2*XCH2O W(60) = RF(60)*XH*XCH2OH - RB(60)*XOH*XCH3 W(61) = RF(61)*XH*XCH2OH - RB(61)*XH2O*XCH2S W(62) = RF(62)*XH*XCH3O - RB(62)*XCH3OH W(63) = RF(63)*XCH3O*XH - RB(63)*XCH2OH*XH W(64) = RF(64)*XH*XCH3O - RB(64)*XH2*XCH2O W(65) = RF(65)*XH*XCH3O - RB(65)*XOH*XCH3 W(66) = RF(66)*XH*XCH3O - RB(66)*XH2O*XCH2S W(67) = RF(67)*XH*XCH3OH - RB(67)*XH2*XCH2OH W(68) = RF(68)*XH*XCH3OH - RB(68)*XH2*XCH3O W(69) = RF(69)*XH*XC2H - RB(69)*XC2H2 W(70) = RF(70)*XH*XC2H2 - RB(70)*XC2H3 W(71) = RF(71)*XH*XC2H3 - RB(71)*XC2H4 W(72) = RF(72)*XH*XC2H3 - RB(72)*XH2*XC2H2 W(73) = RF(73)*XH*XC2H4 - RB(73)*XC2H5 W(74) = RF(74)*XH*XC2H4 - RB(74)*XH2*XC2H3 W(75) = RF(75)*XH*XC2H5 - RB(75)*XC2H6 W(76) = RF(76)*XH*XC2H5 - RB(76)*XH2*XC2H4 W(77) = RF(77)*XH*XC2H6 - RB(77)*XH2*XC2H5 W(78) = RF(78)*XH*XHCCO - RB(78)*XCH2S*XCO W(79) = RF(79)*XH*XCH2CO - RB(79)*XH2*XHCCO W(80) = RF(80)*XH*XCH2CO - RB(80)*XCH3*XCO W(81) = RF(81)*XHCCOH*XH - RB(81)*XCH2CO*XH W(82) = RF(82)*XH2*XCO - RB(82)*XCH2O W(83) = RF(83)*XH2*XOH - RB(83)*XH*XH2O W(84) = RF(84)*XOH*XOH - RB(84)*XH2O2 W(85) = RF(85)*XOH*XOH - RB(85)*XO*XH2O W(86) = RF(86)*XOH*XHO2 - RB(86)*XO2*XH2O W(87) = RF(87)*XOH*XH2O2 - RB(87)*XH2O*XHO2 W(88) = RF(88)*XOH*XH2O2 - RB(88)*XH2O*XHO2 W(89) = RF(89)*XOH*XC - RB(89)*XH*XCO W(90) = RF(90)*XOH*XCH - RB(90)*XH*XHCO W(91) = RF(91)*XOH*XCH2 - RB(91)*XH*XCH2O W(92) = RF(92)*XOH*XCH2 - RB(92)*XH2O*XCH W(93) = RF(93)*XOH*XCH2S - RB(93)*XH*XCH2O W(94) = RF(94)*XOH*XCH3 - RB(94)*XCH3OH W(95) = RF(95)*XOH*XCH3 - RB(95)*XH2O*XCH2 W(96) = RF(96)*XOH*XCH3 - RB(96)*XH2O*XCH2S W(97) = RF(97)*XOH*XCH4 - RB(97)*XH2O*XCH3 W(98) = RF(98)*XOH*XCO - RB(98)*XH*XCO2 W(99) = RF(99)*XOH*XHCO - RB(99)*XH2O*XCO W(100) = RF(100)*XOH*XCH2O - RB(100)*XH2O*XHCO W(101) = RF(101)*XOH*XCH2OH - RB(101)*XH2O*XCH2O W(102) = RF(102)*XOH*XCH3O - RB(102)*XH2O*XCH2O W(103) = RF(103)*XOH*XCH3OH - RB(103)*XH2O*XCH2OH W(104) = RF(104)*XOH*XCH3OH - RB(104)*XH2O*XCH3O W(105) = RF(105)*XOH*XC2H - RB(105)*XH*XHCCO W(106) = RF(106)*XOH*XC2H2 - RB(106)*XH*XCH2CO W(107) = RF(107)*XOH*XC2H2 - RB(107)*XH*XHCCOH W(108) = RF(108)*XOH*XC2H2 - RB(108)*XH2O*XC2H W(109) = RF(109)*XOH*XC2H2 - RB(109)*XCH3*XCO W(110) = RF(110)*XOH*XC2H3 - RB(110)*XH2O*XC2H2 W(111) = RF(111)*XOH*XC2H4 - RB(111)*XH2O*XC2H3 W(112) = RF(112)*XOH*XC2H6 - RB(112)*XH2O*XC2H5 W(113) = RF(113)*XOH*XCH2CO - RB(113)*XH2O*XHCCO W(114) = RF(114)*XHO2*XHO2 - RB(114)*XO2*XH2O2 W(115) = RF(115)*XHO2*XHO2 - RB(115)*XO2*XH2O2 W(116) = RF(116)*XHO2*XCH2 - RB(116)*XOH*XCH2O W(117) = RF(117)*XHO2*XCH3 - RB(117)*XO2*XCH4 W(118) = RF(118)*XHO2*XCH3 - RB(118)*XOH*XCH3O W(119) = RF(119)*XHO2*XCO - RB(119)*XOH*XCO2 W(120) = RF(120)*XHO2*XCH2O - RB(120)*XH2O2*XHCO W(121) = RF(121)*XO2*XC - RB(121)*XO*XCO W(122) = RF(122)*XC*XCH2 - RB(122)*XH*XC2H W(123) = RF(123)*XC*XCH3 - RB(123)*XH*XC2H2 W(124) = RF(124)*XO2*XCH - RB(124)*XO*XHCO W(125) = RF(125)*XH2*XCH - RB(125)*XH*XCH2 W(126) = RF(126)*XH2O*XCH - RB(126)*XH*XCH2O W(127) = RF(127)*XCH*XCH2 - RB(127)*XH*XC2H2 W(128) = RF(128)*XCH*XCH3 - RB(128)*XH*XC2H3 W(129) = RF(129)*XCH*XCH4 - RB(129)*XH*XC2H4 W(130) = RF(130)*XCH*XCO - RB(130)*XHCCO W(131) = RF(131)*XCH*XCO2 - RB(131)*XCO*XHCO W(132) = RF(132)*XCH*XCH2O - RB(132)*XH*XCH2CO W(133) = RF(133)*XCH*XHCCO - RB(133)*XCO*XC2H2 W(134) = RF(134)*XO2*XCH2 - RB(134)*XOH*XHCO W(135) = RF(135)*XH2*XCH2 - RB(135)*XH*XCH3 W(136) = RF(136)*XCH2*XCH2 - RB(136)*XH2*XC2H2 W(137) = RF(137)*XCH2*XCH3 - RB(137)*XH*XC2H4 W(138) = RF(138)*XCH2*XCH4 - RB(138)*XCH3*XCH3 W(139) = RF(139)*XCH2*XCO - RB(139)*XCH2CO W(140) = RF(140)*XCH2*XHCCO - RB(140)*XCO*XC2H3 W(141) = RF(141)*XCH2S*XN2 - RB(141)*XCH2*XN2 W(142) = RF(142)*XO2*XCH2S - RB(142)*XH*XOH*XCO W(143) = RF(143)*XO2*XCH2S - RB(143)*XH2O*XCO W(144) = RF(144)*XH2*XCH2S - RB(144)*XH*XCH3 W(145) = RF(145)*XH2O*XCH2S - RB(145)*XCH3OH W(146) = RF(146)*XCH2S*XH2O - RB(146)*XCH2*XH2O W(147) = RF(147)*XCH2S*XCH3 - RB(147)*XH*XC2H4 W(148) = RF(148)*XCH2S*XCH4 - RB(148)*XCH3*XCH3 W(149) = RF(149)*XCH2S*XCO - RB(149)*XCH2*XCO W(150) = RF(150)*XCH2S*XCO2 - RB(150)*XCH2*XCO2 W(151) = RF(151)*XCH2S*XCO2 - RB(151)*XCO*XCH2O W(152) = RF(152)*XCH2S*XC2H6 - RB(152)*XCH3*XC2H5 W(153) = RF(153)*XO2*XCH3 - RB(153)*XO*XCH3O W(154) = RF(154)*XO2*XCH3 - RB(154)*XOH*XCH2O W(155) = RF(155)*XH2O2*XCH3 - RB(155)*XHO2*XCH4 W(156) = RF(156)*XCH3*XCH3 - RB(156)*XC2H6 W(157) = RF(157)*XCH3*XCH3 - RB(157)*XH*XC2H5 W(158) = RF(158)*XCH3*XHCO - RB(158)*XCH4*XCO W(159) = RF(159)*XCH3*XCH2O - RB(159)*XCH4*XHCO W(160) = RF(160)*XCH3*XCH3OH - RB(160)*XCH4*XCH2OH W(161) = RF(161)*XCH3*XCH3OH - RB(161)*XCH4*XCH3O W(162) = RF(162)*XCH3*XC2H4 - RB(162)*XCH4*XC2H3 W(163) = RF(163)*XCH3*XC2H6 - RB(163)*XCH4*XC2H5 W(164) = RF(164)*XHCO*XH2O - RB(164)*XH*XCO*XH2O W(165) = RF(165)*XHCO*XM(165) - RB(165)*XH*XCO*XM(165) W(166) = RF(166)*XO2*XHCO - RB(166)*XHO2*XCO W(167) = RF(167)*XO2*XCH2OH - RB(167)*XHO2*XCH2O W(168) = RF(168)*XO2*XCH3O - RB(168)*XHO2*XCH2O W(169) = RF(169)*XO2*XC2H - RB(169)*XCO*XHCO W(170) = RF(170)*XH2*XC2H - RB(170)*XH*XC2H2 W(171) = RF(171)*XO2*XC2H3 - RB(171)*XHCO*XCH2O W(172) = RF(172)*XC2H4 - RB(172)*XH2*XC2H2 W(173) = RF(173)*XO2*XC2H5 - RB(173)*XHO2*XC2H4 W(174) = RF(174)*XO2*XHCCO - RB(174)*XOH*XCO*XCO W(175) = RF(175)*XHCCO*XHCCO - RB(175)*XCO*XCO*XC2H2 RETURN END