C                                                                      C
C----------------------------------------------------------------------C
C                                                                      C
C     Automatically generated code
C
C     Tianfeng Lu
C     University of Connecticut
C     191 Auditorium Road U-3139
C     Storrs, CT 06269, USA
C     Email: tlu@engr.uconn.edu
C
C     September 22, 2016
C
C     Copyright © 2018 HyChem   
C
C                                                                      C
C----------------------------------------------------------------------C
C                                                                      C
      SUBROUTINE CKWYP  (P, T, Y, ICKWRK, RCKWRK, WDOT)
      IMPLICIT DOUBLE PRECISION (A-H, O-Z), INTEGER (I-N)
      DIMENSION Y(*),WDOT(*),ICKWRK(*),RCKWRK(*)
      DIMENSION RF(192),RB(192),RKLOW(20),C(29)
      DIMENSION XQ(9)
C
      CALL YTCP(P, T, Y, C)
      CALL RATT(T, RF, RB, RKLOW)
      CALL RATX(T, C, RF, RB, RKLOW)
      CALL QSSA(RF, RB, XQ)
      CALL RDOT(RF, RB, WDOT)
C
      END
C                                                                      C
C----------------------------------------------------------------------C
C                                                                      C
      SUBROUTINE CKWYR  (RHO, T, Y, ICKWRK, RCKWRK, WDOT)
      IMPLICIT DOUBLE PRECISION (A-H, O-Z), INTEGER (I-N)
      DIMENSION Y(*),WDOT(*),ICKWRK(*),RCKWRK(*)
      DIMENSION RF(192),RB(192),RKLOW(20),C(29)
      DIMENSION XQ(9)
C
      CALL YTCR(RHO, T, Y, C)
      CALL RATT(T, RF, RB, RKLOW)
      CALL RATX(T, C, RF, RB, RKLOW)
      CALL QSSA(RF, RB, XQ)
      CALL RDOT(RF, RB, WDOT)
C
      END
C                                                                      C
C----------------------------------------------------------------------C
C                                                                      C
      SUBROUTINE YTCP (P, T, Y, C)
      IMPLICIT DOUBLE PRECISION (A-H, O-Z), INTEGER (I-N)
      DIMENSION Y(*), C(*)
C
      C(1) = Y(1)/1.833601440191269D2
      C(2) = Y(2)/2.805418062210083D1
      C(3) = Y(3)/1.604303026199341D1
      C(4) = Y(4)/4.208127093315125D1
      C(5) = Y(5)/5.610836124420166D1
      C(6) = Y(6)/5.610836124420166D1
      C(7) = Y(7)/2.015939950942993D0
      C(8) = Y(8)/3.007012057304382D1
      C(9) = Y(9)/2.80105504989624D1
      C(10) = Y(10)/7.811472201347351D1
      C(11) = Y(11)/2.603824067115784D1
      C(12) = Y(12)/9.214181232452393D1
      C(13) = Y(13)/1.503506028652191D1
      C(14) = Y(14)/3.199880027770996D1
      C(15) = Y(15)/1.599940013885498D1
      C(16) = Y(16)/1.700737011432648D1
      C(17) = Y(17)/3.300677025318146D1
      C(18) = Y(18)/1.801534008979797D1
      C(19) = Y(19)/3.401474022865295D1
      C(20) = Y(20)/1.007969975471497D0
      C(21) = Y(21)/3.00264904499054D1
      C(22) = Y(22)/4.400995063781738D1
      C(23) = Y(23)/4.203764081001282D1
      C(24) = Y(24)/3.905736100673676D1
      C(25) = Y(25)/4.107330095767975D1
      C(26) = Y(26)/8.008703184127808D1
      C(27) = Y(27)/6.509560167789459D1
      C(28) = Y(28)/9.310615217685699D1
      C(29) = Y(29)/2.801339912414551D1
C
      SUM = 0D0
      DO K = 1, 29
         SUM = SUM + C(K)
      ENDDO
      SUM = P/(SUM*T*8.314510D7)
C
      DO K = 1, 29
         C(K) = C(K)*SUM
      ENDDO
C
      END
C                                                                      C
C----------------------------------------------------------------------C
C                                                                      C
      SUBROUTINE YTCR (RHO, T, Y, C)
      IMPLICIT DOUBLE PRECISION (A-H, O-Z), INTEGER (I-N)
      DIMENSION Y(*), C(*)
      DATA SMALL/1D-50/
C
C     POSF5433
      C(1) = Y(1)/1.833601440191269D2
C     C2H4
      C(2) = Y(2)/2.805418062210083D1
C     CH4
      C(3) = Y(3)/1.604303026199341D1
C     C3H6
      C(4) = Y(4)/4.208127093315125D1
C     iC4H8
      C(5) = Y(5)/5.610836124420166D1
C     C4H81
      C(6) = Y(6)/5.610836124420166D1
C     H2
      C(7) = Y(7)/2.015939950942993D0
C     C2H6
      C(8) = Y(8)/3.007012057304382D1
C     CO
      C(9) = Y(9)/2.80105504989624D1
C     C6H6
      C(10) = Y(10)/7.811472201347351D1
C     C2H2
      C(11) = Y(11)/2.603824067115784D1
C     C6H5CH3
      C(12) = Y(12)/9.214181232452393D1
C     CH3
      C(13) = Y(13)/1.503506028652191D1
C     O2
      C(14) = Y(14)/3.199880027770996D1
C     O
      C(15) = Y(15)/1.599940013885498D1
C     OH
      C(16) = Y(16)/1.700737011432648D1
C     HO2
      C(17) = Y(17)/3.300677025318146D1
C     H2O
      C(18) = Y(18)/1.801534008979797D1
C     H2O2
      C(19) = Y(19)/3.401474022865295D1
C     H
      C(20) = Y(20)/1.007969975471497D0
C     CH2O
      C(21) = Y(21)/3.00264904499054D1
C     CO2
      C(22) = Y(22)/4.400995063781738D1
C     CH2CO
      C(23) = Y(23)/4.203764081001282D1
C     C3H3
      C(24) = Y(24)/3.905736100673676D1
C     aC3H5
      C(25) = Y(25)/4.107330095767975D1
C     C5H4O
      C(26) = Y(26)/8.008703184127808D1
C     C5H5
      C(27) = Y(27)/6.509560167789459D1
C     C6H5O
      C(28) = Y(28)/9.310615217685699D1
C     N2
      C(29) = Y(29)/2.801339912414551D1
C
      DO K = 1, 29
         C(K) = RHO * C(K)
      ENDDO
C
      END
C                                                                      C
C----------------------------------------------------------------------C
C                                                                      C
      SUBROUTINE RATT (T, RF, RB, RKLOW)
      IMPLICIT DOUBLE PRECISION (A-H, O-Z), INTEGER (I-N)
      PARAMETER (RU=8.31451D7, SMALL=1.D-200, PATM=1.01325D6)
      DIMENSION RF(*), RB(*), RKLOW(*)
      DIMENSION SMH(37), EG(37)
C
      ALOGT = LOG(T)
      TI = 1D0/T
      TI2 = TI*TI
C
      CALL RDSMH (T, SMH)
      EG(2) = EXP(SMH(2))
      EG(3) = EXP(SMH(3))
      EG(4) = EXP(SMH(4))
      EG(5) = EXP(SMH(5))
      EG(6) = EXP(SMH(6))
      EG(7) = EXP(SMH(7))
      EG(8) = EXP(SMH(8))
      EG(9) = EXP(SMH(9))
      EG(10) = EXP(SMH(10))
      EG(11) = EXP(SMH(11))
      EG(12) = EXP(SMH(12))
      EG(13) = EXP(SMH(13))
      EG(14) = EXP(SMH(14))
      EG(15) = EXP(SMH(15))
      EG(16) = EXP(SMH(16))
      EG(17) = EXP(SMH(17))
      EG(18) = EXP(SMH(18))
      EG(19) = EXP(SMH(19))
      EG(20) = EXP(SMH(20))
      EG(21) = EXP(SMH(21))
      EG(22) = EXP(SMH(22))
      EG(23) = EXP(SMH(23))
      EG(24) = EXP(SMH(24))
      EG(25) = EXP(SMH(25))
      EG(26) = EXP(SMH(26))
      EG(27) = EXP(SMH(27))
      EG(28) = EXP(SMH(28))
      EG(29) = EXP(SMH(29))
      EG(30) = EXP(SMH(30))
      EG(31) = EXP(SMH(31))
      EG(32) = EXP(SMH(32))
      EG(33) = EXP(SMH(33))
      EG(34) = EXP(SMH(34))
      EG(35) = EXP(SMH(35))
      EG(36) = EXP(SMH(36))
      EG(37) = EXP(SMH(37))
      PFAC1 = PATM / (RU*T)
      PFAC2 = PFAC1*PFAC1
      PFAC3 = PFAC2*PFAC1
C
C
C     R1: POSF5433 = 3.09172C2H4 + 0.803846C3H6 + 0.154586iC4H8 + 0.309172C4H81 + 0.116667C6H6 + 0.05C6H5CH3 + 1.5CH3 + 0.5H
      RF(1) = EXP(5.933277692844006D1 -2.66D0*ALOGT 
     * -4.337576710055816D4*TI)
      RB(1) = 0D0
C     R2: POSF5433 + 0.7H = 3.09911C2H4 + 0.25CH4 + 0.805769C3H6 + 0.154956iC4H8 + 0.309911C4H81 + H2 + 0.175C6H6 + 0.075C6H5CH3 + 0.7CH3
      RF(2) = EXP(-2.569158202235591D0 +4.76D0*ALOGT 
     * -6.51615262926932D2*TI)
      RB(2) = 0D0
C     R3: POSF5433 + 0.3CH3 = 3.09911C2H4 + 1.25CH4 + 0.805769C3H6 + 0.154956iC4H8 + 0.309911C4H81 + 0.175C6H6 + 0.075C6H5CH3 + 0.3H
      RF(3) = EXP(-1.455784798271499D1 +5.95D0*ALOGT 
     * -2.389474024621395D3*TI)
      RB(3) = 0D0
C     R4: POSF5433 + OH = 3.09911C2H4 + 0.25CH4 + 0.805769C3H6 + 0.154956iC4H8 + 0.309911C4H81 + 0.175C6H6 + 0.075C6H5CH3 + 0.7CH3 + H2O + 0.3H
      RF(4) = EXP(2.429279853342778D1 +1.02D0*ALOGT 
     * -1.072857935408309D2*TI)
      RB(4) = 0D0
C     R5: POSF5433 + O2 = 3.09911C2H4 + 0.25CH4 + 0.805769C3H6 + 0.154956iC4H8 + 0.309911C4H81 + 0.175C6H6 + 0.075C6H5CH3 + 0.7CH3 + HO2 + 0.3H
      RF(5) = EXP(3.511538975921468D1 +6.D-2*ALOGT 
     * -2.190632990570157D4*TI)
      RB(5) = 0D0
C     R6: POSF5433 + HO2 = 3.09911C2H4 + 0.25CH4 + 0.805769C3H6 + 0.154956iC4H8 + 0.309911C4H81 + 0.175C6H6 + 0.075C6H5CH3 + 0.7CH3 + H2O2 + 0.3H
      RF(6) = EXP(1.097164063371953D1 +2.94D0*ALOGT 
     * -5.942988844827451D3*TI)
      RB(6) = 0D0
C     R7: POSF5433 + O = 3.09911C2H4 + 0.25CH4 + 0.805769C3H6 + 0.154956iC4H8 + 0.309911C4H81 + 0.175C6H6 + 0.075C6H5CH3 + 0.7CH3 + OH + 0.3H
      RF(7) = EXP(4.493120682179469D0 +3.86D0*ALOGT 
     * -3.850110724113026D2*TI)
      RB(7) = 0D0
C     R8: O2 + H = O + OH
      RF(8) = EXP(3.781365440989417D1 -6.707D-1*ALOGT 
     * -8.575315233251871D3*TI)
      EQK = EG(15)*EG(16)/EG(14)/EG(20)
      RB(8) = RF(8) / MAX(EQK, SMALL)
C     R9: H2 + O = OH + H
      RF(9) = EXP(1.073400250738888D1 +2.7D0*ALOGT 
     * -3.150136339425897D3*TI)
      EQK = EG(16)*EG(20)/EG(7)/EG(15)
      RB(9) = RF(9) / MAX(EQK, SMALL)
C     R10: H2 + OH = H2O + H
      RF(10) = EXP(1.897111162231072D1 +1.51D0*ALOGT 
     * -1.726033170004924D3*TI)
      EQK = EG(18)*EG(20)/EG(7)/EG(16)
      RB(10) = RF(10) / MAX(EQK, SMALL)
C     R11: 2OH = O + H2O
      RF(11) = EXP(1.058986184880864D1 +2.4D0*ALOGT 
     * +1.061787168720231D3*TI)
      EQK = EG(15)*EG(18)/EG(16)/EG(16)
      RB(11) = RF(11) / MAX(EQK, SMALL)
C     R12: 2H = H2
      RF(12) = EXP(4.202314503819682D1 -1.D0*ALOGT)
      EQK = EG(7)/EG(20)/EG(20)/PFAC1
      RB(12) = RF(12) / MAX(EQK, SMALL)
C     R13: 2H = H2
      RF(13) = EXP(3.903858606524095D1 -6.D-1*ALOGT)
      EQK = EG(7)/EG(20)/EG(20)/PFAC1
      RB(13) = RF(13) / MAX(EQK, SMALL)
C     R14: 2H = H2
      RF(14) = EXP(4.547615992139523D1 -1.25D0*ALOGT)
      EQK = EG(7)/EG(20)/EG(20)/PFAC1
      RB(14) = RF(14) / MAX(EQK, SMALL)
C     R15: OH + H = H2O
      RF(15) = EXP(5.213847658679322D1 -2.D0*ALOGT)
      EQK = EG(18)/EG(16)/EG(20)/PFAC1
      RB(15) = RF(15) / MAX(EQK, SMALL)
C     R16: O + H = OH
      RF(16) = EXP(4.369021565896671D1 -1.D0*ALOGT)
      EQK = EG(16)/EG(15)/EG(20)/PFAC1
      RB(16) = RF(16) / MAX(EQK, SMALL)
C     R17: 2O = O2
      RF(17) = EXP(3.932626813769273D1 -1.D0*ALOGT)
      EQK = EG(14)/EG(15)/EG(15)/PFAC1
      RB(17) = RF(17) / MAX(EQK, SMALL)
C     R18: O2 + H = HO2
      RF(18) = EXP(2.926339399964515D1 +4.4D-1*ALOGT)
      EQK = EG(17)/EG(14)/EG(20)/PFAC1
      RB(18) = RF(18) / MAX(EQK, SMALL)
C     R19: H2 + O2 = HO2 + H
      RF(19) = EXP(1.329058600981878D1 +2.433D0*ALOGT 
     * -2.692309815207098D4*TI)
      EQK = EG(17)*EG(20)/EG(7)/EG(14)
      RB(19) = RF(19) / MAX(EQK, SMALL)
C     R20: 2OH = H2O2
      RF(20) = EXP(3.234055131724088D1 -3.7D-1*ALOGT)
      EQK = EG(19)/EG(16)/EG(16)/PFAC1
      RB(20) = RF(20) / MAX(EQK, SMALL)
C     R21: HO2 + H = 2OH
      RF(21) = EXP(3.194650722679419D1 -1.484489169537763D2*TI)
      EQK = EG(16)*EG(16)/EG(17)/EG(20)
      RB(21) = RF(21) / MAX(EQK, SMALL)
C     R22: O + HO2 = O2 + OH
      RF(22) = 4.D13
      EQK = EG(14)*EG(16)/EG(15)/EG(17)
      RB(22) = RF(22) / MAX(EQK, SMALL)
C     R23: 2HO2 = O2 + H2O2
      RF(23) = EXP(2.559080028740199D1 +8.20243168253069D2*TI)
      EQK = EG(14)*EG(19)/EG(17)/EG(17)
      RB(23) = RF(23) / MAX(EQK, SMALL)
C     R24: 2HO2 = O2 + H2O2
      RF(24) = EXP(3.353310785188531D1 -6.038600011679036D3*TI)
      EQK = EG(14)*EG(19)/EG(17)/EG(17)
      RB(24) = RF(24) / MAX(EQK, SMALL)
C     R25: OH + HO2 = O2 + H2O
      RF(25) = EXP(4.17901213782829D1 -1.76D0*ALOGT 
     * -3.019300005839518D1*TI)
      EQK = EG(14)*EG(18)/EG(16)/EG(17)
      RB(25) = RF(25) / MAX(EQK, SMALL)
C     R26: OH + HO2 = O2 + H2O
      RF(26) = EXP(1.62861784418104D2 -1.672D1*ALOGT 
     * -1.655582836535336D4*TI)
      EQK = EG(14)*EG(18)/EG(16)/EG(17)
      RB(26) = RF(26) / MAX(EQK, SMALL)
C     R27: OH + HO2 = O2 + H2O
      RF(27) = EXP(2.855130386907224D1 +2.D0*ALOGT 
     * -2.012866670559679D4*TI)
      EQK = EG(14)*EG(18)/EG(16)/EG(17)
      RB(27) = RF(27) / MAX(EQK, SMALL)
C     R28: H2O2 + H = H2 + HO2
      RF(28) = EXP(1.561556883000702D1 +2.D0*ALOGT 
     * -2.616726671727582D3*TI)
      EQK = EG(7)*EG(17)/EG(19)/EG(20)
      RB(28) = RF(28) / MAX(EQK, SMALL)
C     R29: H2O2 + H = OH + H2O
      RF(29) = EXP(3.081323295642516D1 -1.997770170530481D3*TI)
      EQK = EG(16)*EG(18)/EG(19)/EG(20)
      RB(29) = RF(29) / MAX(EQK, SMALL)
C     R30: O + H2O2 = OH + HO2
      RF(30) = EXP(1.608039378377431D1 +2.D0*ALOGT 
     * -1.997770170530481D3*TI)
      EQK = EG(16)*EG(17)/EG(15)/EG(19)
      RB(30) = RF(30) / MAX(EQK, SMALL)
C     R31: OH + H2O2 = HO2 + H2O
      RF(31) = EXP(2.832416829648849D1 -2.148735170822457D2*TI)
      EQK = EG(17)*EG(18)/EG(16)/EG(19)
      RB(31) = RF(31) / MAX(EQK, SMALL)
C     R32: OH + H2O2 = HO2 + H2O
      RF(32) = EXP(9.538806728516803D1 -7.D0*ALOGT 
     * -1.892094670326098D4*TI)
      EQK = EG(17)*EG(18)/EG(16)/EG(19)
      RB(32) = RF(32) / MAX(EQK, SMALL)
C     R33: CO + O = CO2
      RF(33) = EXP(2.333480513766778D1 -1.199668535653569D3*TI)
      EQK = EG(26)/EG(9)/EG(15)/PFAC1
      RB(33) = RF(33) / MAX(EQK, SMALL)
C     R34: CO + OH = H + CO2
      RF(34) = EXP(1.116280045189523D1 +2.053D0*ALOGT 
     * +1.789790721794902D2*TI)
      EQK = EG(20)*EG(26)/EG(9)/EG(16)
      RB(34) = RF(34) / MAX(EQK, SMALL)
C     R35: CO + OH = H + CO2
      RF(35) = EXP(2.938143762162222D1 -6.64D-1*ALOGT 
     * -1.669823868229545D2*TI)
      EQK = EG(20)*EG(26)/EG(9)/EG(16)
      RB(35) = RF(35) / MAX(EQK, SMALL)
C     R36: CO + O2 = O + CO2
      RF(36) = EXP(2.774345654525834D1 -2.400343504642417D4*TI)
      EQK = EG(15)*EG(26)/EG(9)/EG(14)
      RB(36) = RF(36) / MAX(EQK, SMALL)
C     R37: CO + HO2 = OH + CO2
      RF(37) = EXP(1.196400108433044D1 +2.18D0*ALOGT 
     * -9.0290204129627D3*TI)
      EQK = EG(16)*EG(26)/EG(9)/EG(17)
      RB(37) = RF(37) / MAX(EQK, SMALL)
C     R38: H + HCO = H2 + CO
      RF(38) = 1.2D14
      EQK = EG(7)*EG(9)/EG(20)/EG(23)
      RB(38) = RF(38) / MAX(EQK, SMALL)
C     R39: O + HCO = CO + OH
      RF(39) = 3.D13
      EQK = EG(9)*EG(16)/EG(15)/EG(23)
      RB(39) = RF(39) / MAX(EQK, SMALL)
C     R40: O + HCO = H + CO2
      RF(40) = 3.D13
      EQK = EG(20)*EG(26)/EG(15)/EG(23)
      RB(40) = RF(40) / MAX(EQK, SMALL)
C     R41: OH + HCO = CO + H2O
      RF(41) = 3.02D13
      EQK = EG(9)*EG(18)/EG(16)/EG(23)
      RB(41) = RF(41) / MAX(EQK, SMALL)
C     R42: HCO = CO + H
      RF(42) = EXP(3.976988501176528D1 -1.D0*ALOGT 
     * -8.554683349878635D3*TI)
      EQK = EG(9)*EG(20)/EG(23)*PFAC1
      RB(42) = RF(42) / MAX(EQK, SMALL)
C     R43: HCO = CO + H
      RF(43) = EXP(4.225479166155327D1 -1.D0*ALOGT 
     * -8.554683349878635D3*TI)
      EQK = EG(9)*EG(20)/EG(23)*PFAC1
      RB(43) = RF(43) / MAX(EQK, SMALL)
C     R44: O2 + HCO = CO + HO2
      RF(44) = EXP(2.321150027682708D1 +8.070000000000001D-1*ALOGT 
     * +3.658385173742216D2*TI)
      EQK = EG(9)*EG(17)/EG(14)/EG(23)
      RB(44) = RF(44) / MAX(EQK, SMALL)
C     R45: H2 + CO = CH2O
      RF(45) = EXP(1.757671067365784D1 +1.5D0*ALOGT 
     * -4.005604674413761D4*TI)
      EQK = EG(24)/EG(7)/EG(9)/PFAC1
      RB(45) = RF(45) / MAX(EQK, SMALL)
C     R46: H + HCO = CH2O
      RF(46) = EXP(2.77171988121696D1 +4.8D-1*ALOGT 
     * +1.308363335863791D2*TI)
      EQK = EG(24)/EG(20)/EG(23)/PFAC1
      RB(46) = RF(46) / MAX(EQK, SMALL)
C     R47: H + CH2 = CH3
      RF(47) = EXP(3.775765221977888D1 -8.D-1*ALOGT)
      EQK = EG(13)/EG(20)/EG(21)/PFAC1
      RB(47) = RF(47) / MAX(EQK, SMALL)
C     R48: O + CH2 = H + HCO
      RF(48) = 8.D13
      EQK = EG(20)*EG(23)/EG(15)/EG(21)
      RB(48) = RF(48) / MAX(EQK, SMALL)
C     R49: OH + CH2 = H + CH2O
      RF(49) = 2.D13
      EQK = EG(20)*EG(24)/EG(16)/EG(21)
      RB(49) = RF(49) / MAX(EQK, SMALL)
C     R50: H2 + CH2 = CH3 + H
      RF(50) = EXP(1.312236337740433D1 +2.D0*ALOGT 
     * -3.638256507036619D3*TI)
      EQK = EG(13)*EG(20)/EG(7)/EG(21)
      RB(50) = RF(50) / MAX(EQK, SMALL)
C     R51: O2 + CH2 = OH + HCO
      RF(51) = EXP(2.999187511704657D1 -7.548250014598796D2*TI)
      EQK = EG(16)*EG(23)/EG(14)/EG(21)
      RB(51) = RF(51) / MAX(EQK, SMALL)
