Merge pull request #22 from wirc-sjsu/develop-w21
[WRF-Fire-merge.git] / chem / KPP / mechanisms / cb05_sorg_aq / cb05_sorg_aq.eqn
blob25f4f78fbbb10c8fc4af553fab14dc25c0863d92
1 #EQUATIONS { CB05 MADE/SORG Extension, Implemented by NCSU, Sep 2014}
2 //*********************************************************************
3 //The former rate constant corresponds that for the following reaction in original version of CB05 in CMAQ:
4 //< R34>  HO2  + HO2         = H2O2. But CMAQ uses an older rate constant of (EP3(2.3D-13, 600.0, 1.7D-33, 1000.0, TEMP, C_M).
5 //An updated rate constant of (EP3(2.3D-13, 600.0, 1.7D-33, 1000.0, TEMP, C_M) * ARR2(1.4D-21, 2200.0_dp, TEMP) comes from
6 //Jacobson's book (2005), Table B.4, reaction 33: HO2  + HO2         = H2O2 + O2.
7 //
8 //Additions of Cl equations from CMAQ's CB05 with Cl extension, and also Hg reactions added
9 //*********************************************************************
10 <  R1>  NO2  + hv          = NO        + O            : j(Pj_no2);
11 <  R2>  O   {+ O2    + M}  = O3       {+ M}           : 0.21*C_M*C_M*ARR(6.0D-34, 0.0_dp, -2.4_dp, TEMP);
12 <  R3>  O3   + NO          = NO2      {+ O2}          : ARR2(3.0D-12, 1500.0_dp, TEMP);
13 <  R4>  O    + NO2         = NO       {+ O2}          : ARR2(5.6D-12, -180.0_dp, TEMP);
14 <  R5>  O    + NO2  {+ M}  = NO3      {+ M}           : TROE(2.5D-31, 1.8_dp, 2.2D-11, 0.7_dp, TEMP, C_M);
15 <  R6>  O    + NO   {+ M}  = NO2      {+ M}           : TROE(9.0D-32, 1.5_dp, 3.0D-11, 0.0_dp, TEMP, C_M);
16 <  R7>  NO2  + O3          = NO3      {+ O2}          : ARR2(1.2D-13, 2450.0_dp, TEMP);
17 <  R8>  O3   + hv          = O        {+ O2}          : j(Pj_o33p);
18 <  R9>  O3   + hv          = O1D      {+ O2}          : j(Pj_o31d);
19 < R10>  O1D  + M           = O         + M            : ARR2(2.1D-11, -102.0_dp, TEMP);
20 < R11>  O1D  + H2O         = 2 OH                     : 2.2D-10;
21 < R12>  O3   + OH          = HO2      {+ O2}          : ARR2(1.7D-12, 940.0_dp, TEMP);
22 < R13>  O3   + HO2         = OH       {+ 2 O2}        : ARR2(1.0D-14, 490.0_dp, TEMP);
23 < R14>  NO3  + hv          = NO2       + O            : j(Pj_no3o);
24 < R15>  NO3  + hv          = NO       {+ O2}          : j(Pj_no3o2);
25 < R16>  NO3  + NO          = 2 NO2                    : ARR2(1.5D-11, -170.0_dp, TEMP);
26 < R17>  NO3  + NO2         = NO        + NO2   {+ O2} : ARR2(4.5D-14, 1260.0_dp, TEMP);
27 < R18>  NO3  + NO2  {+ M}  = N2O5     {+ M}           : TROE(2.0D-30, 4.4_dp, 1.4D-12, 0.7_dp, TEMP, C_M);
28 < R19>  N2O5 + H2O         = 2 HNO3                   : 2.5D-22;
29 < R20>  N2O5 + H2O   + H2O = 2 HNO3                   : 1.8D-39;
30 < R21>  N2O5               = NO3       + NO2          : FALL(1.0D-03, 11000.0_dp, -3.5_dp, 9.7D14, 11080.0_dp, 0.1_dp, 0.45_dp, TEMP, C_M);
31 < R22>  NO   + NO {+0.21M} = 2 NO2                    : ARR2(3.3D-39, -530.0_dp, TEMP)*0.21*C_M; 
32 < R23>  NO   + NO2   + H2O = 2 HONO                   : 5.0D-40;
33 < R24>  NO   + OH          = HONO                     : TROE(7.0D-31, 2.6_dp, 3.6D-11, 0.1_dp, TEMP, C_M);
34 < R25>  HONO + hv          = NO        + OH           : j(Pj_hno2);
35 < R26>  OH   + HONO        = NO2      {+ H2O}         : ARR2(1.8D-11, 390.0_dp, TEMP);
36 < R27>  HONO + HONO        = NO        + NO2          : 1.0D-20;
37 < R28>  NO2  + OH          = HNO3                     : TROE(2.0D-30, 3.0_dp, 2.5D-11, 0.0_dp, TEMP, C_M);
38 < R29>  OH   + HNO3        = NO3                      : EP2(2.4D-14, -460.0_dp, 2.7D-17, -2199.0_dp, 6.5D-34, -1335.0_dp, TEMP, C_M);
39 < R30>  HO2  + NO          = OH        + NO2          : ARR2(3.5D-12, -250.0_dp, TEMP);
40 < R31>  HO2  + NO2         = PNA                      : TROE(1.8D-31, 3.2_dp, 4.7D-12, 0.0_dp, TEMP, C_M);
41 < R32>  PNA                = HO2       + NO2          : FALL(4.1D-05, 10650.0_dp, 0.0_dp, 4.8D15, 11170.0_dp, 0.0_dp, 0.6_dp, TEMP, C_M);
42 < R33>  OH   + PNA         = NO2                      : ARR2(1.3D-12, -380.0_dp, TEMP);
43 < R34>  HO2  + HO2         = H2O2                     : EP3(2.3D-13, -600.0_dp, 1.7D-33, -1000.0_dp, TEMP, C_M);
44 < R35>  HO2  + HO2  + H2O  = H2O2                     : EP3(3.22D-34, -2800.0_dp, 2.38D-54, -3200.0_dp, TEMP, C_M);
45 < R36>  H2O2 + hv          = 2 OH                     : j(Pj_h2o2);
46 < R37>  OH   + H2O2        = HO2      {+ H2O}         : ARR2(2.9D-12, 160.0_dp, TEMP);
47 < R38>  O1D  + H2          = OH        + HO2          : 1.1D-10;
48 < R39>  OH   + H2          = HO2                      : ARR2(5.5D-12, 2000.0_dp, TEMP);
49 < R40>  OH   + O           = HO2                      : ARR2(2.2D-11, -120.0_dp, TEMP);
50 < R41>  OH   + OH          = O                        : ARR2(4.2D-12, 240.0_dp, TEMP);
51 < R42>  OH   + OH          = H2O2                     : TROE(6.9D-31, 1.0_dp, 2.6D-11, 0.0_dp, TEMP, C_M);
52 < R43>  OH   + HO2         = H2O      {+ O2}          : ARR2(4.8D-11, -250.0_dp, TEMP);
53 < R44>  HO2  + O           = OH       {+ O2}          : ARR2(3.0D-11, -200.0_dp, TEMP);
54 < R45>  H2O2 + O           = OH        + HO2          : ARR2(1.4D-12, 2000.0_dp, TEMP);
55 < R46>  NO3  + O           = NO2      {+ O2}          : 1.0D-11;
56 < R47>  NO3  + OH          = HO2       + NO2          : 2.2D-11;
57 < R48>  NO3  + HO2         = HNO3     {+ O2}          : 3.5D-12;
58 < R49>  NO3  + O3          = NO2      {+ 2 O2}        : 1.0D-17;
59 < R50>  NO3  + NO3         = 2 NO2    {+ O2}          : ARR2(8.5D-13, 2450.0_dp, TEMP);
60 < R51>  PNA  + hv          = 0.61 HO2  + 0.61 NO2
61                            + 0.39 OH   + 0.39 NO3     : j(Pj_hno4);
62 < R52>  HNO3 + hv          = OH        + NO2          : j(Pj_hno3);
63 < R53>  N2O5 + hv          = NO2       + NO3          : j(Pj_n2o5);
64 < R54>  XO2  + NO          = NO2      {+ XO}          : ARR2(2.6D-12, -365.0_dp, TEMP);
65 < R55>  XO2N + NO          = NTR                      : ARR2(2.6D-12, -365.0_dp, TEMP);
66 < R56>  XO2  + HO2         = ROOH                     : ARR2(7.5D-13, -700.0_dp, TEMP);
67 < R57>  XO2N + HO2         = ROOH                     : ARR2(7.5D-13, -700.0_dp, TEMP);
68 < R58>  XO2  + XO2         = 0.21 M                   : 6.8D-14;
69 < R59>  XO2N + XO2N        = 0.21 M                   : 6.8D-14;
70 < R60>  XO2  + XO2N        = 0.21 M                   : 6.8D-14;
71 < R61>  NTR  + OH          = HNO3      + HO2
72                            + 0.33 FORM + 0.33 ALD2
73                            + 0.33 ALDX - 0.66 PAR     : ARR2(5.9D-13, 360.0_dp, TEMP);
74 < R62>  NTR  + hv          = NO2       + HO2
75                            + 0.33 FORM + 0.33 ALD2
76                            + 0.33 ALDX - 0.66 PAR     : 1.0D-4*j(Pj_no2);
77 < R63>  ROOH + OH          = XO2       + 0.50 ALD2
78                            + 0.50 ALDX                : ARR2(3.01D-12, -190.0_dp, TEMP);
79 < R64>  ROOH + hv          = OH        + HO2
80                            + 0.50 ALD2 + 0.50 ALDX    : 0.7   *j(Pj_h2o2); 
81 < R65>  OH   + CO          = HO2                      : EP3(1.44D-13, 0.0_dp, 3.43D-33, 0.0_dp, TEMP, C_M);
82 < R66>  OH   + CH4         = MEO2                     : ARR2(2.45D-12, 1775.0_dp, TEMP);
83 < R67>  MEO2 + NO          = FORM      + HO2    + NO2 : ARR2(2.8D-12, -300.0_dp, TEMP);
84 < R68>  MEO2 + HO2         = MEPX                     : ARR2(4.1D-13, -750.0_dp, TEMP);
85 < R69>  MEO2 + MEO2        = 1.37 FORM + 0.74 HO2
86                            + 0.63 MEOH                : ARR2(9.5D-14, -390.0_dp, TEMP);
87 < R70>  MEPX + OH          = 0.70 MEO2 + 0.30 XO2
88                            + 0.30 HO2                 : ARR2(3.8D-12, -200.0_dp, TEMP);
89 < R71>  MEPX + hv          = FORM      + HO2    + OH  : 0.7   *j(Pj_h2o2);
90 < R72>  MEOH + OH          = FORM      + HO2          : ARR2(7.3D-12, 620.0_dp, TEMP);
91 < R73>  FORM + OH          = HO2       + CO           : 9.0D-12;
92 < R74>  FORM + hv          = 2.000 HO2 + CO           : j(Pj_ch2or);
93 < R75>  FORM + hv          = CO                       : j(Pj_ch2om);
94 < R76>  FORM + O           = OH        + HO2    + CO  : ARR2(3.4D-11, 1600.0_dp, TEMP);
95 < R77>  FORM + NO3         = HNO3      + HO2    + CO  : 5.8D-16;
96 < R78>  FORM + HO2         = HCO3                     : ARR2(9.7D-15, -625.0_dp, TEMP);
97 < R79>  HCO3               = FORM      + HO2          : ARR2(2.4D+12, 7000.0_dp, TEMP);
98 < R80>  HCO3 + NO          = FACD      + NO2 + HO2    : 5.6D-12;
99 < R81>  HCO3 + HO2         = MEPX                     : ARR2(5.6D-15, -2300.0_dp, TEMP);
100 < R82>  FACD + OH          = HO2                      : 4.0D-13;
101 < R83>  ALD2 + O           = C2O3      + OH           : ARR2(1.8D-11, 1100.0_dp, TEMP);
102 < R84>  ALD2 + OH          = C2O3                     : ARR2(5.6D-12, -270.0_dp, TEMP);
103 < R85>  ALD2 + NO3         = C2O3      + HNO3         : ARR2(1.4D-12, 1900.0_dp, TEMP);
104 < R86>  ALD2 + hv          = MEO2      + CO  + HO2    : 4.6D-4*j(Pj_no2);
105 < R87>  C2O3 + NO          = MEO2      + NO2          : ARR2(8.1D-12, -270.0_dp, TEMP);
106 < R88>  C2O3 + NO2         = PAN          : FALL(2.7D-28, 0.0_dp, -7.1_dp, 1.2D-11, 0.0_dp, -0.9_dp, 0.3_dp, TEMP, C_M);
107 < R89>  PAN                = C2O3      + NO2          : FALL(4.9D-03, 12100.0_dp, 0.0_dp, 5.4D16, 13830.0_dp, 0.0_dp, 0.3_dp, TEMP, C_M);
108 < R90>  PAN  + hv          = C2O3      + NO2          : j(Pj_pan);
109 < R91>  C2O3 + HO2         = 0.80 PACD + 0.20 AACD
110                            + 0.20 O3                  : ARR2(4.3D-13, -1040.0_dp, TEMP);
111 < R92>  C2O3 + MEO2        = 0.90 MEO2 + 0.90 HO2
112                            + FORM      + 0.10 AACD    : ARR2(2.0D-12, -500.0_dp, TEMP);
113 < R93>  C2O3 + XO2         = 0.90 MEO2 + 0.10 AACD    : ARR2(4.4D-13, -1070.0_dp, TEMP);
114 < R94>  C2O3 + C2O3        = 2.00 MEO2                : ARR2(2.9D-12, -500.0_dp, TEMP);
115 < R95>  PACD + OH          = C2O3                     : ARR2(4.0D-13, -200.0_dp, TEMP);
116 < R96>  PACD + hv          = MEO2      + OH           : 0.0*0.7   *j(Pj_h2o2);
117 < R97>  AACD + OH          = MEO2                     : ARR2(4.0D-13, -200.0_dp, TEMP);
118 < R98>  ALDX + O           = CXO3      + OH           : ARR2(1.3D-11, 870.0_dp, TEMP);
119 < R99>  ALDX + OH          = CXO3                     : ARR2(5.1D-12, -405.0_dp, TEMP);
120 <R100>  ALDX + NO3         = CXO3      + HNO3         : 6.5D-15;
121 <R101>  ALDX + hv          = MEO2      + CO     + HO2 : j(Pj_ch3cho);
122 <R102>  CXO3 + NO          = ALD2      + NO2
123                            + HO2       + XO2          : ARR2(6.7D-12, -340.0_dp, TEMP);
124 <R103>  CXO3 + NO2         = PANX                     : FALL(2.7D-28, 0.0_dp, -7.1_dp, 1.2D-11, 0.0_dp, -0.9_dp, 0.3_dp, TEMP, C_M);
125 <R104>  PANX               = CXO3      + NO2          : FALL(4.9D-03, 12100.0_dp, 0.0_dp, 5.4D16, 13830.0_dp, 0.0_dp, 0.3_dp, TEMP, C_M);
126 <R105>  PANX + hv          = CXO3      + NO2          : j(Pj_pan);
127 <R106>  PANX + OH          = ALD2      + NO2          : 3.0D-13;
128 <R107>  CXO3 + HO2         = 0.80 PACD + 0.20 AACD
129                            + 0.20 O3                  : ARR2(4.3D-13, -1040.0_dp, TEMP);
130 <R108>  CXO3 + MEO2        = 0.90 ALD2 + 0.90 XO2
131                            + HO2       + 0.10 AACD
132                            + 0.10 FORM                : ARR2(2.0D-12, -500.0_dp, TEMP);
133 <R109>  CXO3 + XO2         = 0.90 ALD2 + 0.10 AACD    : ARR2(4.4D-13, -1070.0_dp, TEMP);
134 <R110>  CXO3 + CXO3        = 2.00 ALD2 + 2.00 XO2
135                            + 2.00 HO2                 : ARR2(2.9D-12, -500.0_dp, TEMP);
136 <R111>  CXO3 + C2O3        = MEO2      + XO2
137                            + HO2       + ALD2         : ARR2(2.9D-12, -500.0_dp, TEMP);
138 <R112>  PAR  + OH          = 0.87 XO2  + 0.13 XO2N
139                            + 0.11 HO2  + 0.06 ALD2
140                            - 0.11 PAR  + 0.76 ROR
141                            + 0.05 ALDX                : 8.1D-13;
142 <R113>  ROR                = 0.96 XO2  + 0.60 ALD2
143                            + 0.94 HO2  - 2.10 PAR
144                            + 0.04 XO2N + 0.02 ROR
145                            + 0.50 ALDX                : ARR2(1.0D15, 8000.0_dp, TEMP);
146 <R114>  ROR                = HO2                      : 1.6E+03;
147 <R115>  ROR  + NO2         = NTR                      : 1.5D-11;
148 <R116>  O    + OLE         = 0.20 ALD2 + 0.30 ALDX
149                            + 0.30 HO2  + 0.20 XO2
150                            + 0.20 CO   + 0.20 FORM
151                            + 0.01 XO2N + 0.20 PAR
152                            + 0.10 OH                  : ARR2(1.0D-11, 280.0_dp, TEMP);
153 <R117>  OH   + OLE         = 0.80 FORM + 0.33 ALD2
154                            + 0.62 ALDX + 0.80 XO2
155                            + 0.95 HO2  - 0.70 PAR     : 3.2D-11;
156 <R118>  O3   + OLE         = 0.18 ALD2 + 0.74 FORM
157                            + 0.32 ALDX + 0.22 XO2
158                            + 0.10 OH   + 0.33 CO
159                            + 0.44 HO2  - 1.00 PAR     : ARR2(6.5D-15, 1900.0_dp, TEMP);
160 <R119>  NO3  + OLE         = NO2       + FORM
161                            + 0.91 XO2  + 0.09 XO2N
162                            + 0.56 ALDX + 0.35 ALD2
163                            - 1.00 PAR                 : ARR2(7.0D-13, 2160.0_dp, TEMP);
164 <R120>  O    + ETH         = FORM      + 1.70 HO2
165                            + CO        + 0.70 XO2
166                            + 0.30 OH                  : ARR2(1.04D-11, 792.0_dp, TEMP);
167 <R121>  OH   + ETH         = XO2       + 1.56 FORM
168                            + 0.22 ALDX + HO2          : TROE(1.0D-28, 0.8_dp, 8.8D-12, 0.0_dp, TEMP, C_M);
169 <R122>  O3   + ETH         = FORM      + 0.63 CO
170                            + 0.13 HO2  + 0.13 OH
171                            + 0.37 FACD                : ARR2(1.2D-14, 2630.0_dp, TEMP);
172 <R123>  NO3  + ETH         = NO2       + XO2
173                            + 2.0 FORM                 : ARR2(3.3D-12, 2880.0_dp, TEMP);
174 <R124>  IOLE + O           = 1.24 ALD2 + 0.66 ALDX
175                            + 0.10 HO2  + 0.10 XO2
176                            + 0.10 CO   + 0.10 PAR     : 2.3D-11;
177 <R125>  IOLE + OH          = 1.30 ALD2 + 0.70 ALDX
178                            + HO2       + XO2          : ARR2(1.0D-11, -550.0_dp, TEMP);
179 <R126>  IOLE + O3          = 0.65 ALD2 + 0.35 ALDX
180                            + 0.25 FORM + 0.25 CO
181                            + 0.50 O    + 0.50 OH
182                            + 0.50 HO2                 : ARR2(8.4D-15, 1100.0_dp, TEMP);
183 <R127>  IOLE + NO3         = 1.18 ALD2 + 0.64 ALDX
184                            + HO2       + NO2          : ARR2(9.6D-13, 270.0_dp, TEMP);
185 <R128>  TOL  + OH          = 0.44 HO2  + 0.08 XO2
186                            + 0.36 CRES + 0.56 TO2
187                            + 0.071 TOLAER1
188                            + 0.138 TOLAER2            : ARR2(1.8D-12, -355.0_dp, TEMP);
189 <R129>  TO2  + NO          = 0.90 NO2  + 0.90 HO2
190                            + 0.90 OPEN + 0.10 NTR     : 8.1D-12;
191 <R130>  TO2                = CRES      + HO2          : 4.2;
192 <R131>  OH   + CRES        = 0.40 CRO  + 0.60 XO2
193                            + 0.60 HO2  + 0.30 OPEN
194                            + CSLAER                   : 4.1D-11;
195 <R132>  CRES + NO3         = CRO + HNO3 + CSLAER      : 2.2D-11;
196 <R133>  CRO  + NO2         = NTR                      : 1.4D-11;
197 <R134>  CRO  + HO2         = CRES                     : 5.5D-12;
198 <R135>  OPEN + hv          = C2O3       + HO2   + CO  : 9.0    *j(Pj_ch2or);
199 <R136>  OPEN + OH          = XO2        + 2.0 CO
200                            + 2.00 HO2   + C2O3 + FORM : 3.0D-11;
201 <R137>  OPEN + O3          = 0.03 ALDX  + 0.62 C2O3
202                            + 0.70 FORM  + 0.03 XO2
203                            + 0.69 CO    + 0.08 OH
204                            + 0.76 HO2   + 0.20 MGLY   : ARR2(5.4D-17, 500.0_dp, TEMP);
205 <R138>  OH   + XYL         = 0.70 HO2   + 0.50 XO2
206                            + 0.20 CRES  + 0.80 MGLY
207                            + 1.10 PAR   + 0.30 TO2
208                            + 0.038 XYLAER1
209                            + 0.167 XYLAER2            : ARR2(1.7D-11, -116.0_dp, TEMP);
210 <R139>  OH   + MGLY        = XO2        + C2O3        : 1.7D-11;
212 <R140>  MGLY + hv          = C2O3       + HO2   + CO  : 9.64  *j(Pj_ch2or);
213 <R141>  O    + ISOP        = 0.75 ISPD  + 0.50 FORM
214                            + 0.25 XO2   + 0.25 HO2
215                            + 0.25 CXO3
216                            + 0.232 ISOAER1
217                            + 0.0228 ISOAER2           : 3.6D-11;
218 <R142>  OH   + ISOP        = 0.912 ISPD + 0.629 FORM
219                            + 0.991 XO2  + 0.912 HO2
220                            + 0.088 XO2N
221                            + 0.232 ISOAER1
222                            + 0.0288 ISOAER2           : ARR2(2.54D-11, -407.6_dp, TEMP);
223 <R143>  O3   + ISOP        = 0.650 ISPD + 0.600 FORM
224                            + 0.200 XO2  + 0.066 HO2
225                            + 0.266 OH   + 0.200 CXO3
226                            + 0.150 ALDX + 0.350 PAR
227                            + 0.066 CO
228                            + 0.232 ISOAER1
229                            + 0.0288 ISOAER2           : ARR2(7.86D-15, 1912.0_dp, TEMP);
230 <R144>  NO3  + ISOP        = 0.200 ISPD + 0.800 NTR
231                            + XO2        + 0.800 HO2
232                            + 0.200 NO2  + 0.800 ALDX
233                            + 0.232 ISOAER1
234                            + 0.0288 ISOAER2           : ARR2(3.03D-12, 448.0_dp, TEMP);
235 <R145>  OH   + ISPD        = 1.565 PAR  + 0.167 FORM
236                            + 0.713 XO2  + 0.503 HO2
237                            + 0.334 CO   + 0.168 MGLY
238                            + 0.252 ALD2 + 0.210 C2O3
239                            + 0.250 CXO3 + 0.120 ALDX  : 3.36D-11;
240 <R146>  O3   + ISPD        = 0.114 C2O3 + 0.150 FORM
241                            + 0.850 MGLY + 0.154 HO2
242                            + 0.268 OH   + 0.064 XO2
243                            + 0.020 ALD2 + 0.360 PAR
244                            + 0.225 CO                 : 7.1D-18;
245 <R147>  NO3  + ISPD        = 0.357 ALDX + 0.282 FORM
246                            + 1.282 PAR  + 0.925 HO2
247                            + 0.643 CO   + 0.850 NTR
248                            + 0.075 CXO3 + 0.075 XO2
249                            + 0.150 HNO3               : 1.0D-15;
250 <R148>  ISPD + hv          = 0.333 CO   + 0.067 ALD2
251                            + 0.900 FORM + 0.832 PAR
252                            + 1.033 HO2  + 0.700 XO2
253                            + 0.967 C2O3               : 0.0036*0.025 *j(Pj_ch2om);
254 <R149>  TERP + O           = 0.150 OH   + 5.12  PAR
255                            + TERPAER                  : 3.6D-11;
256 <R150>  TERP + OH          = 0.750 HO2  + 1.250 XO2
257                            + 0.250 XO2N + 0.280 FORM
258                            + 1.66  PAR  + 0.470 ALDX
259                            + TERPAER                  : ARR2(1.5D-11, -449.0_dp, TEMP);
260 <R151>  TERP + O3          = 0.570 OH   + 0.070 HO2
261                            + 0.760 XO2  + 0.180 XO2N
262                            + 0.240 FORM + 0.001 CO
263                            + 7.000 PAR  + 0.210 ALDX
264                            + 0.390 CXO3
265                            + TERPAER                  : ARR2(1.2D-15, 821.0_dp, TEMP);
266 <R152>  TERP + NO3         = 0.470 NO2  + 0.280 HO2
267                            + 1.030 XO2  + 0.250 XO2N
268                            + 0.470 ALDX + 0.530 NTR
269                            + TERPAER                  : ARR2(3.7D-12, -175.0_dp, TEMP);
270 <R153>  SO2  + OH          = SULF + HO2 + SULAER      : TROE(3.3D-31, 4.3_dp, 1.6D-12, 0.0_dp, TEMP, C_M);
271 <R154>  OH   + ETOH        = HO2        + 0.900 ALD2
272                            + 0.050 ALDX + 0.100 FORM
273                            + 0.100 XO2                : ARR2(6.9D-12, 230.0_dp, TEMP);
274 <R155>  OH   + ETHA        = 0.991 ALD2 + 0.991 XO2
275                            + 0.009 XO2N + HO2         : ARR2(8.7D-12, 1070.0_dp, TEMP);
276 <R156>  NO2  + ISOP        = 0.200 ISPD + 0.800 NTR
277                            + XO2        + 0.800 HO2
278                            + 0.200 NO   + 0.800 ALDX
279                            + 2.400 PAR                : 1.5D-19;
280 <R157>  CL2  + hv          = 2.000 CL                 : j(Pj_cl2);
281 <R158>  HOCL + hv          = OH        + CL           : j(Pj_hocl);
282 <R159>  CL   + O3          = CLO                      : ARR2(2.3D-11, 200.0_dp, TEMP);
283 <R160>  CLO  + CLO         = 0.300 CL2 + 1.400 CL     : 1.63D-14;
284 <R161>  CLO  + NO          = CL        + NO2          : ARR2(6.4D-12, -290.0_dp, TEMP);
285 <R162>  CLO  + HO2         = HOCL                     : ARR2(2.7D-12, -220.0_dp, TEMP);
286 <R163>  OH   + FMCL        = CL        + CO           : 5.0D-13;
287 <R164>  FMCL + hv          = CL        + CO + HO2     : j(Pj_fmcl);
288 <R165>  CL   + CH4         = HCL       + MEO2         : ARR2(6.6D-12, 1240.0_dp, TEMP);
289 <R166>  CL   + PAR         = HCL
290                              + 0.870 XO2
291                              + 0.130 XO2N
292                              + 0.110 HO2
293                              + 0.060 ALD2
294                              - 0.110 PAR
295                              + 0.760 ROR
296                              + 0.050 ALDX             : 5.0D-11;
297 <R167>  CL   + ETHA       = HCL
298                              + 0.991 ALD2
299                              + 0.991 XO2
300                              + 0.009 XO2N
301                              + HO2                    : ARR2(8.3D-11, 100.0_dp, TEMP);
302 <R168>  CL   + ETH         = FMCL
303                              + 2.000 XO2
304                              + 1.000 HO2
305                              + 1.000 FORM             : 1.07D-10;
306 <R169>  CL   + OLE         = FMCL
307                              + 0.330 ALD2
308                              + 0.670 ALDX
309                              + 2.000 XO2
310                              + 1.000 HO2
311                              - 1.000 PAR              : 2.5D-10;
312 <R170>  CL   + IOLE        = 0.300 HCL
313                              + 0.700 FMCL
314                              + 0.450 ALD2
315                              + 0.550 ALDX
316                              + 0.300 OLE
317                              + 0.300 PAR
318                              + 1.700 XO2
319                              + 1.000 HO2              : 3.5D-10;
320 <R171>  CL   + ISOP        = 0.15 HCL
321                              + 1.000 XO2
322                              + 1.000 HO2
323                              + 0.850 FMCL
324                              + 1.000 ISPD             : 4.3D-10;
325 <R172>  CL   + FORM       = HCL
326                              + 1.000 HO2
327                              + 1.000 CO               : ARR2(8.2D-11, 34.0_dp, TEMP);
328 <R173>  CL   + ALD2       = HCL
329                              + 1.000 C2O3             : 7.9D-11;
330 <R174>  CL   + ALDX       = HCL
331                              + 1.000 CXO3             : 1.3D-10;
332 <R175>  CL   + MEOH       = HCL
333                              + 1.000 HO2
334                              + 1.000 FORM             : 5.5D-11;
335 <R176>  CL   + ETOH       = HCL
336                              + 1.000 HO2
337                              + 1.000 ALD2             : ARR2(8.2D-11, -45.0_dp, TEMP);
338 <R177>  HCL  + OH         = CL                        : ARR(6.58D-13, -58.0_dp, 1.16_dp, TEMP);
339 <R178>  HG0  + O3          = HG2                      : 3.0D-20;
340 <R179>  HG0  + OH          = HG2                      : 8.7D-14;
341 <R180>  HG0  + H2O2        = HG2                      : 8.5D-19;
342 <R181>  HUM  + OH          = HUMAER     + OH          : 2.93D-10;
343 <R182>  LIM  + OH          = 0.239 LIMAER1
344                            + 0.363 LIMAER2
345                            + OH                       : 1.71D-10;
346 <R183>  OCI  + OH          = 0.045 OCIAER1
347                            + 0.149 OCIAER2
348                            + OH                       : 2.52D-10;
349 <R184>  APIN + OH          = 0.038 APINAER1
350                            + 0.326 APINAER2
351                            + OH                       : 5.37D-11;
352 <R185>  APIN + O3          = 0.125 APINAER3
353                            + 0.102 APINAER4
354                            + O3                       : 8.66D-17;
355 <R186>  BPIN + OH          = 0.13  BPINAER1
356                            + 0.0406 BPINAER2
357                            + OH                       : 7.89D-11;
358 <R187>  BPIN + O3          = 0.026 BPINAER3
359                            + 0.485 BPINAER4
360                            + O3                       : 1.36D-17;
361 <R188>  BPIN + NO3         = BPINAER5   + NO3         : 2.31D-12;
362 <R189>  TER  + OH          = 0.091 TERAER1
363                            + 0.367 TERAER2
364                            + OH                       : 2.7D-10;
365 <R190>  ALKH + OH          = 1.173 ALKHAER1
366                            + OH                       : 1.97D-11;
367 <R191>  PAH  + OH          = 0.156 PAHAER1
368                            + 0.777 PAHAER2
369                            + OH                       : 7.7D-11;