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     July 10, 2016
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(200),RB(200),RKLOW(21),C(31)
      DIMENSION XQ(10)
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(200),RB(200),RKLOW(21),C(31)
      DIMENSION XQ(10)
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.392629331350327D2
      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)/4.107330095767975D1
      C(25) = Y(25)/5.308445131778717D1
      C(26) = Y(26)/8.008703184127808D1
      C(27) = Y(27)/6.509560167789459D1
      C(28) = Y(28)/9.113384234905243D1
      C(29) = Y(29)/9.310615217685699D1
      C(30) = Y(30)/1.061252725124359D2
      C(31) = Y(31)/2.801339912414551D1
C
      SUM = 0D0
      DO K = 1, 31
         SUM = SUM + C(K)
      ENDDO
      SUM = P/(SUM*T*8.314510D7)
C
      DO K = 1, 31
         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     POSF12345
      C(1) = Y(1)/1.392629331350327D2
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     aC3H5
      C(24) = Y(24)/4.107330095767975D1
C     c-C4H5
      C(25) = Y(25)/5.308445131778717D1
C     C5H4O
      C(26) = Y(26)/8.008703184127808D1
C     C5H5
      C(27) = Y(27)/6.509560167789459D1
C     C6H5CH2
      C(28) = Y(28)/9.113384234905243D1
C     C6H5O
      C(29) = Y(29)/9.310615217685699D1
C     C6H5CHO
      C(30) = Y(30)/1.061252725124359D2
C     N2
      C(31) = Y(31)/2.801339912414551D1
C
      DO K = 1, 31
         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(40), EG(40)
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))
      EG(38) = EXP(SMH(38))
      EG(39) = EXP(SMH(39))
      EG(40) = EXP(SMH(40))
      PFAC1 = PATM / (RU*T)
      PFAC2 = PFAC1*PFAC1
      PFAC3 = PFAC2*PFAC1
C
C
C     R1: POSF12345 = 1.72699C2H4 + 0.518098C3H6 + 0.0863497iC4H8 + 0.0690798C4H81 + 0.13C6H6 + 0.37C6H5CH3 + CH3 + H
      RF(1) = EXP(6.025247481863583D1 -2.66D0*ALOGT
     * -4.438723260251439D4*TI)
      RB(1) = 0D0
C     R2: POSF12345 + 0.35H = 1.91642C2H4 + 0.01CH4 + 0.574926C3H6 + 0.0958211iC4H8 + 0.0766569C4H81 + H2 + 0.130867C6H6 + 0.372467C6H5CH3 + 0.35CH3
      RF(2) = EXP(-2.569158202235591D0 +4.76D0*ALOGT
     * -6.51615262926932D2*TI)
      RB(2) = 0D0
C     R3: POSF12345 + 0.65CH3 = 1.91642C2H4 + 1.01CH4 + 0.574926C3H6 + 0.0958211iC4H8 + 0.0766569C4H81 + 0.130867C6H6 + 0.372467C6H5CH3 + 0.65H
      RF(3) = EXP(-1.455784798271499D1 +5.95D0*ALOGT
     * -2.389474024621395D3*TI)
      RB(3) = 0D0
C     R4: POSF12345 + OH = 1.91642C2H4 + 0.01CH4 + 0.574926C3H6 + 0.0958211iC4H8 + 0.0766569C4H81 + 0.130867C6H6 + 0.372467C6H5CH3 + 0.35CH3 + H2O + 0.65H
      RF(4) = EXP(2.429279853342778D1 +1.02D0*ALOGT
     * -1.072857935408309D2*TI)
      RB(4) = 0D0
C     R5: POSF12345 + O2 = 1.91642C2H4 + 0.01CH4 + 0.574926C3H6 + 0.0958211iC4H8 + 0.0766569C4H81 + 0.130867C6H6 + 0.372467C6H5CH3 + 0.35CH3 + HO2 + 0.65H
      RF(5) = EXP(3.317228466108698D1 +6.D-2*ALOGT
     * -2.391919657626125D4*TI)
      RB(5) = 0D0
C     R6: POSF12345 + HO2 = 1.91642C2H4 + 0.01CH4 + 0.574926C3H6 + 0.0958211iC4H8 + 0.0766569C4H81 + 0.130867C6H6 + 0.372467C6H5CH3 + 0.35CH3 + H2O2 + 0.65H
      RF(6) = EXP(1.005620863955379D1 +2.94D0*ALOGT
     * -7.20103051392725D3*TI)
      RB(6) = 0D0
C     R7: POSF12345 + O = 1.91642C2H4 + 0.01CH4 + 0.574926C3H6 + 0.0958211iC4H8 + 0.0766569C4H81 + 0.130867C6H6 + 0.372467C6H5CH3 + 0.35CH3 + OH + 0.65H
      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(4.547615992139523D1 -1.25D0*ALOGT)
      EQK = EG(7)/EG(20)/EG(20)/PFAC1
      RB(13) = RF(13) / MAX(EQK, SMALL)
C     R14: OH + H = H2O
      RF(14) = EXP(5.213847658679322D1 -2.D0*ALOGT)
      EQK = EG(18)/EG(16)/EG(20)/PFAC1
      RB(14) = RF(14) / MAX(EQK, SMALL)
C     R15: O + H = OH
      RF(15) = EXP(4.369021565896671D1 -1.D0*ALOGT)
      EQK = EG(16)/EG(15)/EG(20)/PFAC1
      RB(15) = RF(15) / MAX(EQK, SMALL)
C     R16: 2O = O2
      RF(16) = EXP(3.932626813769273D1 -1.D0*ALOGT)
      EQK = EG(14)/EG(15)/EG(15)/PFAC1
      RB(16) = RF(16) / MAX(EQK, SMALL)
C     R17: O2 + H = HO2
      RF(17) = EXP(2.926339399964515D1 +4.4D-1*ALOGT)
      EQK = EG(17)/EG(14)/EG(20)/PFAC1
      RB(17) = RF(17) / MAX(EQK, SMALL)
C     R18: H2 + O2 = HO2 + H
      RF(18) = EXP(1.329058600981878D1 +2.433D0*ALOGT
     * -2.692309815207098D4*TI)
      EQK = EG(17)*EG(20)/EG(7)/EG(14)
      RB(18) = RF(18) / MAX(EQK, SMALL)
C     R19: 2OH = H2O2
      RF(19) = EXP(3.234055131724088D1 -3.7D-1*ALOGT)
      EQK = EG(19)/EG(16)/EG(16)/PFAC1
      RB(19) = RF(19) / MAX(EQK, SMALL)
C     R20: HO2 + H = 2OH
      RF(20) = EXP(3.194650722679419D1 -1.484489169537763D2*TI)
      EQK = EG(16)*EG(16)/EG(17)/EG(20)
      RB(20) = RF(20) / MAX(EQK, SMALL)
C     R21: O + HO2 = O2 + OH
      RF(21) = 4.D13
      EQK = EG(14)*EG(16)/EG(15)/EG(17)
      RB(21) = RF(21) / MAX(EQK, SMALL)
C     R22: 2HO2 = O2 + H2O2
      RF(22) = EXP(2.559080028740199D1 +8.20243168253069D2*TI)
      EQK = EG(14)*EG(19)/EG(17)/EG(17)
      RB(22) = RF(22) / MAX(EQK, SMALL)
C     R23: 2HO2 = O2 + H2O2
      RF(23) = EXP(3.353310785188531D1 -6.038600011679036D3*TI)
      EQK = EG(14)*EG(19)/EG(17)/EG(17)
      RB(23) = RF(23) / MAX(EQK, SMALL)
C     R24: OH + HO2 = O2 + H2O
      RF(24) = EXP(4.17901213782829D1 -1.76D0*ALOGT
     * -3.019300005839518D1*TI)
      EQK = EG(14)*EG(18)/EG(16)/EG(17)
      RB(24) = RF(24) / MAX(EQK, SMALL)