C     R52: O2 + CH2 = 2H + CO2
      RF(52) = EXP(2.860180003308677D1 -7.548250014598796D2*TI)
      EQK = EG(20)*EG(20)*EG(26)/EG(14)/EG(21)*PFAC1
      RB(52) = RF(52) / MAX(EQK, SMALL)
C     R53: CO + CH2 = CH2CO
      RF(53) = EXP(2.74203000846129D1 +5.D-1*ALOGT 
     * -2.269507171056038D3*TI)
      EQK = EG(30)/EG(9)/EG(21)/PFAC1
      RB(53) = RF(53) / MAX(EQK, SMALL)
C     R54: 2CH2 = H2 + C2H2
      RF(54) = 3.2D13
      EQK = EG(7)*EG(11)/EG(21)/EG(21)
      RB(54) = RF(54) / MAX(EQK, SMALL)
C     R55: CH2* = CH2
      RF(55) = EXP(3.033907131703076D1 -3.019300005839518D2*TI)
      EQK = EG(21)/EG(22)
      RB(55) = RF(55) / MAX(EQK, SMALL)
C     R56: O + CH2* = H2 + CO
      RF(56) = 1.5D13
      EQK = EG(7)*EG(9)/EG(15)/EG(22)
      RB(56) = RF(56) / MAX(EQK, SMALL)
C     R57: OH + CH2* = H + CH2O
      RF(57) = 3.D13
      EQK = EG(20)*EG(24)/EG(16)/EG(22)
      RB(57) = RF(57) / MAX(EQK, SMALL)
C     R58: H2 + CH2* = CH3 + H
      RF(58) = 7.D13
      EQK = EG(13)*EG(20)/EG(7)/EG(22)
      RB(58) = RF(58) / MAX(EQK, SMALL)
C     R59: O2 + CH2* = CO + OH + H
      RF(59) = 2.8D13
      EQK = EG(9)*EG(16)*EG(20)/EG(14)/EG(22)*PFAC1
      RB(59) = RF(59) / MAX(EQK, SMALL)
C     R60: O2 + CH2* = CO + H2O
      RF(60) = 1.2D13
      EQK = EG(9)*EG(18)/EG(14)/EG(22)
      RB(60) = RF(60) / MAX(EQK, SMALL)
C     R61: CH2* = CH2
      RF(61) = 3.D13
      EQK = EG(21)/EG(22)
      RB(61) = RF(61) / MAX(EQK, SMALL)
C     R62: CH2* = CH2
      RF(62) = 9.D12
      EQK = EG(21)/EG(22)
      RB(62) = RF(62) / MAX(EQK, SMALL)
C     R63: CH2* = CH2
      RF(63) = 7.D12
      EQK = EG(21)/EG(22)
      RB(63) = RF(63) / MAX(EQK, SMALL)
C     R64: CH2* + CO2 = CO + CH2O
      RF(64) = 1.4D13
      EQK = EG(9)*EG(24)/EG(22)/EG(26)
      RB(64) = RF(64) / MAX(EQK, SMALL)
C     R65: H + CH2O = CH3O
      RF(65) = EXP(2.701483497650473D1 +4.54D-1*ALOGT 
     * -1.308363335863791D3*TI)
      EQK = EG(25)/EG(20)/EG(24)/PFAC1
      RB(65) = RF(65) / MAX(EQK, SMALL)
C     R66: H + CH2O = H2 + HCO
      RF(66) = EXP(2.385876005287556D1 +1.05D0*ALOGT 
     * -1.648034586520737D3*TI)
      EQK = EG(7)*EG(23)/EG(20)/EG(24)
      RB(66) = RF(66) / MAX(EQK, SMALL)
C     R67: O + CH2O = OH + HCO
      RF(67) = EXP(3.129458276205819D1 -1.781387003445316D3*TI)
      EQK = EG(16)*EG(23)/EG(15)/EG(24)
      RB(67) = RF(67) / MAX(EQK, SMALL)
C     R68: OH + CH2O = H2O + HCO
      RF(68) = EXP(2.195582609812426D1 +1.18D0*ALOGT 
     * +2.249378504350441D2*TI)
      EQK = EG(18)*EG(23)/EG(16)/EG(24)
      RB(68) = RF(68) / MAX(EQK, SMALL)
C     R69: O2 + CH2O = HO2 + HCO
      RF(69) = EXP(3.223619130191664D1 -2.012866670559679D4*TI)
      EQK = EG(17)*EG(23)/EG(14)/EG(24)
      RB(69) = RF(69) / MAX(EQK, SMALL)
C     R70: HO2 + CH2O = H2O2 + HCO
      RF(70) = EXP(2.763102111592855D1 -4.025733341119358D3*TI)
      EQK = EG(19)*EG(23)/EG(17)/EG(24)
      RB(70) = RF(70) / MAX(EQK, SMALL)
C     R71: CH3 + H = CH4
      RF(71) = EXP(3.708037838837523D1 -6.3D-1*ALOGT 
     * -1.927319837060892D2*TI)
      EQK = EG(3)/EG(13)/EG(20)/PFAC1
      RB(71) = RF(71) / MAX(EQK, SMALL)
C     R72: CH3 + O = H + CH2O
      RF(72) = 8.43D13
      EQK = EG(20)*EG(24)/EG(13)/EG(15)
      RB(72) = RF(72) / MAX(EQK, SMALL)
C     R73: CH3 + OH = H2O + CH2
      RF(73) = EXP(1.784086224869942D1 +1.6D0*ALOGT 
     * -2.727434338608365D3*TI)
      EQK = EG(18)*EG(21)/EG(13)/EG(16)
      RB(73) = RF(73) / MAX(EQK, SMALL)
C     R74: CH3 + OH = H2O + CH2*
      RF(74) = 2.501D13
      EQK = EG(18)*EG(22)/EG(13)/EG(16)
      RB(74) = RF(74) / MAX(EQK, SMALL)
C     R75: CH3 + O2 = O + CH3O
      RF(75) = EXP(3.105950935782661D1 -1.449264002802969D4*TI)
      EQK = EG(15)*EG(25)/EG(13)/EG(14)
      RB(75) = RF(75) / MAX(EQK, SMALL)
C     R76: CH3 + O2 = OH + CH2O
      RF(76) = EXP(2.430678477540252D1 -4.498757008700882D3*TI)
      EQK = EG(16)*EG(24)/EG(13)/EG(14)
      RB(76) = RF(76) / MAX(EQK, SMALL)
C     R77: CH3 + HO2 = CH4 + O2
      RF(77) = 1.D12
      EQK = EG(3)*EG(14)/EG(13)/EG(17)
      RB(77) = RF(77) / MAX(EQK, SMALL)
C     R78: CH3 + HO2 = OH + CH3O
      RF(78) = 1.34D13
      EQK = EG(16)*EG(25)/EG(13)/EG(17)
      RB(78) = RF(78) / MAX(EQK, SMALL)
C     R79: CH3 + H2O2 = CH4 + HO2
      RF(79) = EXP(1.010642839653282D1 +2.47D0*ALOGT 
     * -2.606662338374784D3*TI)
      EQK = EG(3)*EG(17)/EG(13)/EG(19)
      RB(79) = RF(79) / MAX(EQK, SMALL)
C     R80: CH3 + HCO = CH4 + CO
      RF(80) = 8.48D12
      EQK = EG(3)*EG(9)/EG(13)/EG(23)
      RB(80) = RF(80) / MAX(EQK, SMALL)
C     R81: CH3 + CH2O = CH4 + HCO
      RF(81) = EXP(8.107720061910534D0 +2.81D0*ALOGT 
     * -2.94884967236993D3*TI)
      EQK = EG(3)*EG(23)/EG(13)/EG(24)
      RB(81) = RF(81) / MAX(EQK, SMALL)
C     R82: CH3 + CH2 = C2H4 + H
      RF(82) = 4.D13
      EQK = EG(2)*EG(20)/EG(13)/EG(21)
      RB(82) = RF(82) / MAX(EQK, SMALL)
C     R83: 2CH3 = C2H6
      RF(83) = EXP(3.759277757658865D1 -9.7D-1*ALOGT 
     * -3.119943339367502D2*TI)
      EQK = EG(8)/EG(13)/EG(13)/PFAC1
      RB(83) = RF(83) / MAX(EQK, SMALL)
C     R84: 2CH3 = H + C2H5
      RF(84) = EXP(2.923845702569198D1 +1.D-1*ALOGT 
     * -5.334096676983148D3*TI)
      EQK = EG(20)*EG(28)/EG(13)/EG(13)
      RB(84) = RF(84) / MAX(EQK, SMALL)
C     R85: CH3 + HCCO = C2H4 + CO
      RF(85) = 5.D13
      EQK = EG(2)*EG(9)/EG(13)/EG(29)
      RB(85) = RF(85) / MAX(EQK, SMALL)
C     R86: H + CH3O = H2 + CH2O
      RF(86) = 2.D13
      EQK = EG(7)*EG(24)/EG(20)/EG(25)
      RB(86) = RF(86) / MAX(EQK, SMALL)
C     R87: H + CH3O = CH3 + OH
      RF(87) = 3.2D13
      EQK = EG(13)*EG(16)/EG(20)/EG(25)
      RB(87) = RF(87) / MAX(EQK, SMALL)
C     R88: H + CH3O = H2O + CH2*
      RF(88) = 1.6D13
      EQK = EG(18)*EG(22)/EG(20)/EG(25)
      RB(88) = RF(88) / MAX(EQK, SMALL)
C     R89: OH + CH3O = H2O + CH2O
      RF(89) = 5.D12
      EQK = EG(18)*EG(24)/EG(16)/EG(25)
      RB(89) = RF(89) / MAX(EQK, SMALL)
C     R90: O2 + CH3O = HO2 + CH2O
      RF(90) = EXP(-2.847965319932889D1 +7.6D0*ALOGT 
     * +1.776354836768916D3*TI)
      EQK = EG(17)*EG(24)/EG(14)/EG(25)
      RB(90) = RF(90) / MAX(EQK, SMALL)
C     R91: CH4 + H = H2 + CH3
      RF(91) = EXP(2.030775039298474D1 +1.62D0*ALOGT 
     * -5.454868677216729D3*TI)
      EQK = EG(7)*EG(13)/EG(3)/EG(20)
      RB(91) = RF(91) / MAX(EQK, SMALL)
C     R92: CH4 + O = CH3 + OH
      RF(92) = EXP(2.074306846424259D1 +1.5D0*ALOGT 
     * -4.327663341703309D3*TI)
      EQK = EG(13)*EG(16)/EG(3)/EG(15)
      RB(92) = RF(92) / MAX(EQK, SMALL)
C     R93: CH4 + OH = CH3 + H2O
      RF(93) = EXP(1.842068074395237D1 +1.6D0*ALOGT 
     * -1.570036003036549D3*TI)
      EQK = EG(13)*EG(18)/EG(3)/EG(16)
      RB(93) = RF(93) / MAX(EQK, SMALL)
C     R94: CH4 + CH2 = 2CH3
      RF(94) = EXP(1.471567190790855D1 +2.D0*ALOGT 
     * -4.161601841382136D3*TI)
      EQK = EG(13)*EG(13)/EG(3)/EG(21)
      RB(94) = RF(94) / MAX(EQK, SMALL)
C     R95: CH4 + CH2* = 2CH3
      RF(95) = EXP(3.040360983816833D1 +2.868335005547542D2*TI)
      EQK = EG(13)*EG(13)/EG(3)/EG(22)
      RB(95) = RF(95) / MAX(EQK, SMALL)
C     R96: H + HCCO = CO + CH2*
      RF(96) = 1.D14
      EQK = EG(9)*EG(22)/EG(20)/EG(29)
      RB(96) = RF(96) / MAX(EQK, SMALL)
C     R97: O + HCCO = 2CO + H
      RF(97) = 1.D14
      EQK = EG(9)*EG(9)*EG(20)/EG(15)/EG(29)*PFAC1
      RB(97) = RF(97) / MAX(EQK, SMALL)
C     R98: O2 + HCCO = 2CO + OH
      RF(98) = EXP(2.810102474517428D1 -4.297470341644914D2*TI)
      EQK = EG(9)*EG(9)*EG(16)/EG(14)/EG(29)*PFAC1
      RB(98) = RF(98) / MAX(EQK, SMALL)
C     R99: CH2 + HCCO = CO + C2H3
      RF(99) = 3.D13
      EQK = EG(9)*EG(27)/EG(21)/EG(29)
      RB(99) = RF(99) / MAX(EQK, SMALL)
C     R100: C2H3 = C2H2 + H
      RF(100) = EXP(1.97713479274291D1 +1.62D0*ALOGT 
     * -1.864327174605727D4*TI)
      EQK = EG(11)*EG(20)/EG(27)*PFAC1
      RB(100) = RF(100) / MAX(EQK, SMALL)
C     R101: C2H2 + O = CO + CH2
      RF(101) = EXP(1.522160754638034D1 +2.D0*ALOGT 
     * -9.561116685158474D2*TI)
      EQK = EG(9)*EG(21)/EG(11)/EG(15)
      RB(101) = RF(101) / MAX(EQK, SMALL)
C     R102: C2H2 + O = H + HCCO
      RF(102) = EXP(1.660790190750023D1 +2.D0*ALOGT 
     * -9.561116685158474D2*TI)
      EQK = EG(20)*EG(29)/EG(11)/EG(15)
      RB(102) = RF(102) / MAX(EQK, SMALL)
C     R103: C2H2 + OH = H + CH2CO
      RF(103) = EXP(-8.431015495175185D0 +4.5D0*ALOGT 
     * +5.032166676399197D2*TI)
      EQK = EG(20)*EG(30)/EG(11)/EG(16)
      RB(103) = RF(103) / MAX(EQK, SMALL)
C     R104: C2H2 + OH = CO + CH3
      RF(104) = EXP(-7.635493904311701D0 +4.D0*ALOGT 
     * +1.006433335279839D3*TI)
      EQK = EG(9)*EG(13)/EG(11)/EG(16)
      RB(104) = RF(104) / MAX(EQK, SMALL)
C     R105: C2H2 + HCO = CO + C2H3
      RF(105) = EXP(1.611809565095832D1 +2.D0*ALOGT 
     * -3.019300005839518D3*TI)
      EQK = EG(9)*EG(27)/EG(11)/EG(23)
      RB(105) = RF(105) / MAX(EQK, SMALL)
C     R106: C2H2 + CH2 = H + C3H3
      RF(106) = EXP(3.011592776571655D1 -3.331294339776268D3*TI)
      EQK = EG(20)*EG(32)/EG(11)/EG(21)
      RB(106) = RF(106) / MAX(EQK, SMALL)
C     R107: C2H2 + CH2* = H + C3H3
      RF(107) = 2.D13
      EQK = EG(20)*EG(32)/EG(11)/EG(22)
      RB(107) = RF(107) / MAX(EQK, SMALL)
C     R108: C2H2 + CH3 = aC3H5
      RF(108) = EXP(1.230228267232072D2 -1.282D1*ALOGT 
     * -1.797993153477433D4*TI)
      EQK = EG(33)/EG(11)/EG(13)/PFAC1
      RB(108) = RF(108) / MAX(EQK, SMALL)
C     R109: H + CH2CO = CH2CHO
      RF(109) = EXP(3.343011377038908D1 -6.D-2*ALOGT 
     * -4.277341674939317D3*TI)
      EQK = EG(31)/EG(20)/EG(30)/PFAC1
      RB(109) = RF(109) / MAX(EQK, SMALL)
C     R110: H + CH2CO = H2 + HCCO
      RF(110) = EXP(3.154304412135669D1 -4.025733341119358D3*TI)
      EQK = EG(7)*EG(29)/EG(20)/EG(30)
      RB(110) = RF(110) / MAX(EQK, SMALL)
C     R111: H + CH2CO = CO + CH3
      RF(111) = EXP(2.112873094505457D1 +1.43D0*ALOGT 
     * -1.353652835951384D3*TI)
      EQK = EG(9)*EG(13)/EG(20)/EG(30)
      RB(111) = RF(111) / MAX(EQK, SMALL)
C     R112: O + CH2CO = OH + HCCO
      RF(112) = EXP(2.993360620892259D1 -4.025733341119358D3*TI)
      EQK = EG(16)*EG(29)/EG(15)/EG(30)
      RB(112) = RF(112) / MAX(EQK, SMALL)
C     R113: O + CH2CO = CH2 + CO2
      RF(113) = EXP(2.819063690386397D1 -6.793425013138916D2*TI)
      EQK = EG(21)*EG(26)/EG(15)/EG(30)
      RB(113) = RF(113) / MAX(EQK, SMALL)
C     R114: OH + CH2CO = H2O + HCCO
      RF(114) = EXP(2.964592413647081D1 -1.006433335279839D3*TI)
      EQK = EG(18)*EG(29)/EG(16)/EG(30)
      RB(114) = RF(114) / MAX(EQK, SMALL)
C     R115: H + C2H3 = C2H4
      RF(115) = EXP(2.943602581190662D1 +2.7D-1*ALOGT 
     * -1.409006669391775D2*TI)
      EQK = EG(2)/EG(20)/EG(27)/PFAC1
      RB(115) = RF(115) / MAX(EQK, SMALL)
C     R116: H + C2H3 = H2 + C2H2
      RF(116) = 9.D13
      EQK = EG(7)*EG(11)/EG(20)/EG(27)
      RB(116) = RF(116) / MAX(EQK, SMALL)
C     R117: O + C2H3 = H + CH2CO
      RF(117) = 4.8D13
      EQK = EG(20)*EG(30)/EG(15)/EG(27)
      RB(117) = RF(117) / MAX(EQK, SMALL)
C     R118: O + C2H3 = CO + CH3
      RF(118) = 4.8D13
      EQK = EG(9)*EG(13)/EG(15)/EG(27)
      RB(118) = RF(118) / MAX(EQK, SMALL)
C     R119: OH + C2H3 = C2H2 + H2O
      RF(119) = 3.011D13
      EQK = EG(11)*EG(18)/EG(16)/EG(27)
      RB(119) = RF(119) / MAX(EQK, SMALL)
C     R120: O2 + C2H3 = C2H2 + HO2
      RF(120) = EXP(1.410818017192709D1 +1.61D0*ALOGT 
     * +1.929332703731452D2*TI)
      EQK = EG(11)*EG(17)/EG(14)/EG(27)
      RB(120) = RF(120) / MAX(EQK, SMALL)
C     R121: O2 + C2H3 = O + CH2CHO
      RF(121) = EXP(2.642704831160261D1 +2.9D-1*ALOGT 
     * -5.535383344039117D0*TI)
      EQK = EG(15)*EG(31)/EG(14)/EG(27)
      RB(121) = RF(121) / MAX(EQK, SMALL)
C     R122: O2 + C2H3 = HCO + CH2O
      RF(122) = EXP(3.836741779139978D1 -1.39D0*ALOGT 
     * -5.082488343163189D2*TI)
      EQK = EG(23)*EG(24)/EG(14)/EG(27)
      RB(122) = RF(122) / MAX(EQK, SMALL)
C     R123: CH3 + C2H3 = C3H6
      RF(123) = 2.5D13
      EQK = EG(4)/EG(13)/EG(27)/PFAC1
      RB(123) = RF(123) / MAX(EQK, SMALL)
C     R124: CH3 + C2H3 = H + aC3H5
      RF(124) = EXP(5.566750733996526D1 -2.83D0*ALOGT 
     * -9.368887918120025D3*TI)
      EQK = EG(20)*EG(33)/EG(13)/EG(27)
      RB(124) = RF(124) / MAX(EQK, SMALL)
C     R125: CH2CHO = CO + CH3
      RF(125) = EXP(9.646011254645141D1 -9.147D0*ALOGT 
     * -2.360086171231223D4*TI)
      EQK = EG(9)*EG(13)/EG(31)*PFAC1
      RB(125) = RF(125) / MAX(EQK, SMALL)
C     R126: H + CH2CHO = CH3 + HCO
      RF(126) = 9.D13
      EQK = EG(13)*EG(23)/EG(20)/EG(31)
      RB(126) = RF(126) / MAX(EQK, SMALL)
C     R127: H + CH2CHO = H2 + CH2CO
      RF(127) = EXP(3.062675338948254D1 -2.012866670559679D3*TI)
      EQK = EG(7)*EG(30)/EG(20)/EG(31)
      RB(127) = RF(127) / MAX(EQK, SMALL)
C     R128: O + CH2CHO = OH + CH2CO
      RF(128) = EXP(3.062675338948254D1 -2.012866670559679D3*TI)
      EQK = EG(16)*EG(30)/EG(15)/EG(31)
      RB(128) = RF(128) / MAX(EQK, SMALL)
C     R129: OH + CH2CHO = H2O + CH2CO
      RF(129) = EXP(2.993360620892259D1 -1.006433335279839D3*TI)
      EQK = EG(18)*EG(30)/EG(16)/EG(31)
      RB(129) = RF(129) / MAX(EQK, SMALL)
C     R130: O2 + CH2CHO = HO2 + CH2CO
      RF(130) = 1.4D11
      EQK = EG(17)*EG(30)/EG(14)/EG(31)
      RB(130) = RF(130) / MAX(EQK, SMALL)
C     R131: O2 + CH2CHO = CO + OH + CH2O
      RF(131) = 1.8D10
      EQK = EG(9)*EG(16)*EG(24)/EG(14)/EG(31)*PFAC1
      RB(131) = RF(131) / MAX(EQK, SMALL)
C     R132: C2H4 + H = C2H5
      RF(132) = EXP(2.103807657678644D1 +1.463D0*ALOGT 
     * -6.818585846520912D2*TI)
      EQK = EG(28)/EG(2)/EG(20)/PFAC1
      RB(132) = RF(132) / MAX(EQK, SMALL)
C     R133: C2H4 + H = H2 + C2H3
      RF(133) = EXP(1.774143646856141D1 +1.9D0*ALOGT 
     * -6.51665584593696D3*TI)
      EQK = EG(7)*EG(27)/EG(2)/EG(20)
      RB(133) = RF(133) / MAX(EQK, SMALL)
C     R134: C2H4 + O = OH + C2H3
      RF(134) = EXP(1.653020530178515D1 +1.9D0*ALOGT 
     * -1.8820303369733D3*TI)
      EQK = EG(16)*EG(27)/EG(2)/EG(15)
      RB(134) = RF(134) / MAX(EQK, SMALL)
C     R135: C2H4 + O = CH3 + HCO
      RF(135) = EXP(1.677042083699801D1 +1.83D0*ALOGT 
     * -1.107076668807823D2*TI)
      EQK = EG(13)*EG(23)/EG(2)/EG(15)
      RB(135) = RF(135) / MAX(EQK, SMALL)
C     R136: C2H4 + O = CH2 + CH2O
      RF(136) = EXP(1.285839783156986D1 +1.83D0*ALOGT 
     * -1.107076668807823D2*TI)
      EQK = EG(21)*EG(24)/EG(2)/EG(15)
      RB(136) = RF(136) / MAX(EQK, SMALL)
C     R137: C2H4 + OH = H2O + C2H3
      RF(137) = EXP(1.509644440342634D1 +2.D0*ALOGT 
     * -1.258041669099799D3*TI)
      EQK = EG(18)*EG(27)/EG(2)/EG(16)
      RB(137) = RF(137) / MAX(EQK, SMALL)
C     R138: C2H4 + HCO = CO + C2H5
      RF(138) = EXP(1.611809565095832D1 +2.D0*ALOGT 
     * -4.025733341119358D3*TI)
      EQK = EG(9)*EG(28)/EG(2)/EG(23)
      RB(138) = RF(138) / MAX(EQK, SMALL)
