Merge branch 'master' into devel
[wrffire.git] / wrfv2_fire / chem / KPP / kpp / kpp-2.1 / models / wet.eqn
blob75d22bcd47f443ade38bdf783de946fef9e1ff3d
1 #EQUATIONS
2 { 1.} NO2 + hv =  NO + O              :  8.89E-3*SUN ;
3 { 2.} O {+ O2 + M} = O3               :  ARR2(1.4E+3, 1175.0) ;
4 { 3.} O3 + NO = NO2                   :  ARR2(1.8E-12, -1370.0) ;
5 { 4.} O + NO2 = NO                    :  9.3E-12 ;
6 { 5.} O + NO2   = NO3                 :  ARR2(1.6E-13, 687.0) ;
7 { 6.} O + NO   =  NO2                 :  ARR2(2.2E-13, 602.0) ;
8 { 7.} O3 + NO2 =  NO3                 :  ARR2(1.2E-13, -2450.0) ;
9 { 8.} O3 + hv = O                     :  4.0E-2*RCONST(1) ;
10 { 9.} O3 + hv = O1D                   :  2.8E-3*RCONST(1) ;
11 {10.} O1D   = O                       :  ARR2(1.9E+8, 390.0)  ;
13 {11.} O1D + H2O = 2OH                 :  2.2E-10 ;
14 {12.} O3 + OH = HO2                   :  ARR2(1.6E-12, -940.0) ;
15 {13.} O3 + HO2 = OH                   :  ARR2(1.4E-14, -580.0) ;
16 {14.} NO3 + hv = 0.89 NO2 + 0.89 O 
17                  + 0.11 NO            :  15.5*RCONST(1)  ;
18 {15.} NO3 + NO = 2 NO2                :  ARR2(1.3E-11, 250.0) ;
19 {16.} NO3 + NO2 = NO + NO2            :  ARR2(2.5E-14, -1230.0) ;
20 {17.} NO3 + NO2   =  N2O5             :  ARR2(5.3E-13, 256.0) ;
21 {18.} N2O5 + H2O = 2 HNO3             :  1.3E-21 ;
22 {19.} N2O5   =  NO3 + NO2             :  ARR2(3.5E+14, -10897.0) ;
23 {20.} 2 NO  =  2 NO2                  :  ARR2(1.8E-20, 530.0) ;
26 {21.} NO + NO2 + H2O = 2 HONO         :  4.4E-40 ;
27 {22.} OH + NO   =  HONO               :  ARR2(4.5E-13, 806.0) ;
28 {23.} HONO + hv =  OH + NO            :  0.17*RCONST(1) ;
29 {24.} OH + HONO =  NO2                :  6.6E-12 ;
30 {25.} 2 HONO = NO + NO2               :  1.0E-20 ;
31 {26.} OH + NO2   =  HNO3              :  ARR2(1.0E-12, 713.0) ;
32 {27.} OH + HNO3   =  NO3              :  ARR2(5.1E-15, 1000.0) ;
33 {28.} HO2 + NO = OH + NO2             :  ARR2(3.7E-12, 240.0) ;
34 {29.} HO2 + NO2   =  PNA              :  ARR2(1.2E-13, 749.0) ;
35 {30.} PNA   = HO2 + NO2               :  ARR2(4.8E+13, -10121.0) ;
37 {31.} OH + PNA = NO2                  :  ARR2(1.3E-12, 380.0) ;
38 {32.} 2 HO2 = H2O2                    :  ARR2(5.9E-14, 1150.0)  ;
39 {33.} 2 HO2 + H2O = H2O2              :  ARR2(2.2E-38, 5800.0) ;
40 {34.} H2O2 + hv = 2 OH                :  7.1E-4*RCONST(1) ;
41 {35.} OH + H2O2 = HO2                 :  ARR2(3.1E-12, -187.0) ;
42 {36.} OH + CO  = HO2                  :  2.2E-13 ;
43 {37.} HCHO + OH  =  HO2 + CO    :  1.0E-11 ;
44 {38.} HCHO + hv {+ 2 O2} = 2 HO2 + CO :  3.2E-3*RCONST(1) ;
45 {39.} HCHO +  hv = CO                 :  4.2E-3*RCONST(1) ;
46 {40.} HCHO + O = OH + HO2 + CO        :  ARR2(3.0E-11, -1550.0) ;
48 {41.} HCHO + NO3  = HNO3 
49                         + HO2 + CO    :  6.3E-16 ;
50 {42.} ALD2 + O  =  C2O3 + OH          :  ARR2(1.2E-11, -986.0) ;
51 {43.} ALD2 + OH = C2O3                :  ARR2(7.0E-12, 250.0) ;
52 {44.} ALD2 + NO3  = C2O3 + HNO3 :  2.5E-15  ;
53 {45.} ALD2 + hv {+ 2 O2} = HCHO + XO2 
54                          + CO + 2 HO2 :  4.5E-4*RCONST(1) ;
55 {46.} C2O3 + NO  = HCHO + XO2 
56                        + HO2 + NO2    :  ARR2(5.4E-12, 250.0) ;
57 {47.} C2O3 + NO2 = PAN                :  ARR2(8.0E-20, 5500.0) ;
58 {48.} PAN = C2O3 + NO2                :  ARR2(9.4E+16, -14000.0) ;
59 {49.} 2 C2O3 = 2 HCHO + 2 XO2 + 2 HO2 : 2.0E-12  ;
60 {50.} C2O3 + HO2 = 0.79 HCHO 
61    + 0.79 XO2 + 0.79 HO2 + 0.79 OH    :  6.5E-12 ;
63 {51.} OH = HCHO + XO2 + HO2           :  ARR2(1.1E+2, -1710.0) ;
64 {52.} PAR + OH =  0.87 XO2 + 0.13 XO2N 
65                   + 0.11 HO2 + 0.11 ALD2
66                   + 0.76 ROR - 0.11 PAR        :  8.1E-13 ;
67 {53.} ROR = 1.1 ALD2 + 0.96 XO2 + 0.94 HO2 + 0.04 XO2N
68             + 0.02 ROR - 2.10 PAR     :  ARR2(1.0E+15, -8000.0) ;
69 {54.} ROR = HO2                       :  1.6E+3 ;
70 {55.} ROR + NO2 =  PROD               : 1.5E-11  ;
71 {56.} O + OLE = 0.63 ALD2 + 0.38 HO2 
72                 + 0.28 XO2 + 0.3 CO
73                 + 0.2 HCHO + 0.02 XO2N 
74                 + 0.22 PAR + 0.2 OH   :  ARR2(1.2E-11, -324.0) ;
75 {57.} OH + OLE = HCHO + ALD2 + XO2 
76                  + HO2 - PAR          :  ARR2(5.2E-12, 504.0) ;
77 {58.} O3 + OLE = 0.5 ALD2 + 0.74 HCHO 
78                  + 0.33 CO + 0.44 HO2 
79                  + 0.22 XO2
80                  + 0.1 OH - PAR       :   ARR2(1.4E-14, -2105.0)  ;
81 {59.} NO3 + OLE = 0.91 XO2 + HCHO 
82                   + ALD2 + 0.09 XO2N                          
83                   + NO2 - PAR         :  7.7E-15 ;
84 {60.} O + ETH = HCHO + 0.7 XO2 + CO + 1.7 HO2 + 0.3 OH    :  ARR2(1.0E-11, -792.0) ;
86 {61.} OH + ETH = XO2 + 1.56 HCHO + HO2 + 0.22 ALD2        :  ARR2(2.0E-12, 411.0)  ;
87 {62.} O3 + ETH = HCHO + 0.42 CO + 0.12 HO2                :  ARR2(1.3E-14, -2633.0) ;
88 {63.} OH + TOL = 0.08 XO2 + 0.36 CRES 
89                  + 0.44 HO2 + 0.56 TO2                    :  ARR2(2.1E-12, 322.0) ;
90 {64.} TO2 + NO =  0.9 NO2 + 0.9 OPEN + 0.9 HO2            :  8.1E-12 ;
91 {65.} TO2 = HO2 + CRES                                    :  4.2 ;
92 {66.} OH + CRES = 0.4 CRO + 0.6 XO2 + 0.6 HO2 + 0.3 OPEN  :  4.1E-11  ;
93 {67.} NO3 + CRES = CRO + HNO3                             :  2.2E-11 ;
94 {68.} CRO + NO2 = PROD                                    :  1.4E-11 ;
95 {69.} OH + XYL = 0.7 HO2 + 0.5 XO2 + 0.2 CRES + 0.8 MGLY
96                  + 1.10 PAR + 0.3 TO2                     :  ARR2(1.7E-11, 116.0) ;
97 {70.} OH + OPEN = XO2 + C2O3 + 2 HO2 + 2 CO + HCHO        :  3.0E-11 ;
99 {71.} OPEN + hv = C2O3 + CO + HO2                         :  6.0E-3*RCONST(1) ;
100 {72.} O3 + OPEN = 0.03 ALD2 + 0.62 C2O3 
101                   + 0.7 HCHO + 0.03 XO2 + 0.69 CO 
102                   + 0.08 OH + 0.76 HO2 + 0.2 MGLY         : ARR2(5.4E-17, -500.0)  ;
103 {73.} OH + MGLY =  XO2 + C2O3                             : 1.7E-11 ;
104 {74.} MGLY + hv = C2O3 + CO + HO2                         : 1.86E-2*RCONST(1)  ;
105 {75.} O + ISOP =  0.6 HO2 + 0.8 ALD2 + 0.55 OLE + 0.5 XO2
106                   + 0.5 CO + 0.45 ETH + 0.9 PAR           :   1.8E-11  ;
107 {76.} OH + ISOP = HCHO + XO2 + 0.67 HO2 
108                   + 0.4 MGLY + 0.2 C2O3
109                   + ETH + 0.2 ALD2 + 0.13 XO2N            :   9.6E-11  ;
110 {77.} O3 + ISOP = HCHO + 0.4 ALD2 + 0.55 ETH + 0.2 MGLY 
111                   + 0.06 CO + 0.1 PAR + 0.44 HO2 + 0.1 OH :  1.2E-17 ;
112 {78.} NO3 + ISOP =  XO2N                                  :  3.2E-13 ;
113 {79.} XO2 + NO = NO2                                      :  8.1E-12 ;
114 {80.} 2 XO2 =  PROD                                       :  ARR2(1.7E-14, 1300.0) ;
115 {81.} XO2N + NO =  PROD                                   :  6.8E-13 ;
117 {62/RG71}      SO2  + OH  = HO2  + PROD  :           1.0870e-12 ;
118 {63/RG72}      SO2  = PROD  :                        1.3833e-06 ;
120 {Methane and Methylperoxyl Radical Reactions}
121 {}  CH4 + OH = CH3O2      :                        ARR2(2.4E-12,-1710.0) ;
122 {}  CH3O2 + NO = CH3O + NO2 :                      ARR2(4.15E-12,180.0)  ;
123 {}  CH3O2 + HO2 = PROD :                           ARR2(1.7E-13,1000.0) ;
124 {}  2 CH3O2 = 0.76 CH3O + 0.62 HCHO + PROD :       ARR2(1.62E-13, 220.0) ;
125 {}  CH3O + O2 = HCHO + HO2 :                       ARR2(7.12E-14,-1080.0);
127 {82/RE21}      HO2  = HO2aq :                        3.3000e-01 ;
128 {83/RE22}      HO2aq = HO2 :                         2.2492e+01 ;   
129 {84/RE31}      H2O2  = H2O2aq  :                     3.3000e-01 ;     
130 {85/RE32}      H2O2aq = H2O2  :                      6.0790e-01 ;     
131 {86/RE41}      O3  = O3aq :                          1.0500e-02 ;        
132 {87/RE42}      O3aq = O3  :                          1.3012e+05 ;       
133 {88/RE51}      HCHO  = HCHOaq :                      2.4000e-01 ; 
134 {89/RE52}      HCHOaq = HCHO :                       5.1931e+00 ;  
135 {90/RE61}      HCOOH  = HCOOHaq  :                   2.4000e-01 ;       
136 {91/RE62}      HCOOHaq = HCOOH :                     8.8422e+00 ;        
137 {92/RE71}      ROOH  = CH3O2Haq :                    2.4000e-01 ;    
138 {93/RE72}      CH3O2Haq = ROOH  :                    1.4871e+02 ;   
139 {94/RE81}      XO2  = CH3O2aq :                      2.4000e-01 ;   
140 {95/RE82}      CH3O2aq = XO2  :                      1.6358e+01 ;  
141 {96/RE91}      HNO3  = HNO3aq  :                     2.4000e-01 ;  
142 {97/RE92}      HNO3aq = HNO3  :                      1.5579e-01 ;  
143 {98/RE101}     NO3  = NO3aq  :                       2.4000e-01 ;  
144 {99/RE102}     NO3aq = NO3 :                         2.1811e+03 ;   
145 {100/RE111}    NH3  = NH3aq :                        2.4000e-01 ; 
146 {101/RE112}    NH3aq = NH3  :                        4.3622e+02 ;
147 {102/RE121}    SO2  = SO2aq :                        2.4000e-01 ;         
148 {103/RE122}    SO2aq = SO2 :                         2.6599e+04 ;         
149 {104/RE131}    OH  = OHaq :                          2.4000e-01 ;         
150 {105/RE132}    OHaq = OH  :                          3.6351e+00 ;        
151 {106/RE141}    N2O5  = HNO3aq + HNO3aq :             7.5000e-02 ;
152 {107/RE151}    H2Oaq = HPLUSaq + OHMINaq   :         2.3636e-05 ; 
153 {108/RE152}    HPLUSaq + OHMINaq = H2Oaq  :          7.1958e-04 ; 
154 {109/RE171}    HCOOHaq = HPLUSaq + HCOOMINaq  :      8.6000e+06 ;   
155 {110/RE172}    HPLUSaq + HCOOMINaq = HCOOHaq  :      2.6016e-04 ;  
156 {111/RE181}    HO2aq = HPLUSaq + O2MINaq :           1.6000e+05 ;      
157 {112/RE182}    HPLUSaq + O2MINaq = HO2aq  :          5.5353e-05 ;    
158 {113/RE191}    SO2aq = HSO3MINaq + HPLUSaq :         3.6000e+06 ;  
159 {114/RE192}    HSO3MINaq + HPLUSaq = SO2aq  :        1.6606e-06 ;
160 {115/RE201}    HSO3MINaq = SO32MINaq + HPLUSaq  :    1.0000e+04 ;   
161 {116/RE202}    SO32MINaq + HPLUSaq = HSO3MINaq  :    8.3029e-04 ;  
162 {117/RE211}    NH3aq = NH4aq + OHMINaq  :            6.0000e+05 ;     
163 {118/RE212}    NH4aq + OHMINaq = NH3aq :             1.8820e-04 ;     
164 {119/RE221}    HNO3aq = NO3MINaq + HPLUSaq :         2.2000e+09 ;    
165 {120/RE222}    NO3MINaq + HPLUSaq = HNO3aq  :        7.7494e-07 ;  
167 {121/RW01}      FE3PLUSaq = FE2PLUSaq + OHaq :       0.4583E-03  ; 
168 {122/RW02}      H2O2aq = OHaq + OHaq :               0.7219E-05 ;         
169 {123/RW03 }     HCHOaq + OHaq = HCOOHaq 
170                       + HO2aq + H2Oaq :              0.1107E-04  ;   
171 {124/RW05}      HCOOMINaq + OHaq = O2MINaq + CO2aq : 0.1771E-04  ;   
172 {125/RW06}      HCOOHaq + OHaq = O2MINaq 
173                     + HPLUSaq + CO2aq :              0.7196E-06  ; 
174 {126/RW07}      HSO3MINaq + OHaq = SO42MINaq + HPLUSaq :  0.2491E-04  ;  
175 {127/RW08}      HMSAaq + OHaq = SO42MINaq 
176                       + HCHOaq + HPLUSaq :            0.5535E-05  ;
177 {128/RW09}      FE2PLUSaq + OHaq = FE3PLUSaq
178                       + OHMINaq:                      0.2380E-05  ;
179 {129/RW10}      H2O2aq + OHaq = H2Oaq + HO2aq  :      0.1495E-06  ;      
180 {130/RW11}      O3aq + OHaq = HO2aq  :                0.1661E-04  ;             
182 {131/RW13}      O2MINaq + FE2PLUSaq = FE3PLUSaq
183                      + H2O2aq + OHMINaq + OHMINaq :   0.5535E-07  ;
184 {132/RW14}      O2MINaq + FE3PLUSaq = FE2PLUSaq :     0.8303E-06  ;                      
185 {133/RW15}      O2MINaq + CUPLUSaq = CU2PLUSaq 
186                   + H2O2aq + OHMINaq + OHMINaq :      0.4982E-04 ;
187 {134/RW16}      O2MINaq + CU2PLUSaq = CUPLUSaq  :     0.4428E-04 ;                    
188 {135/RW17}      O2MINaq + O3aq = OHaq + OHMINaq  :    0.8303E-05 ;                    
189 {136/RW18}      HO2aq   + FE2PLUSaq = FE3PLUSaq 
190                     + H2O2aq + OHMINaq   :            0.6642E-08  ;   
191 {137/RW19}      HO2aq + FE3PLUSaq = FE2PLUSaq + HPLUSaq : 0.5535E-10  ;              
192 {138/RW20}      HO2aq + CUPLUSaq = CU2PLUSaq 
193                     + H2O2aq + OHMINaq  :             0.5535E-05 ;   
194 {139/RW21}      HO2aq + CU2PLUSaq = CUPLUSaq + HPLUSaq :  0.2768E-06  ;              
195 {140/RW22}      H2O2aq  + HSO3MINaq = SO42MINaq 
196                     + H2Oaq + HPLUSaq  :              0.5535E-14  ;    
197 {141/RW23}      H2O2aq  + FE2PLUSaq = FE3PLUSaq 
198                    + OHaq + OHMINaq :                 0.5535E-14  ; 
199 {142/RW24}      H2O2aq  + CUPLUSaq = CU2PLUSaq 
200                    + OHaq + OHMINaq  :                0.2214E-08  ;
201 {143/RW25}      CUPLUSaq = CU2PLUSaq + O2MINaq :      0.2500E-03  ;            
202 {144/RW26}      CUPLUSaq + FE3PLUSaq = CU2PLUSaq 
203                    + FE2PLUSaq :                      0.1661E-06 ;
204 {145/RW27}      CUPLUSaq + O3aq = CU2PLUSaq
205                    + OHaq + OHMINaq :                 0.5535E-07  ;
206 {146/RW28}      FE2PLUSaq + O3aq = FE3PLUSaq 
207                    + OHaq + OHMINaq :                 0.4539E-08  ;
208 {147/RW29}      FE3PLUSaq + HSO3MINaq = FE2PLUSaq 
209                       + SO42MINaq + HPLUSaq  :        0.5535E-14  ;
210 {148/RW30}      O3aq + HSO3MINaq = SO42MINaq + HPLUSaq  : 0.1716E-08  ;              
211 {149/RW31}      O3aq + SO32MINaq = SO42MINaq  :       0.5535E-05  ;                 
212 {150/RW32}      O3aq = H2O2aq  :                      0.2355E-04  ;          
213 {151/RW34}      HO2aq + O2MINaq = H2O2aq + OHMINaq :  0.5535E-06  ;            
214 {152/RW35}      HO2aq + HO2aq = H2O2aq  :             0.4594E-08  ;                    
215 {153/RW39}      SO32MINaq + HCHOaq = HMSAaq 
216                      + OHMINaq  :                     0.1384E-06  ;      
217 {154/RW40}      HMSAaq + OHMINaq = SO32MINaq 
218                      + HCHOaq :                       0.1993E-10  ;      
219 {155/RW42}      CH3O2aq + O2MINaq = CH3O2Haq 
220                      + OHMINaq  :                     0.2768E-06  ;   
221 {156/RW43}      CH3O2Haq + OHaq = CH3O2aq + H2Oaq  :  0.1495E-06  ;       
222 {157/RW44}      CH3O2Haq + OHaq = HCHOaq + OHaq  :    0.1052E-06  ; 
224 {CO2 reactions}
225 CO2aq = HPLUSaq + HCO3MINaq :  1.0E-4 ;
226 HPLUSaq + HCO3MINaq = CO2aq :  1.2178E-12 ;
228 {CO2aq + H2Oaq = H2CO3aq :                    0*3.0E-2 ;
229 H2CO3aq = CO2aq + H2Oaq :                     0*1.19E+1 ;
230 H2CO3aq = HPLUSaq + HCO3MINaq :             0*8.0E+6  ;
231 HPLUSaq + HCO3MINaq = H2CO3aq :             0*4.7E+10 ;
232 HPLUSaq + HCO3MINaq  = CO2aq :              0*1.5E+4 ; 
233 CO2aq  + H2Oaq = HPLUSaq + HCO3MINaq :      0*1.0E-2 ;}