C     R25: OH + HO2 = O2 + H2O
      RF(25) = EXP(1.62861784418104D2 -1.672D1*ALOGT
     * -1.655582836535336D4*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(2.855130386907224D1 +2.D0*ALOGT
     * -2.012866670559679D4*TI)
      EQK = EG(14)*EG(18)/EG(16)/EG(17)
      RB(26) = RF(26) / MAX(EQK, SMALL)
C     R27: H2O2 + H = H2 + HO2
      RF(27) = EXP(1.561556883000702D1 +2.D0*ALOGT
     * -2.616726671727582D3*TI)
      EQK = EG(7)*EG(17)/EG(19)/EG(20)
      RB(27) = RF(27) / MAX(EQK, SMALL)
C     R28: H2O2 + H = OH + H2O
      RF(28) = EXP(3.081323295642516D1 -1.997770170530481D3*TI)
      EQK = EG(16)*EG(18)/EG(19)/EG(20)
      RB(28) = RF(28) / MAX(EQK, SMALL)
C     R29: O + H2O2 = OH + HO2
      RF(29) = EXP(1.608039378377431D1 +2.D0*ALOGT
     * -1.997770170530481D3*TI)
      EQK = EG(16)*EG(17)/EG(15)/EG(19)
      RB(29) = RF(29) / MAX(EQK, SMALL)
C     R30: OH + H2O2 = HO2 + H2O
      RF(30) = EXP(9.538806728516803D1 -7.D0*ALOGT
     * -1.892094670326098D4*TI)
      EQK = EG(17)*EG(18)/EG(16)/EG(19)
      RB(30) = RF(30) / MAX(EQK, SMALL)
C     R31: CO + O = CO2
      RF(31) = EXP(2.333480513766778D1 -1.199668535653569D3*TI)
      EQK = EG(26)/EG(9)/EG(15)/PFAC1
      RB(31) = RF(31) / MAX(EQK, SMALL)
C     R32: CO + OH = H + CO2
      RF(32) = EXP(1.116280045189523D1 +2.053D0*ALOGT
     * +1.789790721794902D2*TI)
      EQK = EG(20)*EG(26)/EG(9)/EG(16)
      RB(32) = RF(32) / MAX(EQK, SMALL)
C     R33: CO + OH = H + CO2
      RF(33) = EXP(2.938143762162222D1 -6.64D-1*ALOGT
     * -1.669823868229545D2*TI)
      EQK = EG(20)*EG(26)/EG(9)/EG(16)
      RB(33) = RF(33) / MAX(EQK, SMALL)
C     R34: CO + O2 = O + CO2
      RF(34) = EXP(2.774345654525834D1 -2.400343504642417D4*TI)
      EQK = EG(15)*EG(26)/EG(9)/EG(14)
      RB(34) = RF(34) / MAX(EQK, SMALL)
C     R35: CO + HO2 = OH + CO2
      RF(35) = EXP(1.196400108433044D1 +2.18D0*ALOGT
     * -9.0290204129627D3*TI)
      EQK = EG(16)*EG(26)/EG(9)/EG(17)
      RB(35) = RF(35) / MAX(EQK, SMALL)
C     R36: H + HCO = H2 + CO
      RF(36) = 1.2D14
      EQK = EG(7)*EG(9)/EG(20)/EG(23)
      RB(36) = RF(36) / MAX(EQK, SMALL)
C     R37: O + HCO = CO + OH
      RF(37) = 3.D13
      EQK = EG(9)*EG(16)/EG(15)/EG(23)
      RB(37) = RF(37) / MAX(EQK, SMALL)
C     R38: O + HCO = H + CO2
      RF(38) = 3.D13
      EQK = EG(20)*EG(26)/EG(15)/EG(23)
      RB(38) = RF(38) / MAX(EQK, SMALL)
C     R39: OH + HCO = CO + H2O
      RF(39) = 3.02D13
      EQK = EG(9)*EG(18)/EG(16)/EG(23)
      RB(39) = RF(39) / MAX(EQK, SMALL)
C     R40: HCO = CO + H
      RF(40) = EXP(3.976988501176528D1 -1.D0*ALOGT
     * -8.554683349878635D3*TI)
      EQK = EG(9)*EG(20)/EG(23)*PFAC1
      RB(40) = RF(40) / MAX(EQK, SMALL)
C     R41: HCO = CO + H
      RF(41) = EXP(4.225479166155327D1 -1.D0*ALOGT
     * -8.554683349878635D3*TI)
      EQK = EG(9)*EG(20)/EG(23)*PFAC1
      RB(41) = RF(41) / MAX(EQK, SMALL)
C     R42: O2 + HCO = CO + HO2
      RF(42) = EXP(2.321150027682708D1 +8.070000000000001D-1*ALOGT
     * +3.658385173742216D2*TI)
      EQK = EG(9)*EG(17)/EG(14)/EG(23)
      RB(42) = RF(42) / MAX(EQK, SMALL)
C     R43: H2 + CO = CH2O
      RF(43) = EXP(1.757671067365784D1 +1.5D0*ALOGT
     * -4.005604674413761D4*TI)
      EQK = EG(24)/EG(7)/EG(9)/PFAC1
      RB(43) = RF(43) / MAX(EQK, SMALL)
C     R44: H + HCO = CH2O
      RF(44) = EXP(2.77171988121696D1 +4.8D-1*ALOGT
     * +1.308363335863791D2*TI)
      EQK = EG(24)/EG(20)/EG(23)/PFAC1
      RB(44) = RF(44) / MAX(EQK, SMALL)
C     R45: H + CH2 = CH3
      RF(45) = EXP(3.775765221977888D1 -8.D-1*ALOGT)
      EQK = EG(13)/EG(20)/EG(21)/PFAC1
      RB(45) = RF(45) / MAX(EQK, SMALL)
C     R46: O + CH2 = H + HCO
      RF(46) = 8.D13
      EQK = EG(20)*EG(23)/EG(15)/EG(21)
      RB(46) = RF(46) / MAX(EQK, SMALL)
C     R47: OH + CH2 = H + CH2O
      RF(47) = 2.D13
      EQK = EG(20)*EG(24)/EG(16)/EG(21)
      RB(47) = RF(47) / MAX(EQK, SMALL)
C     R48: H2 + CH2 = CH3 + H
      RF(48) = EXP(1.312236337740433D1 +2.D0*ALOGT
     * -3.638256507036619D3*TI)
      EQK = EG(13)*EG(20)/EG(7)/EG(21)
      RB(48) = RF(48) / MAX(EQK, SMALL)
C     R49: O2 + CH2 = OH + HCO
      RF(49) = EXP(2.999187511704657D1 -7.548250014598796D2*TI)
      EQK = EG(16)*EG(23)/EG(14)/EG(21)
      RB(49) = RF(49) / MAX(EQK, SMALL)
C     R50: O2 + CH2 = 2H + CO2
      RF(50) = EXP(2.860180003308677D1 -7.548250014598796D2*TI)
      EQK = EG(20)*EG(20)*EG(26)/EG(14)/EG(21)*PFAC1
      RB(50) = RF(50) / MAX(EQK, SMALL)
C     R51: CO + CH2 = CH2CO
      RF(51) = EXP(2.74203000846129D1 +5.D-1*ALOGT
     * -2.269507171056038D3*TI)
      EQK = EG(30)/EG(9)/EG(21)/PFAC1
      RB(51) = RF(51) / MAX(EQK, SMALL)
C     R52: CH2* = CH2
      RF(52) = EXP(3.033907131703076D1 -3.019300005839518D2*TI)
      EQK = EG(21)/EG(22)
      RB(52) = RF(52) / MAX(EQK, SMALL)
C     R53: O + CH2* = H2 + CO
      RF(53) = 1.5D13
      EQK = EG(7)*EG(9)/EG(15)/EG(22)
      RB(53) = RF(53) / MAX(EQK, SMALL)
C     R54: OH + CH2* = H + CH2O
      RF(54) = 3.D13
      EQK = EG(20)*EG(24)/EG(16)/EG(22)
      RB(54) = RF(54) / MAX(EQK, SMALL)
C     R55: H2 + CH2* = CH3 + H
      RF(55) = 7.D13
      EQK = EG(13)*EG(20)/EG(7)/EG(22)
      RB(55) = RF(55) / MAX(EQK, SMALL)
C     R56: O2 + CH2* = CO + OH + H
      RF(56) = 2.8D13
      EQK = EG(9)*EG(16)*EG(20)/EG(14)/EG(22)*PFAC1
      RB(56) = RF(56) / MAX(EQK, SMALL)
C     R57: O2 + CH2* = CO + H2O
      RF(57) = 1.2D13
      EQK = EG(9)*EG(18)/EG(14)/EG(22)
      RB(57) = RF(57) / MAX(EQK, SMALL)
C     R58: CH2* = CH2
      RF(58) = 3.D13
      EQK = EG(21)/EG(22)
      RB(58) = RF(58) / MAX(EQK, SMALL)
C     R59: CH2* = CH2
      RF(59) = 9.D12
      EQK = EG(21)/EG(22)
      RB(59) = RF(59) / MAX(EQK, SMALL)
C     R60: CH2* = CH2
      RF(60) = 7.D12
      EQK = EG(21)/EG(22)
      RB(60) = RF(60) / MAX(EQK, SMALL)
C     R61: CH2* + CO2 = CO + CH2O
      RF(61) = 1.4D13
      EQK = EG(9)*EG(24)/EG(22)/EG(26)
      RB(61) = RF(61) / MAX(EQK, SMALL)
C     R62: H + CH2O = CH3O
      RF(62) = EXP(2.701483497650473D1 +4.54D-1*ALOGT
     * -1.308363335863791D3*TI)
      EQK = EG(25)/EG(20)/EG(24)/PFAC1
      RB(62) = RF(62) / MAX(EQK, SMALL)
C     R63: H + CH2O = H2 + HCO
      RF(63) = EXP(2.385876005287556D1 +1.05D0*ALOGT
     * -1.648034586520737D3*TI)
      EQK = EG(7)*EG(23)/EG(20)/EG(24)
      RB(63) = RF(63) / MAX(EQK, SMALL)
C     R64: O + CH2O = OH + HCO
      RF(64) = EXP(3.129458276205819D1 -1.781387003445316D3*TI)
      EQK = EG(16)*EG(23)/EG(15)/EG(24)
      RB(64) = RF(64) / MAX(EQK, SMALL)
C     R65: OH + CH2O = H2O + HCO
      RF(65) = EXP(2.195582609812426D1 +1.18D0*ALOGT
     * +2.249378504350441D2*TI)
      EQK = EG(18)*EG(23)/EG(16)/EG(24)
      RB(65) = RF(65) / MAX(EQK, SMALL)
C     R66: O2 + CH2O = HO2 + HCO
      RF(66) = EXP(3.223619130191664D1 -2.012866670559679D4*TI)
      EQK = EG(17)*EG(23)/EG(14)/EG(24)
      RB(66) = RF(66) / MAX(EQK, SMALL)
C     R67: HO2 + CH2O = H2O2 + HCO
      RF(67) = EXP(2.763102111592855D1 -4.025733341119358D3*TI)
      EQK = EG(19)*EG(23)/EG(17)/EG(24)
      RB(67) = RF(67) / MAX(EQK, SMALL)
C     R68: CH3 + H = CH4
      RF(68) = EXP(3.708037838837523D1 -6.3D-1*ALOGT
     * -1.927319837060892D2*TI)
      EQK = EG(3)/EG(13)/EG(20)/PFAC1
      RB(68) = RF(68) / MAX(EQK, SMALL)
C     R69: CH3 + O = H + CH2O
      RF(69) = 8.43D13
      EQK = EG(20)*EG(24)/EG(13)/EG(15)
      RB(69) = RF(69) / MAX(EQK, SMALL)
C     R70: CH3 + OH = H2O + CH2
      RF(70) = EXP(1.784086224869942D1 +1.6D0*ALOGT
     * -2.727434338608365D3*TI)
      EQK = EG(18)*EG(21)/EG(13)/EG(16)
      RB(70) = RF(70) / MAX(EQK, SMALL)
C     R71: CH3 + OH = H2O + CH2*
      RF(71) = 2.501D13
      EQK = EG(18)*EG(22)/EG(13)/EG(16)
      RB(71) = RF(71) / MAX(EQK, SMALL)
C     R72: CH3 + O2 = O + CH3O
      RF(72) = EXP(3.105950935782661D1 -1.449264002802969D4*TI)
      EQK = EG(15)*EG(25)/EG(13)/EG(14)
      RB(72) = RF(72) / MAX(EQK, SMALL)
C     R73: CH3 + O2 = OH + CH2O
      RF(73) = EXP(2.430678477540252D1 -4.498757008700882D3*TI)
      EQK = EG(16)*EG(24)/EG(13)/EG(14)
      RB(73) = RF(73) / MAX(EQK, SMALL)
C     R74: CH3 + HO2 = CH4 + O2
      RF(74) = 1.D12
      EQK = EG(3)*EG(14)/EG(13)/EG(17)
      RB(74) = RF(74) / MAX(EQK, SMALL)
C     R75: CH3 + HO2 = OH + CH3O
      RF(75) = 1.34D13
      EQK = EG(16)*EG(25)/EG(13)/EG(17)
      RB(75) = RF(75) / MAX(EQK, SMALL)
C     R76: CH3 + H2O2 = CH4 + HO2
      RF(76) = EXP(1.010642839653282D1 +2.47D0*ALOGT
     * -2.606662338374784D3*TI)
      EQK = EG(3)*EG(17)/EG(13)/EG(19)
      RB(76) = RF(76) / MAX(EQK, SMALL)
C     R77: CH3 + HCO = CH4 + CO
      RF(77) = 8.48D12
      EQK = EG(3)*EG(9)/EG(13)/EG(23)
      RB(77) = RF(77) / MAX(EQK, SMALL)
C     R78: CH3 + CH2O = CH4 + HCO
      RF(78) = EXP(8.107720061910534D0 +2.81D0*ALOGT
     * -2.94884967236993D3*TI)
      EQK = EG(3)*EG(23)/EG(13)/EG(24)
      RB(78) = RF(78) / MAX(EQK, SMALL)
C     R79: CH3 + CH2 = C2H4 + H
      RF(79) = 4.D13
      EQK = EG(2)*EG(20)/EG(13)/EG(21)
      RB(79) = RF(79) / MAX(EQK, SMALL)
C     R80: 2CH3 = C2H6
      RF(80) = EXP(3.759277757658865D1 -9.7D-1*ALOGT
     * -3.119943339367502D2*TI)
      EQK = EG(8)/EG(13)/EG(13)/PFAC1
      RB(80) = RF(80) / MAX(EQK, SMALL)
C     R81: 2CH3 = H + C2H5
      RF(81) = EXP(2.923845702569198D1 +1.D-1*ALOGT
     * -5.334096676983148D3*TI)
      EQK = EG(20)*EG(28)/EG(13)/EG(13)
      RB(81) = RF(81) / MAX(EQK, SMALL)
C     R82: CH3 + HCCO = C2H4 + CO
      RF(82) = 5.D13
      EQK = EG(2)*EG(9)/EG(13)/EG(29)
      RB(82) = RF(82) / MAX(EQK, SMALL)
C     R83: H + CH3O = H2 + CH2O
      RF(83) = 2.D13
      EQK = EG(7)*EG(24)/EG(20)/EG(25)
      RB(83) = RF(83) / MAX(EQK, SMALL)
C     R84: H + CH3O = CH3 + OH
      RF(84) = 3.2D13
      EQK = EG(13)*EG(16)/EG(20)/EG(25)
      RB(84) = RF(84) / MAX(EQK, SMALL)
C     R85: H + CH3O = H2O + CH2*
      RF(85) = 1.6D13
      EQK = EG(18)*EG(22)/EG(20)/EG(25)
      RB(85) = RF(85) / MAX(EQK, SMALL)
C     R86: O + CH3O = OH + CH2O
      RF(86) = 1.D13
      EQK = EG(16)*EG(24)/EG(15)/EG(25)
      RB(86) = RF(86) / MAX(EQK, SMALL)
C     R87: OH + CH3O = H2O + CH2O
      RF(87) = 5.D12
      EQK = EG(18)*EG(24)/EG(16)/EG(25)
      RB(87) = RF(87) / MAX(EQK, SMALL)
C     R88: O2 + CH3O = HO2 + CH2O
      RF(88) = EXP(-2.847965319932889D1 +7.6D0*ALOGT
     * +1.776354836768916D3*TI)
      EQK = EG(17)*EG(24)/EG(14)/EG(25)
      RB(88) = RF(88) / MAX(EQK, SMALL)
C     R89: CH4 + H = H2 + CH3
      RF(89) = EXP(2.030775039298474D1 +1.62D0*ALOGT
     * -5.454868677216729D3*TI)
      EQK = EG(7)*EG(13)/EG(3)/EG(20)
      RB(89) = RF(89) / MAX(EQK, SMALL)
C     R90: CH4 + O = CH3 + OH
      RF(90) = EXP(2.074306846424259D1 +1.5D0*ALOGT
     * -4.327663341703309D3*TI)
      EQK = EG(13)*EG(16)/EG(3)/EG(15)
      RB(90) = RF(90) / MAX(EQK, SMALL)
C     R91: CH4 + OH = CH3 + H2O
      RF(91) = EXP(1.842068074395237D1 +1.6D0*ALOGT
     * -1.570036003036549D3*TI)
      EQK = EG(13)*EG(18)/EG(3)/EG(16)
      RB(91) = RF(91) / MAX(EQK, SMALL)
C     R92: CH4 + CH2* = 2CH3
      RF(92) = EXP(3.040360983816833D1 +2.868335005547542D2*TI)
      EQK = EG(13)*EG(13)/EG(3)/EG(22)
      RB(92) = RF(92) / MAX(EQK, SMALL)
C     R93: H + HCCO = CO + CH2*
      RF(93) = 1.D14
      EQK = EG(9)*EG(22)/EG(20)/EG(29)
      RB(93) = RF(93) / MAX(EQK, SMALL)
C     R94: O + HCCO = 2CO + H
      RF(94) = 1.D14
      EQK = EG(9)*EG(9)*EG(20)/EG(15)/EG(29)*PFAC1
      RB(94) = RF(94) / MAX(EQK, SMALL)
C     R95: O2 + HCCO = 2CO + OH
      RF(95) = EXP(2.810102474517428D1 -4.297470341644914D2*TI)
      EQK = EG(9)*EG(9)*EG(16)/EG(14)/EG(29)*PFAC1
      RB(95) = RF(95) / MAX(EQK, SMALL)
C     R96: C2H3 = C2H2 + H
      RF(96) = EXP(1.97713479274291D1 +1.62D0*ALOGT
     * -1.864327174605727D4*TI)
      EQK = EG(11)*EG(20)/EG(27)*PFAC1
      RB(96) = RF(96) / MAX(EQK, SMALL)
C     R97: C2H2 + O = CO + CH2
      RF(97) = EXP(1.522160754638034D1 +2.D0*ALOGT
     * -9.561116685158474D2*TI)
      EQK = EG(9)*EG(21)/EG(11)/EG(15)
      RB(97) = RF(97) / MAX(EQK, SMALL)
C     R98: C2H2 + O = H + HCCO
      RF(98) = EXP(1.660790190750023D1 +2.D0*ALOGT
     * -9.561116685158474D2*TI)
      EQK = EG(20)*EG(29)/EG(11)/EG(15)
      RB(98) = RF(98) / MAX(EQK, SMALL)
C     R99: C2H2 + OH = H + CH2CO
      RF(99) = EXP(-8.431015495175185D0 +4.5D0*ALOGT
     * +5.032166676399197D2*TI)
      EQK = EG(20)*EG(30)/EG(11)/EG(16)
      RB(99) = RF(99) / MAX(EQK, SMALL)
C     R100: C2H2 + OH = CO + CH3
      RF(100) = EXP(-7.635493904311701D0 +4.D0*ALOGT
     * +1.006433335279839D3*TI)
      EQK = EG(9)*EG(13)/EG(11)/EG(16)
      RB(100) = RF(100) / MAX(EQK, SMALL)
C     R101: C2H2 + HCO = CO + C2H3
      RF(101) = EXP(1.611809565095832D1 +2.D0*ALOGT
     * -3.019300005839518D3*TI)
      EQK = EG(9)*EG(27)/EG(11)/EG(23)
      RB(101) = RF(101) / MAX(EQK, SMALL)
C     R102: C2H2 + CH3 = aC3H5
      RF(102) = EXP(1.230228267232072D2 -1.282D1*ALOGT
     * -1.797993153477433D4*TI)
      EQK = EG(32)/EG(11)/EG(13)/PFAC1
      RB(102) = RF(102) / MAX(EQK, SMALL)
C     R103: H + CH2CO = CH2CHO
      RF(103) = EXP(3.343011377038908D1 -6.D-2*ALOGT
     * -4.277341674939317D3*TI)
      EQK = EG(31)/EG(20)/EG(30)/PFAC1
      RB(103) = RF(103) / MAX(EQK, SMALL)
C     R104: H + CH2CO = H2 + HCCO
      RF(104) = EXP(3.154304412135669D1 -4.025733341119358D3*TI)
      EQK = EG(7)*EG(29)/EG(20)/EG(30)
      RB(104) = RF(104) / MAX(EQK, SMALL)
C     R105: H + CH2CO = CO + CH3
      RF(105) = EXP(2.112873094505457D1 +1.43D0*ALOGT
     * -1.353652835951384D3*TI)
      EQK = EG(9)*EG(13)/EG(20)/EG(30)
      RB(105) = RF(105) / MAX(EQK, SMALL)
C     R106: O + CH2CO = OH + HCCO
      RF(106) = EXP(2.993360620892259D1 -4.025733341119358D3*TI)
      EQK = EG(16)*EG(29)/EG(15)/EG(30)
      RB(106) = RF(106) / MAX(EQK, SMALL)
C     R107: O + CH2CO = CH2 + CO2
      RF(107) = EXP(2.819063690386397D1 -6.793425013138916D2*TI)
      EQK = EG(21)*EG(26)/EG(15)/EG(30)
      RB(107) = RF(107) / MAX(EQK, SMALL)
C     R108: OH + CH2CO = H2O + HCCO
      RF(108) = EXP(2.964592413647081D1 -1.006433335279839D3*TI)
      EQK = EG(18)*EG(29)/EG(16)/EG(30)
      RB(108) = RF(108) / MAX(EQK, SMALL)
C     R109: H + C2H3 = C2H4
      RF(109) = EXP(2.943602581190662D1 +2.7D-1*ALOGT
     * -1.409006669391775D2*TI)
      EQK = EG(2)/EG(20)/EG(27)/PFAC1
      RB(109) = RF(109) / MAX(EQK, SMALL)
C     R110: H + C2H3 = H2 + C2H2
      RF(110) = 9.D13
      EQK = EG(7)*EG(11)/EG(20)/EG(27)
      RB(110) = RF(110) / MAX(EQK, SMALL)
C     R111: O + C2H3 = H + CH2CO
      RF(111) = 4.8D13
      EQK = EG(20)*EG(30)/EG(15)/EG(27)
      RB(111) = RF(111) / MAX(EQK, SMALL)
C     R112: O + C2H3 = CO + CH3
      RF(112) = 4.8D13
      EQK = EG(9)*EG(13)/EG(15)/EG(27)
      RB(112) = RF(112) / MAX(EQK, SMALL)
C     R113: OH + C2H3 = C2H2 + H2O
      RF(113) = 3.011D13
      EQK = EG(11)*EG(18)/EG(16)/EG(27)
      RB(113) = RF(113) / MAX(EQK, SMALL)
C     R114: O2 + C2H3 = C2H2 + HO2
      RF(114) = EXP(1.410818017192709D1 +1.61D0*ALOGT
     * +1.929332703731452D2*TI)
      EQK = EG(11)*EG(17)/EG(14)/EG(27)
      RB(114) = RF(114) / MAX(EQK, SMALL)
C     R115: O2 + C2H3 = O + CH2CHO
      RF(115) = EXP(2.642704831160261D1 +2.9D-1*ALOGT
     * -5.535383344039117D0*TI)
      EQK = EG(15)*EG(31)/EG(14)/EG(27)
      RB(115) = RF(115) / MAX(EQK, SMALL)
C     R116: O2 + C2H3 = HCO + CH2O
      RF(116) = EXP(3.836741779139978D1 -1.39D0*ALOGT
     * -5.082488343163189D2*TI)
      EQK = EG(23)*EG(24)/EG(14)/EG(27)
      RB(116) = RF(116) / MAX(EQK, SMALL)
C     R117: CH3 + C2H3 = C3H6
      RF(117) = 2.5D13
      EQK = EG(4)/EG(13)/EG(27)/PFAC1
      RB(117) = RF(117) / MAX(EQK, SMALL)
C     R118: CH3 + C2H3 = H + aC3H5
      RF(118) = EXP(5.566750733996526D1 -2.83D0*ALOGT
     * -9.368887918120025D3*TI)
      EQK = EG(20)*EG(32)/EG(13)/EG(27)
      RB(118) = RF(118) / MAX(EQK, SMALL)
C     R119: CH2CHO = CO + CH3
      RF(119) = EXP(9.646011254645141D1 -9.147D0*ALOGT
     * -2.360086171231223D4*TI)
      EQK = EG(9)*EG(13)/EG(31)*PFAC1
      RB(119) = RF(119) / MAX(EQK, SMALL)
C     R120: H + CH2CHO = CH3 + HCO
      RF(120) = 9.D13
      EQK = EG(13)*EG(23)/EG(20)/EG(31)
      RB(120) = RF(120) / MAX(EQK, SMALL)
C     R121: H + CH2CHO = H2 + CH2CO
      RF(121) = EXP(3.062675338948254D1 -2.012866670559679D3*TI)
      EQK = EG(7)*EG(30)/EG(20)/EG(31)
      RB(121) = RF(121) / MAX(EQK, SMALL)
C     R122: O + CH2CHO = OH + CH2CO
      RF(122) = EXP(3.062675338948254D1 -2.012866670559679D3*TI)
      EQK = EG(16)*EG(30)/EG(15)/EG(31)
      RB(122) = RF(122) / MAX(EQK, SMALL)
C     R123: OH + CH2CHO = H2O + CH2CO
      RF(123) = EXP(2.993360620892259D1 -1.006433335279839D3*TI)
      EQK = EG(18)*EG(30)/EG(16)/EG(31)
      RB(123) = RF(123) / MAX(EQK, SMALL)
C     R124: O2 + CH2CHO = HO2 + CH2CO
      RF(124) = 1.4D11
      EQK = EG(17)*EG(30)/EG(14)/EG(31)
      RB(124) = RF(124) / MAX(EQK, SMALL)
C     R125: C2H4 + H = C2H5
      RF(125) = EXP(2.103807657678644D1 +1.463D0*ALOGT
     * -6.818585846520912D2*TI)
      EQK = EG(28)/EG(2)/EG(20)/PFAC1
      RB(125) = RF(125) / MAX(EQK, SMALL)
C     R126: C2H4 + H = H2 + C2H3
      RF(126) = EXP(1.774143646856141D1 +1.9D0*ALOGT
     * -6.51665584593696D3*TI)
      EQK = EG(7)*EG(27)/EG(2)/EG(20)
      RB(126) = RF(126) / MAX(EQK, SMALL)
C     R127: C2H4 + O = OH + C2H3
      RF(127) = EXP(1.653020530178515D1 +1.9D0*ALOGT
     * -1.8820303369733D3*TI)
      EQK = EG(16)*EG(27)/EG(2)/EG(15)
      RB(127) = RF(127) / MAX(EQK, SMALL)
C     R128: C2H4 + O = CH3 + HCO
      RF(128) = EXP(1.677042083699801D1 +1.83D0*ALOGT
     * -1.107076668807823D2*TI)
      EQK = EG(13)*EG(23)/EG(2)/EG(15)
      RB(128) = RF(128) / MAX(EQK, SMALL)
C     R129: C2H4 + O = CH2 + CH2O
      RF(129) = EXP(1.285839783156986D1 +1.83D0*ALOGT
     * -1.107076668807823D2*TI)
      EQK = EG(21)*EG(24)/EG(2)/EG(15)
      RB(129) = RF(129) / MAX(EQK, SMALL)
C     R130: C2H4 + OH = H2O + C2H3
      RF(130) = EXP(1.509644440342634D1 +2.D0*ALOGT
     * -1.258041669099799D3*TI)
      EQK = EG(18)*EG(27)/EG(2)/EG(16)
      RB(130) = RF(130) / MAX(EQK, SMALL)
C     R131: C2H4 + HCO = CO + C2H5
      RF(131) = EXP(1.611809565095832D1 +2.D0*ALOGT
     * -4.025733341119358D3*TI)
      EQK = EG(9)*EG(28)/EG(2)/EG(23)
      RB(131) = RF(131) / MAX(EQK, SMALL)
C     R132: C2H4 + CH2 = H + aC3H5
      RF(132) = EXP(3.062675338948254D1 -3.019300005839518D3*TI)
      EQK = EG(20)*EG(32)/EG(2)/EG(21)
      RB(132) = RF(132) / MAX(EQK, SMALL)
C     R133: C2H4 + CH2* = H + aC3H5
      RF(133) = 5.D13
      EQK = EG(20)*EG(32)/EG(2)/EG(22)
      RB(133) = RF(133) / MAX(EQK, SMALL)
C     R134: C2H4 + CH3 = CH4 + C2H3
      RF(134) = EXP(1.233270529646354D1 +2.D0*ALOGT
     * -4.629593342287261D3*TI)
      EQK = EG(3)*EG(27)/EG(2)/EG(13)
      RB(134) = RF(134) / MAX(EQK, SMALL)
C     R135: C2H4 + O2 = HO2 + C2H3
      RF(135) = EXP(3.137344133697052D1 -3.059557339250712D4*TI)
      EQK = EG(17)*EG(27)/EG(2)/EG(14)
      RB(135) = RF(135) / MAX(EQK, SMALL)
C     R136: H + C2H5 = C2H6
      RF(136) = EXP(4.079452643666405D1 -9.9D-1*ALOGT
     * -7.95082334871073D2*TI)
      EQK = EG(8)/EG(20)/EG(28)/PFAC1
      RB(136) = RF(136) / MAX(EQK, SMALL)
C     R137: O + C2H5 = CH3 + CH2O
      RF(137) = 1.604D13
      EQK = EG(13)*EG(24)/EG(15)/EG(28)
      RB(137) = RF(137) / MAX(EQK, SMALL)
C     R138: O2 + C2H5 = C2H4 + HO2
      RF(138) = 2.D10
      EQK = EG(2)*EG(17)/EG(14)/EG(28)
      RB(138) = RF(138) / MAX(EQK, SMALL)
C     R139: HO2 + C2H5 = CH3 + OH + CH2O
      RF(139) = 2.4D13
      EQK = EG(13)*EG(16)*EG(24)/EG(17)/EG(28)*PFAC1
      RB(139) = RF(139) / MAX(EQK, SMALL)
C     R140: C2H3 + C2H5 = C4H81
      RF(140) = 1.5D13
      EQK = EG(6)/EG(27)/EG(28)/PFAC1
      RB(140) = RF(140) / MAX(EQK, SMALL)
C     R141: C2H6 + H = H2 + C2H5
      RF(141) = EXP(1.856044268632752D1 +1.9D0*ALOGT
     * -3.789221507328595D3*TI)
      EQK = EG(7)*EG(28)/EG(8)/EG(20)
      RB(141) = RF(141) / MAX(EQK, SMALL)
C     R142: C2H6 + O = OH + C2H5
      RF(142) = EXP(1.831309553327243D1 +1.92D0*ALOGT
     * -2.863302838871143D3*TI)
      EQK = EG(16)*EG(28)/EG(8)/EG(15)
      RB(142) = RF(142) / MAX(EQK, SMALL)
C     R143: C2H6 + OH = H2O + C2H5
      RF(143) = EXP(1.507963728510996D1 +2.12D0*ALOGT
     * -4.377985008467301D2*TI)
      EQK = EG(18)*EG(28)/EG(8)/EG(16)
      RB(143) = RF(143) / MAX(EQK, SMALL)
C     R144: H + aC3H5 = C3H6
      RF(144) = 2.D14
      EQK = EG(4)/EG(20)/EG(32)/PFAC1
      RB(144) = RF(144) / MAX(EQK, SMALL)
C     R145: HO2 + aC3H5 = C3H6 + O2
      RF(145) = 2.66D12
      EQK = EG(4)*EG(14)/EG(17)/EG(32)
      RB(145) = RF(145) / MAX(EQK, SMALL)
C     R146: HO2 + aC3H5 = OH + CH2O + C2H3
      RF(146) = 6.6D12
      EQK = EG(16)*EG(24)*EG(27)/EG(17)/EG(32)*PFAC1
      RB(146) = RF(146) / MAX(EQK, SMALL)
C     R147: HCO + aC3H5 = C3H6 + CO
      RF(147) = 6.D13
      EQK = EG(4)*EG(9)/EG(23)/EG(32)
      RB(147) = RF(147) / MAX(EQK, SMALL)
C     R148: CH3 + aC3H5 = C4H81
      RF(148) = EXP(3.223619130191664D1 -3.2D-1*ALOGT
     * +1.319937319219509D2*TI)
      EQK = EG(6)/EG(13)/EG(32)/PFAC1
      RB(148) = RF(148) / MAX(EQK, SMALL)
C     R149: C3H6 + H = C2H4 + CH3
      RF(149) = EXP(5.043372849455479D1 -2.39D0*ALOGT
     * -5.625962344214302D3*TI)
      EQK = EG(2)*EG(13)/EG(4)/EG(20)
      RB(149) = RF(149) / MAX(EQK, SMALL)
C     R150: C3H6 + H = H2 + aC3H5
      RF(150) = EXP(1.446783574400397D1 +2.309D0*ALOGT
     * -2.352034704548985D3*TI)
      EQK = EG(7)*EG(32)/EG(4)/EG(20)
      RB(150) = RF(150) / MAX(EQK, SMALL)
C     R151: C3H6 + O = CH3 + H + CH2CO
      RF(151) = EXP(1.819753719263815D1 +1.65D0*ALOGT
     * -1.645518503182537D2*TI)
      EQK = EG(13)*EG(20)*EG(30)/EG(4)/EG(15)*PFAC1
      RB(151) = RF(151) / MAX(EQK, SMALL)
C     R152: C3H6 + O = HCO + C2H5
      RF(152) = EXP(1.737085861945369D1 +1.65D0*ALOGT
     * +4.891266009460019D2*TI)
      EQK = EG(23)*EG(28)/EG(4)/EG(15)
      RB(152) = RF(152) / MAX(EQK, SMALL)
C     R153: C3H6 + O = OH + aC3H5
      RF(153) = EXP(2.591622268783662D1 +7.D-1*ALOGT
     * -2.958914005722728D3*TI)
      EQK = EG(16)*EG(32)/EG(4)/EG(15)
      RB(153) = RF(153) / MAX(EQK, SMALL)
C     R154: C3H6 + OH = H2O + aC3H5
      RF(154) = EXP(1.494691266945537D1 +2.D0*ALOGT
     * +1.499585669566961D2*TI)
      EQK = EG(18)*EG(32)/EG(4)/EG(16)
      RB(154) = RF(154) / MAX(EQK, SMALL)
C     R155: C3H6 + CH3 = CH4 + aC3H5
      RF(155) = EXP(3.815512105047302D0 +3.268D0*ALOGT
     * -3.599005606960705D3*TI)
      EQK = EG(3)*EG(32)/EG(4)/EG(13)
      RB(155) = RF(155) / MAX(EQK, SMALL)
C     R156: C4H81 + H = C2H4 + C2H5
      RF(156) = EXP(5.112687567511474D1 -2.39D0*ALOGT
     * -5.625962344214302D3*TI)
      EQK = EG(2)*EG(28)/EG(6)/EG(20)
      RB(156) = RF(156) / MAX(EQK, SMALL)
C     R157: C4H81 + H = C3H6 + CH3
      RF(157) = EXP(5.182002285567469D1 -2.39D0*ALOGT
     * -5.625962344214302D3*TI)
      EQK = EG(4)*EG(13)/EG(6)/EG(20)
      RB(157) = RF(157) / MAX(EQK, SMALL)
C     R158: iC4H8 + H = C3H6 + CH3
      RF(158) = EXP(5.247494882341444D1 -2.538D0*ALOGT
     * -5.851906627984626D3*TI)
      EQK = EG(4)*EG(13)/EG(5)/EG(20)
      RB(158) = RF(158) / MAX(EQK, SMALL)
C     R159: iC4H8 + O = 2CH3 + CH2CO
      RF(159) = EXP(1.929614948130627D1 +1.65D0*ALOGT
     * -1.645518503182537D2*TI)
      EQK = EG(13)*EG(13)*EG(30)/EG(5)/EG(15)*PFAC1
      RB(159) = RF(159) / MAX(EQK, SMALL)
C     R160: CH3 + C6H5 = C6H5CH3
      RF(160) = EXP(3.025568970809171D1 -2.31479667114363D1*TI)
      EQK = EG(12)/EG(13)/EG(36)/PFAC1
      RB(160) = RF(160) / MAX(EQK, SMALL)
C     R161: C6H5CH3 + O2 = HO2 + C6H5CH2
      RF(161) = EXP(3.333480359058475D1 -2.163429097517543D4*TI)
      EQK = EG(17)*EG(37)/EG(12)/EG(14)
      RB(161) = RF(161) / MAX(EQK, SMALL)
C     R162: C6H5CH3 + OH = H2O + C6H5CH2
      RF(162) = EXP(3.041603235816689D1 -1.393910169362578D3*TI)
      EQK = EG(18)*EG(37)/EG(12)/EG(16)
      RB(162) = RF(162) / MAX(EQK, SMALL)
C     R163: C6H5CH3 + H = H2 + C6H5CH2
      RF(163) = EXP(3.246650905697885D1 -4.206388124802089D3*TI)
      EQK = EG(7)*EG(37)/EG(12)/EG(20)
      RB(163) = RF(163) / MAX(EQK, SMALL)
C     R164: C6H5CH3 + H = C6H6 + CH3
      RF(164) = EXP(1.447303056088107D1 +2.17D0*ALOGT
     * -2.094890987384986D3*TI)
      EQK = EG(10)*EG(13)/EG(12)/EG(20)
      RB(164) = RF(164) / MAX(EQK, SMALL)
C     R165: C6H5CH3 + CH3 = CH4 + C6H5CH2
      RF(165) = EXP(2.647900805053332D1 -4.780558342579237D3*TI)
      EQK = EG(3)*EG(37)/EG(12)/EG(13)
      RB(165) = RF(165) / MAX(EQK, SMALL)
C     R166: H + C6H5CH2 = C6H5CH3
      RF(166) = 1.D14
      EQK = EG(12)/EG(20)/EG(37)/PFAC1
      RB(166) = RF(166) / MAX(EQK, SMALL)
C     R167: H + C6H5CH2 = CH3 + C6H5
      RF(167) = EXP(1.523760812457152D2 -1.394D1*ALOGT
     * -3.249773239618601D4*TI)
      EQK = EG(13)*EG(36)/EG(20)/EG(37)
      RB(167) = RF(167) / MAX(EQK, SMALL)
C     R168: O + C6H5CH2 = H + C6H5CHO
      RF(168) = 4.D14
      EQK = EG(20)*EG(40)/EG(15)/EG(37)
      RB(168) = RF(168) / MAX(EQK, SMALL)
C     R169: HO2 + C6H5CH2 = OH + H + C6H5CHO
      RF(169) = 5.D12
      EQK = EG(16)*EG(20)*EG(40)/EG(17)/EG(37)*PFAC1
      RB(169) = RF(169) / MAX(EQK, SMALL)
C     R170: HCO + C6H5 = C6H5CHO
      RF(170) = 1.D13
      EQK = EG(40)/EG(23)/EG(36)/PFAC1
      RB(170) = RF(170) / MAX(EQK, SMALL)
C     R171: C6H5CHO = H + C6H5CO
      RF(171) = EXP(3.592005821420703D1 -4.372952841790902D4*TI)
      EQK = EG(20)*EG(39)/EG(40)*PFAC1
      RB(171) = RF(171) / MAX(EQK, SMALL)
C     R172: O2 + C6H5CHO = HO2 + C6H5CO
      RF(172) = EXP(2.995340883621877D1 -1.960028920457487D4*TI)
      EQK = EG(17)*EG(39)/EG(14)/EG(40)
      RB(172) = RF(172) / MAX(EQK, SMALL)
C     R173: OH + C6H5CHO = H2O + C6H5CO
      RF(173) = EXP(2.388026625809653D1 +7.3D-1*ALOGT
     * +5.585705010803108D2*TI)
      EQK = EG(18)*EG(39)/EG(16)/EG(40)
      RB(173) = RF(173) / MAX(EQK, SMALL)
C     R174: H + C6H5CHO = H2 + C6H5CO
      RF(174) = EXP(2.213425281065667D1 +1.16D0*ALOGT
     * -1.207720002335807D3*TI)
      EQK = EG(7)*EG(39)/EG(20)/EG(40)
      RB(174) = RF(174) / MAX(EQK, SMALL)
C     R175: H + C6H5CHO = C6H6 + HCO
      RF(175) = EXP(1.447303056088107D1 +2.17D0*ALOGT
     * -2.094890987384986D3*TI)
      EQK = EG(10)*EG(23)/EG(20)/EG(40)
      RB(175) = RF(175) / MAX(EQK, SMALL)
C     R176: O + C6H5CHO = OH + C6H5CO
      RF(176) = EXP(2.938887903348092D1 -9.057900017518555D2*TI)
      EQK = EG(16)*EG(39)/EG(15)/EG(40)
      RB(176) = RF(176) / MAX(EQK, SMALL)
C     R177: C6H5CH2 + C6H5CHO = C6H5CH3 + C6H5CO
      RF(177) = EXP(-1.312236337740433D1 +5.6D0*ALOGT
     * -1.237913002394203D3*TI)
      EQK = EG(12)*EG(39)/EG(37)/EG(40)
      RB(177) = RF(177) / MAX(EQK, SMALL)
C     R178: CH3 + C6H5CHO = CH4 + C6H5CO
      RF(178) = EXP(-1.312236337740433D1 +5.6D0*ALOGT
     * -1.237913002394203D3*TI)
      EQK = EG(3)*EG(39)/EG(13)/EG(40)
      RB(178) = RF(178) / MAX(EQK, SMALL)
C     R179: H2O2 + C6H5CO = HO2 + C6H5CHO
      RF(179) = EXP(2.591622268783662D1 -4.13946030800598D3*TI)
      EQK = EG(17)*EG(40)/EG(19)/EG(39)
      RB(179) = RF(179) / MAX(EQK, SMALL)
C     R180: C6H5CO = CO + C6H5
      RF(180) = EXP(3.389822166446991D1 -1.459982517823699D4*TI)
      EQK = EG(9)*EG(36)/EG(39)*PFAC1
      RB(180) = RF(180) / MAX(EQK, SMALL)
C     R181: H + C6H5 = C6H6
      RF(181) = 1.D14
      EQK = EG(10)/EG(20)/EG(36)/PFAC1
      RB(181) = RF(181) / MAX(EQK, SMALL)
C     R182: C6H6 + OH = H2O + C6H5
      RF(182) = EXP(1.289546277721241D1 +2.286D0*ALOGT
     * -5.32403234363035D2*TI)
      EQK = EG(18)*EG(36)/EG(10)/EG(16)
      RB(182) = RF(182) / MAX(EQK, SMALL)
C     R183: C6H6 + O = H + C6H5O
      RF(183) = EXP(3.026290995606519D1 -2.470793838112006D3*TI)
      EQK = EG(20)*EG(38)/EG(10)/EG(15)
      RB(183) = RF(183) / MAX(EQK, SMALL)
C     R184: C6H6 + O = HCO + C5H5
      RF(184) = EXP(3.026290995606519D1 -2.279571504408836D3*TI)
      EQK = EG(23)*EG(35)/EG(10)/EG(15)
      RB(184) = RF(184) / MAX(EQK, SMALL)
C     R185: H2 + C6H5 = C6H6 + H
      RF(185) = EXP(1.095203386353075D1 +2.43D0*ALOGT
     * -3.156678156105216D3*TI)
      EQK = EG(10)*EG(20)/EG(7)/EG(36)
      RB(185) = RF(185) / MAX(EQK, SMALL)
C     R186: O2 + C6H5 = O + C6H5O
      RF(186) = EXP(3.088911765395003D1 -3.079686005956308D3*TI)
      EQK = EG(15)*EG(38)/EG(14)/EG(36)
      RB(186) = RF(186) / MAX(EQK, SMALL)
C     R187: O + C6H5 = CO + C5H5
      RF(187) = 1.D14
      EQK = EG(9)*EG(35)/EG(15)/EG(36)
      RB(187) = RF(187) / MAX(EQK, SMALL)
C     R188: OH + C6H5 = H + C6H5O
      RF(188) = 3.D13
      EQK = EG(20)*EG(38)/EG(16)/EG(36)
      RB(188) = RF(188) / MAX(EQK, SMALL)
C     R189: HO2 + C6H5 = C6H6 + O2
      RF(189) = 1.D12
      EQK = EG(10)*EG(14)/EG(17)/EG(36)
      RB(189) = RF(189) / MAX(EQK, SMALL)
C     R190: CH4 + C6H5 = C6H6 + CH3
      RF(190) = EXP(-5.549346121351782D0 +4.57D0*ALOGT
     * -2.644906805115418D3*TI)
      EQK = EG(10)*EG(13)/EG(3)/EG(36)
      RB(190) = RF(190) / MAX(EQK, SMALL)
C     R191: H + C6H5O = HCO + C5H5
      RF(191) = EXP(2.993360620892259D1 -6.038600011679036D3*TI)
      EQK = EG(23)*EG(35)/EG(20)/EG(38)
      RB(191) = RF(191) / MAX(EQK, SMALL)
C     R192: C6H5O = CO + C5H5
      RF(192) = EXP(1.256640139790803D2 -1.206D1*ALOGT
     * -3.663417340418615D4*TI)
      EQK = EG(9)*EG(35)/EG(38)*PFAC1
      RB(192) = RF(192) / MAX(EQK, SMALL)
C     R193: O + C5H5 = H + C5H4O
      RF(193) = EXP(3.169318677978641D1 -2.D-2*ALOGT
     * -1.006433335279839D1*TI)
      EQK = EG(20)*EG(34)/EG(15)/EG(35)
      RB(193) = RF(193) / MAX(EQK, SMALL)
C     R194: O2 + C5H5 = OH + C5H4O
      RF(194) = EXP(2.821880778083067D1 +8.D-2*ALOGT
     * -9.057900017518554D3*TI)
      EQK = EG(16)*EG(34)/EG(14)/EG(35)
      RB(194) = RF(194) / MAX(EQK, SMALL)
C     R195: C5H4O = CO + 2C2H2
      RF(195) = EXP(9.623053810480693D1 -7.87D0*ALOGT
     * -4.966748509606007D4*TI)
      EQK = EG(9)*EG(11)*EG(11)/EG(34)*PFAC2
      RB(195) = RF(195) / MAX(EQK, SMALL)
C     R196: H + C5H4O = CO + c-C4H5
      RF(196) = EXP(2.218188085964593D1 +1.45D0*ALOGT
     * -1.962545003795687D3*TI)
      EQK = EG(9)*EG(33)/EG(20)/EG(34)
      RB(196) = RF(196) / MAX(EQK, SMALL)
C     R197: H + c-C4H5 = C2H4 + C2H2
      RF(197) = 1.D13
      EQK = EG(2)*EG(11)/EG(20)/EG(33)
      RB(197) = RF(197) / MAX(EQK, SMALL)
C     R198: O + c-C4H5 = C2H2 + CH2CHO
      RF(198) = 1.D14
      EQK = EG(11)*EG(31)/EG(15)/EG(33)
      RB(198) = RF(198) / MAX(EQK, SMALL)
C     R199: O2 + c-C4H5 = CH2CO + CH2CHO
      RF(199) = EXP(2.689705194084835D1 -9.561116685158473D3*TI)
      EQK = EG(30)*EG(31)/EG(14)/EG(33)
      RB(199) = RF(199) / MAX(EQK, SMALL)
C     R200: c-C4H5 = C2H2 + C2H3
      RF(200) = EXP(2.832416829648849D1 -2.918656672311534D4*TI)
      EQK = EG(11)*EG(27)/EG(33)*PFAC1
      RB(200) = RF(200) / 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(2.37261574758191D2 -2.463D1*ALOGT
     * -7.341931180866428D3*TI)
      RKLOW(21) = 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 aC3H5
      IF (T .GT. 1000) THEN
      SMH(32) = -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(32) = 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 c-C4H5
      IF (T .GT. 1000) THEN
      SMH(33) = -1.291288D1 -3.2808359D4*TI
     *         +6.7467155D0*TN(1) +8.6415D-3*TN(2)
     *         -1.086142983333333D-6*TN(3) +8.243131166666666D-11*TN(4)
     *         -1.7302454D-15*TN(5)
      ELSE
      SMH(33) = 3.5701797D1 -3.5373828D4*TI
     *         -2.6397593D0*TN(1) +2.07745785D-2*TN(2)
     *         -3.653492333333333D-6*TN(3) -3.879917833333333D-10*TN(4)
     *         +3.0674445D-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 C6H5CH2
      IF (T .GT. 1000) THEN
      SMH(37) = -5.1665589D1 -1.8564203D4*TI
     *         +1.404398D1*TN(1) +1.17469365D-2*TN(2)
     *         -1.422922783333333D-6*TN(3) +1.157570083333333D-10*TN(4)
     *         -4.180721D-15*TN(5)
      ELSE
      SMH(37) = 2.354882D1 -2.3307027D4*TI
     *         +4.811154D-1*TN(1) +1.9256416D-2*TN(2)
     *         +5.476915333333333D-6*TN(3) -6.414393416666667D-9*TN(4)
     *         +1.7711534D-12*TN(5)
      ENDIF
C C6H5O
      IF (T .GT. 1000) THEN
      SMH(38) = -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(38) = 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 C6H5CO
      IF (T .GT. 1000) THEN
      SMH(39) = -4.4659218D1 -6.9147837D3*TI
     *         +1.3374409D1*TN(1) +1.19996445D-2*TN(2)
     *         -1.744287333333333D-6*TN(3) +1.805760916666667D-10*TN(4)
     *         -9.0035225D-15*TN(5)
      ELSE
      SMH(39) = 3.5778175D1 -1.1255803D4*TI
     *         -2.0251155D0*TN(1) +3.07562705D-2*TN(2)
     *         -5.2672755D-6*TN(3) -5.81038325D-10*TN(4)
     *         +3.991757449999999D-13*TN(5)
      ENDIF
C C6H5CHO
      IF (T .GT. 1000) THEN
      SMH(40) = -4.7965796D1 +1.1019744D4*TI
     *         +1.3650737D1*TN(1) +1.28402095D-2*TN(2)
     *         -1.744454833333333D-6*TN(3) +1.617785833333333D-10*TN(4)
     *         -6.741896D-15*TN(5)
      ELSE
      SMH(40) = 4.0231735D1 +6.1169349D3*TI
     *         -3.1627334D0*TN(1) +3.31846225D-2*TN(2)
     *         -5.8027255D-6*TN(3) -5.249948083333333D-10*TN(4)
     *         +4.29035505D-13*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, 31
         CTOT = CTOT + C(K)
      ENDDO
C
C     R1: POSF12345 = 1.72699C2H4 + 0.518098C3H6 + 0.0863497iC4H8 + 0.0690798C4H81 + 0.13C6H6 + 0.37C6H5CH3 + CH3 + H
      RF(1) = RF(1)*C(1)
C     R2: POSF12345 + 0.35H = 1.91642C2H4 + 0.01CH4 + 0.574926C3H6 + 0.0958211iC4H8 + 0.0766569C4H81 + H2 + 0.130867C6H6 + 0.372467C6H5CH3 + 0.35CH3
      RF(2) = RF(2)*C(1)*C(20)
C     R3: POSF12345 + 0.65CH3 = 1.91642C2H4 + 1.01CH4 + 0.574926C3H6 + 0.0958211iC4H8 + 0.0766569C4H81 + 0.130867C6H6 + 0.372467C6H5CH3 + 0.65H
      RF(3) = RF(3)*C(1)*C(13)
C     R4: POSF12345 + OH = 1.91642C2H4 + 0.01CH4 + 0.574926C3H6 + 0.0958211iC4H8 + 0.0766569C4H81 + 0.130867C6H6 + 0.372467C6H5CH3 + 0.35CH3 + H2O + 0.65H
      RF(4) = RF(4)*C(1)*C(16)
C     R5: POSF12345 + O2 = 1.91642C2H4 + 0.01CH4 + 0.574926C3H6 + 0.0958211iC4H8 + 0.0766569C4H81 + 0.130867C6H6 + 0.372467C6H5CH3 + 0.35CH3 + HO2 + 0.65H
      RF(5) = RF(5)*C(1)*C(14)
C     R6: POSF12345 + HO2 = 1.91642C2H4 + 0.01CH4 + 0.574926C3H6 + 0.0958211iC4H8 + 0.0766569C4H81 + 0.130867C6H6 + 0.372467C6H5CH3 + 0.35CH3 + H2O2 + 0.65H
      RF(6) = RF(6)*C(1)*C(17)
C     R7: POSF12345 + O = 1.91642C2H4 + 0.01CH4 + 0.574926C3H6 + 0.0958211iC4H8 + 0.0766569C4H81 + 0.130867C6H6 + 0.372467C6H5CH3 + 0.35CH3 + OH + 0.65H
      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(18)*C(20)*C(20)
      RB(13) = RB(13)*C(7)*C(18)
C     R14: OH + H = H2O
      CTB = CTOT+C(7)+5.3D0*C(18)
     * +7.5D-1*C(9)+2.6D0*C(22)
      RF(14) = RF(14)*CTB*C(16)*C(20)
      RB(14) = RB(14)*CTB*C(18)
C     R15: O + H = OH
      CTB = CTOT+C(7)+1.1D1*C(18)
     * +7.5D-1*C(9)+2.6D0*C(22)
      RF(15) = RF(15)*CTB*C(15)*C(20)
      RB(15) = RB(15)*CTB*C(16)
C     R16: 2O = O2
      CTB = CTOT+1.4D0*C(7)+1.44D1*C(18)
     * +7.5D-1*C(9)+2.6D0*C(22)
      RF(16) = RF(16)*CTB*C(15)*C(15)
      RB(16) = RB(16)*CTB*C(14)
C     R17: 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(17)
      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(17) = RF(17) * PCOR
      RB(17) = RB(17) * PCOR
      RF(17) = RF(17)*C(14)*C(20)
      RB(17) = RB(17)*C(17)
C     R18: H2 + O2 = HO2 + H
      RF(18) = RF(18)*C(7)*C(14)
      RB(18) = RB(18)*C(17)*C(20)
C     R19: 2OH = H2O2
      CTB = CTOT+C(7)+5.D0*C(18)
     * +7.5D-1*C(9)+2.6D0*C(22)
      PR = RKLOW(2) * CTB / RF(19)
      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(19) = RF(19) * PCOR
      RB(19) = RB(19) * PCOR
      RF(19) = RF(19)*C(16)*C(16)
      RB(19) = RB(19)*C(19)
C     R20: HO2 + H = 2OH
      RF(20) = RF(20)*C(17)*C(20)
      RB(20) = RB(20)*C(16)*C(16)
C     R21: O + HO2 = O2 + OH
      RF(21) = RF(21)*C(15)*C(17)
      RB(21) = RB(21)*C(14)*C(16)
C     R22: 2HO2 = O2 + H2O2
      RF(22) = RF(22)*C(17)*C(17)
      RB(22) = RB(22)*C(14)*C(19)
C     R23: 2HO2 = O2 + H2O2
      RF(23) = RF(23)*C(17)*C(17)
      RB(23) = RB(23)*C(14)*C(19)
C     R24: OH + HO2 = O2 + H2O
      RF(24) = RF(24)*C(16)*C(17)
      RB(24) = RB(24)*C(14)*C(18)
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: H2O2 + H = H2 + HO2
      RF(27) = RF(27)*C(19)*C(20)
      RB(27) = RB(27)*C(7)*C(17)
C     R28: H2O2 + H = OH + H2O
      RF(28) = RF(28)*C(19)*C(20)
      RB(28) = RB(28)*C(16)*C(18)
C     R29: O + H2O2 = OH + HO2
      RF(29) = RF(29)*C(15)*C(19)
      RB(29) = RB(29)*C(16)*C(17)
C     R30: OH + H2O2 = HO2 + H2O
      RF(30) = RF(30)*C(16)*C(19)
      RB(30) = RB(30)*C(17)*C(18)
C     R31: 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(31)
      PCOR = PR / (1.0 + PR)
      RF(31) = RF(31) * PCOR
      RB(31) = RB(31) * PCOR
      RF(31) = RF(31)*C(9)*C(15)
      RB(31) = RB(31)*C(22)
C     R32: CO + OH = H + CO2
      RF(32) = RF(32)*C(9)*C(16)
      RB(32) = RB(32)*C(20)*C(22)
C     R33: CO + OH = H + CO2
      RF(33) = RF(33)*C(9)*C(16)
      RB(33) = RB(33)*C(20)*C(22)
C     R34: CO + O2 = O + CO2
      RF(34) = RF(34)*C(9)*C(14)
      RB(34) = RB(34)*C(15)*C(22)
C     R35: CO + HO2 = OH + CO2
      RF(35) = RF(35)*C(9)*C(17)
      RB(35) = RB(35)*C(16)*C(22)
C     R36: H + HCO = H2 + CO
      RF(36) = RF(36)*C(20)
      RB(36) = RB(36)*C(7)*C(9)
C     R37: O + HCO = CO + OH
      RF(37) = RF(37)*C(15)
      RB(37) = RB(37)*C(9)*C(16)
C     R38: O + HCO = H + CO2
      RF(38) = RF(38)*C(15)
      RB(38) = RB(38)*C(20)*C(22)
C     R39: OH + HCO = CO + H2O
      RF(39) = RF(39)*C(16)
      RB(39) = RB(39)*C(9)*C(18)
C     R40: HCO = CO + H
      CTB = CTOT+C(7)-C(18)
     * +7.5D-1*C(9)+2.6D0*C(22)
      RF(40) = RF(40)*CTB
      RB(40) = RB(40)*CTB*C(9)*C(20)
C     R41: HCO = CO + H
      RF(41) = RF(41)*C(18)
      RB(41) = RB(41)*C(9)*C(18)*C(20)
C     R42: O2 + HCO = CO + HO2
      RF(42) = RF(42)*C(14)
      RB(42) = RB(42)*C(9)*C(17)
C     R43: 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(43)
      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(43) = RF(43) * PCOR
      RB(43) = RB(43) * PCOR
      RF(43) = RF(43)*C(7)*C(9)
      RB(43) = RB(43)*C(21)
C     R44: 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(44)
      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(44) = RF(44) * PCOR
      RB(44) = RB(44) * PCOR
      RF(44) = RF(44)*C(20)
      RB(44) = RB(44)*C(21)
C     R45: 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(45)
      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(45) = RF(45) * PCOR
      RB(45) = RB(45) * PCOR
      RF(45) = RF(45)*C(20)
      RB(45) = RB(45)*C(13)
C     R46: O + CH2 = H + HCO
      RF(46) = RF(46)*C(15)
      RB(46) = RB(46)*C(20)
C     R47: OH + CH2 = H + CH2O
      RF(47) = RF(47)*C(16)
      RB(47) = RB(47)*C(20)*C(21)
C     R48: H2 + CH2 = CH3 + H
      RF(48) = RF(48)*C(7)
      RB(48) = RB(48)*C(13)*C(20)
C     R49: O2 + CH2 = OH + HCO
      RF(49) = RF(49)*C(14)
      RB(49) = RB(49)*C(16)
C     R50: O2 + CH2 = 2H + CO2
      RF(50) = RF(50)*C(14)
      RB(50) = RB(50)*C(20)*C(20)*C(22)
C     R51: 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(51)
      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(51) = RF(51) * PCOR
      RB(51) = RB(51) * PCOR
      RF(51) = RF(51)*C(9)
      RB(51) = RB(51)*C(23)
C     R52: CH2* = CH2
      RF(52) = RF(52)*C(31)
      RB(52) = RB(52)*C(31)
C     R53: O + CH2* = H2 + CO
      RF(53) = RF(53)*C(15)
      RB(53) = RB(53)*C(7)*C(9)
C     R54: OH + CH2* = H + CH2O
      RF(54) = RF(54)*C(16)
      RB(54) = RB(54)*C(20)*C(21)
C     R55: H2 + CH2* = CH3 + H
      RF(55) = RF(55)*C(7)
      RB(55) = RB(55)*C(13)*C(20)
C     R56: O2 + CH2* = CO + OH + H
      RF(56) = RF(56)*C(14)
      RB(56) = RB(56)*C(9)*C(16)*C(20)
C     R57: O2 + CH2* = CO + H2O
      RF(57) = RF(57)*C(14)
      RB(57) = RB(57)*C(9)*C(18)
C     R58: CH2* = CH2
      RF(58) = RF(58)*C(18)
      RB(58) = RB(58)*C(18)
C     R59: CH2* = CH2
      RF(59) = RF(59)*C(9)
      RB(59) = RB(59)*C(9)
C     R60: CH2* = CH2
      RF(60) = RF(60)*C(22)
      RB(60) = RB(60)*C(22)
C     R61: CH2* + CO2 = CO + CH2O
      RF(61) = RF(61)*C(22)
      RB(61) = RB(61)*C(9)*C(21)
C     R62: 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(62)
      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(62) = RF(62) * PCOR
      RB(62) = RB(62) * PCOR
      RF(62) = RF(62)*C(20)*C(21)
C     R63: H + CH2O = H2 + HCO
      RF(63) = RF(63)*C(20)*C(21)
      RB(63) = RB(63)*C(7)
C     R64: O + CH2O = OH + HCO
      RF(64) = RF(64)*C(15)*C(21)
      RB(64) = RB(64)*C(16)
C     R65: OH + CH2O = H2O + HCO
      RF(65) = RF(65)*C(16)*C(21)
      RB(65) = RB(65)*C(18)
C     R66: O2 + CH2O = HO2 + HCO
      RF(66) = RF(66)*C(14)*C(21)
      RB(66) = RB(66)*C(17)
C     R67: HO2 + CH2O = H2O2 + HCO
      RF(67) = RF(67)*C(17)*C(21)
      RB(67) = RB(67)*C(19)
C     R68: 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(68)
      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(68) = RF(68) * PCOR
      RB(68) = RB(68) * PCOR
      RF(68) = RF(68)*C(13)*C(20)
      RB(68) = RB(68)*C(3)
C     R69: CH3 + O = H + CH2O
      RF(69) = RF(69)*C(13)*C(15)
      RB(69) = RB(69)*C(20)*C(21)
C     R70: CH3 + OH = H2O + CH2
      RF(70) = RF(70)*C(13)*C(16)
      RB(70) = RB(70)*C(18)
C     R71: CH3 + OH = H2O + CH2*
      RF(71) = RF(71)*C(13)*C(16)
      RB(71) = RB(71)*C(18)
C     R72: CH3 + O2 = O + CH3O
      RF(72) = RF(72)*C(13)*C(14)
      RB(72) = RB(72)*C(15)
C     R73: CH3 + O2 = OH + CH2O
      RF(73) = RF(73)*C(13)*C(14)
      RB(73) = RB(73)*C(16)*C(21)
C     R74: CH3 + HO2 = CH4 + O2
      RF(74) = RF(74)*C(13)*C(17)
      RB(74) = RB(74)*C(3)*C(14)
C     R75: CH3 + HO2 = OH + CH3O
      RF(75) = RF(75)*C(13)*C(17)
      RB(75) = RB(75)*C(16)
C     R76: CH3 + H2O2 = CH4 + HO2
      RF(76) = RF(76)*C(13)*C(19)
      RB(76) = RB(76)*C(3)*C(17)
C     R77: CH3 + HCO = CH4 + CO
      RF(77) = RF(77)*C(13)
      RB(77) = RB(77)*C(3)*C(9)
C     R78: CH3 + CH2O = CH4 + HCO
      RF(78) = RF(78)*C(13)*C(21)
      RB(78) = RB(78)*C(3)
C     R79: CH3 + CH2 = C2H4 + H
      RF(79) = RF(79)*C(13)
      RB(79) = RB(79)*C(2)*C(20)
C     R80: 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(80)
      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(80) = RF(80) * PCOR
      RB(80) = RB(80) * PCOR
      RF(80) = RF(80)*C(13)*C(13)
      RB(80) = RB(80)*C(8)
C     R81: 2CH3 = H + C2H5
      RF(81) = RF(81)*C(13)*C(13)
      RB(81) = RB(81)*C(20)
C     R82: CH3 + HCCO = C2H4 + CO
      RF(82) = RF(82)*C(13)
      RB(82) = RB(82)*C(2)*C(9)
C     R83: H + CH3O = H2 + CH2O
      RF(83) = RF(83)*C(20)
      RB(83) = RB(83)*C(7)*C(21)
C     R84: H + CH3O = CH3 + OH
      RF(84) = RF(84)*C(20)
      RB(84) = RB(84)*C(13)*C(16)
C     R85: H + CH3O = H2O + CH2*
      RF(85) = RF(85)*C(20)
      RB(85) = RB(85)*C(18)
C     R86: O + CH3O = OH + CH2O
      RF(86) = RF(86)*C(15)
      RB(86) = RB(86)*C(16)*C(21)
C     R87: OH + CH3O = H2O + CH2O
      RF(87) = RF(87)*C(16)
      RB(87) = RB(87)*C(18)*C(21)
C     R88: O2 + CH3O = HO2 + CH2O
      RF(88) = RF(88)*C(14)
      RB(88) = RB(88)*C(17)*C(21)
C     R89: CH4 + H = H2 + CH3
      RF(89) = RF(89)*C(3)*C(20)
      RB(89) = RB(89)*C(7)*C(13)
C     R90: CH4 + O = CH3 + OH
      RF(90) = RF(90)*C(3)*C(15)
      RB(90) = RB(90)*C(13)*C(16)
C     R91: CH4 + OH = CH3 + H2O
      RF(91) = RF(91)*C(3)*C(16)
      RB(91) = RB(91)*C(13)*C(18)
C     R92: CH4 + CH2* = 2CH3
      RF(92) = RF(92)*C(3)
      RB(92) = RB(92)*C(13)*C(13)
C     R93: H + HCCO = CO + CH2*
      RF(93) = RF(93)*C(20)
      RB(93) = RB(93)*C(9)
C     R94: O + HCCO = 2CO + H
      RF(94) = RF(94)*C(15)
      RB(94) = RB(94)*C(9)*C(9)*C(20)
C     R95: O2 + HCCO = 2CO + OH
      RF(95) = RF(95)*C(14)
      RB(95) = RB(95)*C(9)*C(9)*C(16)
C     R96: 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(96)
      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(96) = RF(96) * PCOR
      RB(96) = RB(96) * PCOR
      RB(96) = RB(96)*C(11)*C(20)
C     R97: C2H2 + O = CO + CH2
      RF(97) = RF(97)*C(11)*C(15)
      RB(97) = RB(97)*C(9)
C     R98: C2H2 + O = H + HCCO
      RF(98) = RF(98)*C(11)*C(15)
      RB(98) = RB(98)*C(20)
C     R99: C2H2 + OH = H + CH2CO
      RF(99) = RF(99)*C(11)*C(16)
      RB(99) = RB(99)*C(20)*C(23)
C     R100: C2H2 + OH = CO + CH3
      RF(100) = RF(100)*C(11)*C(16)
      RB(100) = RB(100)*C(9)*C(13)
C     R101: C2H2 + HCO = CO + C2H3
      RF(101) = RF(101)*C(11)
      RB(101) = RB(101)*C(9)
C     R102: C2H2 + CH3 = aC3H5
      RF(102) = RF(102)*C(11)*C(13)
      RB(102) = RB(102)*C(24)
C     R103: 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(103)
      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(103) = RF(103) * PCOR
      RB(103) = RB(103) * PCOR
      RF(103) = RF(103)*C(20)*C(23)
C     R104: H + CH2CO = H2 + HCCO
      RF(104) = RF(104)*C(20)*C(23)
      RB(104) = RB(104)*C(7)
C     R105: H + CH2CO = CO + CH3
      RF(105) = RF(105)*C(20)*C(23)
      RB(105) = RB(105)*C(9)*C(13)
C     R106: O + CH2CO = OH + HCCO
      RF(106) = RF(106)*C(15)*C(23)
      RB(106) = RB(106)*C(16)
C     R107: O + CH2CO = CH2 + CO2
      RF(107) = RF(107)*C(15)*C(23)
      RB(107) = RB(107)*C(22)
C     R108: OH + CH2CO = H2O + HCCO
      RF(108) = RF(108)*C(16)*C(23)
      RB(108) = RB(108)*C(18)
C     R109: 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(109)
      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(109) = RF(109) * PCOR
      RB(109) = RB(109) * PCOR
      RF(109) = RF(109)*C(20)
      RB(109) = RB(109)*C(2)
C     R110: H + C2H3 = H2 + C2H2
      RF(110) = RF(110)*C(20)
      RB(110) = RB(110)*C(7)*C(11)
C     R111: O + C2H3 = H + CH2CO
      RF(111) = RF(111)*C(15)
      RB(111) = RB(111)*C(20)*C(23)
C     R112: O + C2H3 = CO + CH3
      RF(112) = RF(112)*C(15)
      RB(112) = RB(112)*C(9)*C(13)
C     R113: OH + C2H3 = C2H2 + H2O
      RF(113) = RF(113)*C(16)
      RB(113) = RB(113)*C(11)*C(18)
C     R114: O2 + C2H3 = C2H2 + HO2
      RF(114) = RF(114)*C(14)
      RB(114) = RB(114)*C(11)*C(17)
C     R115: O2 + C2H3 = O + CH2CHO
      RF(115) = RF(115)*C(14)
      RB(115) = RB(115)*C(15)
C     R116: O2 + C2H3 = HCO + CH2O
      RF(116) = RF(116)*C(14)
      RB(116) = RB(116)*C(21)
C     R117: 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(117)
      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(117) = RF(117) * PCOR
      RB(117) = RB(117) * PCOR
      RF(117) = RF(117)*C(13)
      RB(117) = RB(117)*C(4)
C     R118: CH3 + C2H3 = H + aC3H5
      RF(118) = RF(118)*C(13)
      RB(118) = RB(118)*C(20)*C(24)
C     R119: CH2CHO = CO + CH3
      RB(119) = RB(119)*C(9)*C(13)
C     R120: H + CH2CHO = CH3 + HCO
      RF(120) = RF(120)*C(20)
      RB(120) = RB(120)*C(13)
C     R121: H + CH2CHO = H2 + CH2CO
      RF(121) = RF(121)*C(20)
      RB(121) = RB(121)*C(7)*C(23)
C     R122: O + CH2CHO = OH + CH2CO
      RF(122) = RF(122)*C(15)
      RB(122) = RB(122)*C(16)*C(23)
C     R123: OH + CH2CHO = H2O + CH2CO
      RF(123) = RF(123)*C(16)
      RB(123) = RB(123)*C(18)*C(23)
C     R124: O2 + CH2CHO = HO2 + CH2CO
      RF(124) = RF(124)*C(14)
      RB(124) = RB(124)*C(17)*C(23)
C     R125: 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(125)
      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(125) = RF(125) * PCOR
      RB(125) = RB(125) * PCOR
      RF(125) = RF(125)*C(2)*C(20)
C     R126: C2H4 + H = H2 + C2H3
      RF(126) = RF(126)*C(2)*C(20)
      RB(126) = RB(126)*C(7)
C     R127: C2H4 + O = OH + C2H3
      RF(127) = RF(127)*C(2)*C(15)
      RB(127) = RB(127)*C(16)
C     R128: C2H4 + O = CH3 + HCO
      RF(128) = RF(128)*C(2)*C(15)
      RB(128) = RB(128)*C(13)
C     R129: C2H4 + O = CH2 + CH2O
      RF(129) = RF(129)*C(2)*C(15)
      RB(129) = RB(129)*C(21)
C     R130: C2H4 + OH = H2O + C2H3
      RF(130) = RF(130)*C(2)*C(16)
      RB(130) = RB(130)*C(18)
C     R131: C2H4 + HCO = CO + C2H5
      RF(131) = RF(131)*C(2)
      RB(131) = RB(131)*C(9)
C     R132: C2H4 + CH2 = H + aC3H5
      RF(132) = RF(132)*C(2)
      RB(132) = RB(132)*C(20)*C(24)
C     R133: C2H4 + CH2* = H + aC3H5
      RF(133) = RF(133)*C(2)
      RB(133) = RB(133)*C(20)*C(24)
C     R134: C2H4 + CH3 = CH4 + C2H3
      RF(134) = RF(134)*C(2)*C(13)
      RB(134) = RB(134)*C(3)
C     R135: C2H4 + O2 = HO2 + C2H3
      RF(135) = RF(135)*C(2)*C(14)
      RB(135) = RB(135)*C(17)
C     R136: 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(136)
      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(136) = RF(136) * PCOR
      RB(136) = RB(136) * PCOR
      RF(136) = RF(136)*C(20)
      RB(136) = RB(136)*C(8)
C     R137: O + C2H5 = CH3 + CH2O
      RF(137) = RF(137)*C(15)
      RB(137) = RB(137)*C(13)*C(21)
C     R138: O2 + C2H5 = C2H4 + HO2
      RF(138) = RF(138)*C(14)
      RB(138) = RB(138)*C(2)*C(17)
C     R139: HO2 + C2H5 = CH3 + OH + CH2O
      RF(139) = RF(139)*C(17)
      RB(139) = RB(139)*C(13)*C(16)*C(21)
C     R140: 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(140)
      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(140) = RF(140) * PCOR
      RB(140) = RB(140) * PCOR
      RB(140) = RB(140)*C(6)
C     R141: C2H6 + H = H2 + C2H5
      RF(141) = RF(141)*C(8)*C(20)
      RB(141) = RB(141)*C(7)
C     R142: C2H6 + O = OH + C2H5
      RF(142) = RF(142)*C(8)*C(15)
      RB(142) = RB(142)*C(16)
C     R143: C2H6 + OH = H2O + C2H5
      RF(143) = RF(143)*C(8)*C(16)
      RB(143) = RB(143)*C(18)
C     R144: 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(144)
      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(144) = RF(144) * PCOR
      RB(144) = RB(144) * PCOR
      RF(144) = RF(144)*C(20)*C(24)
      RB(144) = RB(144)*C(4)
C     R145: HO2 + aC3H5 = C3H6 + O2
      RF(145) = RF(145)*C(17)*C(24)
      RB(145) = RB(145)*C(4)*C(14)
C     R146: HO2 + aC3H5 = OH + CH2O + C2H3
      RF(146) = RF(146)*C(17)*C(24)
      RB(146) = RB(146)*C(16)*C(21)
C     R147: HCO + aC3H5 = C3H6 + CO
      RF(147) = RF(147)*C(24)
      RB(147) = RB(147)*C(4)*C(9)
C     R148: 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(148)
      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(148) = RF(148) * PCOR
      RB(148) = RB(148) * PCOR
      RF(148) = RF(148)*C(13)*C(24)
      RB(148) = RB(148)*C(6)
C     R149: C3H6 + H = C2H4 + CH3
      RF(149) = RF(149)*C(4)*C(20)
      RB(149) = RB(149)*C(2)*C(13)
C     R150: C3H6 + H = H2 + aC3H5
      RF(150) = RF(150)*C(4)*C(20)
      RB(150) = RB(150)*C(7)*C(24)
C     R151: C3H6 + O = CH3 + H + CH2CO
      RF(151) = RF(151)*C(4)*C(15)
      RB(151) = RB(151)*C(13)*C(20)*C(23)
C     R152: C3H6 + O = HCO + C2H5
      RF(152) = RF(152)*C(4)*C(15)
C     R153: C3H6 + O = OH + aC3H5
      RF(153) = RF(153)*C(4)*C(15)
      RB(153) = RB(153)*C(16)*C(24)
C     R154: C3H6 + OH = H2O + aC3H5
      RF(154) = RF(154)*C(4)*C(16)
      RB(154) = RB(154)*C(18)*C(24)
C     R155: C3H6 + CH3 = CH4 + aC3H5
      RF(155) = RF(155)*C(4)*C(13)
      RB(155) = RB(155)*C(3)*C(24)
C     R156: C4H81 + H = C2H4 + C2H5
      RF(156) = RF(156)*C(6)*C(20)
      RB(156) = RB(156)*C(2)
C     R157: C4H81 + H = C3H6 + CH3
      RF(157) = RF(157)*C(6)*C(20)
      RB(157) = RB(157)*C(4)*C(13)
C     R158: iC4H8 + H = C3H6 + CH3
      RF(158) = RF(158)*C(5)*C(20)
      RB(158) = RB(158)*C(4)*C(13)
C     R159: iC4H8 + O = 2CH3 + CH2CO
      RF(159) = RF(159)*C(5)*C(15)
      RB(159) = RB(159)*C(13)*C(13)*C(23)
C     R160: CH3 + C6H5 = C6H5CH3
      RF(160) = RF(160)*C(13)
      RB(160) = RB(160)*C(12)
C     R161: C6H5CH3 + O2 = HO2 + C6H5CH2
      RF(161) = RF(161)*C(12)*C(14)
      RB(161) = RB(161)*C(17)*C(28)
C     R162: C6H5CH3 + OH = H2O + C6H5CH2
      RF(162) = RF(162)*C(12)*C(16)
      RB(162) = RB(162)*C(18)*C(28)
C     R163: C6H5CH3 + H = H2 + C6H5CH2
      RF(163) = RF(163)*C(12)*C(20)
      RB(163) = RB(163)*C(7)*C(28)
C     R164: C6H5CH3 + H = C6H6 + CH3
      RF(164) = RF(164)*C(12)*C(20)
      RB(164) = RB(164)*C(10)*C(13)
C     R165: C6H5CH3 + CH3 = CH4 + C6H5CH2
      RF(165) = RF(165)*C(12)*C(13)
      RB(165) = RB(165)*C(3)*C(28)
C     R166: H + C6H5CH2 = C6H5CH3
      CTB = CTOT+C(7)+5.D0*C(18)
     * +C(3)+5.D-1*C(9)
     * +C(22)+2.D0*C(8)
      PR = RKLOW(20) * CTB / RF(166)
      PCOR = PR / (1.0 + PR)
      PRLOG = LOG10(MAX(PR,SMALL))
      FCENT = 5.69D-1*EXP(-T/3.83D2)
     *      + 4.31D-1*EXP(-T/1.52D2)
     *     + EXP(-4.73D3/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(166) = RF(166) * PCOR
      RB(166) = RB(166) * PCOR
      RF(166) = RF(166)*C(20)*C(28)
      RB(166) = RB(166)*C(12)
C     R167: H + C6H5CH2 = CH3 + C6H5
      RF(167) = RF(167)*C(20)*C(28)
      RB(167) = RB(167)*C(13)
C     R168: O + C6H5CH2 = H + C6H5CHO
      RF(168) = RF(168)*C(15)*C(28)
      RB(168) = RB(168)*C(20)*C(30)
C     R169: HO2 + C6H5CH2 = OH + H + C6H5CHO
      RF(169) = RF(169)*C(17)*C(28)
      RB(169) = RB(169)*C(16)*C(20)*C(30)
C     R170: HCO + C6H5 = C6H5CHO
      RB(170) = RB(170)*C(30)
C     R171: C6H5CHO = H + C6H5CO
      RF(171) = RF(171)*C(30)
      RB(171) = RB(171)*C(20)
C     R172: O2 + C6H5CHO = HO2 + C6H5CO
      RF(172) = RF(172)*C(14)*C(30)
      RB(172) = RB(172)*C(17)
C     R173: OH + C6H5CHO = H2O + C6H5CO
      RF(173) = RF(173)*C(16)*C(30)
      RB(173) = RB(173)*C(18)
C     R174: H + C6H5CHO = H2 + C6H5CO
      RF(174) = RF(174)*C(20)*C(30)
      RB(174) = RB(174)*C(7)
C     R175: H + C6H5CHO = C6H6 + HCO
      RF(175) = RF(175)*C(20)*C(30)
      RB(175) = RB(175)*C(10)
C     R176: O + C6H5CHO = OH + C6H5CO
      RF(176) = RF(176)*C(15)*C(30)
      RB(176) = RB(176)*C(16)
C     R177: C6H5CH2 + C6H5CHO = C6H5CH3 + C6H5CO
      RF(177) = RF(177)*C(28)*C(30)
      RB(177) = RB(177)*C(12)
C     R178: CH3 + C6H5CHO = CH4 + C6H5CO
      RF(178) = RF(178)*C(13)*C(30)
      RB(178) = RB(178)*C(3)
C     R179: H2O2 + C6H5CO = HO2 + C6H5CHO
      RF(179) = RF(179)*C(19)
      RB(179) = RB(179)*C(17)*C(30)
C     R180: C6H5CO = CO + C6H5
      RB(180) = RB(180)*C(9)
C     R181: H + C6H5 = C6H6
      CTB = CTOT+C(7)+5.D0*C(18)
     * +C(3)+5.D-1*C(9)
     * +C(22)
      PR = RKLOW(21) * CTB / RF(181)
      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(181) = RF(181) * PCOR
      RB(181) = RB(181) * PCOR
      RF(181) = RF(181)*C(20)
      RB(181) = RB(181)*C(10)
C     R182: C6H6 + OH = H2O + C6H5
      RF(182) = RF(182)*C(10)*C(16)
      RB(182) = RB(182)*C(18)
C     R183: C6H6 + O = H + C6H5O
      RF(183) = RF(183)*C(10)*C(15)
      RB(183) = RB(183)*C(20)*C(29)
C     R184: C6H6 + O = HCO + C5H5
      RF(184) = RF(184)*C(10)*C(15)
      RB(184) = RB(184)*C(27)
C     R185: H2 + C6H5 = C6H6 + H
      RF(185) = RF(185)*C(7)
      RB(185) = RB(185)*C(10)*C(20)
C     R186: O2 + C6H5 = O + C6H5O
      RF(186) = RF(186)*C(14)
      RB(186) = RB(186)*C(15)*C(29)
C     R187: O + C6H5 = CO + C5H5
      RF(187) = RF(187)*C(15)
      RB(187) = RB(187)*C(9)*C(27)
C     R188: OH + C6H5 = H + C6H5O
      RF(188) = RF(188)*C(16)
      RB(188) = RB(188)*C(20)*C(29)
C     R189: HO2 + C6H5 = C6H6 + O2
      RF(189) = RF(189)*C(17)
      RB(189) = RB(189)*C(10)*C(14)
C     R190: CH4 + C6H5 = C6H6 + CH3
      RF(190) = RF(190)*C(3)
      RB(190) = RB(190)*C(10)*C(13)
C     R191: H + C6H5O = HCO + C5H5
      RF(191) = RF(191)*C(20)*C(29)
      RB(191) = RB(191)*C(27)
C     R192: C6H5O = CO + C5H5
      RF(192) = RF(192)*C(29)
      RB(192) = RB(192)*C(9)*C(27)
C     R193: O + C5H5 = H + C5H4O
      RF(193) = RF(193)*C(15)*C(27)
      RB(193) = RB(193)*C(20)*C(26)
C     R194: O2 + C5H5 = OH + C5H4O
      RF(194) = RF(194)*C(14)*C(27)
      RB(194) = RB(194)*C(16)*C(26)
C     R195: C5H4O = CO + 2C2H2
      RF(195) = RF(195)*C(26)
      RB(195) = RB(195)*C(9)*C(11)*C(11)
C     R196: H + C5H4O = CO + c-C4H5
      RF(196) = RF(196)*C(20)*C(26)
      RB(196) = RB(196)*C(9)*C(25)
C     R197: H + c-C4H5 = C2H4 + C2H2
      RF(197) = RF(197)*C(20)*C(25)
      RB(197) = RB(197)*C(2)*C(11)
C     R198: O + c-C4H5 = C2H2 + CH2CHO
      RF(198) = RF(198)*C(15)*C(25)
      RB(198) = RB(198)*C(11)
C     R199: O2 + c-C4H5 = CH2CO + CH2CHO
      RF(199) = RF(199)*C(14)*C(25)
      RB(199) = RB(199)*C(23)
C     R200: c-C4H5 = C2H2 + C2H3
      RF(200) = RF(200)*C(25)
      RB(200) = RB(200)*C(11)
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, 31
         WDOT(K) = 0D0
      ENDDO
C
      ROP = RF(1)
      WDOT(1) = WDOT(1) -ROP
      WDOT(2) = WDOT(2) +1.72699*ROP
      WDOT(4) = WDOT(4) +0.518098*ROP
      WDOT(5) = WDOT(5) +0.0863497*ROP
      WDOT(6) = WDOT(6) +0.0690798*ROP
      WDOT(10) = WDOT(10) +0.13*ROP
      WDOT(12) = WDOT(12) +0.37*ROP
      WDOT(13) = WDOT(13) +ROP
      WDOT(20) = WDOT(20) +ROP
      ROP = RF(2)
      WDOT(1) = WDOT(1) -ROP
      WDOT(2) = WDOT(2) +1.91642*ROP
      WDOT(3) = WDOT(3) +0.01*ROP
      WDOT(4) = WDOT(4) +0.574926*ROP
      WDOT(5) = WDOT(5) +0.0958211*ROP
      WDOT(6) = WDOT(6) +0.0766569*ROP
      WDOT(7) = WDOT(7) +ROP
      WDOT(10) = WDOT(10) +0.130867*ROP
      WDOT(12) = WDOT(12) +0.372467*ROP
      WDOT(13) = WDOT(13) +0.35*ROP
      WDOT(20) = WDOT(20) -0.35*ROP
      ROP = RF(3)
      WDOT(1) = WDOT(1) -ROP
      WDOT(2) = WDOT(2) +1.91642*ROP
      WDOT(3) = WDOT(3) +1.01*ROP
      WDOT(4) = WDOT(4) +0.574926*ROP
      WDOT(5) = WDOT(5) +0.0958211*ROP
      WDOT(6) = WDOT(6) +0.0766569*ROP
      WDOT(10) = WDOT(10) +0.130867*ROP
      WDOT(12) = WDOT(12) +0.372467*ROP
      WDOT(13) = WDOT(13) -0.65*ROP
      WDOT(20) = WDOT(20) +0.65*ROP
      ROP = RF(4)
      WDOT(1) = WDOT(1) -ROP
      WDOT(2) = WDOT(2) +1.91642*ROP
      WDOT(3) = WDOT(3) +0.01*ROP
      WDOT(4) = WDOT(4) +0.574926*ROP
      WDOT(5) = WDOT(5) +0.0958211*ROP
      WDOT(6) = WDOT(6) +0.0766569*ROP
      WDOT(10) = WDOT(10) +0.130867*ROP
      WDOT(12) = WDOT(12) +0.372467*ROP
      WDOT(13) = WDOT(13) +0.35*ROP
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      WDOT(20) = WDOT(20) +0.65*ROP
      ROP = RF(5)
      WDOT(1) = WDOT(1) -ROP
      WDOT(2) = WDOT(2) +1.91642*ROP
      WDOT(3) = WDOT(3) +0.01*ROP
      WDOT(4) = WDOT(4) +0.574926*ROP
      WDOT(5) = WDOT(5) +0.0958211*ROP
      WDOT(6) = WDOT(6) +0.0766569*ROP
      WDOT(10) = WDOT(10) +0.130867*ROP
      WDOT(12) = WDOT(12) +0.372467*ROP
      WDOT(13) = WDOT(13) +0.35*ROP
      WDOT(14) = WDOT(14) -ROP
      WDOT(17) = WDOT(17) +ROP
      WDOT(20) = WDOT(20) +0.65*ROP
      ROP = RF(6)
      WDOT(1) = WDOT(1) -ROP
      WDOT(2) = WDOT(2) +1.91642*ROP
      WDOT(3) = WDOT(3) +0.01*ROP
      WDOT(4) = WDOT(4) +0.574926*ROP
      WDOT(5) = WDOT(5) +0.0958211*ROP
      WDOT(6) = WDOT(6) +0.0766569*ROP
      WDOT(10) = WDOT(10) +0.130867*ROP
      WDOT(12) = WDOT(12) +0.372467*ROP
      WDOT(13) = WDOT(13) +0.35*ROP
      WDOT(17) = WDOT(17) -ROP
      WDOT(19) = WDOT(19) +ROP
      WDOT(20) = WDOT(20) +0.65*ROP
      ROP = RF(7)
      WDOT(1) = WDOT(1) -ROP
      WDOT(2) = WDOT(2) +1.91642*ROP
      WDOT(3) = WDOT(3) +0.01*ROP
      WDOT(4) = WDOT(4) +0.574926*ROP
      WDOT(5) = WDOT(5) +0.0958211*ROP
      WDOT(6) = WDOT(6) +0.0766569*ROP
      WDOT(10) = WDOT(10) +0.130867*ROP
      WDOT(12) = WDOT(12) +0.372467*ROP
      WDOT(13) = WDOT(13) +0.35*ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(20) = WDOT(20) +0.65*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(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(15)-RB(15)
      WDOT(15) = WDOT(15) -ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(16)-RB(16)
      WDOT(14) = WDOT(14) +ROP
      WDOT(15) = WDOT(15) -2*ROP
      ROP = RF(17)-RB(17)
      WDOT(14) = WDOT(14) -ROP
      WDOT(17) = WDOT(17) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(18)-RB(18)
      WDOT(7) = WDOT(7) -ROP
      WDOT(14) = WDOT(14) -ROP
      WDOT(17) = WDOT(17) +ROP
      WDOT(20) = WDOT(20) +ROP
      ROP = RF(19)-RB(19)
      WDOT(16) = WDOT(16) -2*ROP
      WDOT(19) = WDOT(19) +ROP
      ROP = RF(20)-RB(20)
      WDOT(16) = WDOT(16) +2*ROP
      WDOT(17) = WDOT(17) -ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(21)-RB(21)
      WDOT(14) = WDOT(14) +ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(17) = WDOT(17) -ROP
      ROP = RF(22)-RB(22)
      WDOT(14) = WDOT(14) +ROP
      WDOT(17) = WDOT(17) -2*ROP
      WDOT(19) = WDOT(19) +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(16) = WDOT(16) -ROP
      WDOT(17) = WDOT(17) -ROP
      WDOT(18) = WDOT(18) +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(7) = WDOT(7) +ROP
      WDOT(17) = WDOT(17) +ROP
      WDOT(19) = WDOT(19) -ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(28)-RB(28)
      WDOT(16) = WDOT(16) +ROP
      WDOT(18) = WDOT(18) +ROP
      WDOT(19) = WDOT(19) -ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(29)-RB(29)
      WDOT(15) = WDOT(15) -ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(17) = WDOT(17) +ROP
      WDOT(19) = WDOT(19) -ROP
      ROP = RF(30)-RB(30)
      WDOT(16) = WDOT(16) -ROP
      WDOT(17) = WDOT(17) +ROP
      WDOT(18) = WDOT(18) +ROP
      WDOT(19) = WDOT(19) -ROP
      ROP = RF(31)-RB(31)
      WDOT(9) = WDOT(9) -ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(22) = WDOT(22) +ROP
      ROP = RF(32)-RB(32)
      WDOT(9) = WDOT(9) -ROP
      WDOT(16) = WDOT(16) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(22) = WDOT(22) +ROP
      ROP = RF(33)-RB(33)
      WDOT(9) = WDOT(9) -ROP
      WDOT(16) = WDOT(16) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(22) = WDOT(22) +ROP
      ROP = RF(34)-RB(34)
      WDOT(9) = WDOT(9) -ROP
      WDOT(14) = WDOT(14) -ROP
      WDOT(15) = WDOT(15) +ROP
      WDOT(22) = WDOT(22) +ROP
      ROP = RF(35)-RB(35)
      WDOT(9) = WDOT(9) -ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(17) = WDOT(17) -ROP
      WDOT(22) = WDOT(22) +ROP
      ROP = RF(36)-RB(36)
      WDOT(7) = WDOT(7) +ROP
      WDOT(9) = WDOT(9) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(37)-RB(37)
      WDOT(9) = WDOT(9) +ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(16) = WDOT(16) +ROP
      ROP = RF(38)-RB(38)
      WDOT(15) = WDOT(15) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(22) = WDOT(22) +ROP
      ROP = RF(39)-RB(39)
      WDOT(9) = WDOT(9) +ROP
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      ROP = RF(40)-RB(40)
      WDOT(9) = WDOT(9) +ROP
      WDOT(20) = WDOT(20) +ROP
      ROP = RF(41)-RB(41)
      WDOT(9) = WDOT(9) +ROP
      WDOT(20) = WDOT(20) +ROP
      ROP = RF(42)-RB(42)
      WDOT(9) = WDOT(9) +ROP
      WDOT(14) = WDOT(14) -ROP
      WDOT(17) = WDOT(17) +ROP
      ROP = RF(43)-RB(43)
      WDOT(7) = WDOT(7) -ROP
      WDOT(9) = WDOT(9) -ROP
      WDOT(21) = WDOT(21) +ROP
      ROP = RF(44)-RB(44)
      WDOT(20) = WDOT(20) -ROP
      WDOT(21) = WDOT(21) +ROP
      ROP = RF(45)-RB(45)
      WDOT(13) = WDOT(13) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(46)-RB(46)
      WDOT(15) = WDOT(15) -ROP
      WDOT(20) = WDOT(20) +ROP
      ROP = RF(47)-RB(47)
      WDOT(16) = WDOT(16) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(21) = WDOT(21) +ROP
      ROP = RF(48)-RB(48)
      WDOT(7) = WDOT(7) -ROP
      WDOT(13) = WDOT(13) +ROP
      WDOT(20) = WDOT(20) +ROP
      ROP = RF(49)-RB(49)
      WDOT(14) = WDOT(14) -ROP
      WDOT(16) = WDOT(16) +ROP
      ROP = RF(50)-RB(50)
      WDOT(14) = WDOT(14) -ROP
      WDOT(20) = WDOT(20) +2*ROP
      WDOT(22) = WDOT(22) +ROP
      ROP = RF(51)-RB(51)
      WDOT(9) = WDOT(9) -ROP
      WDOT(23) = WDOT(23) +ROP
      ROP = RF(52)-RB(52)
      ROP = RF(53)-RB(53)
      WDOT(7) = WDOT(7) +ROP
      WDOT(9) = WDOT(9) +ROP
      WDOT(15) = WDOT(15) -ROP
      ROP = RF(54)-RB(54)
      WDOT(16) = WDOT(16) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(21) = WDOT(21) +ROP
      ROP = RF(55)-RB(55)
      WDOT(7) = WDOT(7) -ROP
      WDOT(13) = WDOT(13) +ROP
      WDOT(20) = WDOT(20) +ROP
      ROP = RF(56)-RB(56)
      WDOT(9) = WDOT(9) +ROP
      WDOT(14) = WDOT(14) -ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(20) = WDOT(20) +ROP
      ROP = RF(57)-RB(57)
      WDOT(9) = WDOT(9) +ROP
      WDOT(14) = WDOT(14) -ROP
      WDOT(18) = WDOT(18) +ROP
      ROP = RF(58)-RB(58)
      ROP = RF(59)-RB(59)
      ROP = RF(60)-RB(60)
      ROP = RF(61)-RB(61)
      WDOT(9) = WDOT(9) +ROP
      WDOT(21) = WDOT(21) +ROP
      WDOT(22) = WDOT(22) -ROP
      ROP = RF(62)-RB(62)
      WDOT(20) = WDOT(20) -ROP
      WDOT(21) = WDOT(21) -ROP
      ROP = RF(63)-RB(63)
      WDOT(7) = WDOT(7) +ROP
      WDOT(20) = WDOT(20) -ROP
      WDOT(21) = WDOT(21) -ROP
      ROP = RF(64)-RB(64)
      WDOT(15) = WDOT(15) -ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(21) = WDOT(21) -ROP
      ROP = RF(65)-RB(65)
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      WDOT(21) = WDOT(21) -ROP
      ROP = RF(66)-RB(66)
      WDOT(14) = WDOT(14) -ROP
      WDOT(17) = WDOT(17) +ROP
      WDOT(21) = WDOT(21) -ROP
      ROP = RF(67)-RB(67)
      WDOT(17) = WDOT(17) -ROP
      WDOT(19) = WDOT(19) +ROP
      WDOT(21) = WDOT(21) -ROP
      ROP = RF(68)-RB(68)
      WDOT(3) = WDOT(3) +ROP
      WDOT(13) = WDOT(13) -ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(69)-RB(69)
      WDOT(13) = WDOT(13) -ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(21) = WDOT(21) +ROP
      ROP = RF(70)-RB(70)
      WDOT(13) = WDOT(13) -ROP
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      ROP = RF(71)-RB(71)
      WDOT(13) = WDOT(13) -ROP
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      ROP = RF(72)-RB(72)
      WDOT(13) = WDOT(13) -ROP
      WDOT(14) = WDOT(14) -ROP
      WDOT(15) = WDOT(15) +ROP
      ROP = RF(73)-RB(73)
      WDOT(13) = WDOT(13) -ROP
      WDOT(14) = WDOT(14) -ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(21) = WDOT(21) +ROP
      ROP = RF(74)-RB(74)
      WDOT(3) = WDOT(3) +ROP
      WDOT(13) = WDOT(13) -ROP
      WDOT(14) = WDOT(14) +ROP
      WDOT(17) = WDOT(17) -ROP
      ROP = RF(75)-RB(75)
      WDOT(13) = WDOT(13) -ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(17) = WDOT(17) -ROP
      ROP = RF(76)-RB(76)
      WDOT(3) = WDOT(3) +ROP
      WDOT(13) = WDOT(13) -ROP
      WDOT(17) = WDOT(17) +ROP
      WDOT(19) = WDOT(19) -ROP
      ROP = RF(77)-RB(77)
      WDOT(3) = WDOT(3) +ROP
      WDOT(9) = WDOT(9) +ROP
      WDOT(13) = WDOT(13) -ROP
      ROP = RF(78)-RB(78)
      WDOT(3) = WDOT(3) +ROP
      WDOT(13) = WDOT(13) -ROP
      WDOT(21) = WDOT(21) -ROP
      ROP = RF(79)-RB(79)
      WDOT(2) = WDOT(2) +ROP
      WDOT(13) = WDOT(13) -ROP
      WDOT(20) = WDOT(20) +ROP
      ROP = RF(80)-RB(80)
      WDOT(8) = WDOT(8) +ROP
      WDOT(13) = WDOT(13) -2*ROP
      ROP = RF(81)-RB(81)
      WDOT(13) = WDOT(13) -2*ROP
      WDOT(20) = WDOT(20) +ROP
      ROP = RF(82)-RB(82)
      WDOT(2) = WDOT(2) +ROP
      WDOT(9) = WDOT(9) +ROP
      WDOT(13) = WDOT(13) -ROP
      ROP = RF(83)-RB(83)
      WDOT(7) = WDOT(7) +ROP
      WDOT(20) = WDOT(20) -ROP
      WDOT(21) = WDOT(21) +ROP
      ROP = RF(84)-RB(84)
      WDOT(13) = WDOT(13) +ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(85)-RB(85)
      WDOT(18) = WDOT(18) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(86)-RB(86)
      WDOT(15) = WDOT(15) -ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(21) = WDOT(21) +ROP
      ROP = RF(87)-RB(87)
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      WDOT(21) = WDOT(21) +ROP
      ROP = RF(88)-RB(88)
      WDOT(14) = WDOT(14) -ROP
      WDOT(17) = WDOT(17) +ROP
      WDOT(21) = WDOT(21) +ROP
      ROP = RF(89)-RB(89)
      WDOT(3) = WDOT(3) -ROP
      WDOT(7) = WDOT(7) +ROP
      WDOT(13) = WDOT(13) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(90)-RB(90)
      WDOT(3) = WDOT(3) -ROP
      WDOT(13) = WDOT(13) +ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(16) = WDOT(16) +ROP
      ROP = RF(91)-RB(91)
      WDOT(3) = WDOT(3) -ROP
      WDOT(13) = WDOT(13) +ROP
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      ROP = RF(92)-RB(92)
      WDOT(3) = WDOT(3) -ROP
      WDOT(13) = WDOT(13) +2*ROP
      ROP = RF(93)-RB(93)
      WDOT(9) = WDOT(9) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(94)-RB(94)
      WDOT(9) = WDOT(9) +2*ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(20) = WDOT(20) +ROP
      ROP = RF(95)-RB(95)
      WDOT(9) = WDOT(9) +2*ROP
      WDOT(14) = WDOT(14) -ROP
      WDOT(16) = WDOT(16) +ROP
      ROP = RF(96)-RB(96)
      WDOT(11) = WDOT(11) +ROP
      WDOT(20) = WDOT(20) +ROP
      ROP = RF(97)-RB(97)
      WDOT(9) = WDOT(9) +ROP
      WDOT(11) = WDOT(11) -ROP
      WDOT(15) = WDOT(15) -ROP
      ROP = RF(98)-RB(98)
      WDOT(11) = WDOT(11) -ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(20) = WDOT(20) +ROP
      ROP = RF(99)-RB(99)
      WDOT(11) = WDOT(11) -ROP
      WDOT(16) = WDOT(16) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(23) = WDOT(23) +ROP
      ROP = RF(100)-RB(100)
      WDOT(9) = WDOT(9) +ROP
      WDOT(11) = WDOT(11) -ROP
      WDOT(13) = WDOT(13) +ROP
      WDOT(16) = WDOT(16) -ROP
      ROP = RF(101)-RB(101)
      WDOT(9) = WDOT(9) +ROP
      WDOT(11) = WDOT(11) -ROP
      ROP = RF(102)-RB(102)
      WDOT(11) = WDOT(11) -ROP
      WDOT(13) = WDOT(13) -ROP
      WDOT(24) = WDOT(24) +ROP
      ROP = RF(103)-RB(103)
      WDOT(20) = WDOT(20) -ROP
      WDOT(23) = WDOT(23) -ROP
      ROP = RF(104)-RB(104)
      WDOT(7) = WDOT(7) +ROP
      WDOT(20) = WDOT(20) -ROP
      WDOT(23) = WDOT(23) -ROP
      ROP = RF(105)-RB(105)
      WDOT(9) = WDOT(9) +ROP
      WDOT(13) = WDOT(13) +ROP
      WDOT(20) = WDOT(20) -ROP
      WDOT(23) = WDOT(23) -ROP
      ROP = RF(106)-RB(106)
      WDOT(15) = WDOT(15) -ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(23) = WDOT(23) -ROP
      ROP = RF(107)-RB(107)
      WDOT(15) = WDOT(15) -ROP
      WDOT(22) = WDOT(22) +ROP
      WDOT(23) = WDOT(23) -ROP
      ROP = RF(108)-RB(108)
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      WDOT(23) = WDOT(23) -ROP
      ROP = RF(109)-RB(109)
      WDOT(2) = WDOT(2) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(110)-RB(110)
      WDOT(7) = WDOT(7) +ROP
      WDOT(11) = WDOT(11) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(111)-RB(111)
      WDOT(15) = WDOT(15) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(23) = WDOT(23) +ROP
      ROP = RF(112)-RB(112)
      WDOT(9) = WDOT(9) +ROP
      WDOT(13) = WDOT(13) +ROP
      WDOT(15) = WDOT(15) -ROP
      ROP = RF(113)-RB(113)
      WDOT(11) = WDOT(11) +ROP
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      ROP = RF(114)-RB(114)
      WDOT(11) = WDOT(11) +ROP
      WDOT(14) = WDOT(14) -ROP
      WDOT(17) = WDOT(17) +ROP
      ROP = RF(115)-RB(115)
      WDOT(14) = WDOT(14) -ROP
      WDOT(15) = WDOT(15) +ROP
      ROP = RF(116)-RB(116)
      WDOT(14) = WDOT(14) -ROP
      WDOT(21) = WDOT(21) +ROP
      ROP = RF(117)-RB(117)
      WDOT(4) = WDOT(4) +ROP
      WDOT(13) = WDOT(13) -ROP
      ROP = RF(118)-RB(118)
      WDOT(13) = WDOT(13) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(24) = WDOT(24) +ROP
      ROP = RF(119)-RB(119)
      WDOT(9) = WDOT(9) +ROP
      WDOT(13) = WDOT(13) +ROP
      ROP = RF(120)-RB(120)
      WDOT(13) = WDOT(13) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(121)-RB(121)
      WDOT(7) = WDOT(7) +ROP
      WDOT(20) = WDOT(20) -ROP
      WDOT(23) = WDOT(23) +ROP
      ROP = RF(122)-RB(122)
      WDOT(15) = WDOT(15) -ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(23) = WDOT(23) +ROP
      ROP = RF(123)-RB(123)
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      WDOT(23) = WDOT(23) +ROP
      ROP = RF(124)-RB(124)
      WDOT(14) = WDOT(14) -ROP
      WDOT(17) = WDOT(17) +ROP
      WDOT(23) = WDOT(23) +ROP
      ROP = RF(125)-RB(125)
      WDOT(2) = WDOT(2) -ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(126)-RB(126)
      WDOT(2) = WDOT(2) -ROP
      WDOT(7) = WDOT(7) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(127)-RB(127)
      WDOT(2) = WDOT(2) -ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(16) = WDOT(16) +ROP
      ROP = RF(128)-RB(128)
      WDOT(2) = WDOT(2) -ROP
      WDOT(13) = WDOT(13) +ROP
      WDOT(15) = WDOT(15) -ROP
      ROP = RF(129)-RB(129)
      WDOT(2) = WDOT(2) -ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(21) = WDOT(21) +ROP
      ROP = RF(130)-RB(130)
      WDOT(2) = WDOT(2) -ROP
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      ROP = RF(131)-RB(131)
      WDOT(2) = WDOT(2) -ROP
      WDOT(9) = WDOT(9) +ROP
      ROP = RF(132)-RB(132)
      WDOT(2) = WDOT(2) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(24) = WDOT(24) +ROP
      ROP = RF(133)-RB(133)
      WDOT(2) = WDOT(2) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(24) = WDOT(24) +ROP
      ROP = RF(134)-RB(134)
      WDOT(2) = WDOT(2) -ROP
      WDOT(3) = WDOT(3) +ROP
      WDOT(13) = WDOT(13) -ROP
      ROP = RF(135)-RB(135)
      WDOT(2) = WDOT(2) -ROP
      WDOT(14) = WDOT(14) -ROP
      WDOT(17) = WDOT(17) +ROP
      ROP = RF(136)-RB(136)
      WDOT(8) = WDOT(8) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(137)-RB(137)
      WDOT(13) = WDOT(13) +ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(21) = WDOT(21) +ROP
      ROP = RF(138)-RB(138)
      WDOT(2) = WDOT(2) +ROP
      WDOT(14) = WDOT(14) -ROP
      WDOT(17) = WDOT(17) +ROP
      ROP = RF(139)-RB(139)
      WDOT(13) = WDOT(13) +ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(17) = WDOT(17) -ROP
      WDOT(21) = WDOT(21) +ROP
      ROP = RF(140)-RB(140)
      WDOT(6) = WDOT(6) +ROP
      ROP = RF(141)-RB(141)
      WDOT(7) = WDOT(7) +ROP
      WDOT(8) = WDOT(8) -ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(142)-RB(142)
      WDOT(8) = WDOT(8) -ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(16) = WDOT(16) +ROP
      ROP = RF(143)-RB(143)
      WDOT(8) = WDOT(8) -ROP
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      ROP = RF(144)-RB(144)
      WDOT(4) = WDOT(4) +ROP
      WDOT(20) = WDOT(20) -ROP
      WDOT(24) = WDOT(24) -ROP
      ROP = RF(145)-RB(145)
      WDOT(4) = WDOT(4) +ROP
      WDOT(14) = WDOT(14) +ROP
      WDOT(17) = WDOT(17) -ROP
      WDOT(24) = WDOT(24) -ROP
      ROP = RF(146)-RB(146)
      WDOT(16) = WDOT(16) +ROP
      WDOT(17) = WDOT(17) -ROP
      WDOT(21) = WDOT(21) +ROP
      WDOT(24) = WDOT(24) -ROP
      ROP = RF(147)-RB(147)
      WDOT(4) = WDOT(4) +ROP
      WDOT(9) = WDOT(9) +ROP
      WDOT(24) = WDOT(24) -ROP
      ROP = RF(148)-RB(148)
      WDOT(6) = WDOT(6) +ROP
      WDOT(13) = WDOT(13) -ROP
      WDOT(24) = WDOT(24) -ROP
      ROP = RF(149)-RB(149)
      WDOT(2) = WDOT(2) +ROP
      WDOT(4) = WDOT(4) -ROP
      WDOT(13) = WDOT(13) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(150)-RB(150)
      WDOT(4) = WDOT(4) -ROP
      WDOT(7) = WDOT(7) +ROP
      WDOT(20) = WDOT(20) -ROP
      WDOT(24) = WDOT(24) +ROP
      ROP = RF(151)-RB(151)
      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(152)-RB(152)
      WDOT(4) = WDOT(4) -ROP
      WDOT(15) = WDOT(15) -ROP
      ROP = RF(153)-RB(153)
      WDOT(4) = WDOT(4) -ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(24) = WDOT(24) +ROP
      ROP = RF(154)-RB(154)
      WDOT(4) = WDOT(4) -ROP
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      WDOT(24) = WDOT(24) +ROP
      ROP = RF(155)-RB(155)
      WDOT(3) = WDOT(3) +ROP
      WDOT(4) = WDOT(4) -ROP
      WDOT(13) = WDOT(13) -ROP
      WDOT(24) = WDOT(24) +ROP
      ROP = RF(156)-RB(156)
      WDOT(2) = WDOT(2) +ROP
      WDOT(6) = WDOT(6) -ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(157)-RB(157)
      WDOT(4) = WDOT(4) +ROP
      WDOT(6) = WDOT(6) -ROP
      WDOT(13) = WDOT(13) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(158)-RB(158)
      WDOT(4) = WDOT(4) +ROP
      WDOT(5) = WDOT(5) -ROP
      WDOT(13) = WDOT(13) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(159)-RB(159)
      WDOT(5) = WDOT(5) -ROP
      WDOT(13) = WDOT(13) +2*ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(23) = WDOT(23) +ROP
      ROP = RF(160)-RB(160)
      WDOT(12) = WDOT(12) +ROP
      WDOT(13) = WDOT(13) -ROP
      ROP = RF(161)-RB(161)
      WDOT(12) = WDOT(12) -ROP
      WDOT(14) = WDOT(14) -ROP
      WDOT(17) = WDOT(17) +ROP
      WDOT(28) = WDOT(28) +ROP
      ROP = RF(162)-RB(162)
      WDOT(12) = WDOT(12) -ROP
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      WDOT(28) = WDOT(28) +ROP
      ROP = RF(163)-RB(163)
      WDOT(7) = WDOT(7) +ROP
      WDOT(12) = WDOT(12) -ROP
      WDOT(20) = WDOT(20) -ROP
      WDOT(28) = WDOT(28) +ROP
      ROP = RF(164)-RB(164)
      WDOT(10) = WDOT(10) +ROP
      WDOT(12) = WDOT(12) -ROP
      WDOT(13) = WDOT(13) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(165)-RB(165)
      WDOT(3) = WDOT(3) +ROP
      WDOT(12) = WDOT(12) -ROP
      WDOT(13) = WDOT(13) -ROP
      WDOT(28) = WDOT(28) +ROP
      ROP = RF(166)-RB(166)
      WDOT(12) = WDOT(12) +ROP
      WDOT(20) = WDOT(20) -ROP
      WDOT(28) = WDOT(28) -ROP
      ROP = RF(167)-RB(167)
      WDOT(13) = WDOT(13) +ROP
      WDOT(20) = WDOT(20) -ROP
      WDOT(28) = WDOT(28) -ROP
      ROP = RF(168)-RB(168)
      WDOT(15) = WDOT(15) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(28) = WDOT(28) -ROP
      WDOT(30) = WDOT(30) +ROP
      ROP = RF(169)-RB(169)
      WDOT(16) = WDOT(16) +ROP
      WDOT(17) = WDOT(17) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(28) = WDOT(28) -ROP
      WDOT(30) = WDOT(30) +ROP
      ROP = RF(170)-RB(170)
      WDOT(30) = WDOT(30) +ROP
      ROP = RF(171)-RB(171)
      WDOT(20) = WDOT(20) +ROP
      WDOT(30) = WDOT(30) -ROP
      ROP = RF(172)-RB(172)
      WDOT(14) = WDOT(14) -ROP
      WDOT(17) = WDOT(17) +ROP
      WDOT(30) = WDOT(30) -ROP
      ROP = RF(173)-RB(173)
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      WDOT(30) = WDOT(30) -ROP
      ROP = RF(174)-RB(174)
      WDOT(7) = WDOT(7) +ROP
      WDOT(20) = WDOT(20) -ROP
      WDOT(30) = WDOT(30) -ROP
      ROP = RF(175)-RB(175)
      WDOT(10) = WDOT(10) +ROP
      WDOT(20) = WDOT(20) -ROP
      WDOT(30) = WDOT(30) -ROP
      ROP = RF(176)-RB(176)
      WDOT(15) = WDOT(15) -ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(30) = WDOT(30) -ROP
      ROP = RF(177)-RB(177)
      WDOT(12) = WDOT(12) +ROP
      WDOT(28) = WDOT(28) -ROP
      WDOT(30) = WDOT(30) -ROP
      ROP = RF(178)-RB(178)
      WDOT(3) = WDOT(3) +ROP
      WDOT(13) = WDOT(13) -ROP
      WDOT(30) = WDOT(30) -ROP
      ROP = RF(179)-RB(179)
      WDOT(17) = WDOT(17) +ROP
      WDOT(19) = WDOT(19) -ROP
      WDOT(30) = WDOT(30) +ROP
      ROP = RF(180)-RB(180)
      WDOT(9) = WDOT(9) +ROP
      ROP = RF(181)-RB(181)
      WDOT(10) = WDOT(10) +ROP
      WDOT(20) = WDOT(20) -ROP
      ROP = RF(182)-RB(182)
      WDOT(10) = WDOT(10) -ROP
      WDOT(16) = WDOT(16) -ROP
      WDOT(18) = WDOT(18) +ROP
      ROP = RF(183)-RB(183)
      WDOT(10) = WDOT(10) -ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(29) = WDOT(29) +ROP
      ROP = RF(184)-RB(184)
      WDOT(10) = WDOT(10) -ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(27) = WDOT(27) +ROP
      ROP = RF(185)-RB(185)
      WDOT(7) = WDOT(7) -ROP
      WDOT(10) = WDOT(10) +ROP
      WDOT(20) = WDOT(20) +ROP
      ROP = RF(186)-RB(186)
      WDOT(14) = WDOT(14) -ROP
      WDOT(15) = WDOT(15) +ROP
      WDOT(29) = WDOT(29) +ROP
      ROP = RF(187)-RB(187)
      WDOT(9) = WDOT(9) +ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(27) = WDOT(27) +ROP
      ROP = RF(188)-RB(188)
      WDOT(16) = WDOT(16) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(29) = WDOT(29) +ROP
      ROP = RF(189)-RB(189)
      WDOT(10) = WDOT(10) +ROP
      WDOT(14) = WDOT(14) +ROP
      WDOT(17) = WDOT(17) -ROP
      ROP = RF(190)-RB(190)
      WDOT(3) = WDOT(3) -ROP
      WDOT(10) = WDOT(10) +ROP
      WDOT(13) = WDOT(13) +ROP
      ROP = RF(191)-RB(191)
      WDOT(20) = WDOT(20) -ROP
      WDOT(27) = WDOT(27) +ROP
      WDOT(29) = WDOT(29) -ROP
      ROP = RF(192)-RB(192)
      WDOT(9) = WDOT(9) +ROP
      WDOT(27) = WDOT(27) +ROP
      WDOT(29) = WDOT(29) -ROP
      ROP = RF(193)-RB(193)
      WDOT(15) = WDOT(15) -ROP
      WDOT(20) = WDOT(20) +ROP
      WDOT(26) = WDOT(26) +ROP
      WDOT(27) = WDOT(27) -ROP
      ROP = RF(194)-RB(194)
      WDOT(14) = WDOT(14) -ROP
      WDOT(16) = WDOT(16) +ROP
      WDOT(26) = WDOT(26) +ROP
      WDOT(27) = WDOT(27) -ROP
      ROP = RF(195)-RB(195)
      WDOT(9) = WDOT(9) +ROP
      WDOT(11) = WDOT(11) +2*ROP
      WDOT(26) = WDOT(26) -ROP
      ROP = RF(196)-RB(196)
      WDOT(9) = WDOT(9) +ROP
      WDOT(20) = WDOT(20) -ROP
      WDOT(25) = WDOT(25) +ROP
      WDOT(26) = WDOT(26) -ROP
      ROP = RF(197)-RB(197)
      WDOT(2) = WDOT(2) +ROP
      WDOT(11) = WDOT(11) +ROP
      WDOT(20) = WDOT(20) -ROP
      WDOT(25) = WDOT(25) -ROP
      ROP = RF(198)-RB(198)
      WDOT(11) = WDOT(11) +ROP
      WDOT(15) = WDOT(15) -ROP
      WDOT(25) = WDOT(25) -ROP
      ROP = RF(199)-RB(199)
      WDOT(14) = WDOT(14) -ROP
      WDOT(23) = WDOT(23) +ROP
      WDOT(25) = WDOT(25) -ROP
      ROP = RF(200)-RB(200)
      WDOT(11) = WDOT(11) +ROP
      WDOT(25) = WDOT(25) -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(140) = 0.D0
      RB(152) = 0.D0
      RF(170) = 0.D0
C
C     CH2
      DEN = +RF( 45) +RF( 46) +RF( 47) +RF( 48) +RF( 49)
     *  +RF( 50) +RF( 51) +RF( 79) +RF(132) +RB( 52) +RB( 58)
     *  +RB( 59) +RB( 60) +RB( 70) +RB( 97) +RB(107) +RB(129)
      A1_0 = ( +RB( 45) +RB( 47) +RB( 48) +RB( 50) +RB( 51)
     *  +RF( 70) +RB( 79) +RF( 97) +RF(107) +RF(129) +RB(132) )/DEN
      A1_2 = ( +RF( 52) +RF( 58) +RF( 59) +RF( 60) )/DEN
      A1_3 = ( +RB( 46) +RB( 49) )/DEN
C     CH2*
      DEN = +RF( 52) +RF( 53) +RF( 54) +RF( 55) +RF( 56)
     *  +RF( 57) +RF( 58) +RF( 59) +RF( 60) +RF( 61) +RF( 92)
     *  +RF(133) +RB( 71) +RB( 85) +RB( 93)
      A2_0 = ( +RB( 53) +RB( 54) +RB( 55) +RB( 56) +RB( 57)
     *  +RB( 61) +RF( 71) +RB( 92) +RB(133) )/DEN
      A2_1 = ( +RB( 52) +RB( 58) +RB( 59) +RB( 60) )/DEN
      A2_4 = ( +RF( 85) )/DEN
      A2_7 = ( +RF( 93) )/DEN
C     HCO
      DEN = +RF( 36) +RF( 37) +RF( 38) +RF( 39) +RF( 40)
     *  +RF( 41) +RF( 42) +RF( 44) +RF( 77) +RF(101) +RF(131)
     *  +RF(147) +RB( 46) +RB( 49) +RB( 63) +RB( 64) +RB( 65)
     *  +RB( 66) +RB( 67) +RB( 78) +RB(116) +RB(120) +RB(128)
     *  +RB(175) +RB(184) +RB(191)
      A3_0 = ( +RB( 36) +RB( 37) +RB( 38) +RB( 39) +RB( 40)
     *  +RB( 41) +RB( 42) +RB( 44) +RF( 63) +RF( 64) +RF( 65)
     *  +RF( 66) +RF( 67) +RB( 77) +RF( 78) +RF(128) +RB(147)
     *  +RF(152) +RB(170) +RF(175) +RF(184) +RF(191) )/DEN
      A3_1 = ( +RF( 46) +RF( 49) )/DEN
      A3_5 = ( +RB(101) +RF(116) )/DEN
      A3_6 = ( +RB(131) )/DEN
      A3_8 = ( +RF(120) )/DEN
C     CH3O
      DEN = +RF( 83) +RF( 84) +RF( 85) +RF( 86) +RF( 87)
     *  +RF( 88) +RB( 62) +RB( 72) +RB( 75)
      A4_0 = ( +RF( 62) +RF( 72) +RF( 75) +RB( 83) +RB( 84)
     *  +RB( 86) +RB( 87) +RB( 88) )/DEN
      A4_2 = ( +RB( 85) )/DEN
C     C2H3
      DEN = +RF( 96) +RF(109) +RF(110) +RF(111) +RF(112)
     *  +RF(113) +RF(114) +RF(115) +RF(116) +RF(117) +RF(118)
     *  +RB(101) +RB(126) +RB(127) +RB(130) +RB(134) +RB(135)
     *  +RB(146) +RB(200)
      A5_0 = ( +RB( 96) +RB(109) +RB(110) +RB(111) +RB(112)
     *  +RB(113) +RB(114) +RB(117) +RB(118) +RF(126) +RF(127)
     *  +RF(130) +RF(134) +RF(135) +RB(140) +RF(146) +RF(200) )/DEN
      A5_3 = ( +RF(101) +RB(116) )/DEN
      A5_8 = ( +RB(115) )/DEN
C     C2H5
      DEN = +RF(136) +RF(137) +RF(138) +RF(139) +RB( 81)
     *  +RB(125) +RB(131) +RB(141) +RB(142) +RB(143) +RB(156)
      A6_0 = ( +RF( 81) +RF(125) +RB(136) +RB(137) +RB(138)
     *  +RB(139) +RB(140) +RF(141) +RF(142) +RF(143) +RF(152)
     *  +RF(156) )/DEN
      A6_3 = ( +RF(131) )/DEN
C     HCCO
      DEN = +RF( 82) +RF( 93) +RF( 94) +RF( 95) +RB( 98)
     *  +RB(104) +RB(106) +RB(108)
      A7_0 = ( +RB( 82) +RB( 94) +RB( 95) +RF( 98) +RF(104)
     *  +RF(106) +RF(108) )/DEN
      A7_2 = ( +RB( 93) )/DEN
C     CH2CHO
      DEN = +RF(119) +RF(120) +RF(121) +RF(122) +RF(123)
     *  +RF(124) +RB(103) +RB(115) +RB(198) +RB(199)
      A8_0 = ( +RF(103) +RB(119) +RB(121) +RB(122) +RB(123)
     *  +RB(124) +RF(198) +RF(199) )/DEN
      A8_3 = ( +RB(120) )/DEN
      A8_5 = ( +RF(115) )/DEN
C     C6H5
      DEN = +RF(160) +RF(181) +RF(185) +RF(186) +RF(187)
     *  +RF(188) +RF(189) +RF(190) +RB(167) +RB(180) +RB(182)
      A9_0 = ( +RB(160) +RF(167) +RB(170) +RB(181) +RF(182)
     *  +RB(185) +RB(186) +RB(187) +RB(188) +RB(189) +RB(190) )/DEN
      A9_10 = ( +RF(180) )/DEN
C     C6H5CO
      DEN = +RF(179) +RF(180) +RB(171) +RB(172) +RB(173)
     *  +RB(174) +RB(176) +RB(177) +RB(178)
      A10_0 = ( +RF(171) +RF(172) +RF(173) +RF(174) +RF(176)
     *  +RF(177) +RF(178) +RB(179) )/DEN
      A10_9 = ( +RB(180) )/DEN
C
      A3_0 = A3_0 + A3_6*A6_0
      DEN = 1 -A3_6*A6_3
      A3_0 = A3_0/DEN
      A3_8 = A3_8/DEN
      A3_5 = A3_5/DEN
      A3_1 = A3_1/DEN
      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
      DEN = 1 -A2_7*A7_2
      A2_0 = A2_0/DEN
      A2_1 = A2_1/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
      A3_8 = A3_8/DEN
      A3_5 = A3_5/DEN
      A2_0 = A2_0 + A2_1*A1_0
      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_5*A5_0
      A3_8 = A3_8 + A3_5*A5_8
      DEN = 1 -A3_5*A5_3
      A3_0 = A3_0/DEN
      A3_2 = A3_2/DEN
      A3_8 = A3_8/DEN
      A8_0 = A8_0 + A8_5*A5_0
      A8_3 = A8_3 + A8_5*A5_3
      DEN = 1 -A8_5*A5_8
      A8_0 = A8_0/DEN
      A8_3 = A8_3/DEN
      A3_0 = A3_0 + A3_8*A8_0
      DEN = 1 -A3_8*A8_3
      A3_0 = A3_0/DEN
      A3_2 = A3_2/DEN
      A3_0 = A3_0 + A3_2*A2_0
      DEN = 1 -A3_2*A2_3
      A3_0 = A3_0/DEN
      XQ(3) = A3_0
      XQ(2) = A2_0 +A2_3*XQ(3)
      XQ(8) = A8_0 +A8_3*XQ(3)
      XQ(5) = A5_0 +A5_3*XQ(3) +A5_8*XQ(8)
      XQ(1) = A1_0 +A1_3*XQ(3) +A1_2*XQ(2)
      XQ(7) = A7_0 +A7_2*XQ(2)
      XQ(4) = A4_0 +A4_2*XQ(2)
      XQ(6) = A6_0 +A6_3*XQ(3)
      A9_0 = A9_0 + A9_10*A10_0
      DEN = 1 -A9_10*A10_9
      A9_0 = A9_0/DEN
      XQ(9) = A9_0
      XQ(10) = A10_0 +A10_9*XQ(9)
C
      RF( 36) = RF( 36)*XQ( 3)
      RF( 37) = RF( 37)*XQ( 3)
      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( 44) = RF( 44)*XQ( 3)
      RF( 45) = RF( 45)*XQ( 1)
      RF( 46) = RF( 46)*XQ( 1)
      RB( 46) = RB( 46)*XQ( 3)
      RF( 47) = RF( 47)*XQ( 1)
      RF( 48) = RF( 48)*XQ( 1)
      RF( 49) = RF( 49)*XQ( 1)
      RB( 49) = RB( 49)*XQ( 3)
      RF( 50) = RF( 50)*XQ( 1)
      RF( 51) = RF( 51)*XQ( 1)
      RF( 52) = RF( 52)*XQ( 2)
      RB( 52) = RB( 52)*XQ( 1)
      RF( 53) = RF( 53)*XQ( 2)
      RF( 54) = RF( 54)*XQ( 2)
      RF( 55) = RF( 55)*XQ( 2)
      RF( 56) = RF( 56)*XQ( 2)
      RF( 57) = RF( 57)*XQ( 2)
      RF( 58) = RF( 58)*XQ( 2)
      RB( 58) = RB( 58)*XQ( 1)
      RF( 59) = RF( 59)*XQ( 2)
      RB( 59) = RB( 59)*XQ( 1)
      RF( 60) = RF( 60)*XQ( 2)
      RB( 60) = RB( 60)*XQ( 1)
      RF( 61) = RF( 61)*XQ( 2)
      RB( 62) = RB( 62)*XQ( 4)
      RB( 63) = RB( 63)*XQ( 3)
      RB( 64) = RB( 64)*XQ( 3)
      RB( 65) = RB( 65)*XQ( 3)
      RB( 66) = RB( 66)*XQ( 3)
      RB( 67) = RB( 67)*XQ( 3)
      RB( 70) = RB( 70)*XQ( 1)
      RB( 71) = RB( 71)*XQ( 2)
      RB( 72) = RB( 72)*XQ( 4)
      RB( 75) = RB( 75)*XQ( 4)
      RF( 77) = RF( 77)*XQ( 3)
      RB( 78) = RB( 78)*XQ( 3)
      RF( 79) = RF( 79)*XQ( 1)
      RB( 81) = RB( 81)*XQ( 6)
      RF( 82) = RF( 82)*XQ( 7)
      RF( 83) = RF( 83)*XQ( 4)
      RF( 84) = RF( 84)*XQ( 4)
      RF( 85) = RF( 85)*XQ( 4)
      RB( 85) = RB( 85)*XQ( 2)
      RF( 86) = RF( 86)*XQ( 4)
      RF( 87) = RF( 87)*XQ( 4)
      RF( 88) = RF( 88)*XQ( 4)
      RF( 92) = RF( 92)*XQ( 2)
      RF( 93) = RF( 93)*XQ( 7)
      RB( 93) = RB( 93)*XQ( 2)
      RF( 94) = RF( 94)*XQ( 7)
      RF( 95) = RF( 95)*XQ( 7)
      RF( 96) = RF( 96)*XQ( 5)
      RB( 97) = RB( 97)*XQ( 1)
      RB( 98) = RB( 98)*XQ( 7)
      RF(101) = RF(101)*XQ( 3)
      RB(101) = RB(101)*XQ( 5)
      RB(103) = RB(103)*XQ( 8)
      RB(104) = RB(104)*XQ( 7)
      RB(106) = RB(106)*XQ( 7)
      RB(107) = RB(107)*XQ( 1)
      RB(108) = RB(108)*XQ( 7)
      RF(109) = RF(109)*XQ( 5)
      RF(110) = RF(110)*XQ( 5)
      RF(111) = RF(111)*XQ( 5)
      RF(112) = RF(112)*XQ( 5)
      RF(113) = RF(113)*XQ( 5)
      RF(114) = RF(114)*XQ( 5)
      RF(115) = RF(115)*XQ( 5)
      RB(115) = RB(115)*XQ( 8)
      RF(116) = RF(116)*XQ( 5)
      RB(116) = RB(116)*XQ( 3)
      RF(117) = RF(117)*XQ( 5)
      RF(118) = RF(118)*XQ( 5)
      RF(119) = RF(119)*XQ( 8)
      RF(120) = RF(120)*XQ( 8)
      RB(120) = RB(120)*XQ( 3)
      RF(121) = RF(121)*XQ( 8)
      RF(122) = RF(122)*XQ( 8)
      RF(123) = RF(123)*XQ( 8)
      RF(124) = RF(124)*XQ( 8)
      RB(125) = RB(125)*XQ( 6)
      RB(126) = RB(126)*XQ( 5)
      RB(127) = RB(127)*XQ( 5)
      RB(128) = RB(128)*XQ( 3)
      RB(129) = RB(129)*XQ( 1)
      RB(130) = RB(130)*XQ( 5)
      RF(131) = RF(131)*XQ( 3)
      RB(131) = RB(131)*XQ( 6)
      RF(132) = RF(132)*XQ( 1)
      RF(133) = RF(133)*XQ( 2)
      RB(134) = RB(134)*XQ( 5)
      RB(135) = RB(135)*XQ( 5)
      RF(136) = RF(136)*XQ( 6)
      RF(137) = RF(137)*XQ( 6)
      RF(138) = RF(138)*XQ( 6)
      RF(139) = RF(139)*XQ( 6)
      RB(141) = RB(141)*XQ( 6)
      RB(142) = RB(142)*XQ( 6)
      RB(143) = RB(143)*XQ( 6)
      RB(146) = RB(146)*XQ( 5)
      RF(147) = RF(147)*XQ( 3)
      RB(156) = RB(156)*XQ( 6)
      RF(160) = RF(160)*XQ( 9)
      RB(167) = RB(167)*XQ( 9)
      RB(171) = RB(171)*XQ(10)
      RB(172) = RB(172)*XQ(10)
      RB(173) = RB(173)*XQ(10)
      RB(174) = RB(174)*XQ(10)
      RB(175) = RB(175)*XQ( 3)
      RB(176) = RB(176)*XQ(10)
      RB(177) = RB(177)*XQ(10)
      RB(178) = RB(178)*XQ(10)
      RF(179) = RF(179)*XQ(10)
      RF(180) = RF(180)*XQ(10)
      RB(180) = RB(180)*XQ( 9)
      RF(181) = RF(181)*XQ( 9)
      RB(182) = RB(182)*XQ( 9)
      RB(184) = RB(184)*XQ( 3)
      RF(185) = RF(185)*XQ( 9)
      RF(186) = RF(186)*XQ( 9)
      RF(187) = RF(187)*XQ( 9)
      RF(188) = RF(188)*XQ( 9)
      RF(189) = RF(189)*XQ( 9)
      RF(190) = RF(190)*XQ( 9)
      RB(191) = RB(191)*XQ( 3)
      RB(198) = RB(198)*XQ( 8)
      RB(199) = RB(199)*XQ( 8)
      RB(200) = RB(200)*XQ( 5)
C
      END