C     R139: C2H4 + CH2 = H + aC3H5
      RF(139) = EXP(3.062675338948254D1 -3.019300005839518D3*TI)
      EQK = EG(20)*EG(33)/EG(2)/EG(21)
      RB(139) = RF(139) / MAX(EQK, SMALL)
C     R140: C2H4 + CH2* = H + aC3H5
      RF(140) = 5.D13
      EQK = EG(20)*EG(33)/EG(2)/EG(22)
      RB(140) = RF(140) / MAX(EQK, SMALL)
C     R141: C2H4 + CH3 = CH4 + C2H3
      RF(141) = EXP(1.233270529646354D1 +2.D0*ALOGT 
     * -4.629593342287261D3*TI)
      EQK = EG(3)*EG(27)/EG(2)/EG(13)
      RB(141) = RF(141) / MAX(EQK, SMALL)
C     R142: C2H4 + O2 = HO2 + C2H3
      RF(142) = EXP(3.137344133697052D1 -3.059557339250712D4*TI)
      EQK = EG(17)*EG(27)/EG(2)/EG(14)
      RB(142) = RF(142) / MAX(EQK, SMALL)
C     R143: H + C2H5 = C2H6
      RF(143) = EXP(4.079452643666405D1 -9.9D-1*ALOGT 
     * -7.95082334871073D2*TI)
      EQK = EG(8)/EG(20)/EG(28)/PFAC1
      RB(143) = RF(143) / MAX(EQK, SMALL)
C     R144: H + C2H5 = C2H4 + H2
      RF(144) = 2.D12
      EQK = EG(2)*EG(7)/EG(20)/EG(28)
      RB(144) = RF(144) / MAX(EQK, SMALL)
C     R145: O + C2H5 = CH3 + CH2O
      RF(145) = 1.604D13
      EQK = EG(13)*EG(24)/EG(15)/EG(28)
      RB(145) = RF(145) / MAX(EQK, SMALL)
C     R146: O2 + C2H5 = C2H4 + HO2
      RF(146) = 2.D10
      EQK = EG(2)*EG(17)/EG(14)/EG(28)
      RB(146) = RF(146) / MAX(EQK, SMALL)
C     R147: HO2 + C2H5 = CH3 + OH + CH2O
      RF(147) = 2.4D13
      EQK = EG(13)*EG(16)*EG(24)/EG(17)/EG(28)*PFAC1
      RB(147) = RF(147) / MAX(EQK, SMALL)
C     R148: C2H3 + C2H5 = C4H81
      RF(148) = 1.5D13
      EQK = EG(6)/EG(27)/EG(28)/PFAC1
      RB(148) = RF(148) / MAX(EQK, SMALL)
C     R149: C2H3 + C2H5 = CH3 + aC3H5
      RF(149) = EXP(7.504369952894507D1 -5.22D0*ALOGT 
     * -9.937019535885494D3*TI)
      EQK = EG(13)*EG(33)/EG(27)/EG(28)
      RB(149) = RF(149) / MAX(EQK, SMALL)
C     R150: C2H6 + H = H2 + C2H5
      RF(150) = EXP(1.856044268632752D1 +1.9D0*ALOGT 
     * -3.789221507328595D3*TI)
      EQK = EG(7)*EG(28)/EG(8)/EG(20)
      RB(150) = RF(150) / MAX(EQK, SMALL)
C     R151: C2H6 + O = OH + C2H5
      RF(151) = EXP(1.831309553327243D1 +1.92D0*ALOGT 
     * -2.863302838871143D3*TI)
      EQK = EG(16)*EG(28)/EG(8)/EG(15)
      RB(151) = RF(151) / MAX(EQK, SMALL)
C     R152: C2H6 + OH = H2O + C2H5
      RF(152) = EXP(1.507963728510996D1 +2.12D0*ALOGT 
     * -4.377985008467301D2*TI)
      EQK = EG(18)*EG(28)/EG(8)/EG(16)
      RB(152) = RF(152) / MAX(EQK, SMALL)
C     R153: C2H6 + CH3 = CH4 + C2H5
      RF(153) = EXP(1.563033530012333D1 +1.74D0*ALOGT 
     * -5.258614176837161D3*TI)
      EQK = EG(3)*EG(28)/EG(8)/EG(13)
      RB(153) = RF(153) / MAX(EQK, SMALL)
C     R154: O2 + C3H3 = HCO + CH2CO
      RF(154) = EXP(2.412446321860857D1 -1.44322540279129D3*TI)
      EQK = EG(23)*EG(30)/EG(14)/EG(32)
      RB(154) = RF(154) / MAX(EQK, SMALL)
C     R155: HO2 + C3H3 = CO + OH + C2H3
      RF(155) = 8.D11
      EQK = EG(9)*EG(16)*EG(27)/EG(17)/EG(32)*PFAC1
      RB(155) = RF(155) / MAX(EQK, SMALL)
C     R156: C2H2 + C3H3 = C5H5
      RF(156) = EXP(1.285693442209068D2 -1.25D1*ALOGT 
     * -2.114768045756763D4*TI)
      EQK = EG(35)/EG(11)/EG(32)/PFAC1
      RB(156) = RF(156) / MAX(EQK, SMALL)
C     R157: 2C3H3 = H + C6H5
      RF(157) = 5.D12
      RB(157) = 0D0
C     R158: 2C3H3 = C6H6
      RF(158) = 2.D12
      RB(158) = 0D0
C     R159: H + aC3H5 = C3H6
      RF(159) = 2.D14
      EQK = EG(4)/EG(20)/EG(33)/PFAC1
      RB(159) = RF(159) / MAX(EQK, SMALL)
C     R160: HO2 + aC3H5 = C3H6 + O2
      RF(160) = 2.66D12
      EQK = EG(4)*EG(14)/EG(17)/EG(33)
      RB(160) = RF(160) / MAX(EQK, SMALL)
C     R161: HO2 + aC3H5 = OH + CH2O + C2H3
      RF(161) = 6.6D12
      EQK = EG(16)*EG(24)*EG(27)/EG(17)/EG(33)*PFAC1
      RB(161) = RF(161) / MAX(EQK, SMALL)
C     R162: HCO + aC3H5 = C3H6 + CO
      RF(162) = 6.D13
      EQK = EG(4)*EG(9)/EG(23)/EG(33)
      RB(162) = RF(162) / MAX(EQK, SMALL)
C     R163: CH3 + aC3H5 = C4H81
      RF(163) = EXP(3.223619130191664D1 -3.2D-1*ALOGT 
     * +1.319937319219509D2*TI)
      EQK = EG(6)/EG(13)/EG(33)/PFAC1
      RB(163) = RF(163) / MAX(EQK, SMALL)
C     R164: C3H6 + H = C2H4 + CH3
      RF(164) = EXP(5.043372849455479D1 -2.39D0*ALOGT 
     * -5.625962344214302D3*TI)
      EQK = EG(2)*EG(13)/EG(4)/EG(20)
      RB(164) = RF(164) / MAX(EQK, SMALL)
C     R165: C3H6 + H = H2 + aC3H5
      RF(165) = EXP(1.446783574400397D1 +2.309D0*ALOGT 
     * -2.352034704548985D3*TI)
      EQK = EG(7)*EG(33)/EG(4)/EG(20)
      RB(165) = RF(165) / MAX(EQK, SMALL)
C     R166: C3H6 + O = CH3 + H + CH2CO
      RF(166) = EXP(1.819753719263815D1 +1.65D0*ALOGT 
     * -1.645518503182537D2*TI)
      EQK = EG(13)*EG(20)*EG(30)/EG(4)/EG(15)*PFAC1
      RB(166) = RF(166) / MAX(EQK, SMALL)
C     R167: C3H6 + O = HCO + C2H5
      RF(167) = EXP(1.737085861945369D1 +1.65D0*ALOGT 
     * +4.891266009460019D2*TI)
      EQK = EG(23)*EG(28)/EG(4)/EG(15)
      RB(167) = RF(167) / MAX(EQK, SMALL)
C     R168: C3H6 + O = OH + aC3H5
      RF(168) = EXP(2.591622268783662D1 +7.D-1*ALOGT 
     * -2.958914005722728D3*TI)
      EQK = EG(16)*EG(33)/EG(4)/EG(15)
      RB(168) = RF(168) / MAX(EQK, SMALL)
C     R169: C3H6 + OH = H2O + aC3H5
      RF(169) = EXP(1.494691266945537D1 +2.D0*ALOGT 
     * +1.499585669566961D2*TI)
      EQK = EG(18)*EG(33)/EG(4)/EG(16)
      RB(169) = RF(169) / MAX(EQK, SMALL)
C     R170: C3H6 + CH3 = CH4 + aC3H5
      RF(170) = EXP(3.815512105047302D0 +3.268D0*ALOGT 
     * -3.599005606960705D3*TI)
      EQK = EG(3)*EG(33)/EG(4)/EG(13)
      RB(170) = RF(170) / MAX(EQK, SMALL)
C     R171: C4H81 + H = C2H4 + C2H5
      RF(171) = EXP(5.112687567511474D1 -2.39D0*ALOGT 
     * -5.625962344214302D3*TI)
      EQK = EG(2)*EG(28)/EG(6)/EG(20)
      RB(171) = RF(171) / MAX(EQK, SMALL)
C     R172: C4H81 + H = C3H6 + CH3
      RF(172) = EXP(5.182002285567469D1 -2.39D0*ALOGT 
     * -5.625962344214302D3*TI)
      EQK = EG(4)*EG(13)/EG(6)/EG(20)
      RB(172) = RF(172) / MAX(EQK, SMALL)
C     R173: iC4H8 + H = C3H6 + CH3
      RF(173) = EXP(5.247494882341444D1 -2.538D0*ALOGT 
     * -5.851906627984626D3*TI)
      EQK = EG(4)*EG(13)/EG(5)/EG(20)
      RB(173) = RF(173) / MAX(EQK, SMALL)
C     R174: iC4H8 + O = 2CH3 + CH2CO
      RF(174) = EXP(1.929614948130627D1 +1.65D0*ALOGT 
     * -1.645518503182537D2*TI)
      EQK = EG(13)*EG(13)*EG(30)/EG(5)/EG(15)*PFAC1
      RB(174) = RF(174) / MAX(EQK, SMALL)
C     R175: iC4H8 + HO2 = C3H6 + OH + CH2O
      RF(175) = EXP(2.841034599272955D1 -7.145676680486859D3*TI)
      EQK = EG(4)*EG(16)*EG(24)/EG(5)/EG(17)*PFAC1
      RB(175) = RF(175) / MAX(EQK, SMALL)
C     R176: CH3 + C6H5 = C6H5CH3
      RF(176) = EXP(3.025568970809171D1 -2.31479667114363D1*TI)
      EQK = EG(12)/EG(13)/EG(36)/PFAC1
      RB(176) = RF(176) / MAX(EQK, SMALL)
C     R177: C6H5CH3 + H = C6H6 + CH3
      RF(177) = EXP(1.447303056088107D1 +2.17D0*ALOGT 
     * -2.094890987384986D3*TI)
      EQK = EG(10)*EG(13)/EG(12)/EG(20)
      RB(177) = RF(177) / MAX(EQK, SMALL)
C     R178: H + C6H5 = C6H6
      RF(178) = 1.D14
      EQK = EG(10)/EG(20)/EG(36)/PFAC1
      RB(178) = RF(178) / MAX(EQK, SMALL)
C     R179: C6H6 + OH = H2O + C6H5
      RF(179) = EXP(1.289546277721241D1 +2.286D0*ALOGT 
     * -5.32403234363035D2*TI)
      EQK = EG(18)*EG(36)/EG(10)/EG(16)
      RB(179) = RF(179) / MAX(EQK, SMALL)
C     R180: C6H6 + O = H + C6H5O
      RF(180) = EXP(3.026290995606519D1 -2.470793838112006D3*TI)
      EQK = EG(20)*EG(37)/EG(10)/EG(15)
      RB(180) = RF(180) / MAX(EQK, SMALL)
C     R181: C6H6 + O = HCO + C5H5
      RF(181) = EXP(3.026290995606519D1 -2.279571504408836D3*TI)
      EQK = EG(23)*EG(35)/EG(10)/EG(15)
      RB(181) = RF(181) / MAX(EQK, SMALL)
C     R182: H2 + C6H5 = C6H6 + H
      RF(182) = EXP(1.095203386353075D1 +2.43D0*ALOGT 
     * -3.156678156105216D3*TI)
      EQK = EG(10)*EG(20)/EG(7)/EG(36)
      RB(182) = RF(182) / MAX(EQK, SMALL)
C     R183: O2 + C6H5 = O + C6H5O
      RF(183) = EXP(3.088911765395003D1 -3.079686005956308D3*TI)
      EQK = EG(15)*EG(37)/EG(14)/EG(36)
      RB(183) = RF(183) / MAX(EQK, SMALL)
C     R184: O + C6H5 = CO + C5H5
      RF(184) = 1.D14
      EQK = EG(9)*EG(35)/EG(15)/EG(36)
      RB(184) = RF(184) / MAX(EQK, SMALL)
C     R185: OH + C6H5 = H + C6H5O
      RF(185) = 3.D13
      EQK = EG(20)*EG(37)/EG(16)/EG(36)
      RB(185) = RF(185) / MAX(EQK, SMALL)
C     R186: HO2 + C6H5 = C6H6 + O2
      RF(186) = 1.D12
      EQK = EG(10)*EG(14)/EG(17)/EG(36)
      RB(186) = RF(186) / MAX(EQK, SMALL)
C     R187: CH4 + C6H5 = C6H6 + CH3
      RF(187) = EXP(-5.549346121351782D0 +4.57D0*ALOGT 
     * -2.644906805115418D3*TI)
      EQK = EG(10)*EG(13)/EG(3)/EG(36)
      RB(187) = RF(187) / MAX(EQK, SMALL)
C     R188: C6H5O = CO + C5H5
      RF(188) = EXP(1.256640139790803D2 -1.206D1*ALOGT 
     * -3.663417340418615D4*TI)
      EQK = EG(9)*EG(35)/EG(37)*PFAC1
      RB(188) = RF(188) / MAX(EQK, SMALL)
C     R189: O + C5H5 = H + C5H4O
      RF(189) = EXP(3.169318677978641D1 -2.D-2*ALOGT 
     * -1.006433335279839D1*TI)
      EQK = EG(20)*EG(34)/EG(15)/EG(35)
      RB(189) = RF(189) / MAX(EQK, SMALL)
C     R190: O2 + C5H5 = OH + C5H4O
      RF(190) = EXP(2.821880778083067D1 +8.D-2*ALOGT 
     * -9.057900017518554D3*TI)
      EQK = EG(16)*EG(34)/EG(14)/EG(35)
      RB(190) = RF(190) / MAX(EQK, SMALL)
C     R191: C5H4O = CO + 2C2H2
      RF(191) = EXP(9.623053810480693D1 -7.87D0*ALOGT 
     * -4.966748509606007D4*TI)
      EQK = EG(9)*EG(11)*EG(11)/EG(34)*PFAC2
      RB(191) = RF(191) / MAX(EQK, SMALL)
C     R192: O + C5H4O = CO + HCO + C3H3
      RF(192) = EXP(2.024523003600341D1 +1.45D0*ALOGT 
     * +4.317599008350511D2*TI)
      EQK = EG(9)*EG(23)*EG(32)/EG(15)/EG(34)*PFAC1
      RB(192) = RF(192) / MAX(EQK, SMALL)
C
      RKLOW(1) = EXP(4.559410099735222D1 -1.4D0*ALOGT)
      RKLOW(2) = EXP(3.984208130296976D1 -5.84D-1*ALOGT 
     * +1.153875818898336D3*TI)
      RKLOW(3) = EXP(5.542160680152843D1 -2.79D0*ALOGT 
     * -2.108981054078904D3*TI)
      RKLOW(4) = EXP(6.379313832844232D1 -3.42D0*ALOGT 
     * -4.244632591542722D4*TI)
      RKLOW(5) = EXP(5.556214682430743D1 -2.57D0*ALOGT 
     * -7.170837513868855D2*TI)
      RKLOW(6) = EXP(6.333294832064492D1 -3.14D0*ALOGT 
     * -6.189565011971013D2*TI)
      RKLOW(7) = EXP(7.697484926241725D1 -5.11D0*ALOGT 
     * -3.57032225690523D3*TI)
      RKLOW(8) = EXP(6.986601015018564D1 -4.8D0*ALOGT 
     * -2.797884672077953D3*TI)
      RKLOW(9) = EXP(7.689235621931073D1 -4.76D0*ALOGT 
     * -1.227848669041404D3*TI)
      RKLOW(10) = EXP(1.15700234196288D2 -9.67D0*ALOGT 
     * -3.1300076727203D3*TI)
      RKLOW(11) = EXP(6.311175598946197D1 -3.4D0*ALOGT 
     * -1.801451258417455D4*TI)
      RKLOW(12) = EXP(9.574098987948821D1 -7.64D0*ALOGT 
     * -5.988278344915044D3*TI)
      RKLOW(13) = EXP(6.941402502644259D1 -3.86D0*ALOGT 
     * -1.670679336564533D3*TI)
      RKLOW(14) = EXP(1.350015492208952D2 -1.194D1*ALOGT 
     * -4.916326199508487D3*TI)
      RKLOW(15) = EXP(9.086552936376022D1 -6.642D0*ALOGT 
     * -2.903056955614697D3*TI)
      RKLOW(16) = EXP(9.509412345149228D1 -7.08D0*ALOGT 
     * -3.364003423172863D3*TI)
      RKLOW(17) = EXP(1.293830201385977D2 -1.179D1*ALOGT 
     * -4.521150150410858D3*TI)
      RKLOW(18) = EXP(1.384402845218764D2 -1.2D1*ALOGT 
     * -3.003096429141513D3*TI)
      RKLOW(19) = EXP(1.3951864295364D2 -1.281D1*ALOGT 
     * -3.145104172749498D3*TI)
      RKLOW(20) = EXP(1.745809516235858D2 -1.63D1*ALOGT 
     * -3.522516673479438D3*TI)
C
      END
C                                                                      C
C----------------------------------------------------------------------C
C                                                                      C
      SUBROUTINE RDSMH  (T, SMH)
      IMPLICIT DOUBLE PRECISION (A-H, O-Z), INTEGER (I-N)
      DIMENSION SMH(*), TN(5)
C
      TLOG = LOG(T)
      TI = 1D0/T
C
      TN(1) = TLOG - 1D0
      TN(2) = T
      TN(3) = TN(2)*T
      TN(4) = TN(3)*T
      TN(5) = TN(4)*T
C C2H4
      IF (T .GT. 1000) THEN
      SMH(2) = 1.03053693D1 -4.93988614D3*TI 
     *         +2.03611116D0*TN(1) +7.32270755D-3*TN(2) 
     *         -1.118463191666667D-6*TN(3) +1.226857691666667D-10*TN(4) 
     *         -6.28530305D-15*TN(5) 
      ELSE
      SMH(2) = 4.09733096D0 -5.08977593D3*TI 
     *         +3.95920148D0*TN(1) -3.785261235D-3*TN(2) 
     *         +9.516504866666667D-6*TN(3) -5.763239608333333D-9*TN(4) 
     *         +1.349421865D-12*TN(5) 
      ENDIF
C CH4
      IF (T .GT. 1000) THEN
      SMH(3) = 1.8437318D1 +9.468344590000001D3*TI 
     *         +7.4851495D-2*TN(1) +6.69547335D-3*TN(2) 
     *         -9.554763483333333D-7*TN(3) +1.019104458333333D-10*TN(4) 
     *         -5.0907615D-15*TN(5) 
      ELSE
      SMH(3) = -4.64130376D0 +1.02466476D4*TI 
     *         +5.14987613D0*TN(1) -6.8354894D-3*TN(2) 
     *         +8.19667665D-6*TN(3) -4.039525216666667D-9*TN(4) 
     *         +8.3346978D-13*TN(5) 
      ENDIF
C C3H6
      IF (T .GT. 1000) THEN
      SMH(4) = -1.331335D1 +9.235703D2*TI 
     *         +6.732257D0*TN(1) +7.45417D-3*TN(2) 
     *         -8.249831666666666D-7*TN(3) +6.010018333333334D-11*TN(4) 
     *         -1.883102D-15*TN(5) 
      ELSE
      SMH(4) = 1.614534D1 -1.074826D3*TI 
     *         +1.493307D0*TN(1) +1.046259D-2*TN(2) 
     *         +7.47799D-7*TN(3) -1.39076D-9*TN(4) 
     *         +3.579073D-13*TN(5) 
      ENDIF
C iC4H8
      IF (T .GT. 1000) THEN
      SMH(5) = 1.0671549D0 +5.0066758D3*TI 
     *         +4.460947D0*TN(1) +1.48057435D-2*TN(2) 
     *         -2.1795215D-6*TN(3) +2.214327833333333D-10*TN(4) 
     *         -1.00673565D-14*TN(5) 
      ELSE
      SMH(5) = 1.2676388D1 +4.0373069D3*TI 
     *         +2.6471405D0*TN(1) +1.29514785D-2*TN(2) 
     *         +1.366422566666667D-6*TN(3) -1.84943825D-9*TN(4) 
     *         +4.447929D-13*TN(5) 
      ENDIF
C C4H81
      IF (T .GT. 1000) THEN
      SMH(6) = 1.5543201D1 +2.1397231D3*TI 
     *         +2.0535841D0*TN(1) +1.71752535D-2*TN(2) 
     *         -2.6471995D-6*TN(3) +2.757471833333334D-10*TN(4) 
     *         -1.26805225D-14*TN(5) 
      ELSE
      SMH(6) = 2.1062469D1 +1.7904004D3*TI 
     *         +1.181138D0*TN(1) +1.542669D-2*TN(2) 
     *         +8.477541166666667D-7*TN(3) -2.054574D-9*TN(4) 
     *         +5.555096499999999D-13*TN(5) 
      ENDIF
C H2
      IF (T .GT. 1000) THEN
      SMH(7) = -3.20502331D0 +9.50158922D2*TI 
     *         +3.3372792D0*TN(1) -2.470123655D-5*TN(2) 
     *         +8.324279633333333D-8*TN(3) -1.496386616666667D-11*TN(4) 
     *         +1.00127688D-15*TN(5) 
      ELSE
      SMH(7) = 6.83010238D-1 +9.17935173D2*TI 
     *         +2.34433112D0*TN(1) +3.990260375D-3*TN(2) 
     *         -3.2463585D-6*TN(3) +1.67976745D-9*TN(4) 
     *         -3.688058805D-13*TN(5) 
      ENDIF
C C2H6
      IF (T .GT. 1000) THEN
      SMH(8) = 1.51156107D1 +1.14263932D4*TI 
     *         +1.0718815D0*TN(1) +1.084263385D-2*TN(2) 
     *         -1.67093445D-6*TN(3) +1.845100008333333D-10*TN(4) 
     *         -9.5001445D-15*TN(5) 
      ELSE
      SMH(8) = 2.66682316D0 +1.15222055D4*TI 
     *         +4.29142492D0*TN(1) -2.75077135D-3*TN(2) 
     *         +9.990638133333334D-6*TN(3) -5.903885708333334D-9*TN(4) 
     *         +1.343428855D-12*TN(5) 
      ENDIF
C CO
      IF (T .GT. 1000) THEN
      SMH(9) = 7.81868772D0 +1.41518724D4*TI 
     *         +2.71518561D0*TN(1) +1.031263715D-3*TN(2) 
     *         -1.664709618333334D-7*TN(3) +1.9171084D-11*TN(4) 
     *         -1.01823858D-15*TN(5) 
      ELSE
      SMH(9) = 3.50840928D0 +1.4344086D4*TI 
     *         +3.57953347D0*TN(1) -3.0517684D-4*TN(2) 
     *         +1.69469055D-7*TN(3) +7.558382366666667D-11*TN(4) 
     *         -4.522122495D-14*TN(5) 
      ENDIF
C C6H6
      IF (T .GT. 1000) THEN
      SMH(10) = -2.9115665D1 -5.2043462D3*TI 
     *         +9.1381245D0*TN(1) +1.19272165D-2*TN(2) 
     *         -1.468795433333333D-6*TN(3) +1.00825175D-10*TN(4) 
     *         -9.110751500000001D-16*TN(5) 
      ELSE
      SMH(10) = 4.3889832D1 -9.1817773D3*TI 
     *         -4.8437734D0*TN(1) +2.92138065D-2*TN(2) 
     *         -4.914309166666667D-6*TN(3) -5.782536666666667D-10*TN(4) 
     *         +4.10626265D-13*TN(5) 
      ENDIF
C C2H2
      IF (T .GT. 1000) THEN
      SMH(11) = -1.23028121D0 -2.59359992D4*TI 
     *         +4.14756964D0*TN(1) +2.98083332D-3*TN(2) 
     *         -3.9549142D-7*TN(3) +3.895101425D-11*TN(4) 
     *         -1.806176065D-15*TN(5) 
      ELSE
      SMH(11) = 1.39397051D1 -2.64289807D4*TI 
     *         +8.08681094D-1*TN(1) +1.168078145D-2*TN(2) 
     *         -5.91953025D-6*TN(3) +2.334603641666667D-9*TN(4) 
     *         -4.25036487D-13*TN(5) 
      ENDIF
C C6H5CH3
      IF (T .GT. 1000) THEN
      SMH(12) = -4.6728785D1 +6.9764908D2*TI 
     *         +1.2940034D1*TN(1) +1.33456435D-2*TN(2) 
     *         -1.613975083333333D-6*TN(3) +1.311552416666667D-10*TN(4) 
     *         -4.73318005D-15*TN(5) 
      ELSE
      SMH(12) = 2.028221D1 -4.07563D3*TI 
     *         +1.6152663D0*TN(1) +1.0549719D-2*TN(2) 
     *         +1.422766966666667D-5*TN(3) -1.105088833333333D-8*TN(4) 
     *         +2.7978302D-12*TN(5) 
      ENDIF
C CH3
      IF (T .GT. 1000) THEN
      SMH(13) = 8.48007179D0 -1.67755843D4*TI 
     *         +2.28571772D0*TN(1) +3.619950185D-3*TN(2) 
     *         -4.978572466666667D-7*TN(3) +4.9640387D-11*TN(4) 
     *         -2.33577197D-15*TN(5) 
      ELSE
      SMH(13) = 1.60456433D0 -1.64449988D4*TI 
     *         +3.6735904D0*TN(1) +1.005475875D-3*TN(2) 
     *         +9.550364266666668D-7*TN(3) -5.725978541666666D-10*TN(4) 
     *         +1.27192867D-13*TN(5) 
      ENDIF
C O2
      IF (T .GT. 1000) THEN
      SMH(14) = 5.45323129D0 +1.08845772D3*TI 
     *         +3.28253784D0*TN(1) +7.4154377D-4*TN(2) 
     *         -1.263277781666667D-7*TN(3) +1.745587958333333D-11*TN(4) 
     *         -1.08358897D-15*TN(5) 
      ELSE
      SMH(14) = 3.65767573D0 +1.06394356D3*TI 
     *         +3.78245636D0*TN(1) -1.49836708D-3*TN(2) 
     *         +1.641217001666667D-6*TN(3) -8.067745908333334D-10*TN(4) 
     *         +1.621864185D-13*TN(5) 
      ENDIF
C O
      IF (T .GT. 1000) THEN
      SMH(15) = 4.78433864D0 -2.92175791D4*TI 
     *         +2.56942078D0*TN(1) -4.298705685D-5*TN(2) 
     *         +6.991409816666667D-9*TN(3) -8.348149916666666D-13*TN(4) 
     *         +6.141684549999999D-17*TN(5) 
      ELSE
      SMH(15) = 2.05193346D0 -2.91222592D4*TI 
     *         +3.1682671D0*TN(1) -1.63965942D-3*TN(2) 
     *         +1.107177326666667D-6*TN(3) -5.106721866666666D-10*TN(4) 
     *         +1.056329855D-13*TN(5) 
      ENDIF
C OH
      IF (T .GT. 1000) THEN
      SMH(16) = 5.70164073D0 -3.71885774D3*TI 
     *         +2.86472886D0*TN(1) +5.2825224D-4*TN(2) 
     *         -4.318045966666667D-8*TN(3) +2.54348895D-12*TN(4) 
     *         -6.6597938D-17*TN(5) 
      ELSE
      SMH(16) = -6.9043296D-1 -3.38153812D3*TI 
     *         +4.12530561D0*TN(1) -1.612724695D-3*TN(2) 
     *         +1.087941151666667D-6*TN(3) -4.832113691666666D-10*TN(4) 
     *         +1.031186895D-13*TN(5) 
      ENDIF
C HO2
      IF (T .GT. 1000) THEN
      SMH(17) = 3.78510215D0 -1.11856713D2*TI 
     *         +4.0172109D0*TN(1) +1.119910065D-3*TN(2) 
     *         -1.056096916666667D-7*TN(3) +9.520530833333334D-12*TN(4) 
     *         -5.39542675D-16*TN(5) 
      ELSE
      SMH(17) = 3.71666245D0 -2.9480804D2*TI 
     *         +4.30179801D0*TN(1) -2.374560255D-3*TN(2) 
     *         +3.526381516666666D-6*TN(3) -2.02303245D-9*TN(4) 
     *         +4.646125620000001D-13*TN(5) 
      ENDIF
C H2O
      IF (T .GT. 1000) THEN
      SMH(18) = 4.9667701D0 +3.00042971D4*TI 
     *         +3.03399249D0*TN(1) +1.08845902D-3*TN(2) 
     *         -2.734541966666666D-8*TN(3) -8.08683225D-12*TN(4) 
     *         +8.4100496D-16*TN(5) 
      ELSE
      SMH(18) = -8.49032208D-1 +3.02937267D4*TI 
     *         +4.19864056D0*TN(1) -1.01821705D-3*TN(2) 
     *         +1.086733685D-6*TN(3) -4.57330885D-10*TN(4) 
     *         +8.85989085D-14*TN(5) 
      ENDIF
C H2O2
      IF (T .GT. 1000) THEN
      SMH(19) = 2.91615662D0 +1.78617877D4*TI 
     *         +4.16500285D0*TN(1) +2.45415847D-3*TN(2) 
     *         -3.168987083333333D-7*TN(3) +3.09321655D-11*TN(4) 
     *         -1.439541525D-15*TN(5) 
      ELSE
      SMH(19) = 3.43505074D0 +1.77025821D4*TI 
     *         +4.27611269D0*TN(1) -2.714112085D-4*TN(2) 
     *         +2.78892835D-6*TN(3) -1.798090108333333D-9*TN(4) 
     *         +4.312271815D-13*TN(5) 
      ENDIF
C H
      IF (T .GT. 1000) THEN
      SMH(20) = -4.46682914D-1 -2.54736599D4*TI 
     *         +2.50000001D0*TN(1) -1.154214865D-11*TN(2) 
     *         +2.692699133333334D-15*TN(3) -3.945960291666667D-19*TN(4) 
     *         +2.490986785D-23*TN(5) 
      ELSE
      SMH(20) = -4.46682853D-1 -2.54736599D4*TI 
     *         +2.5D0*TN(1) +3.526664095D-13*TN(2) 
     *         -3.326532733333333D-16*TN(3) +1.917346933333333D-19*TN(4) 
     *         -4.63866166D-23*TN(5) 
      ENDIF
C CH2
      IF (T .GT. 1000) THEN
      SMH(21) = 6.17119324D0 -4.6263604D4*TI 
     *         +2.87410113D0*TN(1) +1.82819646D-3*TN(2) 
     *         -2.348243283333333D-7*TN(3) +2.168162908333333D-11*TN(4) 
     *         -9.38637835D-16*TN(5) 
      ELSE
      SMH(21) = 1.56253185D0 -4.60040401D4*TI 
     *         +3.76267867D0*TN(1) +4.844360715D-4*TN(2) 
     *         +4.658164016666667D-7*TN(3) -3.209092941666667D-10*TN(4) 
     *         +8.43708595D-14*TN(5) 
      ENDIF
C CH2*
      IF (T .GT. 1000) THEN
      SMH(22) = 8.62650169D0 -5.09259997D4*TI 
     *         +2.29203842D0*TN(1) +2.327943185D-3*TN(2) 
     *         -3.353199116666667D-7*TN(3) +3.48255D-11*TN(4) 
     *         -1.698581825D-15*TN(5) 
      ELSE
      SMH(22) = -7.69118967D-1 -5.04968163D4*TI 
     *         +4.19860411D0*TN(1) -1.183307095D-3*TN(2) 
     *         +1.372160366666667D-6*TN(3) -5.573466508333334D-10*TN(4) 
     *         +9.71573685D-14*TN(5) 
      ENDIF
C HCO
      IF (T .GT. 1000) THEN
      SMH(23) = 9.79834492D0 -4.01191815D3*TI 
     *         +2.77217438D0*TN(1) +2.47847763D-3*TN(2) 
     *         -4.140760216666667D-7*TN(3) +4.909681483333334D-11*TN(4) 
     *         -2.667543555D-15*TN(5) 
      ELSE
      SMH(23) = 3.39437243D0 -3.83956496D3*TI 
     *         +4.22118584D0*TN(1) -1.62196266D-3*TN(2) 
     *         +2.296657433333333D-6*TN(3) -1.109534108333333D-9*TN(4) 
     *         +2.168844325D-13*TN(5) 
      ENDIF
C CH2O
      IF (T .GT. 1000) THEN
      SMH(24) = 1.3656323D1 +1.39958323D4*TI 
     *         +1.76069008D0*TN(1) +4.60000041D-3*TN(2) 
     *         -7.370980216666666D-7*TN(3) +8.386767666666666D-11*TN(4) 
     *         -4.419278200000001D-15*TN(5) 
      ELSE
      SMH(24) = 6.028129D-1 +1.43089567D4*TI 
     *         +4.79372315D0*TN(1) -4.954166845D-3*TN(2) 
     *         +6.220333466666666D-6*TN(3) -3.160710508333333D-9*TN(4) 
     *         +6.5886326D-13*TN(5) 
      ENDIF
C CH3O
      IF (T .GT. 1000) THEN
      SMH(25) = -1.96680028D0 -3.7811194D2*TI 
     *         +4.75779238D0*TN(1) +3.72071237D-3*TN(2) 
     *         -4.495086266666666D-7*TN(3) +3.6507542D-11*TN(4) 
     *         -1.31768549D-15*TN(5) 
      ELSE
      SMH(25) = 6.57240864D0 -1.2956976D3*TI 
     *         +3.71180502D0*TN(1) -1.40231653D-3*TN(2) 
     *         +6.275849516666667D-6*TN(3) -3.942267408333333D-9*TN(4) 
     *         +9.329421000000001D-13*TN(5) 
      ENDIF
C CO2
      IF (T .GT. 1000) THEN
      SMH(26) = 2.27163806D0 +4.8759166D4*TI 
     *         +3.85746029D0*TN(1) +2.20718513D-3*TN(2) 
     *         -3.691356733333334D-7*TN(3) +4.362418233333334D-11*TN(4) 
     *         -2.36042082D-15*TN(5) 
      ELSE
      SMH(26) = 9.90105222D0 +4.83719697D4*TI 
     *         +2.35677352D0*TN(1) +4.492298385D-3*TN(2) 
     *         -1.187260448333333D-6*TN(3) +2.049325183333333D-10*TN(4) 
     *         -7.184977399999999D-15*TN(5) 
      ENDIF
C C2H3
      IF (T .GT. 1000) THEN
      SMH(27) = 7.78732378D0 -3.46128739D4*TI 
     *         +3.016724D0*TN(1) +5.1651146D-3*TN(2) 
     *         -7.801372483333333D-7*TN(3) +8.480274D-11*TN(4) 
     *         -4.313035205D-15*TN(5) 
      ELSE
      SMH(27) = 8.51054025D0 -3.48598468D4*TI 
     *         +3.21246645D0*TN(1) +7.5739581D-4*TN(2) 
     *         +4.320156866666666D-6*TN(3) -2.980482058333333D-9*TN(4) 
     *         +7.35754365D-13*TN(5) 
      ENDIF
C C2H5
      IF (T .GT. 1000) THEN
      SMH(28) = 1.34624343D1 -1.285752D4*TI 
     *         +1.95465642D0*TN(1) +8.698636100000001D-3*TN(2) 
     *         -1.330344446666667D-6*TN(3) +1.460147408333333D-10*TN(4) 
     *         -7.4820788D-15*TN(5) 
      ELSE
      SMH(28) = 4.70720924D0 -1.28416265D4*TI 
     *         +4.30646568D0*TN(1) -2.09329446D-3*TN(2) 
     *         +8.28571345D-6*TN(3) -4.992721716666666D-9*TN(4) 
     *         +1.15254502D-12*TN(5) 
      ENDIF
C HCCO
      IF (T .GT. 1000) THEN
      SMH(29) = -3.9302595D0 -1.9327215D4*TI 
     *         +5.6282058D0*TN(1) +2.04267005D-3*TN(2) 
     *         -2.655757833333333D-7*TN(3) +2.385504333333333D-11*TN(4) 
     *         -9.703915999999999D-16*TN(5) 
      ELSE
      SMH(29) = 1.2490417D1 -2.0059449D4*TI 
     *         +2.2517214D0*TN(1) +8.8275105D-3*TN(2) 
     *         -3.954850166666666D-6*TN(3) +1.439646583333334D-9*TN(4) 
     *         -2.53324055D-13*TN(5) 
      ENDIF
C CH2CO
      IF (T .GT. 1000) THEN
      SMH(30) = 6.32247205D-1 +7.7785D3*TI 
     *         +4.51129732D0*TN(1) +4.501798725D-3*TN(2) 
     *         -6.948993916666666D-7*TN(3) +7.694549016666667D-11*TN(4) 
     *         -3.974191005D-15*TN(5) 
      ELSE
      SMH(30) = 1.2215648D1 +7.27D3*TI 
     *         +2.1358363D0*TN(1) +9.05943605D-3*TN(2) 
     *         -2.899124566666666D-6*TN(3) +7.7866464D-10*TN(4) 
     *         -1.007288075D-13*TN(5) 
      ENDIF
C CH2CHO
      IF (T .GT. 1000) THEN
      SMH(31) = -5.0320879D0 +9.695D2*TI 
     *         +5.9756699D0*TN(1) +4.0652957D-3*TN(2) 
     *         -4.5727075D-7*TN(3) +3.391920083333333D-11*TN(4) 
     *         -1.08800855D-15*TN(5) 
      ELSE
      SMH(31) = 9.571453500000001D0 -6.2D1*TI 
     *         +3.4090624D0*TN(1) +5.369287D-3*TN(2) 
     *         +3.1524875D-7*TN(3) -5.965485916666667D-10*TN(4) 
     *         +1.43369255D-13*TN(5) 
      ENDIF
C C3H3
      IF (T .GT. 1000) THEN
      SMH(32) = -1.25848436D1 -3.89087427D4*TI 
     *         +7.1422188D0*TN(1) +3.809510025D-3*TN(2) 
     *         -4.457665833333333D-7*TN(3) +3.540956675D-11*TN(4) 
     *         -1.257377075D-15*TN(5) 
      ELSE
      SMH(32) = 1.52058924D1 -4.01057783D4*TI 
     *         +1.35110927D0*TN(1) +1.637056115D-2*TN(2) 
     *         -7.897118916666667D-6*TN(3) +3.135915066666667D-9*TN(4) 
     *         -5.92704615D-13*TN(5) 
      ENDIF
C aC3H5
      IF (T .GT. 1000) THEN
      SMH(33) = -1.124305D1 -1.7482449D4*TI 
     *         +6.5007877D0*TN(1) +7.1623655D-3*TN(2) 
     *         -9.463605333333332D-7*TN(3) +9.234000833333333D-11*TN(4) 
     *         -4.518194349999999D-15*TN(5) 
      ELSE
      SMH(33) = 1.7173214D1 -1.9245629D4*TI 
     *         +1.3631835D0*TN(1) +9.906910499999999D-3*TN(2) 
     *         +2.082843333333333D-6*TN(3) -2.779629583333333D-9*TN(4) 
     *         +7.9232855D-13*TN(5) 
      ENDIF
C C5H4O
      IF (T .GT. 1000) THEN
      SMH(34) = -2.94521623D1 -1.94364771D3*TI 
     *         +1.00806824D1*TN(1) +8.057173250000001D-3*TN(2) 
     *         -9.721908483333333D-7*TN(3) +7.889661D-11*TN(4) 
     *         -2.84486103D-15*TN(5) 
      ELSE
      SMH(34) = 2.35409513D1 -5.11159287D3*TI 
     *         +2.64576497D-1*TN(1) +1.674369135D-2*TN(2) 
     *         +2.795641166666666D-7*TN(3) -2.468395458333334D-9*TN(4) 
     *         +7.721573800000001D-13*TN(5) 
      ENDIF
C C5H5
      IF (T .GT. 1000) THEN
      SMH(35) = -1.6133D1 -2.8086D4*TI 
     *         +7.4743938D0*TN(1) +8.006366500000001D-3*TN(2) 
     *         -1.0803858D-9*TN(3) -2.98497525D-10*TN(4) 
     *         +4.618253550000001D-14*TN(5) 
      ELSE
      SMH(35) = 1.6585519D1 -2.9626D4*TI 
     *         +9.8349822D-1*TN(1) +1.6825738D-2*TN(2) 
     *         -1.8423635D-8*TN(3) -3.0619495D-9*TN(4) 
     *         +1.1570592D-12*TN(5) 
      ENDIF
C C6H5
      IF (T .GT. 1000) THEN
      SMH(36) = -2.2954643D1 -3.6261047D4*TI 
     *         +8.597310999999999D0*TN(1) +1.1120815D-2*TN(2) 
     *         -1.453332966666666D-6*TN(3) +1.149065416666667D-10*TN(4) 
     *         -2.6573028D-15*TN(5) 
      ELSE
      SMH(36) = 4.1332535D1 -3.977959D4*TI 
     *         -3.6931453D0*TN(1) +2.6089484D-2*TN(2) 
     *         -4.259737833333333D-6*TN(3) -5.88842675D-10*TN(4) 
     *         +3.79169875D-13*TN(5) 
      ENDIF
C C6H5O
      IF (T .GT. 1000) THEN
      SMH(37) = -4.8818168D1 -2.87274751D2*TI 
     *         +1.3722172D1*TN(1) +8.734438549999999D-3*TN(2) 
     *         -1.0591742D-6*TN(3) +8.624359000000001D-11*TN(4) 
     *         -3.11705252D-15*TN(5) 
      ELSE
      SMH(37) = 2.76990274D1 -4.77858391D3*TI 
     *         -4.66204455D-1*TN(1) +2.067219875D-2*TN(2) 
     *         +2.206883183333333D-6*TN(3) -4.773939741666667D-9*TN(4) 
     *         +1.448818535D-12*TN(5) 
      ENDIF
C
      END
C                                                                      C
C----------------------------------------------------------------------C
C                                                                      C
      SUBROUTINE RATX (T, C, RF, RB, RKLOW)
      IMPLICIT DOUBLE PRECISION (A-H, O-Z), INTEGER (I-N)
      PARAMETER (SMALL = 1D-200)
      DIMENSION C(*), RF(*), RB(*), RKLOW(*)
C
      ALOGT = LOG(T)
      CTOT = 0.0
      DO K = 1, 29
         CTOT = CTOT + C(K)
      ENDDO
C
C     R1: POSF5433 = 3.09172C2H4 + 0.803846C3H6 + 0.154586iC4H8 + 0.309172C4H81 + 0.116667C6H6 + 0.05C6H5CH3 + 1.5CH3 + 0.5H
      RF(1) = RF(1)*C(1)
C     R2: POSF5433 + 0.7H = 3.09911C2H4 + 0.25CH4 + 0.805769C3H6 + 0.154956iC4H8 + 0.309911C4H81 + H2 + 0.175C6H6 + 0.075C6H5CH3 + 0.7CH3
      RF(2) = RF(2)*C(1)*C(20)
C     R3: POSF5433 + 0.3CH3 = 3.09911C2H4 + 1.25CH4 + 0.805769C3H6 + 0.154956iC4H8 + 0.309911C4H81 + 0.175C6H6 + 0.075C6H5CH3 + 0.3H
      RF(3) = RF(3)*C(1)*C(13)
C     R4: POSF5433 + OH = 3.09911C2H4 + 0.25CH4 + 0.805769C3H6 + 0.154956iC4H8 + 0.309911C4H81 + 0.175C6H6 + 0.075C6H5CH3 + 0.7CH3 + H2O + 0.3H
      RF(4) = RF(4)*C(1)*C(16)
C     R5: POSF5433 + O2 = 3.09911C2H4 + 0.25CH4 + 0.805769C3H6 + 0.154956iC4H8 + 0.309911C4H81 + 0.175C6H6 + 0.075C6H5CH3 + 0.7CH3 + HO2 + 0.3H
      RF(5) = RF(5)*C(1)*C(14)
C     R6: POSF5433 + HO2 = 3.09911C2H4 + 0.25CH4 + 0.805769C3H6 + 0.154956iC4H8 + 0.309911C4H81 + 0.175C6H6 + 0.075C6H5CH3 + 0.7CH3 + H2O2 + 0.3H
      RF(6) = RF(6)*C(1)*C(17)
C     R7: POSF5433 + O = 3.09911C2H4 + 0.25CH4 + 0.805769C3H6 + 0.154956iC4H8 + 0.309911C4H81 + 0.175C6H6 + 0.075C6H5CH3 + 0.7CH3 + OH + 0.3H
      RF(7) = RF(7)*C(1)*C(15)
C     R8: O2 + H = O + OH
      RF(8) = RF(8)*C(14)*C(20)
      RB(8) = RB(8)*C(15)*C(16)
C     R9: H2 + O = OH + H
      RF(9) = RF(9)*C(7)*C(15)
      RB(9) = RB(9)*C(16)*C(20)
C     R10: H2 + OH = H2O + H
      RF(10) = RF(10)*C(7)*C(16)
      RB(10) = RB(10)*C(18)*C(20)
C     R11: 2OH = O + H2O
      RF(11) = RF(11)*C(16)*C(16)
      RB(11) = RB(11)*C(15)*C(18)
C     R12: 2H = H2
      CTB = CTOT-C(7)-C(18)
     * -C(22)
      RF(12) = RF(12)*CTB*C(20)*C(20)
      RB(12) = RB(12)*CTB*C(7)
C     R13: 2H = H2
      RF(13) = RF(13)*C(7)*C(20)*C(20)
      RB(13) = RB(13)*C(7)*C(7)
C     R14: 2H = H2
      RF(14) = RF(14)*C(18)*C(20)*C(20)
      RB(14) = RB(14)*C(7)*C(18)
C     R15: OH + H = H2O
      CTB = CTOT+C(7)+5.3D0*C(18)
     * +7.5D-1*C(9)+2.6D0*C(22)
      RF(15) = RF(15)*CTB*C(16)*C(20)
      RB(15) = RB(15)*CTB*C(18)
C     R16: O + H = OH
      CTB = CTOT+C(7)+1.1D1*C(18)
     * +7.5D-1*C(9)+2.6D0*C(22)
      RF(16) = RF(16)*CTB*C(15)*C(20)
      RB(16) = RB(16)*CTB*C(16)
C     R17: 2O = O2
      CTB = CTOT+1.4D0*C(7)+1.44D1*C(18)
     * +7.5D-1*C(9)+2.6D0*C(22)
      RF(17) = RF(17)*CTB*C(15)*C(15)
      RB(17) = RB(17)*CTB*C(14)
C     R18: O2 + H = HO2
      CTB = CTOT-1.5D-1*C(14)+1.089D1*C(18)
     * +9.000000000000008D-2*C(9)+1.18D0*C(22)
      PR = RKLOW(1) * CTB / RF(18)
      PCOR = PR / (1.0 + PR)
      PRLOG = LOG10(MAX(PR,SMALL))
      FC = (PRLOG -1.983099029051326D-1)
     *     /(1.132308094493256D0 -0.14D0*(PRLOG -1.983099029051326D-1))
      FC = EXP(-6.931471805599453D-1 /(1.0D0 + FC*FC))
      PCOR = FC * PCOR
      RF(18) = RF(18) * PCOR
      RB(18) = RB(18) * PCOR
      RF(18) = RF(18)*C(14)*C(20)
      RB(18) = RB(18)*C(17)
C     R19: H2 + O2 = HO2 + H
      RF(19) = RF(19)*C(7)*C(14)
      RB(19) = RB(19)*C(17)*C(20)
C     R20: 2OH = H2O2
      CTB = CTOT+C(7)+5.D0*C(18)
     * +7.5D-1*C(9)+2.6D0*C(22)
      PR = RKLOW(2) * CTB / RF(20)
      PCOR = PR / (1.0 + PR)
      PRLOG = LOG10(MAX(PR,SMALL))
      FCENT = 2.654D-1*EXP(-T/9.4D1)
     *      + 7.346D-1*EXP(-T/1.756D3)
     *     + EXP(-5.182D3/T)
      FCLOG = LOG10(MAX(FCENT,SMALL))
      XN    = 0.75 - 1.27*FCLOG
      CPRLOG= PRLOG - (0.4 + 0.67*FCLOG)
      FLOG = FCLOG/(1.0 + (CPRLOG/(XN-0.14*CPRLOG))**2)
      FC = 10.0**FLOG
      PCOR = FC * PCOR
      RF(20) = RF(20) * PCOR
      RB(20) = RB(20) * PCOR
      RF(20) = RF(20)*C(16)*C(16)
      RB(20) = RB(20)*C(19)
C     R21: HO2 + H = 2OH
      RF(21) = RF(21)*C(17)*C(20)
      RB(21) = RB(21)*C(16)*C(16)
C     R22: O + HO2 = O2 + OH
      RF(22) = RF(22)*C(15)*C(17)
      RB(22) = RB(22)*C(14)*C(16)
C     R23: 2HO2 = O2 + H2O2
      RF(23) = RF(23)*C(17)*C(17)
      RB(23) = RB(23)*C(14)*C(19)
C     R24: 2HO2 = O2 + H2O2
      RF(24) = RF(24)*C(17)*C(17)
      RB(24) = RB(24)*C(14)*C(19)
C     R25: OH + HO2 = O2 + H2O
      RF(25) = RF(25)*C(16)*C(17)
      RB(25) = RB(25)*C(14)*C(18)
C     R26: OH + HO2 = O2 + H2O
      RF(26) = RF(26)*C(16)*C(17)
      RB(26) = RB(26)*C(14)*C(18)
C     R27: OH + HO2 = O2 + H2O
      RF(27) = RF(27)*C(16)*C(17)
      RB(27) = RB(27)*C(14)*C(18)
C     R28: H2O2 + H = H2 + HO2
      RF(28) = RF(28)*C(19)*C(20)
      RB(28) = RB(28)*C(7)*C(17)
C     R29: H2O2 + H = OH + H2O
      RF(29) = RF(29)*C(19)*C(20)
      RB(29) = RB(29)*C(16)*C(18)
C     R30: O + H2O2 = OH + HO2
      RF(30) = RF(30)*C(15)*C(19)
      RB(30) = RB(30)*C(16)*C(17)
C     R31: OH + H2O2 = HO2 + H2O
      RF(31) = RF(31)*C(16)*C(19)
      RB(31) = RB(31)*C(17)*C(18)
C     R32: OH + H2O2 = HO2 + H2O
      RF(32) = RF(32)*C(16)*C(19)
      RB(32) = RB(32)*C(17)*C(18)
C     R33: CO + O = CO2
      CTB = CTOT+C(7)+1.1D1*C(18)
     * +7.5D-1*C(9)+2.6D0*C(22)
      PR = RKLOW(3) * CTB / RF(33)
      PCOR = PR / (1.0 + PR)
      RF(33) = RF(33) * PCOR
      RB(33) = RB(33) * PCOR
      RF(33) = RF(33)*C(9)*C(15)
      RB(33) = RB(33)*C(22)
C     R34: CO + OH = H + CO2
      RF(34) = RF(34)*C(9)*C(16)
      RB(34) = RB(34)*C(20)*C(22)
C     R35: CO + OH = H + CO2
      RF(35) = RF(35)*C(9)*C(16)
      RB(35) = RB(35)*C(20)*C(22)
C     R36: CO + O2 = O + CO2
      RF(36) = RF(36)*C(9)*C(14)
      RB(36) = RB(36)*C(15)*C(22)
C     R37: CO + HO2 = OH + CO2
      RF(37) = RF(37)*C(9)*C(17)
      RB(37) = RB(37)*C(16)*C(22)
C     R38: H + HCO = H2 + CO
      RF(38) = RF(38)*C(20)
      RB(38) = RB(38)*C(7)*C(9)
C     R39: O + HCO = CO + OH
      RF(39) = RF(39)*C(15)
      RB(39) = RB(39)*C(9)*C(16)
C     R40: O + HCO = H + CO2
      RF(40) = RF(40)*C(15)
      RB(40) = RB(40)*C(20)*C(22)
C     R41: OH + HCO = CO + H2O
      RF(41) = RF(41)*C(16)
      RB(41) = RB(41)*C(9)*C(18)
C     R42: HCO = CO + H
      CTB = CTOT+C(7)-C(18)
     * +7.5D-1*C(9)+2.6D0*C(22)
      RF(42) = RF(42)*CTB
      RB(42) = RB(42)*CTB*C(9)*C(20)
C     R43: HCO = CO + H
      RF(43) = RF(43)*C(18)
      RB(43) = RB(43)*C(9)*C(18)*C(20)
C     R44: O2 + HCO = CO + HO2
      RF(44) = RF(44)*C(14)
      RB(44) = RB(44)*C(9)*C(17)
C     R45: H2 + CO = CH2O
      CTB = CTOT+C(7)+5.D0*C(18)
     * +C(3)+5.D-1*C(9)
     * +C(22)+2.D0*C(8)
      PR = RKLOW(4) * CTB / RF(45)
      PCOR = PR / (1.0 + PR)
      PRLOG = LOG10(MAX(PR,SMALL))
      FCENT = 6.799999999999995D-2*EXP(-T/1.97D2)
     *      + 9.320000000000001D-1*EXP(-T/1.54D3)
     *     + EXP(-1.03D4/T)
      FCLOG = LOG10(MAX(FCENT,SMALL))
      XN    = 0.75 - 1.27*FCLOG
      CPRLOG= PRLOG - (0.4 + 0.67*FCLOG)
      FLOG = FCLOG/(1.0 + (CPRLOG/(XN-0.14*CPRLOG))**2)
      FC = 10.0**FLOG
      PCOR = FC * PCOR
      RF(45) = RF(45) * PCOR
      RB(45) = RB(45) * PCOR
      RF(45) = RF(45)*C(7)*C(9)
      RB(45) = RB(45)*C(21)
C     R46: H + HCO = CH2O
      CTB = CTOT+C(7)+5.D0*C(18)
     * +C(3)+5.D-1*C(9)
     * +C(22)+2.D0*C(8)
      PR = RKLOW(5) * CTB / RF(46)
      PCOR = PR / (1.0 + PR)
      PRLOG = LOG10(MAX(PR,SMALL))
      FCENT = 2.176D-1*EXP(-T/2.71D2)
     *      + 7.824D-1*EXP(-T/2.755D3)
     *     + EXP(-6.57D3/T)
      FCLOG = LOG10(MAX(FCENT,SMALL))
      XN    = 0.75 - 1.27*FCLOG
      CPRLOG= PRLOG - (0.4 + 0.67*FCLOG)
      FLOG = FCLOG/(1.0 + (CPRLOG/(XN-0.14*CPRLOG))**2)
      FC = 10.0**FLOG
      PCOR = FC * PCOR
      RF(46) = RF(46) * PCOR
      RB(46) = RB(46) * PCOR
      RF(46) = RF(46)*C(20)
      RB(46) = RB(46)*C(21)
C     R47: H + CH2 = CH3
      CTB = CTOT+C(7)+5.D0*C(18)
     * +C(3)+5.D-1*C(9)
     * +C(22)+2.D0*C(8)
      PR = RKLOW(6) * CTB / RF(47)
      PCOR = PR / (1.0 + PR)
      PRLOG = LOG10(MAX(PR,SMALL))
      FCENT = 3.2D-1*EXP(-T/7.8D1)
     *      + 6.8D-1*EXP(-T/1.995D3)
     *     + EXP(-5.59D3/T)
      FCLOG = LOG10(MAX(FCENT,SMALL))
      XN    = 0.75 - 1.27*FCLOG
      CPRLOG= PRLOG - (0.4 + 0.67*FCLOG)
      FLOG = FCLOG/(1.0 + (CPRLOG/(XN-0.14*CPRLOG))**2)
      FC = 10.0**FLOG
      PCOR = FC * PCOR
      RF(47) = RF(47) * PCOR
      RB(47) = RB(47) * PCOR
      RF(47) = RF(47)*C(20)
      RB(47) = RB(47)*C(13)
C     R48: O + CH2 = H + HCO
      RF(48) = RF(48)*C(15)
      RB(48) = RB(48)*C(20)
C     R49: OH + CH2 = H + CH2O
      RF(49) = RF(49)*C(16)
      RB(49) = RB(49)*C(20)*C(21)
C     R50: H2 + CH2 = CH3 + H
      RF(50) = RF(50)*C(7)
      RB(50) = RB(50)*C(13)*C(20)
C     R51: O2 + CH2 = OH + HCO
      RF(51) = RF(51)*C(14)
      RB(51) = RB(51)*C(16)
C     R52: O2 + CH2 = 2H + CO2
      RF(52) = RF(52)*C(14)
      RB(52) = RB(52)*C(20)*C(20)*C(22)
C     R53: CO + CH2 = CH2CO
      CTB = CTOT+C(7)+5.D0*C(18)
     * +C(3)+5.D-1*C(9)
     * +C(22)+2.D0*C(8)
      PR = RKLOW(7) * CTB / RF(53)
      PCOR = PR / (1.0 + PR)
      PRLOG = LOG10(MAX(PR,SMALL))
      FCENT = 4.093D-1*EXP(-T/2.75D2)
     *      + 5.907D-1*EXP(-T/1.226D3)
     *     + EXP(-5.185D3/T)
      FCLOG = LOG10(MAX(FCENT,SMALL))
      XN    = 0.75 - 1.27*FCLOG
      CPRLOG= PRLOG - (0.4 + 0.67*FCLOG)
      FLOG = FCLOG/(1.0 + (CPRLOG/(XN-0.14*CPRLOG))**2)
      FC = 10.0**FLOG
      PCOR = FC * PCOR
      RF(53) = RF(53) * PCOR
      RB(53) = RB(53) * PCOR
      RF(53) = RF(53)*C(9)
      RB(53) = RB(53)*C(23)
C     R54: 2CH2 = H2 + C2H2
      RB(54) = RB(54)*C(7)*C(11)
C     R55: CH2* = CH2
      RF(55) = RF(55)*C(29)
      RB(55) = RB(55)*C(29)
C     R56: O + CH2* = H2 + CO
      RF(56) = RF(56)*C(15)
      RB(56) = RB(56)*C(7)*C(9)
C     R57: OH + CH2* = H + CH2O
      RF(57) = RF(57)*C(16)
      RB(57) = RB(57)*C(20)*C(21)
C     R58: H2 + CH2* = CH3 + H
      RF(58) = RF(58)*C(7)
      RB(58) = RB(58)*C(13)*C(20)
C     R59: O2 + CH2* = CO + OH + H
      RF(59) = RF(59)*C(14)
      RB(59) = RB(59)*C(9)*C(16)*C(20)
C     R60: O2 + CH2* = CO + H2O
      RF(60) = RF(60)*C(14)
      RB(60) = RB(60)*C(9)*C(18)
C     R61: CH2* = CH2
      RF(61) = RF(61)*C(18)
      RB(61) = RB(61)*C(18)
C     R62: CH2* = CH2
      RF(62) = RF(62)*C(9)
      RB(62) = RB(62)*C(9)
C     R63: CH2* = CH2
      RF(63) = RF(63)*C(22)
      RB(63) = RB(63)*C(22)
C     R64: CH2* + CO2 = CO + CH2O
      RF(64) = RF(64)*C(22)
      RB(64) = RB(64)*C(9)*C(21)
C     R65: H + CH2O = CH3O
      CTB = CTOT+C(7)+5.D0*C(18)
     * +C(3)+5.D-1*C(9)
     * +C(22)+2.D0*C(8)
      PR = RKLOW(8) * CTB / RF(65)
      PCOR = PR / (1.0 + PR)
      PRLOG = LOG10(MAX(PR,SMALL))
      FCENT = 2.42D-1*EXP(-T/9.4D1)
     *      + 7.58D-1*EXP(-T/1.555D3)
     *     + EXP(-4.2D3/T)
      FCLOG = LOG10(MAX(FCENT,SMALL))
      XN    = 0.75 - 1.27*FCLOG
      CPRLOG= PRLOG - (0.4 + 0.67*FCLOG)
      FLOG = FCLOG/(1.0 + (CPRLOG/(XN-0.14*CPRLOG))**2)
      FC = 10.0**FLOG
      PCOR = FC * PCOR
      RF(65) = RF(65) * PCOR
      RB(65) = RB(65) * PCOR
      RF(65) = RF(65)*C(20)*C(21)
C     R66: H + CH2O = H2 + HCO
      RF(66) = RF(66)*C(20)*C(21)
      RB(66) = RB(66)*C(7)
C     R67: O + CH2O = OH + HCO
      RF(67) = RF(67)*C(15)*C(21)
      RB(67) = RB(67)*C(16)
C     R68: OH + CH2O = H2O + HCO
      RF(68) = RF(68)*C(16)*C(21)
      RB(68) = RB(68)*C(18)
C     R69: O2 + CH2O = HO2 + HCO
      RF(69) = RF(69)*C(14)*C(21)
      RB(69) = RB(69)*C(17)
C     R70: HO2 + CH2O = H2O2 + HCO
      RF(70) = RF(70)*C(17)*C(21)
      RB(70) = RB(70)*C(19)
C     R71: CH3 + H = CH4
      CTB = CTOT+C(7)+5.D0*C(18)
     * +C(3)+5.D-1*C(9)
     * +C(22)+2.D0*C(8)
      PR = RKLOW(9) * CTB / RF(71)
      PCOR = PR / (1.0 + PR)
      PRLOG = LOG10(MAX(PR,SMALL))
      FCENT = 2.17D-1*EXP(-T/7.4D1)
     *      + 7.83D-1*EXP(-T/2.941D3)
     *     + EXP(-6.964D3/T)
      FCLOG = LOG10(MAX(FCENT,SMALL))
      XN    = 0.75 - 1.27*FCLOG
      CPRLOG= PRLOG - (0.4 + 0.67*FCLOG)
      FLOG = FCLOG/(1.0 + (CPRLOG/(XN-0.14*CPRLOG))**2)
      FC = 10.0**FLOG
      PCOR = FC * PCOR
      RF(71) = RF(71) * PCOR
      RB(71) = RB(71) * PCOR
      RF(71) = RF(71)*C(13)*C(20)
      RB(71) = RB(71)*C(3)
C     R72: CH3 + O = H + CH2O
      RF(72) = RF(72)*C(13)*C(15)
      RB(72) = RB(72)*C(20)*C(21)
C     R73: CH3 + OH = H2O + CH2
      RF(73) = RF(73)*C(13)*C(16)
      RB(73) = RB(73)*C(18)
C     R74: CH3 + OH = H2O + CH2*
      RF(74) = RF(74)*C(13)*C(16)
      RB(74) = RB(74)*C(18)
C     R75: CH3 + O2 = O + CH3O
      RF(75) = RF(75)*C(13)*C(14)
      RB(75) = RB(75)*C(15)
C     R76: CH3 + O2 = OH + CH2O
      RF(76) = RF(76)*C(13)*C(14)
      RB(76) = RB(76)*C(16)*C(21)
C     R77: CH3 + HO2 = CH4 + O2
      RF(77) = RF(77)*C(13)*C(17)
      RB(77) = RB(77)*C(3)*C(14)
C     R78: CH3 + HO2 = OH + CH3O
      RF(78) = RF(78)*C(13)*C(17)
      RB(78) = RB(78)*C(16)
C     R79: CH3 + H2O2 = CH4 + HO2
      RF(79) = RF(79)*C(13)*C(19)
      RB(79) = RB(79)*C(3)*C(17)
C     R80: CH3 + HCO = CH4 + CO
      RF(80) = RF(80)*C(13)
      RB(80) = RB(80)*C(3)*C(9)
C     R81: CH3 + CH2O = CH4 + HCO
      RF(81) = RF(81)*C(13)*C(21)
      RB(81) = RB(81)*C(3)
C     R82: CH3 + CH2 = C2H4 + H
      RF(82) = RF(82)*C(13)
      RB(82) = RB(82)*C(2)*C(20)
C     R83: 2CH3 = C2H6
      CTB = CTOT+C(7)+5.D0*C(18)
     * +C(3)+5.D-1*C(9)
     * +C(22)+2.D0*C(8)
      PR = RKLOW(10) * CTB / RF(83)
      PCOR = PR / (1.0 + PR)
      PRLOG = LOG10(MAX(PR,SMALL))
      FCENT = 4.675D-1*EXP(-T/1.51D2)
     *      + 5.325D-1*EXP(-T/1.038D3)
     *     + EXP(-4.97D3/T)
      FCLOG = LOG10(MAX(FCENT,SMALL))
      XN    = 0.75 - 1.27*FCLOG
      CPRLOG= PRLOG - (0.4 + 0.67*FCLOG)
      FLOG = FCLOG/(1.0 + (CPRLOG/(XN-0.14*CPRLOG))**2)
      FC = 10.0**FLOG
      PCOR = FC * PCOR
      RF(83) = RF(83) * PCOR
      RB(83) = RB(83) * PCOR
      RF(83) = RF(83)*C(13)*C(13)
      RB(83) = RB(83)*C(8)
C     R84: 2CH3 = H + C2H5
      RF(84) = RF(84)*C(13)*C(13)
      RB(84) = RB(84)*C(20)
C     R85: CH3 + HCCO = C2H4 + CO
      RF(85) = RF(85)*C(13)
      RB(85) = RB(85)*C(2)*C(9)
C     R86: H + CH3O = H2 + CH2O
      RF(86) = RF(86)*C(20)
      RB(86) = RB(86)*C(7)*C(21)
C     R87: H + CH3O = CH3 + OH
      RF(87) = RF(87)*C(20)
      RB(87) = RB(87)*C(13)*C(16)
C     R88: H + CH3O = H2O + CH2*
      RF(88) = RF(88)*C(20)
      RB(88) = RB(88)*C(18)
C     R89: OH + CH3O = H2O + CH2O
      RF(89) = RF(89)*C(16)
      RB(89) = RB(89)*C(18)*C(21)
C     R90: O2 + CH3O = HO2 + CH2O
      RF(90) = RF(90)*C(14)
      RB(90) = RB(90)*C(17)*C(21)
C     R91: CH4 + H = H2 + CH3
      RF(91) = RF(91)*C(3)*C(20)
      RB(91) = RB(91)*C(7)*C(13)
C     R92: CH4 + O = CH3 + OH
      RF(92) = RF(92)*C(3)*C(15)
      RB(92) = RB(92)*C(13)*C(16)
C     R93: CH4 + OH = CH3 + H2O
      RF(93) = RF(93)*C(3)*C(16)
      RB(93) = RB(93)*C(13)*C(18)
C     R94: CH4 + CH2 = 2CH3
      RF(94) = RF(94)*C(3)
      RB(94) = RB(94)*C(13)*C(13)
C     R95: CH4 + CH2* = 2CH3
      RF(95) = RF(95)*C(3)
      RB(95) = RB(95)*C(13)*C(13)
C     R96: H + HCCO = CO + CH2*
      RF(96) = RF(96)*C(20)
      RB(96) = RB(96)*C(9)
C     R97: O + HCCO = 2CO + H
      RF(97) = RF(97)*C(15)
      RB(97) = RB(97)*C(9)*C(9)*C(20)
C     R98: O2 + HCCO = 2CO + OH
      RF(98) = RF(98)*C(14)
      RB(98) = RB(98)*C(9)*C(9)*C(16)
C     R99: CH2 + HCCO = CO + C2H3
      RB(99) = RB(99)*C(9)
C     R100: C2H3 = C2H2 + H
      CTB = CTOT+C(7)+5.D0*C(18)
     * +C(3)+5.D-1*C(9)
     * +C(22)+2.D0*C(8)
     * +2.D0*C(11)+2.D0*C(2)
      PR = RKLOW(11) * CTB / RF(100)
      PCOR = PR / (1.0 + PR)
      PRLOG = LOG10(MAX(PR,SMALL))
      FCENT = -9.816D-1*EXP(-T/5.3837D3)
     *      + 1.9816D0*EXP(-T/4.2932D0)
     *     + EXP(7.95D-2/T)
      FCLOG = LOG10(MAX(FCENT,SMALL))
      XN    = 0.75 - 1.27*FCLOG
      CPRLOG= PRLOG - (0.4 + 0.67*FCLOG)
      FLOG = FCLOG/(1.0 + (CPRLOG/(XN-0.14*CPRLOG))**2)
      FC = 10.0**FLOG
      PCOR = FC * PCOR
      RF(100) = RF(100) * PCOR
      RB(100) = RB(100) * PCOR
      RB(100) = RB(100)*C(11)*C(20)
C     R101: C2H2 + O = CO + CH2
      RF(101) = RF(101)*C(11)*C(15)
      RB(101) = RB(101)*C(9)
C     R102: C2H2 + O = H + HCCO
      RF(102) = RF(102)*C(11)*C(15)
      RB(102) = RB(102)*C(20)
C     R103: C2H2 + OH = H + CH2CO
      RF(103) = RF(103)*C(11)*C(16)
      RB(103) = RB(103)*C(20)*C(23)
C     R104: C2H2 + OH = CO + CH3
      RF(104) = RF(104)*C(11)*C(16)
      RB(104) = RB(104)*C(9)*C(13)
C     R105: C2H2 + HCO = CO + C2H3
      RF(105) = RF(105)*C(11)
      RB(105) = RB(105)*C(9)
C     R106: C2H2 + CH2 = H + C3H3
      RF(106) = RF(106)*C(11)
      RB(106) = RB(106)*C(20)*C(24)
C     R107: C2H2 + CH2* = H + C3H3
      RF(107) = RF(107)*C(11)
      RB(107) = RB(107)*C(20)*C(24)
C     R108: C2H2 + CH3 = aC3H5
      RF(108) = RF(108)*C(11)*C(13)
      RB(108) = RB(108)*C(25)
C     R109: H + CH2CO = CH2CHO
      CTB = CTOT+C(7)+5.D0*C(18)
     * +C(3)+5.D-1*C(9)
     * +C(22)+2.D0*C(8)
     * +2.D0*C(11)+2.D0*C(2)
      PR = RKLOW(12) * CTB / RF(109)
      PCOR = PR / (1.0 + PR)
      PRLOG = LOG10(MAX(PR,SMALL))
      FCENT = 6.63D-1*EXP(-T/1.707D3)
     *      + 3.37D-1*EXP(-T/3.2D3)
     *     + EXP(-4.131D3/T)
      FCLOG = LOG10(MAX(FCENT,SMALL))
      XN    = 0.75 - 1.27*FCLOG
      CPRLOG= PRLOG - (0.4 + 0.67*FCLOG)
      FLOG = FCLOG/(1.0 + (CPRLOG/(XN-0.14*CPRLOG))**2)
      FC = 10.0**FLOG
      PCOR = FC * PCOR
      RF(109) = RF(109) * PCOR
      RB(109) = RB(109) * PCOR
      RF(109) = RF(109)*C(20)*C(23)
C     R110: H + CH2CO = H2 + HCCO
      RF(110) = RF(110)*C(20)*C(23)
      RB(110) = RB(110)*C(7)
C     R111: H + CH2CO = CO + CH3
      RF(111) = RF(111)*C(20)*C(23)
      RB(111) = RB(111)*C(9)*C(13)
C     R112: O + CH2CO = OH + HCCO
      RF(112) = RF(112)*C(15)*C(23)
      RB(112) = RB(112)*C(16)
C     R113: O + CH2CO = CH2 + CO2
      RF(113) = RF(113)*C(15)*C(23)
      RB(113) = RB(113)*C(22)
C     R114: OH + CH2CO = H2O + HCCO
      RF(114) = RF(114)*C(16)*C(23)
      RB(114) = RB(114)*C(18)
C     R115: H + C2H3 = C2H4
      CTB = CTOT+C(7)+5.D0*C(18)
     * +C(3)+5.D-1*C(9)
     * +C(22)+2.D0*C(8)
     * +2.D0*C(11)+2.D0*C(2)
      PR = RKLOW(13) * CTB / RF(115)
      PCOR = PR / (1.0 + PR)
      PRLOG = LOG10(MAX(PR,SMALL))
      FCENT = 2.18D-1*EXP(-T/2.075D2)
     *      + 7.82D-1*EXP(-T/2.663D3)
     *     + EXP(-6.095D3/T)
      FCLOG = LOG10(MAX(FCENT,SMALL))
      XN    = 0.75 - 1.27*FCLOG
      CPRLOG= PRLOG - (0.4 + 0.67*FCLOG)
      FLOG = FCLOG/(1.0 + (CPRLOG/(XN-0.14*CPRLOG))**2)
      FC = 10.0**FLOG
      PCOR = FC * PCOR
      RF(115) = RF(115) * PCOR
      RB(115) = RB(115) * PCOR
      RF(115) = RF(115)*C(20)
      RB(115) = RB(115)*C(2)
C     R116: H + C2H3 = H2 + C2H2
      RF(116) = RF(116)*C(20)
      RB(116) = RB(116)*C(7)*C(11)
C     R117: O + C2H3 = H + CH2CO
      RF(117) = RF(117)*C(15)
      RB(117) = RB(117)*C(20)*C(23)
C     R118: O + C2H3 = CO + CH3
      RF(118) = RF(118)*C(15)
      RB(118) = RB(118)*C(9)*C(13)
C     R119: OH + C2H3 = C2H2 + H2O
      RF(119) = RF(119)*C(16)
      RB(119) = RB(119)*C(11)*C(18)
C     R120: O2 + C2H3 = C2H2 + HO2
      RF(120) = RF(120)*C(14)
      RB(120) = RB(120)*C(11)*C(17)
C     R121: O2 + C2H3 = O + CH2CHO
      RF(121) = RF(121)*C(14)
      RB(121) = RB(121)*C(15)
C     R122: O2 + C2H3 = HCO + CH2O
      RF(122) = RF(122)*C(14)
      RB(122) = RB(122)*C(21)
C     R123: CH3 + C2H3 = C3H6
      CTB = CTOT+C(7)+5.D0*C(18)
     * +C(3)+5.D-1*C(9)
     * +C(22)+2.D0*C(8)
     * +2.D0*C(2)
      PR = RKLOW(14) * CTB / RF(123)
      PCOR = PR / (1.0 + PR)
      PRLOG = LOG10(MAX(PR,SMALL))
      FCENT = 8.25D-1*EXP(-T/1.3406D3)
     *      + 1.75D-1*EXP(-T/6.D4)
     *     + EXP(-1.01398D4/T)
      FCLOG = LOG10(MAX(FCENT,SMALL))
      XN    = 0.75 - 1.27*FCLOG
      CPRLOG= PRLOG - (0.4 + 0.67*FCLOG)
      FLOG = FCLOG/(1.0 + (CPRLOG/(XN-0.14*CPRLOG))**2)
      FC = 10.0**FLOG
      PCOR = FC * PCOR
      RF(123) = RF(123) * PCOR
      RB(123) = RB(123) * PCOR
      RF(123) = RF(123)*C(13)
      RB(123) = RB(123)*C(4)
C     R124: CH3 + C2H3 = H + aC3H5
      RF(124) = RF(124)*C(13)
      RB(124) = RB(124)*C(20)*C(25)
C     R125: CH2CHO = CO + CH3
      RB(125) = RB(125)*C(9)*C(13)
C     R126: H + CH2CHO = CH3 + HCO
      RF(126) = RF(126)*C(20)
      RB(126) = RB(126)*C(13)
C     R127: H + CH2CHO = H2 + CH2CO
      RF(127) = RF(127)*C(20)
      RB(127) = RB(127)*C(7)*C(23)
C     R128: O + CH2CHO = OH + CH2CO
      RF(128) = RF(128)*C(15)
      RB(128) = RB(128)*C(16)*C(23)
C     R129: OH + CH2CHO = H2O + CH2CO
      RF(129) = RF(129)*C(16)
      RB(129) = RB(129)*C(18)*C(23)
C     R130: O2 + CH2CHO = HO2 + CH2CO
      RF(130) = RF(130)*C(14)
      RB(130) = RB(130)*C(17)*C(23)
C     R131: O2 + CH2CHO = CO + OH + CH2O
      RF(131) = RF(131)*C(14)
      RB(131) = RB(131)*C(9)*C(16)*C(21)
C     R132: C2H4 + H = C2H5
      CTB = CTOT+C(7)+5.D0*C(18)
     * +C(3)+5.D-1*C(9)
     * +C(22)+2.D0*C(8)
      PR = RKLOW(15) * CTB / RF(132)
      PCOR = PR / (1.0 + PR)
      PRLOG = LOG10(MAX(PR,SMALL))
      FCENT = -5.69D-1*EXP(T/9.147D3)
     *      + 1.569D0*EXP(-T/2.99D2)
     *     + EXP(-1.524D2/T)
      FCLOG = LOG10(MAX(FCENT,SMALL))
      XN    = 0.75 - 1.27*FCLOG
      CPRLOG= PRLOG - (0.4 + 0.67*FCLOG)
      FLOG = FCLOG/(1.0 + (CPRLOG/(XN-0.14*CPRLOG))**2)
      FC = 10.0**FLOG
      PCOR = FC * PCOR
      RF(132) = RF(132) * PCOR
      RB(132) = RB(132) * PCOR
      RF(132) = RF(132)*C(2)*C(20)
C     R133: C2H4 + H = H2 + C2H3
      RF(133) = RF(133)*C(2)*C(20)
      RB(133) = RB(133)*C(7)
C     R134: C2H4 + O = OH + C2H3
      RF(134) = RF(134)*C(2)*C(15)
      RB(134) = RB(134)*C(16)
C     R135: C2H4 + O = CH3 + HCO
      RF(135) = RF(135)*C(2)*C(15)
      RB(135) = RB(135)*C(13)
C     R136: C2H4 + O = CH2 + CH2O
      RF(136) = RF(136)*C(2)*C(15)
      RB(136) = RB(136)*C(21)
C     R137: C2H4 + OH = H2O + C2H3
      RF(137) = RF(137)*C(2)*C(16)
      RB(137) = RB(137)*C(18)
C     R138: C2H4 + HCO = CO + C2H5
      RF(138) = RF(138)*C(2)
      RB(138) = RB(138)*C(9)
C     R139: C2H4 + CH2 = H + aC3H5
      RF(139) = RF(139)*C(2)
      RB(139) = RB(139)*C(20)*C(25)
C     R140: C2H4 + CH2* = H + aC3H5
      RF(140) = RF(140)*C(2)
      RB(140) = RB(140)*C(20)*C(25)
C     R141: C2H4 + CH3 = CH4 + C2H3
      RF(141) = RF(141)*C(2)*C(13)
      RB(141) = RB(141)*C(3)
C     R142: C2H4 + O2 = HO2 + C2H3
      RF(142) = RF(142)*C(2)*C(14)
      RB(142) = RB(142)*C(17)
C     R143: H + C2H5 = C2H6
      CTB = CTOT+C(7)+5.D0*C(18)
     * +C(3)+5.D-1*C(9)
     * +C(22)+2.D0*C(8)
      PR = RKLOW(16) * CTB / RF(143)
      PCOR = PR / (1.0 + PR)
      PRLOG = LOG10(MAX(PR,SMALL))
      FCENT = 1.578000000000001D-1*EXP(-T/1.25D2)
     *      + 8.421999999999999D-1*EXP(-T/2.219D3)
     *     + EXP(-6.882D3/T)
      FCLOG = LOG10(MAX(FCENT,SMALL))
      XN    = 0.75 - 1.27*FCLOG
      CPRLOG= PRLOG - (0.4 + 0.67*FCLOG)
      FLOG = FCLOG/(1.0 + (CPRLOG/(XN-0.14*CPRLOG))**2)
      FC = 10.0**FLOG
      PCOR = FC * PCOR
      RF(143) = RF(143) * PCOR
      RB(143) = RB(143) * PCOR
      RF(143) = RF(143)*C(20)
      RB(143) = RB(143)*C(8)
C     R144: H + C2H5 = C2H4 + H2
      RF(144) = RF(144)*C(20)
      RB(144) = RB(144)*C(2)*C(7)
C     R145: O + C2H5 = CH3 + CH2O
      RF(145) = RF(145)*C(15)
      RB(145) = RB(145)*C(13)*C(21)
C     R146: O2 + C2H5 = C2H4 + HO2
      RF(146) = RF(146)*C(14)
      RB(146) = RB(146)*C(2)*C(17)
C     R147: HO2 + C2H5 = CH3 + OH + CH2O
      RF(147) = RF(147)*C(17)
      RB(147) = RB(147)*C(13)*C(16)*C(21)
C     R148: C2H3 + C2H5 = C4H81
      CTB = CTOT+C(7)+5.D0*C(18)
     * +C(3)+5.D-1*C(9)
     * +C(22)+2.D0*C(8)
      PR = RKLOW(17) * CTB / RF(148)
      PCOR = PR / (1.0 + PR)
      PRLOG = LOG10(MAX(PR,SMALL))
      FCENT = 8.02D-1*EXP(-T/2.2779D3)
     *      + 1.98D-1*EXP(-T/6.D4)
     *     + EXP(-5.7232D3/T)
      FCLOG = LOG10(MAX(FCENT,SMALL))
      XN    = 0.75 - 1.27*FCLOG
      CPRLOG= PRLOG - (0.4 + 0.67*FCLOG)
      FLOG = FCLOG/(1.0 + (CPRLOG/(XN-0.14*CPRLOG))**2)
      FC = 10.0**FLOG
      PCOR = FC * PCOR
      RF(148) = RF(148) * PCOR
      RB(148) = RB(148) * PCOR
      RB(148) = RB(148)*C(6)
C     R149: C2H3 + C2H5 = CH3 + aC3H5
      RB(149) = RB(149)*C(13)*C(25)
C     R150: C2H6 + H = H2 + C2H5
      RF(150) = RF(150)*C(8)*C(20)
      RB(150) = RB(150)*C(7)
C     R151: C2H6 + O = OH + C2H5
      RF(151) = RF(151)*C(8)*C(15)
      RB(151) = RB(151)*C(16)
C     R152: C2H6 + OH = H2O + C2H5
      RF(152) = RF(152)*C(8)*C(16)
      RB(152) = RB(152)*C(18)
C     R153: C2H6 + CH3 = CH4 + C2H5
      RF(153) = RF(153)*C(8)*C(13)
      RB(153) = RB(153)*C(3)
C     R154: O2 + C3H3 = HCO + CH2CO
      RF(154) = RF(154)*C(14)*C(24)
      RB(154) = RB(154)*C(23)
C     R155: HO2 + C3H3 = CO + OH + C2H3
      RF(155) = RF(155)*C(17)*C(24)
      RB(155) = RB(155)*C(9)*C(16)
C     R156: C2H2 + C3H3 = C5H5
      RF(156) = RF(156)*C(11)*C(24)
      RB(156) = RB(156)*C(27)
C     R157: 2C3H3 = H + C6H5
      RF(157) = RF(157)*C(24)*C(24)
C     R158: 2C3H3 = C6H6
      RF(158) = RF(158)*C(24)*C(24)
C     R159: H + aC3H5 = C3H6
      CTB = CTOT+C(7)+5.D0*C(18)
     * +C(3)+5.D-1*C(9)
     * +C(22)+2.D0*C(8)
      PR = RKLOW(18) * CTB / RF(159)
      PCOR = PR / (1.0 + PR)
      PRLOG = LOG10(MAX(PR,SMALL))
      FCENT = 9.8D-1*EXP(-T/1.0966D3)
     *      + 2.D-2*EXP(-T/1.0966D3)
     *     + EXP(-6.8595D3/T)
      FCLOG = LOG10(MAX(FCENT,SMALL))
      XN    = 0.75 - 1.27*FCLOG
      CPRLOG= PRLOG - (0.4 + 0.67*FCLOG)
      FLOG = FCLOG/(1.0 + (CPRLOG/(XN-0.14*CPRLOG))**2)
      FC = 10.0**FLOG
      PCOR = FC * PCOR
      RF(159) = RF(159) * PCOR
      RB(159) = RB(159) * PCOR
      RF(159) = RF(159)*C(20)*C(25)
      RB(159) = RB(159)*C(4)
C     R160: HO2 + aC3H5 = C3H6 + O2
      RF(160) = RF(160)*C(17)*C(25)
      RB(160) = RB(160)*C(4)*C(14)
C     R161: HO2 + aC3H5 = OH + CH2O + C2H3
      RF(161) = RF(161)*C(17)*C(25)
      RB(161) = RB(161)*C(16)*C(21)
C     R162: HCO + aC3H5 = C3H6 + CO
      RF(162) = RF(162)*C(25)
      RB(162) = RB(162)*C(4)*C(9)
C     R163: CH3 + aC3H5 = C4H81
      CTB = CTOT+C(7)+5.D0*C(18)
     * +C(3)+5.D-1*C(9)
     * +C(22)+2.D0*C(8)
      PR = RKLOW(19) * CTB / RF(163)
      PCOR = PR / (1.0 + PR)
      PRLOG = LOG10(MAX(PR,SMALL))
      FCENT = 8.96D-1*EXP(-T/1.606D3)
     *      + 1.04D-1*EXP(-T/6.D4)
     *     + EXP(-6.1184D3/T)
      FCLOG = LOG10(MAX(FCENT,SMALL))
      XN    = 0.75 - 1.27*FCLOG
      CPRLOG= PRLOG - (0.4 + 0.67*FCLOG)
      FLOG = FCLOG/(1.0 + (CPRLOG/(XN-0.14*CPRLOG))**2)
      FC = 10.0**FLOG
      PCOR = FC * PCOR
      RF(163) = RF(163) * PCOR
      RB(163) = RB(163) * PCOR
      RF(163) = RF(163)*C(13)*C(25)
      RB(163) = RB(163)*C(6)
C     R164: C3H6 + H = C2H4 + CH3
      RF(164) = RF(164)*C(4)*C(20)
      RB(164) = RB(164)*C(2)*C(13)
C     R165: C3H6 + H = H2 + aC3H5
      RF(165) = RF(165)*C(4)*C(20)
      RB(165) = RB(165)*C(7)*C(25)
C     R166: C3H6 + O = CH3 + H + CH2CO
      RF(166) = RF(166)*C(4)*C(15)
      RB(166) = RB(166)*C(13)*C(20)*C(23)
C     R167: C3H6 + O = HCO + C2H5
      RF(167) = RF(167)*C(4)*C(15)
C     R168: C3H6 + O = OH + aC3H5
      RF(168) = RF(168)*C(4)*C(15)
      RB(168) = RB(168)*C(16)*C(25)
C     R169: C3H6 + OH = H2O + aC3H5
      RF(169) = RF(169)*C(4)*C(16)
      RB(169) = RB(169)*C(18)*C(25)
C     R170: C3H6 + CH3 = CH4 + aC3H5
      RF(170) = RF(170)*C(4)*C(13)
      RB(170) = RB(170)*C(3)*C(25)
C     R171: C4H81 + H = C2H4 + C2H5
      RF(171) = RF(171)*C(6)*C(20)
      RB(171) = RB(171)*C(2)
C     R172: C4H81 + H = C3H6 + CH3
      RF(172) = RF(172)*C(6)*C(20)
      RB(172) = RB(172)*C(4)*C(13)
C     R173: iC4H8 + H = C3H6 + CH3
      RF(173) = RF(173)*C(5)*C(20)
      RB(173) = RB(173)*C(4)*C(13)
C     R174: iC4H8 + O = 2CH3 + CH2CO
      RF(174) = RF(174)*C(5)*C(15)
      RB(174) = RB(174)*C(13)*C(13)*C(23)
C     R175: iC4H8 + HO2 = C3H6 + OH + CH2O
      RF(175) = RF(175)*C(5)*C(17)
      RB(175) = RB(175)*C(4)*C(16)*C(21)
C     R176: CH3 + C6H5 = C6H5CH3
      RF(176) = RF(176)*C(13)
      RB(176) = RB(176)*C(12)
C     R177: C6H5CH3 + H = C6H6 + CH3
      RF(177) = RF(177)*C(12)*C(20)
      RB(177) = RB(177)*C(10)*C(13)
C     R178: H + C6H5 = C6H6
      CTB = CTOT+C(7)+5.D0*C(18)
     * +C(3)+5.D-1*C(9)
     * +C(22)
      PR = RKLOW(20) * CTB / RF(178)
      PCOR = PR / (1.0 + PR)
      PRLOG = LOG10(MAX(PR,SMALL))
      FCENT = 0.D0*EXP(-T/1.D-1)
     *      + 1.D0*EXP(-T/5.849D2)
     *     + EXP(-6.113D3/T)
      FCLOG = LOG10(MAX(FCENT,SMALL))
      XN    = 0.75 - 1.27*FCLOG
      CPRLOG= PRLOG - (0.4 + 0.67*FCLOG)
      FLOG = FCLOG/(1.0 + (CPRLOG/(XN-0.14*CPRLOG))**2)
      FC = 10.0**FLOG
      PCOR = FC * PCOR
      RF(178) = RF(178) * PCOR
      RB(178) = RB(178) * PCOR
      RF(178) = RF(178)*C(20)
      RB(178) = RB(178)*C(10)
C     R179: C6H6 + OH = H2O + C6H5
      RF(179) = RF(179)*C(10)*C(16)
      RB(179) = RB(179)*C(18)
C     R180: C6H6 + O = H + C6H5O
      RF(180) = RF(180)*C(10)*C(15)
      RB(180) = RB(180)*C(20)*C(28)
C     R181: C6H6 + O = HCO + C5H5
      RF(181) = RF(181)*C(10)*C(15)
      RB(181) = RB(181)*C(27)
C     R182: H2 + C6H5 = C6H6 + H
      RF(182) = RF(182)*C(7)
      RB(182) = RB(182)*C(10)*C(20)
C     R183: O2 + C6H5 = O + C6H5O
      RF(183) = RF(183)*C(14)
      RB(183) = RB(183)*C(15)*C(28)
C     R184: O + C6H5 = CO + C5H5
      RF(184) = RF(184)*C(15)
      RB(184) = RB(184)*C(9)*C(27)
C     R185: OH + C6H5 = H + C6H5O
      RF(185) = RF(185)*C(16)
      RB(185) = RB(185)*C(20)*C(28)
C     R186: HO2 + C6H5 = C6H6 + O2
      RF(186) = RF(186)*C(17)
      RB(186) = RB(186)*C(10)*C(14)
C     R187: CH4 + C6H5 = C6H6 + CH3
      RF(187) = RF(187)*C(3)
      RB(187) = RB(187)*C(10)*C(13)
C     R188: C6H5O = CO + C5H5
      RF(188) = RF(188)*C(28)
      RB(188) = RB(188)*C(9)*C(27)
C     R189: O + C5H5 = H + C5H4O
      RF(189) = RF(189)*C(15)*C(27)
      RB(189) = RB(189)*C(20)*C(26)
C     R190: O2 + C5H5 = OH + C5H4O
      RF(190) = RF(190)*C(14)*C(27)
      RB(190) = RB(190)*C(16)*C(26)
C     R191: C5H4O = CO + 2C2H2
      RF(191) = RF(191)*C(26)
      RB(191) = RB(191)*C(9)*C(11)*C(11)
C     R192: O + C5H4O = CO + HCO + C3H3
      RF(192) = RF(192)*C(15)*C(26)
      RB(192) = RB(192)*C(9)*C(24)
C
      END
C                                                                      C
C----------------------------------------------------------------------C
C                                                                      C
      SUBROUTINE RDOT(RF, RB, WDOT)
      IMPLICIT DOUBLE PRECISION (A-H, O-Z), INTEGER (I-N)
      DIMENSION RF(*), RB(*), WDOT(*)
C
      DO K = 1, 29
         WDOT(K) = 0D0
      ENDDO
C
      ROP = RF(1)
      WDOT(1) = WDOT(1) -ROP
      WDOT(2) = WDOT(2) +3.09172*ROP
      WDOT(4) = WDOT(4) +0.803846*ROP
      WDOT(5) = WDOT(5) +0.154586*ROP
      WDOT(6) = WDOT(6) +0.309172*ROP
      WDOT(10) = WDOT(10) +0.116667*ROP
      WDOT(12) = WDOT(12) +0.05*ROP
      WDOT(13) = WDOT(13) +1.5*ROP
      WDOT(20) = WDOT(20) +0.5*ROP
      ROP = RF(2)
      WDOT(1) = WDOT(1) -ROP
      WDOT(2) = WDOT(2) +3.09911*ROP
      WDOT(3) = WDOT(3) +0.25*ROP
      WDOT(4) = WDOT(4) +0.805769*ROP
      WDOT(5) = WDOT(5) +0.154956*ROP
      WDOT(6) = WDOT(6) +0.309911*ROP
      WDOT(7) = WDOT(7) +ROP
      WDOT(10) = WDOT(10) +0.175*ROP
      WDOT(12) = WDOT(12) +0.075*ROP
      WDOT(13) = WDOT(13) +0.7*ROP
      WDOT(20) = WDOT(20) -0.7*ROP
      ROP = RF(3)
      WDOT(1) = WDOT(1) -ROP
      WDOT(2) = WDOT(2) +3.09911*ROP
      WDOT(3) = WDOT(3) +1.25*ROP
      WDOT(4) = WDOT(4) +0.805769*ROP
      WDOT(5) = WDOT(5) +0.154956*ROP
      WDOT(6) = WDOT(6) +0.309911*ROP
      WDOT(10) = WDOT(10) +0.175*ROP
      WDOT(12) = WDOT(12) +0.075*ROP
      WDOT(13) = WDOT(13) -0.3*ROP
      WDOT(20) = WDOT(20) +0.3*ROP
      ROP = RF(4)
      WDOT(1) = WDOT(1) -ROP
      WDOT(2) = WDOT(2) +3.09911*ROP
      WDOT(3) = WDOT(3) +0.25*ROP
      WDOT(4) = WDOT(4) +0.805769*ROP
      WDOT(5) = WDOT(5) +0.154956*ROP
      WDOT(6) = WDOT(6) +0.309911*ROP
      WDOT(10) = WDOT(10) +0.175*ROP
      WDOT(12) = WDOT(12) +0.075*ROP
      WDOT(13) = WDOT(13) +0.7*ROP
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      WDOT(20) = WDOT(20) +0.3*ROP
      ROP = RF(5)
      WDOT(1) = WDOT(1) -ROP
      WDOT(2) = WDOT(2) +3.09911*ROP
      WDOT(3) = WDOT(3) +0.25*ROP
      WDOT(4) = WDOT(4) +0.805769*ROP
      WDOT(5) = WDOT(5) +0.154956*ROP
      WDOT(6) = WDOT(6) +0.309911*ROP
      WDOT(10) = WDOT(10) +0.175*ROP
      WDOT(12) = WDOT(12) +0.075*ROP
      WDOT(13) = WDOT(13) +0.7*ROP
      WDOT(14) = WDOT(14) -ROP
      WDOT(17) = WDOT(17) +ROP
      WDOT(20) = WDOT(20) +0.3*ROP
      ROP = RF(6)
      WDOT(1) = WDOT(1) -ROP
      WDOT(2) = WDOT(2) +3.09911*ROP
      WDOT(3) = WDOT(3) +0.25*ROP
      WDOT(4) = WDOT(4) +0.805769*ROP
      WDOT(5) = WDOT(5) +0.154956*ROP
      WDOT(6) = WDOT(6) +0.309911*ROP
      WDOT(10) = WDOT(10) +0.175*ROP
      WDOT(12) = WDOT(12) +0.075*ROP
      WDOT(13) = WDOT(13) +0.7*ROP
      WDOT(17) = WDOT(17) -ROP
      WDOT(19) = WDOT(19) +ROP
      WDOT(20) = WDOT(20) +0.3*ROP
      ROP = RF(7)
      WDOT(1) = WDOT(1) -ROP
      WDOT(2) = WDOT(2) +3.09911*ROP
      WDOT(3) = WDOT(3) +0.25*ROP
      WDOT(4) = WDOT(4) +0.805769*ROP
      WDOT(5) = WDOT(5) +0.154956*ROP
      WDOT(6) = WDOT(6) +0.309911*ROP
      WDOT(10) = WDOT(10) +0.175*ROP
      WDOT(12) = WDOT(12) +0.075*ROP
      WDOT(13) = WDOT(13) +0.7*ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(20) = WDOT(20) +0.3*ROP
      ROP = RF(8)-RB(8)
      WDOT(14) = WDOT(14) -ROP
      WDOT(15) = WDOT(15) +ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(9)-RB(9)
      WDOT(7) = WDOT(7) -ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(20) = WDOT(20) +ROP
      ROP = RF(10)-RB(10)
      WDOT(7) = WDOT(7) -ROP
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      WDOT(20) = WDOT(20) +ROP
      ROP = RF(11)-RB(11)
      WDOT(15) = WDOT(15) +ROP
      WDOT(16) = WDOT(16) -2*ROP
      WDOT(18) = WDOT(18) +ROP
      ROP = RF(12)-RB(12)
      WDOT(7) = WDOT(7) +ROP
      WDOT(20) = WDOT(20) -2*ROP
      ROP = RF(13)-RB(13)
      WDOT(7) = WDOT(7) +ROP
      WDOT(20) = WDOT(20) -2*ROP
      ROP = RF(14)-RB(14)
      WDOT(7) = WDOT(7) +ROP
      WDOT(20) = WDOT(20) -2*ROP
      ROP = RF(15)-RB(15)
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(16)-RB(16)
      WDOT(15) = WDOT(15) -ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(17)-RB(17)
      WDOT(14) = WDOT(14) +ROP
      WDOT(15) = WDOT(15) -2*ROP
      ROP = RF(18)-RB(18)
      WDOT(14) = WDOT(14) -ROP
      WDOT(17) = WDOT(17) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(19)-RB(19)
      WDOT(7) = WDOT(7) -ROP
      WDOT(14) = WDOT(14) -ROP
      WDOT(17) = WDOT(17) +ROP
      WDOT(20) = WDOT(20) +ROP
      ROP = RF(20)-RB(20)
      WDOT(16) = WDOT(16) -2*ROP
      WDOT(19) = WDOT(19) +ROP
      ROP = RF(21)-RB(21)
      WDOT(16) = WDOT(16) +2*ROP
      WDOT(17) = WDOT(17) -ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(22)-RB(22)
      WDOT(14) = WDOT(14) +ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(17) = WDOT(17) -ROP
      ROP = RF(23)-RB(23)
      WDOT(14) = WDOT(14) +ROP
      WDOT(17) = WDOT(17) -2*ROP
      WDOT(19) = WDOT(19) +ROP
      ROP = RF(24)-RB(24)
      WDOT(14) = WDOT(14) +ROP
      WDOT(17) = WDOT(17) -2*ROP
      WDOT(19) = WDOT(19) +ROP
      ROP = RF(25)-RB(25)
      WDOT(14) = WDOT(14) +ROP
      WDOT(16) = WDOT(16) -ROP
      WDOT(17) = WDOT(17) -ROP
      WDOT(18) = WDOT(18) +ROP
      ROP = RF(26)-RB(26)
      WDOT(14) = WDOT(14) +ROP
      WDOT(16) = WDOT(16) -ROP
      WDOT(17) = WDOT(17) -ROP
      WDOT(18) = WDOT(18) +ROP
      ROP = RF(27)-RB(27)
      WDOT(14) = WDOT(14) +ROP
      WDOT(16) = WDOT(16) -ROP
      WDOT(17) = WDOT(17) -ROP
      WDOT(18) = WDOT(18) +ROP
      ROP = RF(28)-RB(28)
      WDOT(7) = WDOT(7) +ROP
      WDOT(17) = WDOT(17) +ROP
      WDOT(19) = WDOT(19) -ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(29)-RB(29)
      WDOT(16) = WDOT(16) +ROP
      WDOT(18) = WDOT(18) +ROP
      WDOT(19) = WDOT(19) -ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(30)-RB(30)
      WDOT(15) = WDOT(15) -ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(17) = WDOT(17) +ROP
      WDOT(19) = WDOT(19) -ROP
      ROP = RF(31)-RB(31)
      WDOT(16) = WDOT(16) -ROP
      WDOT(17) = WDOT(17) +ROP
      WDOT(18) = WDOT(18) +ROP
      WDOT(19) = WDOT(19) -ROP
      ROP = RF(32)-RB(32)
      WDOT(16) = WDOT(16) -ROP
      WDOT(17) = WDOT(17) +ROP
      WDOT(18) = WDOT(18) +ROP
      WDOT(19) = WDOT(19) -ROP
      ROP = RF(33)-RB(33)
      WDOT(9) = WDOT(9) -ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(22) = WDOT(22) +ROP
      ROP = RF(34)-RB(34)
      WDOT(9) = WDOT(9) -ROP
      WDOT(16) = WDOT(16) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(22) = WDOT(22) +ROP
      ROP = RF(35)-RB(35)
      WDOT(9) = WDOT(9) -ROP
      WDOT(16) = WDOT(16) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(22) = WDOT(22) +ROP
      ROP = RF(36)-RB(36)
      WDOT(9) = WDOT(9) -ROP
      WDOT(14) = WDOT(14) -ROP
      WDOT(15) = WDOT(15) +ROP
      WDOT(22) = WDOT(22) +ROP
      ROP = RF(37)-RB(37)
      WDOT(9) = WDOT(9) -ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(17) = WDOT(17) -ROP
      WDOT(22) = WDOT(22) +ROP
      ROP = RF(38)-RB(38)
      WDOT(7) = WDOT(7) +ROP
      WDOT(9) = WDOT(9) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(39)-RB(39)
      WDOT(9) = WDOT(9) +ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(16) = WDOT(16) +ROP
      ROP = RF(40)-RB(40)
      WDOT(15) = WDOT(15) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(22) = WDOT(22) +ROP
      ROP = RF(41)-RB(41)
      WDOT(9) = WDOT(9) +ROP
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      ROP = RF(42)-RB(42)
      WDOT(9) = WDOT(9) +ROP
      WDOT(20) = WDOT(20) +ROP
      ROP = RF(43)-RB(43)
      WDOT(9) = WDOT(9) +ROP
      WDOT(20) = WDOT(20) +ROP
      ROP = RF(44)-RB(44)
      WDOT(9) = WDOT(9) +ROP
      WDOT(14) = WDOT(14) -ROP
      WDOT(17) = WDOT(17) +ROP
      ROP = RF(45)-RB(45)
      WDOT(7) = WDOT(7) -ROP
      WDOT(9) = WDOT(9) -ROP
      WDOT(21) = WDOT(21) +ROP
      ROP = RF(46)-RB(46)
      WDOT(20) = WDOT(20) -ROP
      WDOT(21) = WDOT(21) +ROP
      ROP = RF(47)-RB(47)
      WDOT(13) = WDOT(13) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(48)-RB(48)
      WDOT(15) = WDOT(15) -ROP
      WDOT(20) = WDOT(20) +ROP
      ROP = RF(49)-RB(49)
      WDOT(16) = WDOT(16) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(21) = WDOT(21) +ROP
      ROP = RF(50)-RB(50)
      WDOT(7) = WDOT(7) -ROP
      WDOT(13) = WDOT(13) +ROP
      WDOT(20) = WDOT(20) +ROP
      ROP = RF(51)-RB(51)
      WDOT(14) = WDOT(14) -ROP
      WDOT(16) = WDOT(16) +ROP
      ROP = RF(52)-RB(52)
      WDOT(14) = WDOT(14) -ROP
      WDOT(20) = WDOT(20) +2*ROP
      WDOT(22) = WDOT(22) +ROP
      ROP = RF(53)-RB(53)
      WDOT(9) = WDOT(9) -ROP
      WDOT(23) = WDOT(23) +ROP
      ROP = RF(54)-RB(54)
      WDOT(7) = WDOT(7) +ROP
      WDOT(11) = WDOT(11) +ROP
      ROP = RF(55)-RB(55)
      ROP = RF(56)-RB(56)
      WDOT(7) = WDOT(7) +ROP
      WDOT(9) = WDOT(9) +ROP
      WDOT(15) = WDOT(15) -ROP
      ROP = RF(57)-RB(57)
      WDOT(16) = WDOT(16) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(21) = WDOT(21) +ROP
      ROP = RF(58)-RB(58)
      WDOT(7) = WDOT(7) -ROP
      WDOT(13) = WDOT(13) +ROP
      WDOT(20) = WDOT(20) +ROP
      ROP = RF(59)-RB(59)
      WDOT(9) = WDOT(9) +ROP
      WDOT(14) = WDOT(14) -ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(20) = WDOT(20) +ROP
      ROP = RF(60)-RB(60)
      WDOT(9) = WDOT(9) +ROP
      WDOT(14) = WDOT(14) -ROP
      WDOT(18) = WDOT(18) +ROP
      ROP = RF(61)-RB(61)
      ROP = RF(62)-RB(62)
      ROP = RF(63)-RB(63)
      ROP = RF(64)-RB(64)
      WDOT(9) = WDOT(9) +ROP
      WDOT(21) = WDOT(21) +ROP
      WDOT(22) = WDOT(22) -ROP
      ROP = RF(65)-RB(65)
      WDOT(20) = WDOT(20) -ROP
      WDOT(21) = WDOT(21) -ROP
      ROP = RF(66)-RB(66)
      WDOT(7) = WDOT(7) +ROP
      WDOT(20) = WDOT(20) -ROP
      WDOT(21) = WDOT(21) -ROP
      ROP = RF(67)-RB(67)
      WDOT(15) = WDOT(15) -ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(21) = WDOT(21) -ROP
      ROP = RF(68)-RB(68)
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      WDOT(21) = WDOT(21) -ROP
      ROP = RF(69)-RB(69)
      WDOT(14) = WDOT(14) -ROP
      WDOT(17) = WDOT(17) +ROP
      WDOT(21) = WDOT(21) -ROP
      ROP = RF(70)-RB(70)
      WDOT(17) = WDOT(17) -ROP
      WDOT(19) = WDOT(19) +ROP
      WDOT(21) = WDOT(21) -ROP
      ROP = RF(71)-RB(71)
      WDOT(3) = WDOT(3) +ROP
      WDOT(13) = WDOT(13) -ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(72)-RB(72)
      WDOT(13) = WDOT(13) -ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(21) = WDOT(21) +ROP
      ROP = RF(73)-RB(73)
      WDOT(13) = WDOT(13) -ROP
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      ROP = RF(74)-RB(74)
      WDOT(13) = WDOT(13) -ROP
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      ROP = RF(75)-RB(75)
      WDOT(13) = WDOT(13) -ROP
      WDOT(14) = WDOT(14) -ROP
      WDOT(15) = WDOT(15) +ROP
      ROP = RF(76)-RB(76)
      WDOT(13) = WDOT(13) -ROP
      WDOT(14) = WDOT(14) -ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(21) = WDOT(21) +ROP
      ROP = RF(77)-RB(77)
      WDOT(3) = WDOT(3) +ROP
      WDOT(13) = WDOT(13) -ROP
      WDOT(14) = WDOT(14) +ROP
      WDOT(17) = WDOT(17) -ROP
      ROP = RF(78)-RB(78)
      WDOT(13) = WDOT(13) -ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(17) = WDOT(17) -ROP
      ROP = RF(79)-RB(79)
      WDOT(3) = WDOT(3) +ROP
      WDOT(13) = WDOT(13) -ROP
      WDOT(17) = WDOT(17) +ROP
      WDOT(19) = WDOT(19) -ROP
      ROP = RF(80)-RB(80)
      WDOT(3) = WDOT(3) +ROP
      WDOT(9) = WDOT(9) +ROP
      WDOT(13) = WDOT(13) -ROP
      ROP = RF(81)-RB(81)
      WDOT(3) = WDOT(3) +ROP
      WDOT(13) = WDOT(13) -ROP
      WDOT(21) = WDOT(21) -ROP
      ROP = RF(82)-RB(82)
      WDOT(2) = WDOT(2) +ROP
      WDOT(13) = WDOT(13) -ROP
      WDOT(20) = WDOT(20) +ROP
      ROP = RF(83)-RB(83)
      WDOT(8) = WDOT(8) +ROP
      WDOT(13) = WDOT(13) -2*ROP
      ROP = RF(84)-RB(84)
      WDOT(13) = WDOT(13) -2*ROP
      WDOT(20) = WDOT(20) +ROP
      ROP = RF(85)-RB(85)
      WDOT(2) = WDOT(2) +ROP
      WDOT(9) = WDOT(9) +ROP
      WDOT(13) = WDOT(13) -ROP
      ROP = RF(86)-RB(86)
      WDOT(7) = WDOT(7) +ROP
      WDOT(20) = WDOT(20) -ROP
      WDOT(21) = WDOT(21) +ROP
      ROP = RF(87)-RB(87)
      WDOT(13) = WDOT(13) +ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(88)-RB(88)
      WDOT(18) = WDOT(18) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(89)-RB(89)
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      WDOT(21) = WDOT(21) +ROP
      ROP = RF(90)-RB(90)
      WDOT(14) = WDOT(14) -ROP
      WDOT(17) = WDOT(17) +ROP
      WDOT(21) = WDOT(21) +ROP
      ROP = RF(91)-RB(91)
      WDOT(3) = WDOT(3) -ROP
      WDOT(7) = WDOT(7) +ROP
      WDOT(13) = WDOT(13) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(92)-RB(92)
      WDOT(3) = WDOT(3) -ROP
      WDOT(13) = WDOT(13) +ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(16) = WDOT(16) +ROP
      ROP = RF(93)-RB(93)
      WDOT(3) = WDOT(3) -ROP
      WDOT(13) = WDOT(13) +ROP
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      ROP = RF(94)-RB(94)
      WDOT(3) = WDOT(3) -ROP
      WDOT(13) = WDOT(13) +2*ROP
      ROP = RF(95)-RB(95)
      WDOT(3) = WDOT(3) -ROP
      WDOT(13) = WDOT(13) +2*ROP
      ROP = RF(96)-RB(96)
      WDOT(9) = WDOT(9) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(97)-RB(97)
      WDOT(9) = WDOT(9) +2*ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(20) = WDOT(20) +ROP
      ROP = RF(98)-RB(98)
      WDOT(9) = WDOT(9) +2*ROP
      WDOT(14) = WDOT(14) -ROP
      WDOT(16) = WDOT(16) +ROP
      ROP = RF(99)-RB(99)
      WDOT(9) = WDOT(9) +ROP
      ROP = RF(100)-RB(100)
      WDOT(11) = WDOT(11) +ROP
      WDOT(20) = WDOT(20) +ROP
      ROP = RF(101)-RB(101)
      WDOT(9) = WDOT(9) +ROP
      WDOT(11) = WDOT(11) -ROP
      WDOT(15) = WDOT(15) -ROP
      ROP = RF(102)-RB(102)
      WDOT(11) = WDOT(11) -ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(20) = WDOT(20) +ROP
      ROP = RF(103)-RB(103)
      WDOT(11) = WDOT(11) -ROP
      WDOT(16) = WDOT(16) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(23) = WDOT(23) +ROP
      ROP = RF(104)-RB(104)
      WDOT(9) = WDOT(9) +ROP
      WDOT(11) = WDOT(11) -ROP
      WDOT(13) = WDOT(13) +ROP
      WDOT(16) = WDOT(16) -ROP
      ROP = RF(105)-RB(105)
      WDOT(9) = WDOT(9) +ROP
      WDOT(11) = WDOT(11) -ROP
      ROP = RF(106)-RB(106)
      WDOT(11) = WDOT(11) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(24) = WDOT(24) +ROP
      ROP = RF(107)-RB(107)
      WDOT(11) = WDOT(11) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(24) = WDOT(24) +ROP
      ROP = RF(108)-RB(108)
      WDOT(11) = WDOT(11) -ROP
      WDOT(13) = WDOT(13) -ROP
      WDOT(25) = WDOT(25) +ROP
      ROP = RF(109)-RB(109)
      WDOT(20) = WDOT(20) -ROP
      WDOT(23) = WDOT(23) -ROP
      ROP = RF(110)-RB(110)
      WDOT(7) = WDOT(7) +ROP
      WDOT(20) = WDOT(20) -ROP
      WDOT(23) = WDOT(23) -ROP
      ROP = RF(111)-RB(111)
      WDOT(9) = WDOT(9) +ROP
      WDOT(13) = WDOT(13) +ROP
      WDOT(20) = WDOT(20) -ROP
      WDOT(23) = WDOT(23) -ROP
      ROP = RF(112)-RB(112)
      WDOT(15) = WDOT(15) -ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(23) = WDOT(23) -ROP
      ROP = RF(113)-RB(113)
      WDOT(15) = WDOT(15) -ROP
      WDOT(22) = WDOT(22) +ROP
      WDOT(23) = WDOT(23) -ROP
      ROP = RF(114)-RB(114)
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      WDOT(23) = WDOT(23) -ROP
      ROP = RF(115)-RB(115)
      WDOT(2) = WDOT(2) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(116)-RB(116)
      WDOT(7) = WDOT(7) +ROP
      WDOT(11) = WDOT(11) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(117)-RB(117)
      WDOT(15) = WDOT(15) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(23) = WDOT(23) +ROP
      ROP = RF(118)-RB(118)
      WDOT(9) = WDOT(9) +ROP
      WDOT(13) = WDOT(13) +ROP
      WDOT(15) = WDOT(15) -ROP
      ROP = RF(119)-RB(119)
      WDOT(11) = WDOT(11) +ROP
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      ROP = RF(120)-RB(120)
      WDOT(11) = WDOT(11) +ROP
      WDOT(14) = WDOT(14) -ROP
      WDOT(17) = WDOT(17) +ROP
      ROP = RF(121)-RB(121)
      WDOT(14) = WDOT(14) -ROP
      WDOT(15) = WDOT(15) +ROP
      ROP = RF(122)-RB(122)
      WDOT(14) = WDOT(14) -ROP
      WDOT(21) = WDOT(21) +ROP
      ROP = RF(123)-RB(123)
      WDOT(4) = WDOT(4) +ROP
      WDOT(13) = WDOT(13) -ROP
      ROP = RF(124)-RB(124)
      WDOT(13) = WDOT(13) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(25) = WDOT(25) +ROP
      ROP = RF(125)-RB(125)
      WDOT(9) = WDOT(9) +ROP
      WDOT(13) = WDOT(13) +ROP
      ROP = RF(126)-RB(126)
      WDOT(13) = WDOT(13) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(127)-RB(127)
      WDOT(7) = WDOT(7) +ROP
      WDOT(20) = WDOT(20) -ROP
      WDOT(23) = WDOT(23) +ROP
      ROP = RF(128)-RB(128)
      WDOT(15) = WDOT(15) -ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(23) = WDOT(23) +ROP
      ROP = RF(129)-RB(129)
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      WDOT(23) = WDOT(23) +ROP
      ROP = RF(130)-RB(130)
      WDOT(14) = WDOT(14) -ROP
      WDOT(17) = WDOT(17) +ROP
      WDOT(23) = WDOT(23) +ROP
      ROP = RF(131)-RB(131)
      WDOT(9) = WDOT(9) +ROP
      WDOT(14) = WDOT(14) -ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(21) = WDOT(21) +ROP
      ROP = RF(132)-RB(132)
      WDOT(2) = WDOT(2) -ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(133)-RB(133)
      WDOT(2) = WDOT(2) -ROP
      WDOT(7) = WDOT(7) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(134)-RB(134)
      WDOT(2) = WDOT(2) -ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(16) = WDOT(16) +ROP
      ROP = RF(135)-RB(135)
      WDOT(2) = WDOT(2) -ROP
      WDOT(13) = WDOT(13) +ROP
      WDOT(15) = WDOT(15) -ROP
      ROP = RF(136)-RB(136)
      WDOT(2) = WDOT(2) -ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(21) = WDOT(21) +ROP
      ROP = RF(137)-RB(137)
      WDOT(2) = WDOT(2) -ROP
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      ROP = RF(138)-RB(138)
      WDOT(2) = WDOT(2) -ROP
      WDOT(9) = WDOT(9) +ROP
      ROP = RF(139)-RB(139)
      WDOT(2) = WDOT(2) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(25) = WDOT(25) +ROP
      ROP = RF(140)-RB(140)
      WDOT(2) = WDOT(2) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(25) = WDOT(25) +ROP
      ROP = RF(141)-RB(141)
      WDOT(2) = WDOT(2) -ROP
      WDOT(3) = WDOT(3) +ROP
      WDOT(13) = WDOT(13) -ROP
      ROP = RF(142)-RB(142)
      WDOT(2) = WDOT(2) -ROP
      WDOT(14) = WDOT(14) -ROP
      WDOT(17) = WDOT(17) +ROP
      ROP = RF(143)-RB(143)
      WDOT(8) = WDOT(8) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(144)-RB(144)
      WDOT(2) = WDOT(2) +ROP
      WDOT(7) = WDOT(7) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(145)-RB(145)
      WDOT(13) = WDOT(13) +ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(21) = WDOT(21) +ROP
      ROP = RF(146)-RB(146)
      WDOT(2) = WDOT(2) +ROP
      WDOT(14) = WDOT(14) -ROP
      WDOT(17) = WDOT(17) +ROP
      ROP = RF(147)-RB(147)
      WDOT(13) = WDOT(13) +ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(17) = WDOT(17) -ROP
      WDOT(21) = WDOT(21) +ROP
      ROP = RF(148)-RB(148)
      WDOT(6) = WDOT(6) +ROP
      ROP = RF(149)-RB(149)
      WDOT(13) = WDOT(13) +ROP
      WDOT(25) = WDOT(25) +ROP
      ROP = RF(150)-RB(150)
      WDOT(7) = WDOT(7) +ROP
      WDOT(8) = WDOT(8) -ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(151)-RB(151)
      WDOT(8) = WDOT(8) -ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(16) = WDOT(16) +ROP
      ROP = RF(152)-RB(152)
      WDOT(8) = WDOT(8) -ROP
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      ROP = RF(153)-RB(153)
      WDOT(3) = WDOT(3) +ROP
      WDOT(8) = WDOT(8) -ROP
      WDOT(13) = WDOT(13) -ROP
      ROP = RF(154)-RB(154)
      WDOT(14) = WDOT(14) -ROP
      WDOT(23) = WDOT(23) +ROP
      WDOT(24) = WDOT(24) -ROP
      ROP = RF(155)-RB(155)
      WDOT(9) = WDOT(9) +ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(17) = WDOT(17) -ROP
      WDOT(24) = WDOT(24) -ROP
      ROP = RF(156)-RB(156)
      WDOT(11) = WDOT(11) -ROP
      WDOT(24) = WDOT(24) -ROP
      WDOT(27) = WDOT(27) +ROP
      ROP = RF(157)
      WDOT(20) = WDOT(20) +ROP
      WDOT(24) = WDOT(24) -2*ROP
      ROP = RF(158)
      WDOT(10) = WDOT(10) +ROP
      WDOT(24) = WDOT(24) -2*ROP
      ROP = RF(159)-RB(159)
      WDOT(4) = WDOT(4) +ROP
      WDOT(20) = WDOT(20) -ROP
      WDOT(25) = WDOT(25) -ROP
      ROP = RF(160)-RB(160)
      WDOT(4) = WDOT(4) +ROP
      WDOT(14) = WDOT(14) +ROP
      WDOT(17) = WDOT(17) -ROP
      WDOT(25) = WDOT(25) -ROP
      ROP = RF(161)-RB(161)
      WDOT(16) = WDOT(16) +ROP
      WDOT(17) = WDOT(17) -ROP
      WDOT(21) = WDOT(21) +ROP
      WDOT(25) = WDOT(25) -ROP
      ROP = RF(162)-RB(162)
      WDOT(4) = WDOT(4) +ROP
      WDOT(9) = WDOT(9) +ROP
      WDOT(25) = WDOT(25) -ROP
      ROP = RF(163)-RB(163)
      WDOT(6) = WDOT(6) +ROP
      WDOT(13) = WDOT(13) -ROP
      WDOT(25) = WDOT(25) -ROP
      ROP = RF(164)-RB(164)
      WDOT(2) = WDOT(2) +ROP
      WDOT(4) = WDOT(4) -ROP
      WDOT(13) = WDOT(13) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(165)-RB(165)
      WDOT(4) = WDOT(4) -ROP
      WDOT(7) = WDOT(7) +ROP
      WDOT(20) = WDOT(20) -ROP
      WDOT(25) = WDOT(25) +ROP
      ROP = RF(166)-RB(166)
      WDOT(4) = WDOT(4) -ROP
      WDOT(13) = WDOT(13) +ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(23) = WDOT(23) +ROP
      ROP = RF(167)-RB(167)
      WDOT(4) = WDOT(4) -ROP
      WDOT(15) = WDOT(15) -ROP
      ROP = RF(168)-RB(168)
      WDOT(4) = WDOT(4) -ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(25) = WDOT(25) +ROP
      ROP = RF(169)-RB(169)
      WDOT(4) = WDOT(4) -ROP
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      WDOT(25) = WDOT(25) +ROP
      ROP = RF(170)-RB(170)
      WDOT(3) = WDOT(3) +ROP
      WDOT(4) = WDOT(4) -ROP
      WDOT(13) = WDOT(13) -ROP
      WDOT(25) = WDOT(25) +ROP
      ROP = RF(171)-RB(171)
      WDOT(2) = WDOT(2) +ROP
      WDOT(6) = WDOT(6) -ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(172)-RB(172)
      WDOT(4) = WDOT(4) +ROP
      WDOT(6) = WDOT(6) -ROP
      WDOT(13) = WDOT(13) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(173)-RB(173)
      WDOT(4) = WDOT(4) +ROP
      WDOT(5) = WDOT(5) -ROP
      WDOT(13) = WDOT(13) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(174)-RB(174)
      WDOT(5) = WDOT(5) -ROP
      WDOT(13) = WDOT(13) +2*ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(23) = WDOT(23) +ROP
      ROP = RF(175)-RB(175)
      WDOT(4) = WDOT(4) +ROP
      WDOT(5) = WDOT(5) -ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(17) = WDOT(17) -ROP
      WDOT(21) = WDOT(21) +ROP
      ROP = RF(176)-RB(176)
      WDOT(12) = WDOT(12) +ROP
      WDOT(13) = WDOT(13) -ROP
      ROP = RF(177)-RB(177)
      WDOT(10) = WDOT(10) +ROP
      WDOT(12) = WDOT(12) -ROP
      WDOT(13) = WDOT(13) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(178)-RB(178)
      WDOT(10) = WDOT(10) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(179)-RB(179)
      WDOT(10) = WDOT(10) -ROP
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      ROP = RF(180)-RB(180)
      WDOT(10) = WDOT(10) -ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(28) = WDOT(28) +ROP
      ROP = RF(181)-RB(181)
      WDOT(10) = WDOT(10) -ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(27) = WDOT(27) +ROP
      ROP = RF(182)-RB(182)
      WDOT(7) = WDOT(7) -ROP
      WDOT(10) = WDOT(10) +ROP
      WDOT(20) = WDOT(20) +ROP
      ROP = RF(183)-RB(183)
      WDOT(14) = WDOT(14) -ROP
      WDOT(15) = WDOT(15) +ROP
      WDOT(28) = WDOT(28) +ROP
      ROP = RF(184)-RB(184)
      WDOT(9) = WDOT(9) +ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(27) = WDOT(27) +ROP
      ROP = RF(185)-RB(185)
      WDOT(16) = WDOT(16) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(28) = WDOT(28) +ROP
      ROP = RF(186)-RB(186)
      WDOT(10) = WDOT(10) +ROP
      WDOT(14) = WDOT(14) +ROP
      WDOT(17) = WDOT(17) -ROP
      ROP = RF(187)-RB(187)
      WDOT(3) = WDOT(3) -ROP
      WDOT(10) = WDOT(10) +ROP
      WDOT(13) = WDOT(13) +ROP
      ROP = RF(188)-RB(188)
      WDOT(9) = WDOT(9) +ROP
      WDOT(27) = WDOT(27) +ROP
      WDOT(28) = WDOT(28) -ROP
      ROP = RF(189)-RB(189)
      WDOT(15) = WDOT(15) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(26) = WDOT(26) +ROP
      WDOT(27) = WDOT(27) -ROP
      ROP = RF(190)-RB(190)
      WDOT(14) = WDOT(14) -ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(26) = WDOT(26) +ROP
      WDOT(27) = WDOT(27) -ROP
      ROP = RF(191)-RB(191)
      WDOT(9) = WDOT(9) +ROP
      WDOT(11) = WDOT(11) +2*ROP
      WDOT(26) = WDOT(26) -ROP
      ROP = RF(192)-RB(192)
      WDOT(9) = WDOT(9) +ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(24) = WDOT(24) +ROP
      WDOT(26) = WDOT(26) -ROP
C
      END
C                                                                      C
C----------------------------------------------------------------------C
C                                                                      C
      SUBROUTINE QSSA(RF, RB, XQ)
      IMPLICIT DOUBLE PRECISION (A-H, O-Z), INTEGER (I-N)
      DIMENSION RF(*), RB(*), XQ(*)
C
      RF(54) = 0.D0
      RF(99) = 0.D0
      RF(148) = 0.D0
      RF(149) = 0.D0
      RB(167) = 0.D0
C
C     CH2
      DEN = +RF( 47) +RF( 48) +RF( 49) +RF( 50) +RF( 51) 
     *  +RF( 52) +RF( 53) +RF( 82) +RF( 94) +RF(106) +RF(139) 
     *  +RB( 55) +RB( 61) +RB( 62) +RB( 63) +RB( 73) +RB(101) 
     *  +RB(113) +RB(136) 
      A1_0 = ( +RB( 47) +RB( 49) +RB( 50) +RB( 52) +RB( 53) 
     *  +2*RB( 54) +RF( 73) +RB( 82) +RB( 94) +RF(101) +RB(106) 
     *  +RF(113) +RF(136) +RB(139) )/DEN
      A1_2 = ( +RF( 55) +RF( 61) +RF( 62) +RF( 63) )/DEN
      A1_3 = ( +RB( 48) +RB( 51) )/DEN
      A1_5 = ( +RB( 99) )/DEN
C     CH2*
      DEN = +RF( 55) +RF( 56) +RF( 57) +RF( 58) +RF( 59) 
     *  +RF( 60) +RF( 61) +RF( 62) +RF( 63) +RF( 64) +RF( 95) 
     *  +RF(107) +RF(140) +RB( 74) +RB( 88) +RB( 96) 
      A2_0 = ( +RB( 56) +RB( 57) +RB( 58) +RB( 59) +RB( 60) 
     *  +RB( 64) +RF( 74) +RB( 95) +RB(107) +RB(140) )/DEN
      A2_1 = ( +RB( 55) +RB( 61) +RB( 62) +RB( 63) )/DEN
      A2_4 = ( +RF( 88) )/DEN
      A2_7 = ( +RF( 96) )/DEN
C     HCO
      DEN = +RF( 38) +RF( 39) +RF( 40) +RF( 41) +RF( 42) 
     *  +RF( 43) +RF( 44) +RF( 46) +RF( 80) +RF(105) +RF(138) 
     *  +RF(162) +RB( 48) +RB( 51) +RB( 66) +RB( 67) +RB( 68) 
     *  +RB( 69) +RB( 70) +RB( 81) +RB(122) +RB(126) +RB(135) 
     *  +RB(154) +RB(181) +RB(192) 
      A3_0 = ( +RB( 38) +RB( 39) +RB( 40) +RB( 41) +RB( 42) 
     *  +RB( 43) +RB( 44) +RB( 46) +RF( 66) +RF( 67) +RF( 68) 
     *  +RF( 69) +RF( 70) +RB( 80) +RF( 81) +RF(135) +RF(154) 
     *  +RB(162) +RF(167) +RF(181) +RF(192) )/DEN
      A3_1 = ( +RF( 48) +RF( 51) )/DEN
      A3_5 = ( +RB(105) +RF(122) )/DEN
      A3_6 = ( +RB(138) )/DEN
      A3_8 = ( +RF(126) )/DEN
C     CH3O
      DEN = +RF( 86) +RF( 87) +RF( 88) +RF( 89) +RF( 90) 
     *  +RB( 65) +RB( 75) +RB( 78) 
      A4_0 = ( +RF( 65) +RF( 75) +RF( 78) +RB( 86) +RB( 87) 
     *  +RB( 89) +RB( 90) )/DEN
      A4_2 = ( +RB( 88) )/DEN
C     C2H3
      DEN = +RF(100) +RF(115) +RF(116) +RF(117) +RF(118) 
     *  +RF(119) +RF(120) +RF(121) +RF(122) +RF(123) +RF(124) 
     *  +RB( 99) +RB(105) +RB(133) +RB(134) +RB(137) +RB(141) 
     *  +RB(142) +RB(155) +RB(161) 
      A5_0 = ( +RB(100) +RB(115) +RB(116) +RB(117) +RB(118) 
     *  +RB(119) +RB(120) +RB(123) +RB(124) +RF(133) +RF(134) 
     *  +RF(137) +RF(141) +RF(142) +RB(148) +RB(149) +RF(155) 
     *  +RF(161) )/DEN
      A5_3 = ( +RF(105) +RB(122) )/DEN
      A5_8 = ( +RB(121) )/DEN
C     C2H5
      DEN = +RF(143) +RF(144) +RF(145) +RF(146) +RF(147) 
     *  +RB( 84) +RB(132) +RB(138) +RB(150) +RB(151) +RB(152) 
     *  +RB(153) +RB(171) 
      A6_0 = ( +RF( 84) +RF(132) +RB(143) +RB(144) +RB(145) 
     *  +RB(146) +RB(147) +RB(148) +RB(149) +RF(150) +RF(151) 
     *  +RF(152) +RF(153) +RF(167) +RF(171) )/DEN
      A6_3 = ( +RF(138) )/DEN
C     HCCO
      DEN = +RF( 85) +RF( 96) +RF( 97) +RF( 98) +RB(102) 
     *  +RB(110) +RB(112) +RB(114) 
      A7_0 = ( +RB( 85) +RB( 97) +RB( 98) +RF(102) +RF(110) 
     *  +RF(112) +RF(114) )/DEN
      A7_2 = ( +RB( 96) )/DEN
      A7_5 = ( +RB( 99) )/DEN
C     CH2CHO
      DEN = +RF(125) +RF(126) +RF(127) +RF(128) +RF(129) 
     *  +RF(130) +RF(131) +RB(109) +RB(121) 
      A8_0 = ( +RF(109) +RB(125) +RB(127) +RB(128) +RB(129) 
     *  +RB(130) +RB(131) )/DEN
      A8_3 = ( +RB(126) )/DEN
      A8_5 = ( +RF(121) )/DEN
C     C6H5
      DEN = +RF(176) +RF(178) +RF(182) +RF(183) +RF(184) 
     *  +RF(185) +RF(186) +RF(187) +RB(157) +RB(179) 
      A9_0 = ( +RF(157) +RB(176) +RB(178) +RF(179) +RB(182) 
     *  +RB(183) +RB(184) +RB(185) +RB(186) +RB(187) )/DEN
C
      A3_0 = A3_0 + A3_6*A6_0
      DEN = 1 -A3_6*A6_3
      A3_0 = A3_0/DEN
      A3_1 = A3_1/DEN
      A3_5 = A3_5/DEN
      A3_8 = A3_8/DEN
      A3_0 = A3_0 + A3_8*A8_0
      A3_5 = A3_5 + A3_8*A8_5
      DEN = 1 -A3_8*A8_3
      A3_0 = A3_0/DEN
      A3_1 = A3_1/DEN
      A3_5 = A3_5/DEN
      A5_0 = A5_0 + A5_8*A8_0
      A5_3 = A5_3 + A5_8*A8_3
      DEN = 1 -A5_8*A8_5
      A5_0 = A5_0/DEN
      A5_3 = A5_3/DEN
      A3_0 = A3_0 + A3_5*A5_0
      DEN = 1 -A3_5*A5_3
      A3_0 = A3_0/DEN
      A3_1 = A3_1/DEN
      A1_0 = A1_0 + A1_5*A5_0
      A1_3 = A1_3 + A1_5*A5_3
      A7_0 = A7_0 + A7_5*A5_0
      A7_3 = A7_5*A5_3
      A2_0 = A2_0 + A2_4*A4_0
      DEN = 1 -A2_4*A4_2
      A2_0 = A2_0/DEN
      A2_1 = A2_1/DEN
      A2_7 = A2_7/DEN
      A2_0 = A2_0 + A2_7*A7_0
      A2_3 = A2_7*A7_3
      DEN = 1 -A2_7*A7_2
      A2_0 = A2_0/DEN
      A2_3 = A2_3/DEN
      A2_1 = A2_1/DEN
      A2_0 = A2_0 + A2_1*A1_0
      A2_3 = A2_3 + A2_1*A1_3
      DEN = 1 -A2_1*A1_2
      A2_0 = A2_0/DEN
      A2_3 = A2_3/DEN
      A3_0 = A3_0 + A3_1*A1_0
      A3_2 = A3_1*A1_2
      DEN = 1 -A3_1*A1_3
      A3_0 = A3_0/DEN
      A3_2 = A3_2/DEN
      A2_0 = A2_0 + A2_3*A3_0
      DEN = 1 -A2_3*A3_2
      A2_0 = A2_0/DEN
      XQ(2) = A2_0
      XQ(3) = A3_0 +A3_2*XQ(2)
      XQ(1) = A1_0 +A1_2*XQ(2) +A1_3*XQ(3)
      XQ(7) = A7_0 +A7_2*XQ(2) +A7_3*XQ(3)
      XQ(4) = A4_0 +A4_2*XQ(2)
      XQ(5) = A5_0 +A5_3*XQ(3)
      XQ(8) = A8_0 +A8_3*XQ(3) +A8_5*XQ(5)
      XQ(6) = A6_0 +A6_3*XQ(3)
      XQ(9) = A9_0
C
      RF( 38) = RF( 38)*XQ( 3)
      RF( 39) = RF( 39)*XQ( 3)
      RF( 40) = RF( 40)*XQ( 3)
      RF( 41) = RF( 41)*XQ( 3)
      RF( 42) = RF( 42)*XQ( 3)
      RF( 43) = RF( 43)*XQ( 3)
      RF( 44) = RF( 44)*XQ( 3)
      RF( 46) = RF( 46)*XQ( 3)
      RF( 47) = RF( 47)*XQ( 1)
      RF( 48) = RF( 48)*XQ( 1)
      RB( 48) = RB( 48)*XQ( 3)
      RF( 49) = RF( 49)*XQ( 1)
      RF( 50) = RF( 50)*XQ( 1)
      RF( 51) = RF( 51)*XQ( 1)
      RB( 51) = RB( 51)*XQ( 3)
      RF( 52) = RF( 52)*XQ( 1)
      RF( 53) = RF( 53)*XQ( 1)
      RF( 55) = RF( 55)*XQ( 2)
      RB( 55) = RB( 55)*XQ( 1)
      RF( 56) = RF( 56)*XQ( 2)
      RF( 57) = RF( 57)*XQ( 2)
      RF( 58) = RF( 58)*XQ( 2)
      RF( 59) = RF( 59)*XQ( 2)
      RF( 60) = RF( 60)*XQ( 2)
      RF( 61) = RF( 61)*XQ( 2)
      RB( 61) = RB( 61)*XQ( 1)
      RF( 62) = RF( 62)*XQ( 2)
      RB( 62) = RB( 62)*XQ( 1)
      RF( 63) = RF( 63)*XQ( 2)
      RB( 63) = RB( 63)*XQ( 1)
      RF( 64) = RF( 64)*XQ( 2)
      RB( 65) = RB( 65)*XQ( 4)
      RB( 66) = RB( 66)*XQ( 3)
      RB( 67) = RB( 67)*XQ( 3)
      RB( 68) = RB( 68)*XQ( 3)
      RB( 69) = RB( 69)*XQ( 3)
      RB( 70) = RB( 70)*XQ( 3)
      RB( 73) = RB( 73)*XQ( 1)
      RB( 74) = RB( 74)*XQ( 2)
      RB( 75) = RB( 75)*XQ( 4)
      RB( 78) = RB( 78)*XQ( 4)
      RF( 80) = RF( 80)*XQ( 3)
      RB( 81) = RB( 81)*XQ( 3)
      RF( 82) = RF( 82)*XQ( 1)
      RB( 84) = RB( 84)*XQ( 6)
      RF( 85) = RF( 85)*XQ( 7)
      RF( 86) = RF( 86)*XQ( 4)
      RF( 87) = RF( 87)*XQ( 4)
      RF( 88) = RF( 88)*XQ( 4)
      RB( 88) = RB( 88)*XQ( 2)
      RF( 89) = RF( 89)*XQ( 4)
      RF( 90) = RF( 90)*XQ( 4)
      RF( 94) = RF( 94)*XQ( 1)
      RF( 95) = RF( 95)*XQ( 2)
      RF( 96) = RF( 96)*XQ( 7)
      RB( 96) = RB( 96)*XQ( 2)
      RF( 97) = RF( 97)*XQ( 7)
      RF( 98) = RF( 98)*XQ( 7)
      RB( 99) = RB( 99)*XQ( 5)
      RF(100) = RF(100)*XQ( 5)
      RB(101) = RB(101)*XQ( 1)
      RB(102) = RB(102)*XQ( 7)
      RF(105) = RF(105)*XQ( 3)
      RB(105) = RB(105)*XQ( 5)
      RF(106) = RF(106)*XQ( 1)
      RF(107) = RF(107)*XQ( 2)
      RB(109) = RB(109)*XQ( 8)
      RB(110) = RB(110)*XQ( 7)
      RB(112) = RB(112)*XQ( 7)
      RB(113) = RB(113)*XQ( 1)
      RB(114) = RB(114)*XQ( 7)
      RF(115) = RF(115)*XQ( 5)
      RF(116) = RF(116)*XQ( 5)
      RF(117) = RF(117)*XQ( 5)
      RF(118) = RF(118)*XQ( 5)
      RF(119) = RF(119)*XQ( 5)
      RF(120) = RF(120)*XQ( 5)
      RF(121) = RF(121)*XQ( 5)
      RB(121) = RB(121)*XQ( 8)
      RF(122) = RF(122)*XQ( 5)
      RB(122) = RB(122)*XQ( 3)
      RF(123) = RF(123)*XQ( 5)
      RF(124) = RF(124)*XQ( 5)
      RF(125) = RF(125)*XQ( 8)
      RF(126) = RF(126)*XQ( 8)
      RB(126) = RB(126)*XQ( 3)
      RF(127) = RF(127)*XQ( 8)
      RF(128) = RF(128)*XQ( 8)
      RF(129) = RF(129)*XQ( 8)
      RF(130) = RF(130)*XQ( 8)
      RF(131) = RF(131)*XQ( 8)
      RB(132) = RB(132)*XQ( 6)
      RB(133) = RB(133)*XQ( 5)
      RB(134) = RB(134)*XQ( 5)
      RB(135) = RB(135)*XQ( 3)
      RB(136) = RB(136)*XQ( 1)
      RB(137) = RB(137)*XQ( 5)
      RF(138) = RF(138)*XQ( 3)
      RB(138) = RB(138)*XQ( 6)
      RF(139) = RF(139)*XQ( 1)
      RF(140) = RF(140)*XQ( 2)
      RB(141) = RB(141)*XQ( 5)
      RB(142) = RB(142)*XQ( 5)
      RF(143) = RF(143)*XQ( 6)
      RF(144) = RF(144)*XQ( 6)
      RF(145) = RF(145)*XQ( 6)
      RF(146) = RF(146)*XQ( 6)
      RF(147) = RF(147)*XQ( 6)
      RB(150) = RB(150)*XQ( 6)
      RB(151) = RB(151)*XQ( 6)
      RB(152) = RB(152)*XQ( 6)
      RB(153) = RB(153)*XQ( 6)
      RB(154) = RB(154)*XQ( 3)
      RB(155) = RB(155)*XQ( 5)
      RB(157) = RB(157)*XQ( 9)
      RB(161) = RB(161)*XQ( 5)
      RF(162) = RF(162)*XQ( 3)
      RB(171) = RB(171)*XQ( 6)
      RF(176) = RF(176)*XQ( 9)
      RF(178) = RF(178)*XQ( 9)
      RB(179) = RB(179)*XQ( 9)
      RB(181) = RB(181)*XQ( 3)
      RF(182) = RF(182)*XQ( 9)
      RF(183) = RF(183)*XQ( 9)
      RF(184) = RF(184)*XQ( 9)
      RF(185) = RF(185)*XQ( 9)
      RF(186) = RF(186)*XQ( 9)
      RF(187) = RF(187)*XQ( 9)
      RB(192) = RB(192)*XQ( 3)
C
      END
