Line data Source code
1 :
2 : module mo_setrxt
3 :
4 : use shr_kind_mod, only : r8 => shr_kind_r8
5 :
6 : private
7 : public :: setrxt
8 : public :: setrxt_hrates
9 :
10 : contains
11 :
12 72960 : subroutine setrxt( rate, temp, m, ncol )
13 :
14 : use ppgrid, only : pcols, pver
15 :
16 :
17 : use chem_mods, only : rxntot
18 : use mo_jpl, only : jpl
19 :
20 : implicit none
21 :
22 : !-------------------------------------------------------
23 : ! ... dummy arguments
24 : !-------------------------------------------------------
25 : integer, intent(in) :: ncol
26 : real(r8), intent(in) :: temp(pcols,pver)
27 : real(r8), intent(in) :: m(ncol*pver)
28 : real(r8), intent(inout) :: rate(ncol*pver,max(1,rxntot))
29 :
30 : !-------------------------------------------------------
31 : ! ... local variables
32 : !-------------------------------------------------------
33 : integer :: n
34 : integer :: offset
35 145920 : real(r8) :: itemp(ncol*pver)
36 145920 : real(r8) :: exp_fac(ncol*pver)
37 145920 : real(r8) :: ko(ncol*pver)
38 145920 : real(r8) :: kinf(ncol*pver)
39 :
40 58907904 : rate(:,168) = 1.2e-10_r8
41 58907904 : rate(:,172) = 1.2e-10_r8
42 58907904 : rate(:,178) = 6.9e-12_r8
43 58907904 : rate(:,179) = 7.2e-11_r8
44 58907904 : rate(:,180) = 1.6e-12_r8
45 58907904 : rate(:,186) = 1.8e-12_r8
46 58907904 : rate(:,190) = 1.8e-12_r8
47 58907904 : rate(:,202) = 3.5e-12_r8
48 58907904 : rate(:,204) = 1.3e-11_r8
49 58907904 : rate(:,205) = 2.2e-11_r8
50 58907904 : rate(:,206) = 5e-11_r8
51 58907904 : rate(:,241) = 1.7e-13_r8
52 58907904 : rate(:,243) = 2.607e-10_r8
53 58907904 : rate(:,244) = 9.75e-11_r8
54 58907904 : rate(:,245) = 2.07e-10_r8
55 58907904 : rate(:,246) = 2.088e-10_r8
56 58907904 : rate(:,247) = 1.17e-10_r8
57 58907904 : rate(:,248) = 4.644e-11_r8
58 58907904 : rate(:,249) = 1.204e-10_r8
59 58907904 : rate(:,250) = 9.9e-11_r8
60 58907904 : rate(:,251) = 3.3e-12_r8
61 58907904 : rate(:,270) = 4.5e-11_r8
62 58907904 : rate(:,271) = 4.62e-10_r8
63 58907904 : rate(:,272) = 1.2e-10_r8
64 58907904 : rate(:,273) = 9e-11_r8
65 58907904 : rate(:,274) = 3e-11_r8
66 58907904 : rate(:,279) = 2.14e-11_r8
67 58907904 : rate(:,280) = 1.9e-10_r8
68 58907904 : rate(:,293) = 2.57e-10_r8
69 58907904 : rate(:,294) = 1.8e-10_r8
70 58907904 : rate(:,295) = 1.794e-10_r8
71 58907904 : rate(:,296) = 1.3e-10_r8
72 58907904 : rate(:,297) = 7.65e-11_r8
73 58907904 : rate(:,310) = 4e-13_r8
74 58907904 : rate(:,315) = 1.31e-10_r8
75 58907904 : rate(:,316) = 3.5e-11_r8
76 58907904 : rate(:,317) = 9e-12_r8
77 58907904 : rate(:,324) = 6.8e-14_r8
78 58907904 : rate(:,325) = 2e-13_r8
79 58907904 : rate(:,340) = 1e-12_r8
80 58907904 : rate(:,344) = 1e-14_r8
81 58907904 : rate(:,345) = 1e-11_r8
82 58907904 : rate(:,346) = 1.15e-11_r8
83 58907904 : rate(:,347) = 3.3e-11_r8
84 58907904 : rate(:,348) = 3.4e-12_r8
85 58907904 : rate(:,349) = 4e-14_r8
86 58907904 : rate(:,362) = 3e-12_r8
87 58907904 : rate(:,363) = 1.2e-11_r8
88 58907904 : rate(:,364) = 6.7e-13_r8
89 58907904 : rate(:,374) = 3.5e-13_r8
90 58907904 : rate(:,375) = 5.4e-11_r8
91 58907904 : rate(:,376) = 3.77e-11_r8
92 58907904 : rate(:,379) = 2e-12_r8
93 58907904 : rate(:,380) = 1.29e-11_r8
94 58907904 : rate(:,382) = 4.5e-14_r8
95 58907904 : rate(:,387) = 3.77e-11_r8
96 58907904 : rate(:,393) = 4e-12_r8
97 58907904 : rate(:,399) = 1.78e-12_r8
98 58907904 : rate(:,401) = 6.1e-13_r8
99 58907904 : rate(:,405) = 4.8e-11_r8
100 58907904 : rate(:,408) = 1.6e-12_r8
101 58907904 : rate(:,410) = 6.7e-12_r8
102 58907904 : rate(:,413) = 3.5e-12_r8
103 58907904 : rate(:,418) = 6.42e-11_r8
104 58907904 : rate(:,425) = 1.6e-13_r8
105 58907904 : rate(:,431) = 1.4e-12_r8
106 58907904 : rate(:,436) = 7.5e-13_r8
107 58907904 : rate(:,437) = 1.4e-13_r8
108 58907904 : rate(:,438) = 7.5e-13_r8
109 58907904 : rate(:,439) = 3.6e-13_r8
110 58907904 : rate(:,440) = 6.5e-13_r8
111 58907904 : rate(:,441) = 2.1e-13_r8
112 58907904 : rate(:,442) = 6.5e-13_r8
113 58907904 : rate(:,443) = 4.9e-13_r8
114 58907904 : rate(:,445) = 1.2e-12_r8
115 58907904 : rate(:,449) = 9.8e-13_r8
116 58907904 : rate(:,452) = 1.85e-11_r8
117 58907904 : rate(:,453) = 1.63e-12_r8
118 58907904 : rate(:,454) = 2.5e-11_r8
119 58907904 : rate(:,455) = 1.1e-11_r8
120 58907904 : rate(:,456) = 3.3e-11_r8
121 58907904 : rate(:,459) = 2.8e-17_r8
122 58907904 : rate(:,460) = 8e-11_r8
123 58907904 : rate(:,463) = 3e-11_r8
124 58907904 : rate(:,466) = 4.2e-11_r8
125 58907904 : rate(:,469) = 2.8e-17_r8
126 58907904 : rate(:,470) = 1.1e-10_r8
127 58907904 : rate(:,472) = 3.9e-11_r8
128 58907904 : rate(:,475) = 1.3e-12_r8
129 58907904 : rate(:,477) = 5e-12_r8
130 58907904 : rate(:,478) = 2.3e-12_r8
131 58907904 : rate(:,481) = 3.9e-11_r8
132 58907904 : rate(:,484) = 2.8e-17_r8
133 58907904 : rate(:,485) = 9.2e-11_r8
134 58907904 : rate(:,488) = 3.85e-11_r8
135 58907904 : rate(:,492) = 1.2e-12_r8
136 58907904 : rate(:,496) = 9.8e-13_r8
137 58907904 : rate(:,501) = 4.4e-18_r8
138 58907904 : rate(:,502) = 3.6e-11_r8
139 58907904 : rate(:,554) = 4.7e-11_r8
140 58907904 : rate(:,567) = 2.1e-12_r8
141 58907904 : rate(:,568) = 2.8e-13_r8
142 58907904 : rate(:,576) = 1.7e-11_r8
143 58907904 : rate(:,582) = 8.4e-11_r8
144 58907904 : rate(:,585) = 5.3e-13_r8
145 58907904 : rate(:,587) = 2e-12_r8
146 58907904 : rate(:,590) = 2.3e-12_r8
147 58907904 : rate(:,595) = 2e-12_r8
148 58907904 : rate(:,598) = 2.3e-12_r8
149 58907904 : rate(:,604) = 1.9e-11_r8
150 58907904 : rate(:,605) = 5.3e-13_r8
151 58907904 : rate(:,607) = 2e-12_r8
152 58907904 : rate(:,610) = 2.3e-12_r8
153 58907904 : rate(:,615) = 2e-12_r8
154 58907904 : rate(:,618) = 2.3e-12_r8
155 58907904 : rate(:,622) = 1.2e-14_r8
156 58907904 : rate(:,623) = 2e-10_r8
157 58907904 : rate(:,624) = 2.5e-12_r8
158 58907904 : rate(:,625) = 5.3e-13_r8
159 58907904 : rate(:,627) = 2e-12_r8
160 58907904 : rate(:,630) = 2.3e-12_r8
161 58907904 : rate(:,635) = 2e-12_r8
162 58907904 : rate(:,638) = 2.3e-12_r8
163 58907904 : rate(:,644) = 1.2e-11_r8
164 58907904 : rate(:,646) = 2e-12_r8
165 58907904 : rate(:,648) = 5.3e-13_r8
166 58907904 : rate(:,650) = 2.3e-12_r8
167 58907904 : rate(:,655) = 2e-12_r8
168 58907904 : rate(:,658) = 2.3e-12_r8
169 58907904 : rate(:,664) = 1.1e-11_r8
170 58907904 : rate(:,666) = 2e-12_r8
171 58907904 : rate(:,668) = 5.3e-13_r8
172 58907904 : rate(:,670) = 2.3e-12_r8
173 58907904 : rate(:,675) = 2e-12_r8
174 58907904 : rate(:,678) = 2.3e-12_r8
175 58907904 : rate(:,683) = 2.1e-10_r8
176 58907904 : rate(:,689) = 8.9e-11_r8
177 58907904 : rate(:,690) = 8.9e-11_r8
178 58907904 : rate(:,694) = 2e-12_r8
179 58907904 : rate(:,697) = 2.3e-12_r8
180 58907904 : rate(:,705) = 4e-12_r8
181 58907904 : rate(:,708) = 2e-14_r8
182 58907904 : rate(:,710) = 2e-12_r8
183 58907904 : rate(:,713) = 2.3e-12_r8
184 58907904 : rate(:,718) = 2.52e-11_r8
185 58907904 : rate(:,723) = 4e-12_r8
186 58907904 : rate(:,727) = 2e-14_r8
187 58907904 : rate(:,729) = 2e-12_r8
188 58907904 : rate(:,732) = 2.3e-12_r8
189 58907904 : rate(:,737) = 1.92e-11_r8
190 58907904 : rate(:,739) = 2e-12_r8
191 58907904 : rate(:,742) = 2.3e-12_r8
192 58907904 : rate(:,746) = 8.8e-12_r8
193 58907904 : rate(:,747) = 8.8e-12_r8
194 58907904 : rate(:,748) = 8.8e-12_r8
195 58907904 : rate(:,753) = 4e-12_r8
196 58907904 : rate(:,755) = 2e-14_r8
197 58907904 : rate(:,757) = 3.66e-12_r8
198 58907904 : rate(:,758) = 2.8e-11_r8
199 58907904 : rate(:,759) = 2.6e-13_r8
200 58907904 : rate(:,762) = 8.3e-18_r8
201 58907904 : rate(:,763) = 1.1e-10_r8
202 58907904 : rate(:,767) = 1.1e-16_r8
203 58907904 : rate(:,769) = 3.64e-12_r8
204 58907904 : rate(:,770) = 2.8e-11_r8
205 58907904 : rate(:,771) = 1.7e-11_r8
206 58907904 : rate(:,774) = 1.1e-10_r8
207 58907904 : rate(:,775) = 9.58e-12_r8
208 58907904 : rate(:,778) = 1.1e-10_r8
209 58907904 : rate(:,779) = 1.23e-11_r8
210 58907904 : rate(:,782) = 1.1e-10_r8
211 58907904 : rate(:,783) = 3.64e-12_r8
212 58907904 : rate(:,786) = 1.1e-10_r8
213 58907904 : rate(:,787) = 5.5e-12_r8
214 58907904 : rate(:,788) = 4.65e-11_r8
215 58907904 : rate(:,789) = 2.8e-11_r8
216 58907904 : rate(:,797) = 2.3e-12_r8
217 58907904 : rate(:,799) = 1.2e-11_r8
218 58907904 : rate(:,800) = 5.7e-11_r8
219 58907904 : rate(:,801) = 2.8e-11_r8
220 58907904 : rate(:,802) = 6.6e-11_r8
221 58907904 : rate(:,803) = 1.4e-11_r8
222 58907904 : rate(:,806) = 1.9e-12_r8
223 58907904 : rate(:,829) = 6.34e-08_r8
224 58907904 : rate(:,846) = 1.9e-11_r8
225 58907904 : rate(:,849) = 1.2e-14_r8
226 58907904 : rate(:,850) = 2e-10_r8
227 58907904 : rate(:,854) = 2.5e-12_r8
228 58907904 : rate(:,866) = 1.34e-11_r8
229 58907904 : rate(:,867) = 1.2e-11_r8
230 58907904 : rate(:,872) = 1.1e-11_r8
231 58907904 : rate(:,876) = 2.1e-10_r8
232 58907904 : rate(:,877) = 1.34e-11_r8
233 58907904 : rate(:,881) = 1.7e-11_r8
234 58907904 : rate(:,901) = 1.29e-07_r8
235 58907904 : rate(:,902) = 2.31e-07_r8
236 58907904 : rate(:,903) = 2.31e-06_r8
237 58907904 : rate(:,904) = 4.63e-07_r8
238 :
239 4158720 : do n = 1,pver
240 4085760 : offset = (n-1)*ncol
241 62993664 : itemp(offset+1:offset+ncol) = 1._r8 / temp(:ncol,n)
242 : end do
243 :
244 58907904 : rate(:,169) = 1.63e-10_r8 * exp( 60._r8 * itemp(:) )
245 58907904 : rate(:,170) = 2.15e-11_r8 * exp( 110._r8 * itemp(:) )
246 58907904 : rate(:,171) = 3.3e-11_r8 * exp( 55._r8 * itemp(:) )
247 58907904 : rate(:,173) = 8e-12_r8 * exp( -2060._r8 * itemp(:) )
248 58907904 : rate(:,176) = 1.6e-11_r8 * exp( -4570._r8 * itemp(:) )
249 58907904 : rate(:,177) = 1.4e-12_r8 * exp( -2000._r8 * itemp(:) )
250 58907904 : exp_fac(:) = exp( 200._r8 * itemp(:) )
251 58907904 : rate(:,182) = 3e-11_r8 * exp_fac(:)
252 58907904 : rate(:,268) = 5.5e-12_r8 * exp_fac(:)
253 58907904 : rate(:,307) = 3.8e-12_r8 * exp_fac(:)
254 58907904 : rate(:,329) = 3.8e-12_r8 * exp_fac(:)
255 58907904 : rate(:,358) = 3.8e-12_r8 * exp_fac(:)
256 58907904 : rate(:,367) = 3.8e-12_r8 * exp_fac(:)
257 58907904 : rate(:,371) = 3.8e-12_r8 * exp_fac(:)
258 58907904 : rate(:,397) = 3.8e-12_r8 * exp_fac(:)
259 58907904 : rate(:,412) = 3.8e-12_r8 * exp_fac(:)
260 58907904 : rate(:,489) = 5.53e-12_r8 * exp_fac(:)
261 58907904 : rate(:,546) = 3.8e-12_r8 * exp_fac(:)
262 58907904 : rate(:,549) = 3.8e-12_r8 * exp_fac(:)
263 58907904 : rate(:,553) = 3.8e-12_r8 * exp_fac(:)
264 58907904 : rate(:,569) = 3.8e-12_r8 * exp_fac(:)
265 58907904 : rate(:,573) = 3.8e-12_r8 * exp_fac(:)
266 58907904 : rate(:,579) = 3.8e-12_r8 * exp_fac(:)
267 58907904 : rate(:,583) = 3.8e-12_r8 * exp_fac(:)
268 58907904 : rate(:,183) = 1e-14_r8 * exp( -490._r8 * itemp(:) )
269 58907904 : rate(:,184) = 1.4e-10_r8 * exp( -470._r8 * itemp(:) )
270 58907904 : rate(:,185) = 2.8e-12_r8 * exp( -1800._r8 * itemp(:) )
271 58907904 : exp_fac(:) = exp( 250._r8 * itemp(:) )
272 58907904 : rate(:,187) = 4.8e-11_r8 * exp_fac(:)
273 58907904 : rate(:,266) = 1.7e-11_r8 * exp_fac(:)
274 58907904 : exp_fac(:) = exp( 180._r8 * itemp(:) )
275 58907904 : rate(:,188) = 1.8e-11_r8 * exp_fac(:)
276 58907904 : rate(:,342) = 4.2e-12_r8 * exp_fac(:)
277 58907904 : rate(:,357) = 4.2e-12_r8 * exp_fac(:)
278 58907904 : rate(:,366) = 4.2e-12_r8 * exp_fac(:)
279 58907904 : rate(:,395) = 4.2e-12_r8 * exp_fac(:)
280 58907904 : rate(:,189) = 1.7e-12_r8 * exp( -940._r8 * itemp(:) )
281 58907904 : rate(:,193) = 4.5e-13_r8 * exp( 610._r8 * itemp(:) )
282 58907904 : rate(:,194) = 2.1e-11_r8 * exp( 100._r8 * itemp(:) )
283 58907904 : exp_fac(:) = exp( 220._r8 * itemp(:) )
284 58907904 : rate(:,195) = 2.9e-12_r8 * exp_fac(:)
285 58907904 : rate(:,196) = 1.45e-12_r8 * exp_fac(:)
286 58907904 : rate(:,197) = 1.45e-12_r8 * exp_fac(:)
287 58907904 : rate(:,198) = 3.3e-12_r8 * exp( -3150._r8 * itemp(:) )
288 58907904 : rate(:,199) = 5.1e-12_r8 * exp( 210._r8 * itemp(:) )
289 58907904 : exp_fac(:) = exp( -2450._r8 * itemp(:) )
290 58907904 : rate(:,200) = 1.2e-13_r8 * exp_fac(:)
291 58907904 : rate(:,226) = 3e-11_r8 * exp_fac(:)
292 58907904 : exp_fac(:) = exp( 125._r8 * itemp(:) )
293 58907904 : rate(:,203) = 1.7e-11_r8 * exp_fac(:)
294 58907904 : rate(:,301) = 5.5e-12_r8 * exp_fac(:)
295 58907904 : exp_fac(:) = exp( 260._r8 * itemp(:) )
296 58907904 : rate(:,207) = 3.44e-12_r8 * exp_fac(:)
297 58907904 : rate(:,259) = 2.3e-12_r8 * exp_fac(:)
298 58907904 : rate(:,262) = 8.8e-12_r8 * exp_fac(:)
299 58907904 : exp_fac(:) = exp( -1500._r8 * itemp(:) )
300 58907904 : rate(:,208) = 3e-12_r8 * exp_fac(:)
301 58907904 : rate(:,267) = 5.8e-12_r8 * exp_fac(:)
302 58907904 : exp_fac(:) = exp( 20._r8 * itemp(:) )
303 58907904 : rate(:,210) = 7.26e-11_r8 * exp_fac(:)
304 58907904 : rate(:,211) = 4.64e-11_r8 * exp_fac(:)
305 58907904 : rate(:,218) = 8.1e-11_r8 * exp( -30._r8 * itemp(:) )
306 58907904 : exp_fac(:) = exp( -1270._r8 * itemp(:) )
307 58907904 : rate(:,219) = 7.1e-12_r8 * exp_fac(:)
308 58907904 : rate(:,642) = 1.35e-15_r8 * exp_fac(:)
309 58907904 : rate(:,857) = 1.35e-15_r8 * exp_fac(:)
310 58907904 : rate(:,220) = 3.05e-11_r8 * exp( -2270._r8 * itemp(:) )
311 58907904 : rate(:,221) = 1.1e-11_r8 * exp( -980._r8 * itemp(:) )
312 58907904 : exp_fac(:) = exp( 270._r8 * itemp(:) )
313 58907904 : rate(:,222) = 1.4e-11_r8 * exp_fac(:)
314 58907904 : rate(:,236) = 7.4e-12_r8 * exp_fac(:)
315 58907904 : rate(:,338) = 8.1e-12_r8 * exp_fac(:)
316 58907904 : rate(:,392) = 8.1e-12_r8 * exp_fac(:)
317 58907904 : rate(:,704) = 8.1e-12_r8 * exp_fac(:)
318 58907904 : rate(:,722) = 8.1e-12_r8 * exp_fac(:)
319 58907904 : rate(:,752) = 8.1e-12_r8 * exp_fac(:)
320 58907904 : rate(:,223) = 3.6e-11_r8 * exp( -375._r8 * itemp(:) )
321 58907904 : rate(:,224) = 2.3e-11_r8 * exp( -200._r8 * itemp(:) )
322 58907904 : rate(:,225) = 3.3e-12_r8 * exp( -115._r8 * itemp(:) )
323 58907904 : rate(:,227) = 1e-12_r8 * exp( -1590._r8 * itemp(:) )
324 58907904 : rate(:,228) = 3.5e-13_r8 * exp( -1370._r8 * itemp(:) )
325 58907904 : exp_fac(:) = exp( 290._r8 * itemp(:) )
326 58907904 : rate(:,229) = 2.6e-12_r8 * exp_fac(:)
327 58907904 : rate(:,230) = 6.4e-12_r8 * exp_fac(:)
328 58907904 : rate(:,260) = 4.1e-13_r8 * exp_fac(:)
329 58907904 : rate(:,542) = 7.5e-12_r8 * exp_fac(:)
330 58907904 : rate(:,556) = 7.5e-12_r8 * exp_fac(:)
331 58907904 : rate(:,559) = 7.5e-12_r8 * exp_fac(:)
332 58907904 : rate(:,562) = 7.5e-12_r8 * exp_fac(:)
333 58907904 : rate(:,231) = 6.5e-12_r8 * exp( 135._r8 * itemp(:) )
334 58907904 : exp_fac(:) = exp( -840._r8 * itemp(:) )
335 58907904 : rate(:,233) = 3.6e-12_r8 * exp_fac(:)
336 58907904 : rate(:,282) = 2e-12_r8 * exp_fac(:)
337 58907904 : rate(:,234) = 1.2e-12_r8 * exp( -330._r8 * itemp(:) )
338 58907904 : rate(:,235) = 2.8e-11_r8 * exp( 85._r8 * itemp(:) )
339 58907904 : exp_fac(:) = exp( 230._r8 * itemp(:) )
340 58907904 : rate(:,237) = 6e-13_r8 * exp_fac(:)
341 58907904 : rate(:,257) = 1.5e-12_r8 * exp_fac(:)
342 58907904 : rate(:,265) = 1.9e-11_r8 * exp_fac(:)
343 58907904 : rate(:,238) = 1e-11_r8 * exp( -3300._r8 * itemp(:) )
344 58907904 : rate(:,239) = 1.8e-12_r8 * exp( -250._r8 * itemp(:) )
345 58907904 : rate(:,240) = 3.4e-12_r8 * exp( -130._r8 * itemp(:) )
346 58907904 : exp_fac(:) = exp( -500._r8 * itemp(:) )
347 58907904 : rate(:,242) = 3e-12_r8 * exp_fac(:)
348 58907904 : rate(:,276) = 1.4e-10_r8 * exp_fac(:)
349 58907904 : exp_fac(:) = exp( -800._r8 * itemp(:) )
350 58907904 : rate(:,254) = 1.7e-11_r8 * exp_fac(:)
351 58907904 : rate(:,281) = 6.3e-12_r8 * exp_fac(:)
352 58907904 : rate(:,255) = 4.8e-12_r8 * exp( -310._r8 * itemp(:) )
353 58907904 : rate(:,256) = 1.6e-11_r8 * exp( -780._r8 * itemp(:) )
354 58907904 : rate(:,258) = 9.5e-13_r8 * exp( 550._r8 * itemp(:) )
355 58907904 : exp_fac(:) = exp( 460._r8 * itemp(:) )
356 58907904 : rate(:,261) = 4.5e-12_r8 * exp_fac(:)
357 58907904 : rate(:,643) = 1.62e-11_r8 * exp_fac(:)
358 58907904 : rate(:,858) = 1.62e-11_r8 * exp_fac(:)
359 58907904 : rate(:,264) = 1.9e-11_r8 * exp( 215._r8 * itemp(:) )
360 58907904 : rate(:,269) = 1.2e-10_r8 * exp( -430._r8 * itemp(:) )
361 58907904 : rate(:,275) = 1.6e-10_r8 * exp( -260._r8 * itemp(:) )
362 58907904 : exp_fac(:) = exp( 0._r8 * itemp(:) )
363 58907904 : rate(:,277) = 1.4e-11_r8 * exp_fac(:)
364 58907904 : rate(:,279) = 2.14e-11_r8 * exp_fac(:)
365 58907904 : rate(:,280) = 1.9e-10_r8 * exp_fac(:)
366 58907904 : rate(:,293) = 2.57e-10_r8 * exp_fac(:)
367 58907904 : rate(:,294) = 1.8e-10_r8 * exp_fac(:)
368 58907904 : rate(:,295) = 1.794e-10_r8 * exp_fac(:)
369 58907904 : rate(:,296) = 1.3e-10_r8 * exp_fac(:)
370 58907904 : rate(:,297) = 7.65e-11_r8 * exp_fac(:)
371 58907904 : rate(:,310) = 4e-13_r8 * exp_fac(:)
372 58907904 : rate(:,315) = 1.31e-10_r8 * exp_fac(:)
373 58907904 : rate(:,316) = 3.5e-11_r8 * exp_fac(:)
374 58907904 : rate(:,317) = 9e-12_r8 * exp_fac(:)
375 58907904 : rate(:,324) = 6.8e-14_r8 * exp_fac(:)
376 58907904 : rate(:,325) = 2e-13_r8 * exp_fac(:)
377 58907904 : rate(:,340) = 1e-12_r8 * exp_fac(:)
378 58907904 : rate(:,344) = 1e-14_r8 * exp_fac(:)
379 58907904 : rate(:,345) = 1e-11_r8 * exp_fac(:)
380 58907904 : rate(:,346) = 1.15e-11_r8 * exp_fac(:)
381 58907904 : rate(:,347) = 3.3e-11_r8 * exp_fac(:)
382 58907904 : rate(:,348) = 3.4e-12_r8 * exp_fac(:)
383 58907904 : rate(:,349) = 4e-14_r8 * exp_fac(:)
384 58907904 : rate(:,362) = 3e-12_r8 * exp_fac(:)
385 58907904 : rate(:,363) = 1.2e-11_r8 * exp_fac(:)
386 58907904 : rate(:,364) = 6.7e-13_r8 * exp_fac(:)
387 58907904 : rate(:,374) = 3.5e-13_r8 * exp_fac(:)
388 58907904 : rate(:,375) = 5.4e-11_r8 * exp_fac(:)
389 58907904 : rate(:,376) = 3.77e-11_r8 * exp_fac(:)
390 58907904 : rate(:,379) = 2e-12_r8 * exp_fac(:)
391 58907904 : rate(:,380) = 1.29e-11_r8 * exp_fac(:)
392 58907904 : rate(:,382) = 4.5e-14_r8 * exp_fac(:)
393 58907904 : rate(:,387) = 3.77e-11_r8 * exp_fac(:)
394 58907904 : rate(:,393) = 4e-12_r8 * exp_fac(:)
395 58907904 : rate(:,399) = 1.78e-12_r8 * exp_fac(:)
396 58907904 : rate(:,401) = 6.1e-13_r8 * exp_fac(:)
397 58907904 : rate(:,405) = 4.8e-11_r8 * exp_fac(:)
398 58907904 : rate(:,408) = 1.6e-12_r8 * exp_fac(:)
399 58907904 : rate(:,410) = 6.7e-12_r8 * exp_fac(:)
400 58907904 : rate(:,413) = 3.5e-12_r8 * exp_fac(:)
401 58907904 : rate(:,418) = 6.42e-11_r8 * exp_fac(:)
402 58907904 : rate(:,425) = 1.6e-13_r8 * exp_fac(:)
403 58907904 : rate(:,431) = 1.4e-12_r8 * exp_fac(:)
404 58907904 : rate(:,436) = 7.5e-13_r8 * exp_fac(:)
405 58907904 : rate(:,437) = 1.4e-13_r8 * exp_fac(:)
406 58907904 : rate(:,438) = 7.5e-13_r8 * exp_fac(:)
407 58907904 : rate(:,439) = 3.6e-13_r8 * exp_fac(:)
408 58907904 : rate(:,440) = 6.5e-13_r8 * exp_fac(:)
409 58907904 : rate(:,441) = 2.1e-13_r8 * exp_fac(:)
410 58907904 : rate(:,442) = 6.5e-13_r8 * exp_fac(:)
411 58907904 : rate(:,443) = 4.9e-13_r8 * exp_fac(:)
412 58907904 : rate(:,445) = 1.2e-12_r8 * exp_fac(:)
413 58907904 : rate(:,449) = 9.8e-13_r8 * exp_fac(:)
414 58907904 : rate(:,452) = 1.85e-11_r8 * exp_fac(:)
415 58907904 : rate(:,453) = 1.63e-12_r8 * exp_fac(:)
416 58907904 : rate(:,454) = 2.5e-11_r8 * exp_fac(:)
417 58907904 : rate(:,455) = 1.1e-11_r8 * exp_fac(:)
418 58907904 : rate(:,456) = 3.3e-11_r8 * exp_fac(:)
419 58907904 : rate(:,459) = 2.8e-17_r8 * exp_fac(:)
420 58907904 : rate(:,460) = 8e-11_r8 * exp_fac(:)
421 58907904 : rate(:,463) = 3e-11_r8 * exp_fac(:)
422 58907904 : rate(:,466) = 4.2e-11_r8 * exp_fac(:)
423 58907904 : rate(:,469) = 2.8e-17_r8 * exp_fac(:)
424 58907904 : rate(:,470) = 1.1e-10_r8 * exp_fac(:)
425 58907904 : rate(:,472) = 3.9e-11_r8 * exp_fac(:)
426 58907904 : rate(:,475) = 1.3e-12_r8 * exp_fac(:)
427 58907904 : rate(:,477) = 5e-12_r8 * exp_fac(:)
428 58907904 : rate(:,478) = 2.3e-12_r8 * exp_fac(:)
429 58907904 : rate(:,481) = 3.9e-11_r8 * exp_fac(:)
430 58907904 : rate(:,484) = 2.8e-17_r8 * exp_fac(:)
431 58907904 : rate(:,485) = 9.2e-11_r8 * exp_fac(:)
432 58907904 : rate(:,488) = 3.85e-11_r8 * exp_fac(:)
433 58907904 : rate(:,492) = 1.2e-12_r8 * exp_fac(:)
434 58907904 : rate(:,496) = 9.8e-13_r8 * exp_fac(:)
435 58907904 : rate(:,501) = 4.4e-18_r8 * exp_fac(:)
436 58907904 : rate(:,502) = 3.6e-11_r8 * exp_fac(:)
437 58907904 : rate(:,554) = 4.7e-11_r8 * exp_fac(:)
438 58907904 : rate(:,567) = 2.1e-12_r8 * exp_fac(:)
439 58907904 : rate(:,568) = 2.8e-13_r8 * exp_fac(:)
440 58907904 : rate(:,576) = 1.7e-11_r8 * exp_fac(:)
441 58907904 : rate(:,582) = 8.4e-11_r8 * exp_fac(:)
442 58907904 : rate(:,585) = 5.3e-13_r8 * exp_fac(:)
443 58907904 : rate(:,587) = 2e-12_r8 * exp_fac(:)
444 58907904 : rate(:,590) = 2.3e-12_r8 * exp_fac(:)
445 58907904 : rate(:,595) = 2e-12_r8 * exp_fac(:)
446 58907904 : rate(:,598) = 2.3e-12_r8 * exp_fac(:)
447 58907904 : rate(:,604) = 1.9e-11_r8 * exp_fac(:)
448 58907904 : rate(:,605) = 5.3e-13_r8 * exp_fac(:)
449 58907904 : rate(:,607) = 2e-12_r8 * exp_fac(:)
450 58907904 : rate(:,610) = 2.3e-12_r8 * exp_fac(:)
451 58907904 : rate(:,615) = 2e-12_r8 * exp_fac(:)
452 58907904 : rate(:,618) = 2.3e-12_r8 * exp_fac(:)
453 58907904 : rate(:,622) = 1.2e-14_r8 * exp_fac(:)
454 58907904 : rate(:,623) = 2e-10_r8 * exp_fac(:)
455 58907904 : rate(:,624) = 2.5e-12_r8 * exp_fac(:)
456 58907904 : rate(:,625) = 5.3e-13_r8 * exp_fac(:)
457 58907904 : rate(:,627) = 2e-12_r8 * exp_fac(:)
458 58907904 : rate(:,630) = 2.3e-12_r8 * exp_fac(:)
459 58907904 : rate(:,635) = 2e-12_r8 * exp_fac(:)
460 58907904 : rate(:,638) = 2.3e-12_r8 * exp_fac(:)
461 58907904 : rate(:,644) = 1.2e-11_r8 * exp_fac(:)
462 58907904 : rate(:,646) = 2e-12_r8 * exp_fac(:)
463 58907904 : rate(:,648) = 5.3e-13_r8 * exp_fac(:)
464 58907904 : rate(:,650) = 2.3e-12_r8 * exp_fac(:)
465 58907904 : rate(:,655) = 2e-12_r8 * exp_fac(:)
466 58907904 : rate(:,658) = 2.3e-12_r8 * exp_fac(:)
467 58907904 : rate(:,664) = 1.1e-11_r8 * exp_fac(:)
468 58907904 : rate(:,666) = 2e-12_r8 * exp_fac(:)
469 58907904 : rate(:,668) = 5.3e-13_r8 * exp_fac(:)
470 58907904 : rate(:,670) = 2.3e-12_r8 * exp_fac(:)
471 58907904 : rate(:,675) = 2e-12_r8 * exp_fac(:)
472 58907904 : rate(:,678) = 2.3e-12_r8 * exp_fac(:)
473 58907904 : rate(:,683) = 2.1e-10_r8 * exp_fac(:)
474 58907904 : rate(:,689) = 8.9e-11_r8 * exp_fac(:)
475 58907904 : rate(:,690) = 8.9e-11_r8 * exp_fac(:)
476 58907904 : rate(:,694) = 2e-12_r8 * exp_fac(:)
477 58907904 : rate(:,697) = 2.3e-12_r8 * exp_fac(:)
478 58907904 : rate(:,705) = 4e-12_r8 * exp_fac(:)
479 58907904 : rate(:,708) = 2e-14_r8 * exp_fac(:)
480 58907904 : rate(:,710) = 2e-12_r8 * exp_fac(:)
481 58907904 : rate(:,713) = 2.3e-12_r8 * exp_fac(:)
482 58907904 : rate(:,718) = 2.52e-11_r8 * exp_fac(:)
483 58907904 : rate(:,723) = 4e-12_r8 * exp_fac(:)
484 58907904 : rate(:,727) = 2e-14_r8 * exp_fac(:)
485 58907904 : rate(:,729) = 2e-12_r8 * exp_fac(:)
486 58907904 : rate(:,732) = 2.3e-12_r8 * exp_fac(:)
487 58907904 : rate(:,737) = 1.92e-11_r8 * exp_fac(:)
488 58907904 : rate(:,739) = 2e-12_r8 * exp_fac(:)
489 58907904 : rate(:,742) = 2.3e-12_r8 * exp_fac(:)
490 58907904 : rate(:,746) = 8.8e-12_r8 * exp_fac(:)
491 58907904 : rate(:,747) = 8.8e-12_r8 * exp_fac(:)
492 58907904 : rate(:,748) = 8.8e-12_r8 * exp_fac(:)
493 58907904 : rate(:,753) = 4e-12_r8 * exp_fac(:)
494 58907904 : rate(:,755) = 2e-14_r8 * exp_fac(:)
495 58907904 : rate(:,757) = 3.66e-12_r8 * exp_fac(:)
496 58907904 : rate(:,758) = 2.8e-11_r8 * exp_fac(:)
497 58907904 : rate(:,759) = 2.6e-13_r8 * exp_fac(:)
498 58907904 : rate(:,762) = 8.3e-18_r8 * exp_fac(:)
499 58907904 : rate(:,763) = 1.1e-10_r8 * exp_fac(:)
500 58907904 : rate(:,767) = 1.1e-16_r8 * exp_fac(:)
501 58907904 : rate(:,769) = 3.64e-12_r8 * exp_fac(:)
502 58907904 : rate(:,770) = 2.8e-11_r8 * exp_fac(:)
503 58907904 : rate(:,771) = 1.7e-11_r8 * exp_fac(:)
504 58907904 : rate(:,774) = 1.1e-10_r8 * exp_fac(:)
505 58907904 : rate(:,775) = 9.58e-12_r8 * exp_fac(:)
506 58907904 : rate(:,778) = 1.1e-10_r8 * exp_fac(:)
507 58907904 : rate(:,779) = 1.23e-11_r8 * exp_fac(:)
508 58907904 : rate(:,782) = 1.1e-10_r8 * exp_fac(:)
509 58907904 : rate(:,783) = 3.64e-12_r8 * exp_fac(:)
510 58907904 : rate(:,786) = 1.1e-10_r8 * exp_fac(:)
511 58907904 : rate(:,787) = 5.5e-12_r8 * exp_fac(:)
512 58907904 : rate(:,788) = 4.65e-11_r8 * exp_fac(:)
513 58907904 : rate(:,789) = 2.8e-11_r8 * exp_fac(:)
514 58907904 : rate(:,797) = 2.3e-12_r8 * exp_fac(:)
515 58907904 : rate(:,799) = 1.2e-11_r8 * exp_fac(:)
516 58907904 : rate(:,800) = 5.7e-11_r8 * exp_fac(:)
517 58907904 : rate(:,801) = 2.8e-11_r8 * exp_fac(:)
518 58907904 : rate(:,802) = 6.6e-11_r8 * exp_fac(:)
519 58907904 : rate(:,803) = 1.4e-11_r8 * exp_fac(:)
520 58907904 : rate(:,806) = 1.9e-12_r8 * exp_fac(:)
521 58907904 : rate(:,829) = 6.34e-08_r8 * exp_fac(:)
522 58907904 : rate(:,846) = 1.9e-11_r8 * exp_fac(:)
523 58907904 : rate(:,849) = 1.2e-14_r8 * exp_fac(:)
524 58907904 : rate(:,850) = 2e-10_r8 * exp_fac(:)
525 58907904 : rate(:,854) = 2.5e-12_r8 * exp_fac(:)
526 58907904 : rate(:,866) = 1.34e-11_r8 * exp_fac(:)
527 58907904 : rate(:,867) = 1.2e-11_r8 * exp_fac(:)
528 58907904 : rate(:,872) = 1.1e-11_r8 * exp_fac(:)
529 58907904 : rate(:,876) = 2.1e-10_r8 * exp_fac(:)
530 58907904 : rate(:,877) = 1.34e-11_r8 * exp_fac(:)
531 58907904 : rate(:,881) = 1.7e-11_r8 * exp_fac(:)
532 58907904 : rate(:,901) = 1.29e-07_r8 * exp_fac(:)
533 58907904 : rate(:,902) = 2.31e-07_r8 * exp_fac(:)
534 58907904 : rate(:,903) = 2.31e-06_r8 * exp_fac(:)
535 58907904 : rate(:,904) = 4.63e-07_r8 * exp_fac(:)
536 58907904 : rate(:,278) = 6e-12_r8 * exp( 400._r8 * itemp(:) )
537 58907904 : rate(:,283) = 1.46e-11_r8 * exp( -1040._r8 * itemp(:) )
538 58907904 : rate(:,284) = 1.42e-12_r8 * exp( -1150._r8 * itemp(:) )
539 58907904 : exp_fac(:) = exp( -1520._r8 * itemp(:) )
540 58907904 : rate(:,285) = 1.64e-12_r8 * exp_fac(:)
541 58907904 : rate(:,403) = 8.5e-16_r8 * exp_fac(:)
542 58907904 : exp_fac(:) = exp( -1100._r8 * itemp(:) )
543 58907904 : rate(:,286) = 2.03e-11_r8 * exp_fac(:)
544 58907904 : rate(:,805) = 3.4e-12_r8 * exp_fac(:)
545 58907904 : rate(:,287) = 1.96e-12_r8 * exp( -1200._r8 * itemp(:) )
546 58907904 : rate(:,288) = 4.85e-12_r8 * exp( -850._r8 * itemp(:) )
547 58907904 : rate(:,289) = 9e-13_r8 * exp( -360._r8 * itemp(:) )
548 58907904 : exp_fac(:) = exp( -1600._r8 * itemp(:) )
549 58907904 : rate(:,290) = 1.25e-12_r8 * exp_fac(:)
550 58907904 : rate(:,300) = 3.4e-11_r8 * exp_fac(:)
551 58907904 : rate(:,291) = 1.3e-12_r8 * exp( -1770._r8 * itemp(:) )
552 58907904 : rate(:,292) = 9.2e-13_r8 * exp( -1560._r8 * itemp(:) )
553 58907904 : rate(:,298) = 9.7e-15_r8 * exp( 625._r8 * itemp(:) )
554 58907904 : rate(:,299) = 6e-13_r8 * exp( -2058._r8 * itemp(:) )
555 58907904 : rate(:,302) = 5e-13_r8 * exp( -424._r8 * itemp(:) )
556 58907904 : rate(:,303) = 1.9e-14_r8 * exp( 706._r8 * itemp(:) )
557 58907904 : rate(:,304) = 4.1e-13_r8 * exp( 750._r8 * itemp(:) )
558 58907904 : exp_fac(:) = exp( 300._r8 * itemp(:) )
559 58907904 : rate(:,305) = 2.8e-12_r8 * exp_fac(:)
560 58907904 : rate(:,370) = 2.9e-12_r8 * exp_fac(:)
561 58907904 : rate(:,306) = 2.9e-12_r8 * exp( -345._r8 * itemp(:) )
562 58907904 : rate(:,308) = 2.45e-12_r8 * exp( -1775._r8 * itemp(:) )
563 58907904 : exp_fac(:) = exp( 500._r8 * itemp(:) )
564 58907904 : rate(:,311) = 1.3e-12_r8 * exp_fac(:)
565 58907904 : rate(:,335) = 2.9e-12_r8 * exp_fac(:)
566 58907904 : rate(:,336) = 2e-12_r8 * exp_fac(:)
567 58907904 : rate(:,368) = 7.1e-13_r8 * exp_fac(:)
568 58907904 : rate(:,381) = 2e-12_r8 * exp_fac(:)
569 58907904 : rate(:,388) = 2.9e-12_r8 * exp_fac(:)
570 58907904 : rate(:,389) = 2e-12_r8 * exp_fac(:)
571 58907904 : rate(:,391) = 2.9e-12_r8 * exp_fac(:)
572 58907904 : rate(:,400) = 2e-12_r8 * exp_fac(:)
573 58907904 : rate(:,424) = 2e-12_r8 * exp_fac(:)
574 58907904 : rate(:,430) = 2e-12_r8 * exp_fac(:)
575 58907904 : rate(:,444) = 2e-12_r8 * exp_fac(:)
576 58907904 : rate(:,448) = 2e-12_r8 * exp_fac(:)
577 58907904 : rate(:,474) = 2e-12_r8 * exp_fac(:)
578 58907904 : rate(:,491) = 2e-12_r8 * exp_fac(:)
579 58907904 : rate(:,495) = 2e-12_r8 * exp_fac(:)
580 58907904 : rate(:,586) = 2e-12_r8 * exp_fac(:)
581 58907904 : rate(:,591) = 2e-12_r8 * exp_fac(:)
582 58907904 : rate(:,592) = 2e-12_r8 * exp_fac(:)
583 58907904 : rate(:,593) = 2e-12_r8 * exp_fac(:)
584 58907904 : rate(:,594) = 2e-12_r8 * exp_fac(:)
585 58907904 : rate(:,599) = 2e-12_r8 * exp_fac(:)
586 58907904 : rate(:,600) = 2e-12_r8 * exp_fac(:)
587 58907904 : rate(:,601) = 2e-12_r8 * exp_fac(:)
588 58907904 : rate(:,606) = 2e-12_r8 * exp_fac(:)
589 58907904 : rate(:,611) = 2e-12_r8 * exp_fac(:)
590 58907904 : rate(:,612) = 2e-12_r8 * exp_fac(:)
591 58907904 : rate(:,613) = 2e-12_r8 * exp_fac(:)
592 58907904 : rate(:,614) = 2e-12_r8 * exp_fac(:)
593 58907904 : rate(:,619) = 2e-12_r8 * exp_fac(:)
594 58907904 : rate(:,620) = 2e-12_r8 * exp_fac(:)
595 58907904 : rate(:,621) = 2e-12_r8 * exp_fac(:)
596 58907904 : rate(:,626) = 2e-12_r8 * exp_fac(:)
597 58907904 : rate(:,631) = 2e-12_r8 * exp_fac(:)
598 58907904 : rate(:,632) = 2e-12_r8 * exp_fac(:)
599 58907904 : rate(:,633) = 2e-12_r8 * exp_fac(:)
600 58907904 : rate(:,634) = 2e-12_r8 * exp_fac(:)
601 58907904 : rate(:,639) = 2e-12_r8 * exp_fac(:)
602 58907904 : rate(:,640) = 2e-12_r8 * exp_fac(:)
603 58907904 : rate(:,641) = 2e-12_r8 * exp_fac(:)
604 58907904 : rate(:,645) = 2e-12_r8 * exp_fac(:)
605 58907904 : rate(:,651) = 2e-12_r8 * exp_fac(:)
606 58907904 : rate(:,652) = 2e-12_r8 * exp_fac(:)
607 58907904 : rate(:,653) = 2e-12_r8 * exp_fac(:)
608 58907904 : rate(:,654) = 2e-12_r8 * exp_fac(:)
609 58907904 : rate(:,659) = 2e-12_r8 * exp_fac(:)
610 58907904 : rate(:,660) = 2e-12_r8 * exp_fac(:)
611 58907904 : rate(:,661) = 2e-12_r8 * exp_fac(:)
612 58907904 : rate(:,665) = 2e-12_r8 * exp_fac(:)
613 58907904 : rate(:,671) = 2e-12_r8 * exp_fac(:)
614 58907904 : rate(:,672) = 2e-12_r8 * exp_fac(:)
615 58907904 : rate(:,673) = 2e-12_r8 * exp_fac(:)
616 58907904 : rate(:,674) = 2e-12_r8 * exp_fac(:)
617 58907904 : rate(:,679) = 2e-12_r8 * exp_fac(:)
618 58907904 : rate(:,680) = 2e-12_r8 * exp_fac(:)
619 58907904 : rate(:,681) = 2e-12_r8 * exp_fac(:)
620 58907904 : rate(:,693) = 2e-12_r8 * exp_fac(:)
621 58907904 : rate(:,698) = 2e-12_r8 * exp_fac(:)
622 58907904 : rate(:,699) = 2e-12_r8 * exp_fac(:)
623 58907904 : rate(:,700) = 2e-12_r8 * exp_fac(:)
624 58907904 : rate(:,701) = 2.9e-12_r8 * exp_fac(:)
625 58907904 : rate(:,702) = 2e-12_r8 * exp_fac(:)
626 58907904 : rate(:,706) = 2.9e-12_r8 * exp_fac(:)
627 58907904 : rate(:,707) = 2.9e-12_r8 * exp_fac(:)
628 58907904 : rate(:,709) = 2e-12_r8 * exp_fac(:)
629 58907904 : rate(:,714) = 2e-12_r8 * exp_fac(:)
630 58907904 : rate(:,715) = 2e-12_r8 * exp_fac(:)
631 58907904 : rate(:,716) = 2e-12_r8 * exp_fac(:)
632 58907904 : rate(:,719) = 2.9e-12_r8 * exp_fac(:)
633 58907904 : rate(:,720) = 2e-12_r8 * exp_fac(:)
634 58907904 : rate(:,724) = 2.9e-12_r8 * exp_fac(:)
635 58907904 : rate(:,725) = 2.9e-12_r8 * exp_fac(:)
636 58907904 : rate(:,726) = 2.9e-12_r8 * exp_fac(:)
637 58907904 : rate(:,728) = 2e-12_r8 * exp_fac(:)
638 58907904 : rate(:,733) = 2e-12_r8 * exp_fac(:)
639 58907904 : rate(:,734) = 2e-12_r8 * exp_fac(:)
640 58907904 : rate(:,735) = 2e-12_r8 * exp_fac(:)
641 58907904 : rate(:,738) = 2e-12_r8 * exp_fac(:)
642 58907904 : rate(:,743) = 2e-12_r8 * exp_fac(:)
643 58907904 : rate(:,744) = 2e-12_r8 * exp_fac(:)
644 58907904 : rate(:,745) = 2e-12_r8 * exp_fac(:)
645 58907904 : rate(:,749) = 2.9e-12_r8 * exp_fac(:)
646 58907904 : rate(:,750) = 2e-12_r8 * exp_fac(:)
647 58907904 : rate(:,754) = 2.9e-12_r8 * exp_fac(:)
648 58907904 : rate(:,312) = 5.6e-15_r8 * exp( 2300._r8 * itemp(:) )
649 58907904 : rate(:,313) = 2.4e+12_r8 * exp( -7000._r8 * itemp(:) )
650 58907904 : rate(:,314) = 2.6e-12_r8 * exp( 265._r8 * itemp(:) )
651 58907904 : rate(:,318) = 1.08e-10_r8 * exp( 105._r8 * itemp(:) )
652 58907904 : rate(:,323) = 1.2e-14_r8 * exp( -2630._r8 * itemp(:) )
653 58907904 : exp_fac(:) = exp( 700._r8 * itemp(:) )
654 58907904 : rate(:,326) = 7.5e-13_r8 * exp_fac(:)
655 58907904 : rate(:,341) = 7.5e-13_r8 * exp_fac(:)
656 58907904 : rate(:,356) = 7.5e-13_r8 * exp_fac(:)
657 58907904 : rate(:,365) = 7.5e-13_r8 * exp_fac(:)
658 58907904 : rate(:,369) = 8.6e-13_r8 * exp_fac(:)
659 58907904 : rate(:,394) = 7.5e-13_r8 * exp_fac(:)
660 58907904 : rate(:,409) = 7.5e-13_r8 * exp_fac(:)
661 58907904 : rate(:,544) = 7.5e-13_r8 * exp_fac(:)
662 58907904 : rate(:,548) = 7.5e-13_r8 * exp_fac(:)
663 58907904 : rate(:,551) = 7.5e-13_r8 * exp_fac(:)
664 58907904 : rate(:,564) = 7.5e-13_r8 * exp_fac(:)
665 58907904 : rate(:,571) = 7.5e-13_r8 * exp_fac(:)
666 58907904 : rate(:,577) = 7.5e-13_r8 * exp_fac(:)
667 58907904 : rate(:,580) = 7.5e-13_r8 * exp_fac(:)
668 58907904 : rate(:,852) = 7.5e-13_r8 * exp_fac(:)
669 58907904 : rate(:,864) = 7.5e-13_r8 * exp_fac(:)
670 58907904 : rate(:,879) = 7.5e-13_r8 * exp_fac(:)
671 58907904 : rate(:,882) = 7.5e-13_r8 * exp_fac(:)
672 58907904 : exp_fac(:) = exp( 365._r8 * itemp(:) )
673 58907904 : rate(:,327) = 2.6e-12_r8 * exp_fac(:)
674 58907904 : rate(:,545) = 2.6e-12_r8 * exp_fac(:)
675 58907904 : rate(:,550) = 2.6e-12_r8 * exp_fac(:)
676 58907904 : rate(:,552) = 2.6e-12_r8 * exp_fac(:)
677 58907904 : rate(:,565) = 2.6e-12_r8 * exp_fac(:)
678 58907904 : rate(:,572) = 2.6e-12_r8 * exp_fac(:)
679 58907904 : rate(:,578) = 2.6e-12_r8 * exp_fac(:)
680 58907904 : rate(:,581) = 2.6e-12_r8 * exp_fac(:)
681 58907904 : rate(:,853) = 2.6e-12_r8 * exp_fac(:)
682 58907904 : rate(:,865) = 2.6e-12_r8 * exp_fac(:)
683 58907904 : rate(:,880) = 2.6e-12_r8 * exp_fac(:)
684 58907904 : rate(:,883) = 2.6e-12_r8 * exp_fac(:)
685 58907904 : rate(:,328) = 6.9e-12_r8 * exp( -230._r8 * itemp(:) )
686 58907904 : rate(:,330) = 7.2e-11_r8 * exp( -70._r8 * itemp(:) )
687 58907904 : rate(:,331) = 7.66e-12_r8 * exp( -1020._r8 * itemp(:) )
688 58907904 : exp_fac(:) = exp( -1900._r8 * itemp(:) )
689 58907904 : rate(:,332) = 1.4e-12_r8 * exp_fac(:)
690 58907904 : rate(:,354) = 6.5e-15_r8 * exp_fac(:)
691 58907904 : rate(:,333) = 4.63e-12_r8 * exp( 350._r8 * itemp(:) )
692 58907904 : rate(:,334) = 7.8e-13_r8 * exp( -1050._r8 * itemp(:) )
693 58907904 : exp_fac(:) = exp( 1040._r8 * itemp(:) )
694 58907904 : rate(:,337) = 4.3e-13_r8 * exp_fac(:)
695 58907904 : rate(:,390) = 4.3e-13_r8 * exp_fac(:)
696 58907904 : rate(:,541) = 4.3e-13_r8 * exp_fac(:)
697 58907904 : rate(:,555) = 4.3e-13_r8 * exp_fac(:)
698 58907904 : rate(:,558) = 4.3e-13_r8 * exp_fac(:)
699 58907904 : rate(:,561) = 4.3e-13_r8 * exp_fac(:)
700 58907904 : rate(:,703) = 4.3e-13_r8 * exp_fac(:)
701 58907904 : rate(:,721) = 4.3e-13_r8 * exp_fac(:)
702 58907904 : rate(:,751) = 4.3e-13_r8 * exp_fac(:)
703 58907904 : rate(:,339) = 3.15e-14_r8 * exp( 920._r8 * itemp(:) )
704 58907904 : rate(:,343) = 1.6e+11_r8 * exp( -4150._r8 * itemp(:) )
705 58907904 : rate(:,353) = 4.6e-13_r8 * exp( -1156._r8 * itemp(:) )
706 58907904 : rate(:,355) = 3.75e-13_r8 * exp( -40._r8 * itemp(:) )
707 58907904 : rate(:,359) = 9.19e-12_r8 * exp( -630._r8 * itemp(:) )
708 58907904 : rate(:,360) = 1.4e-12_r8 * exp( -1860._r8 * itemp(:) )
709 58907904 : rate(:,361) = 8.4e-13_r8 * exp( 830._r8 * itemp(:) )
710 58907904 : rate(:,377) = 4.8e-12_r8 * exp( 120._r8 * itemp(:) )
711 58907904 : rate(:,378) = 5.1e-14_r8 * exp( 693._r8 * itemp(:) )
712 58907904 : exp_fac(:) = exp( 1300._r8 * itemp(:) )
713 58907904 : rate(:,383) = 2.11e-13_r8 * exp_fac(:)
714 58907904 : rate(:,402) = 2.11e-13_r8 * exp_fac(:)
715 58907904 : rate(:,421) = 2.38e-13_r8 * exp_fac(:)
716 58907904 : rate(:,426) = 2.12e-13_r8 * exp_fac(:)
717 58907904 : rate(:,432) = 2.12e-13_r8 * exp_fac(:)
718 58907904 : rate(:,446) = 2.12e-13_r8 * exp_fac(:)
719 58907904 : rate(:,450) = 2.12e-13_r8 * exp_fac(:)
720 58907904 : rate(:,457) = 2.6e-13_r8 * exp_fac(:)
721 58907904 : rate(:,461) = 2.6e-13_r8 * exp_fac(:)
722 58907904 : rate(:,464) = 2.6e-13_r8 * exp_fac(:)
723 58907904 : rate(:,467) = 2.6e-13_r8 * exp_fac(:)
724 58907904 : rate(:,471) = 2.6e-13_r8 * exp_fac(:)
725 58907904 : rate(:,476) = 2.47e-13_r8 * exp_fac(:)
726 58907904 : rate(:,479) = 2.64e-13_r8 * exp_fac(:)
727 58907904 : rate(:,482) = 2.64e-13_r8 * exp_fac(:)
728 58907904 : rate(:,493) = 2.12e-13_r8 * exp_fac(:)
729 58907904 : rate(:,497) = 2.12e-13_r8 * exp_fac(:)
730 58907904 : rate(:,499) = 2.6e-13_r8 * exp_fac(:)
731 58907904 : rate(:,588) = 2.71e-13_r8 * exp_fac(:)
732 58907904 : rate(:,596) = 2.6e-13_r8 * exp_fac(:)
733 58907904 : rate(:,608) = 2.78e-13_r8 * exp_fac(:)
734 58907904 : rate(:,616) = 2.75e-13_r8 * exp_fac(:)
735 58907904 : rate(:,628) = 2.71e-13_r8 * exp_fac(:)
736 58907904 : rate(:,636) = 2.6e-13_r8 * exp_fac(:)
737 58907904 : rate(:,647) = 2.71e-13_r8 * exp_fac(:)
738 58907904 : rate(:,656) = 2.6e-13_r8 * exp_fac(:)
739 58907904 : rate(:,667) = 2.71e-13_r8 * exp_fac(:)
740 58907904 : rate(:,676) = 2.6e-13_r8 * exp_fac(:)
741 58907904 : rate(:,687) = 2.71e-13_r8 * exp_fac(:)
742 58907904 : rate(:,691) = 2.71e-13_r8 * exp_fac(:)
743 58907904 : rate(:,695) = 2.54e-13_r8 * exp_fac(:)
744 58907904 : rate(:,711) = 2.62e-13_r8 * exp_fac(:)
745 58907904 : rate(:,730) = 2.66e-13_r8 * exp_fac(:)
746 58907904 : rate(:,740) = 2.51e-13_r8 * exp_fac(:)
747 58907904 : rate(:,760) = 2.68e-13_r8 * exp_fac(:)
748 58907904 : rate(:,765) = 2.47e-13_r8 * exp_fac(:)
749 58907904 : rate(:,772) = 2.76e-13_r8 * exp_fac(:)
750 58907904 : rate(:,776) = 2.76e-13_r8 * exp_fac(:)
751 58907904 : rate(:,780) = 2.75e-13_r8 * exp_fac(:)
752 58907904 : rate(:,784) = 2.75e-13_r8 * exp_fac(:)
753 58907904 : rate(:,842) = 2.6e-13_r8 * exp_fac(:)
754 58907904 : rate(:,847) = 2.75e-13_r8 * exp_fac(:)
755 58907904 : rate(:,855) = 2.6e-13_r8 * exp_fac(:)
756 58907904 : rate(:,860) = 2.12e-13_r8 * exp_fac(:)
757 58907904 : rate(:,868) = 2.6e-13_r8 * exp_fac(:)
758 58907904 : rate(:,873) = 2.6e-13_r8 * exp_fac(:)
759 58907904 : rate(:,384) = 2.9e+07_r8 * exp( -5297._r8 * itemp(:) )
760 58907904 : rate(:,385) = 1.5e-15_r8 * exp( -2100._r8 * itemp(:) )
761 58907904 : exp_fac(:) = exp( 360._r8 * itemp(:) )
762 58907904 : rate(:,386) = 9.6e-12_r8 * exp_fac(:)
763 58907904 : rate(:,589) = 2.7e-12_r8 * exp_fac(:)
764 58907904 : rate(:,597) = 2.7e-12_r8 * exp_fac(:)
765 58907904 : rate(:,609) = 2.7e-12_r8 * exp_fac(:)
766 58907904 : rate(:,617) = 2.7e-12_r8 * exp_fac(:)
767 58907904 : rate(:,629) = 2.7e-12_r8 * exp_fac(:)
768 58907904 : rate(:,637) = 2.7e-12_r8 * exp_fac(:)
769 58907904 : rate(:,649) = 2.7e-12_r8 * exp_fac(:)
770 58907904 : rate(:,657) = 2.7e-12_r8 * exp_fac(:)
771 58907904 : rate(:,669) = 2.7e-12_r8 * exp_fac(:)
772 58907904 : rate(:,677) = 2.7e-12_r8 * exp_fac(:)
773 58907904 : rate(:,688) = 2.7e-12_r8 * exp_fac(:)
774 58907904 : rate(:,692) = 2.7e-12_r8 * exp_fac(:)
775 58907904 : rate(:,696) = 2.7e-12_r8 * exp_fac(:)
776 58907904 : rate(:,712) = 2.7e-12_r8 * exp_fac(:)
777 58907904 : rate(:,731) = 2.7e-12_r8 * exp_fac(:)
778 58907904 : rate(:,741) = 2.7e-12_r8 * exp_fac(:)
779 58907904 : rate(:,761) = 2.7e-12_r8 * exp_fac(:)
780 58907904 : rate(:,766) = 2.7e-12_r8 * exp_fac(:)
781 58907904 : rate(:,773) = 2.7e-12_r8 * exp_fac(:)
782 58907904 : rate(:,777) = 2.7e-12_r8 * exp_fac(:)
783 58907904 : rate(:,781) = 2.7e-12_r8 * exp_fac(:)
784 58907904 : rate(:,785) = 2.7e-12_r8 * exp_fac(:)
785 58907904 : rate(:,843) = 2.7e-12_r8 * exp_fac(:)
786 58907904 : rate(:,848) = 2.7e-12_r8 * exp_fac(:)
787 58907904 : rate(:,856) = 2.7e-12_r8 * exp_fac(:)
788 58907904 : rate(:,861) = 2.7e-12_r8 * exp_fac(:)
789 58907904 : rate(:,869) = 2.7e-12_r8 * exp_fac(:)
790 58907904 : rate(:,874) = 2.7e-12_r8 * exp_fac(:)
791 58907904 : rate(:,396) = 2.3e-12_r8 * exp( -170._r8 * itemp(:) )
792 58907904 : rate(:,404) = 2.7e-12_r8 * exp( 580._r8 * itemp(:) )
793 58907904 : rate(:,411) = 5.4e-14_r8 * exp( 870._r8 * itemp(:) )
794 58907904 : exp_fac(:) = exp( 450._r8 * itemp(:) )
795 58907904 : rate(:,414) = 1.17e-11_r8 * exp_fac(:)
796 58907904 : rate(:,415) = 1.17e-11_r8 * exp_fac(:)
797 58907904 : exp_fac(:) = exp( 390._r8 * itemp(:) )
798 58907904 : rate(:,416) = 2.2e-11_r8 * exp_fac(:)
799 58907904 : rate(:,417) = 3.5e-11_r8 * exp_fac(:)
800 58907904 : rate(:,487) = 2.7e-11_r8 * exp_fac(:)
801 58907904 : rate(:,490) = 2.08e-11_r8 * exp_fac(:)
802 58907904 : rate(:,768) = 2.7e-11_r8 * exp_fac(:)
803 58907904 : rate(:,863) = 2.7e-11_r8 * exp_fac(:)
804 58907904 : exp_fac(:) = exp( 410._r8 * itemp(:) )
805 58907904 : rate(:,419) = 9.85e-12_r8 * exp_fac(:)
806 58907904 : rate(:,603) = 1.34e-11_r8 * exp_fac(:)
807 58907904 : rate(:,845) = 1.34e-11_r8 * exp_fac(:)
808 58907904 : exp_fac(:) = exp( -400._r8 * itemp(:) )
809 58907904 : rate(:,420) = 4.43e-11_r8 * exp_fac(:)
810 58907904 : rate(:,422) = 4.43e-11_r8 * exp_fac(:)
811 58907904 : rate(:,423) = 3.22e-11_r8 * exp_fac(:)
812 58907904 : rate(:,427) = 1.04e+11_r8 * exp( -9746._r8 * itemp(:) )
813 58907904 : rate(:,428) = 2.24e+15_r8 * exp( -10865._r8 * itemp(:) )
814 58907904 : rate(:,429) = 2.22e+15_r8 * exp( -10355._r8 * itemp(:) )
815 58907904 : rate(:,433) = 1.88e+11_r8 * exp( -9752._r8 * itemp(:) )
816 58907904 : rate(:,434) = 2.49e+15_r8 * exp( -11112._r8 * itemp(:) )
817 58907904 : rate(:,435) = 2.49e+15_r8 * exp( -10890._r8 * itemp(:) )
818 58907904 : rate(:,447) = 1.83e+14_r8 * exp( -8930._r8 * itemp(:) )
819 58907904 : rate(:,451) = 2.08e+14_r8 * exp( -9400._r8 * itemp(:) )
820 58907904 : exp_fac(:) = exp( -10000._r8 * itemp(:) )
821 58907904 : rate(:,458) = 1.256e+13_r8 * exp_fac(:)
822 58907904 : rate(:,462) = 1.875e+13_r8 * exp_fac(:)
823 58907904 : rate(:,465) = 1.875e+13_r8 * exp_fac(:)
824 58907904 : rate(:,468) = 5.092e+12_r8 * exp_fac(:)
825 58907904 : rate(:,480) = 8.72e+12_r8 * exp_fac(:)
826 58907904 : rate(:,483) = 6.55e+12_r8 * exp_fac(:)
827 58907904 : exp_fac(:) = exp( -450._r8 * itemp(:) )
828 58907904 : rate(:,473) = 2.95e-12_r8 * exp_fac(:)
829 58907904 : rate(:,764) = 2.95e-12_r8 * exp_fac(:)
830 58907904 : rate(:,859) = 2.95e-12_r8 * exp_fac(:)
831 58907904 : exp_fac(:) = exp( -1995._r8 * itemp(:) )
832 58907904 : rate(:,486) = 1.03e-14_r8 * exp_fac(:)
833 58907904 : rate(:,862) = 1.03e-14_r8 * exp_fac(:)
834 58907904 : rate(:,494) = 1.79e+14_r8 * exp( -8830._r8 * itemp(:) )
835 58907904 : rate(:,498) = 1.75e+14_r8 * exp( -9054._r8 * itemp(:) )
836 58907904 : rate(:,500) = 1e+07_r8 * exp( -5000._r8 * itemp(:) )
837 58907904 : exp_fac(:) = exp( -193._r8 * itemp(:) )
838 58907904 : rate(:,543) = 2.3e-12_r8 * exp_fac(:)
839 58907904 : rate(:,851) = 2.3e-12_r8 * exp_fac(:)
840 58907904 : rate(:,547) = 5.9e-12_r8 * exp( 225._r8 * itemp(:) )
841 58907904 : rate(:,566) = 4.7e-13_r8 * exp( 1220._r8 * itemp(:) )
842 58907904 : exp_fac(:) = exp( 352._r8 * itemp(:) )
843 58907904 : rate(:,574) = 1.7e-12_r8 * exp_fac(:)
844 58907904 : rate(:,878) = 1.7e-12_r8 * exp_fac(:)
845 58907904 : exp_fac(:) = exp( 490._r8 * itemp(:) )
846 58907904 : rate(:,584) = 1.2e-12_r8 * exp_fac(:)
847 58907904 : rate(:,841) = 1.2e-12_r8 * exp_fac(:)
848 58907904 : exp_fac(:) = exp( -640._r8 * itemp(:) )
849 58907904 : rate(:,602) = 8.05e-16_r8 * exp_fac(:)
850 58907904 : rate(:,844) = 8.05e-16_r8 * exp_fac(:)
851 58907904 : exp_fac(:) = exp( -770._r8 * itemp(:) )
852 58907904 : rate(:,662) = 2.8e-15_r8 * exp_fac(:)
853 58907904 : rate(:,870) = 2.8e-15_r8 * exp_fac(:)
854 58907904 : exp_fac(:) = exp( 470._r8 * itemp(:) )
855 58907904 : rate(:,663) = 3.41e-11_r8 * exp_fac(:)
856 58907904 : rate(:,871) = 3.41e-11_r8 * exp_fac(:)
857 58907904 : exp_fac(:) = exp( -520._r8 * itemp(:) )
858 58907904 : rate(:,682) = 2.65e-15_r8 * exp_fac(:)
859 58907904 : rate(:,875) = 2.65e-15_r8 * exp_fac(:)
860 58907904 : exp_fac(:) = exp( 600._r8 * itemp(:) )
861 58907904 : rate(:,717) = 5.2e-12_r8 * exp_fac(:)
862 58907904 : rate(:,736) = 5.2e-12_r8 * exp_fac(:)
863 58907904 : rate(:,756) = 5.2e-12_r8 * exp_fac(:)
864 58907904 : rate(:,793) = 1.9e-13_r8 * exp( 520._r8 * itemp(:) )
865 58907904 : rate(:,794) = 1.1e-11_r8 * exp( -280._r8 * itemp(:) )
866 58907904 : rate(:,795) = 2.1e-11_r8 * exp( -2200._r8 * itemp(:) )
867 58907904 : rate(:,796) = 7.2e-14_r8 * exp( -1070._r8 * itemp(:) )
868 58907904 : rate(:,804) = 1.6e-13_r8 * exp( -2280._r8 * itemp(:) )
869 58907904 : rate(:,807) = 2.6e-11_r8 * exp( 330._r8 * itemp(:) )
870 58907904 : rate(:,810) = 1.7e-12_r8 * exp( -710._r8 * itemp(:) )
871 :
872 58907904 : itemp(:) = 300._r8 * itemp(:)
873 :
874 72960 : n = ncol*pver
875 :
876 58907904 : ko(:) = 5.3e-32_r8 * itemp(:)**1.8_r8
877 58907904 : kinf(:) = 9.5e-11_r8 * itemp(:)**(-0.4_r8)
878 72960 : call jpl( rate(:,181), m, 0.6_r8, ko, kinf, n )
879 :
880 58907904 : ko(:) = 6.9e-31_r8 * itemp(:)**1._r8
881 58907904 : kinf(:) = 2.6e-11_r8
882 72960 : call jpl( rate(:,191), m, 0.6_r8, ko, kinf, n )
883 :
884 58907904 : ko(:) = 2.5e-31_r8 * itemp(:)**1.8_r8
885 58907904 : kinf(:) = 2.2e-11_r8 * itemp(:)**0.7_r8
886 72960 : call jpl( rate(:,201), m, 0.6_r8, ko, kinf, n )
887 :
888 58907904 : ko(:) = 9e-32_r8 * itemp(:)**1.5_r8
889 58907904 : kinf(:) = 3e-11_r8
890 72960 : call jpl( rate(:,209), m, 0.6_r8, ko, kinf, n )
891 :
892 58907904 : ko(:) = 1.9e-31_r8 * itemp(:)**3.4_r8
893 58907904 : kinf(:) = 4e-12_r8 * itemp(:)**0.3_r8
894 72960 : call jpl( rate(:,212), m, 0.6_r8, ko, kinf, n )
895 :
896 58907904 : ko(:) = 2.4e-30_r8 * itemp(:)**3._r8
897 58907904 : kinf(:) = 1.6e-12_r8 * itemp(:)**(-0.1_r8)
898 72960 : call jpl( rate(:,213), m, 0.6_r8, ko, kinf, n )
899 :
900 58907904 : ko(:) = 1.8e-30_r8 * itemp(:)**3._r8
901 58907904 : kinf(:) = 2.8e-11_r8
902 72960 : call jpl( rate(:,214), m, 0.6_r8, ko, kinf, n )
903 :
904 58907904 : ko(:) = 1.8e-31_r8 * itemp(:)**3.4_r8
905 58907904 : kinf(:) = 1.5e-11_r8 * itemp(:)**1.9_r8
906 72960 : call jpl( rate(:,232), m, 0.6_r8, ko, kinf, n )
907 :
908 58907904 : ko(:) = 1.9e-32_r8 * itemp(:)**3.6_r8
909 58907904 : kinf(:) = 3.7e-12_r8 * itemp(:)**1.6_r8
910 72960 : call jpl( rate(:,252), m, 0.6_r8, ko, kinf, n )
911 :
912 58907904 : ko(:) = 5.2e-31_r8 * itemp(:)**3.2_r8
913 58907904 : kinf(:) = 6.9e-12_r8 * itemp(:)**2.9_r8
914 72960 : call jpl( rate(:,263), m, 0.6_r8, ko, kinf, n )
915 :
916 58907904 : ko(:) = 6.1e-33_r8 * itemp(:)**1.5_r8
917 58907904 : kinf(:) = 9.8e-15_r8 * itemp(:)**(-4.6_r8)
918 72960 : call jpl( rate(:,309), m, 0.8_r8, ko, kinf, n )
919 :
920 58907904 : ko(:) = 5.2e-30_r8 * itemp(:)**2.4_r8
921 58907904 : kinf(:) = 2.2e-10_r8 * itemp(:)**0.7_r8
922 72960 : call jpl( rate(:,320), m, 0.6_r8, ko, kinf, n )
923 :
924 58907904 : ko(:) = 5.5e-30_r8
925 58907904 : kinf(:) = 8.3e-13_r8 * itemp(:)**(-2._r8)
926 72960 : call jpl( rate(:,321), m, 0.6_r8, ko, kinf, n )
927 :
928 58907904 : ko(:) = 1.6e-29_r8 * itemp(:)**3.3_r8
929 58907904 : kinf(:) = 3.1e-10_r8 * itemp(:)
930 72960 : call jpl( rate(:,322), m, 0.6_r8, ko, kinf, n )
931 :
932 58907904 : ko(:) = 8.6e-29_r8 * itemp(:)**3.1_r8
933 58907904 : kinf(:) = 9e-12_r8 * itemp(:)**0.85_r8
934 72960 : call jpl( rate(:,350), m, 0.48_r8, ko, kinf, n )
935 :
936 58907904 : ko(:) = 7.3e-29_r8 * itemp(:)**4.1_r8
937 58907904 : kinf(:) = 9.5e-12_r8 * itemp(:)**1.6_r8
938 72960 : call jpl( rate(:,351), m, 0.6_r8, ko, kinf, n )
939 :
940 58907904 : ko(:) = 8e-27_r8 * itemp(:)**3.5_r8
941 58907904 : kinf(:) = 3e-11_r8
942 72960 : call jpl( rate(:,372), m, 0.5_r8, ko, kinf, n )
943 :
944 58907904 : ko(:) = 8e-27_r8 * itemp(:)**3.5_r8
945 58907904 : kinf(:) = 3e-11_r8
946 72960 : call jpl( rate(:,398), m, 0.5_r8, ko, kinf, n )
947 :
948 58907904 : ko(:) = 9.7e-29_r8 * itemp(:)**5.6_r8
949 58907904 : kinf(:) = 9.3e-12_r8 * itemp(:)**1.5_r8
950 72960 : call jpl( rate(:,406), m, 0.6_r8, ko, kinf, n )
951 :
952 58907904 : ko(:) = 9.7e-29_r8 * itemp(:)**5.6_r8
953 58907904 : kinf(:) = 9.3e-12_r8 * itemp(:)**1.5_r8
954 72960 : call jpl( rate(:,557), m, 0.6_r8, ko, kinf, n )
955 :
956 58907904 : ko(:) = 9.7e-29_r8 * itemp(:)**5.6_r8
957 58907904 : kinf(:) = 9.3e-12_r8 * itemp(:)**1.5_r8
958 72960 : call jpl( rate(:,560), m, 0.6_r8, ko, kinf, n )
959 :
960 58907904 : ko(:) = 9.7e-29_r8 * itemp(:)**5.6_r8
961 58907904 : kinf(:) = 9.3e-12_r8 * itemp(:)**1.5_r8
962 72960 : call jpl( rate(:,563), m, 0.6_r8, ko, kinf, n )
963 :
964 58907904 : ko(:) = 9.7e-29_r8 * itemp(:)**5.6_r8
965 58907904 : kinf(:) = 9.3e-12_r8 * itemp(:)**1.5_r8
966 72960 : call jpl( rate(:,570), m, 0.6_r8, ko, kinf, n )
967 :
968 58907904 : ko(:) = 9.7e-29_r8 * itemp(:)**5.6_r8
969 58907904 : kinf(:) = 9.3e-12_r8 * itemp(:)**1.5_r8
970 72960 : call jpl( rate(:,684), m, 0.6_r8, ko, kinf, n )
971 :
972 58907904 : ko(:) = 9.7e-29_r8 * itemp(:)**5.6_r8
973 58907904 : kinf(:) = 9.3e-12_r8 * itemp(:)**1.5_r8
974 72960 : call jpl( rate(:,685), m, 0.6_r8, ko, kinf, n )
975 :
976 58907904 : ko(:) = 9.7e-29_r8 * itemp(:)**5.6_r8
977 58907904 : kinf(:) = 9.3e-12_r8 * itemp(:)**1.5_r8
978 72960 : call jpl( rate(:,686), m, 0.6_r8, ko, kinf, n )
979 :
980 58907904 : ko(:) = 2.9e-31_r8 * itemp(:)**4.1_r8
981 58907904 : kinf(:) = 1.7e-12_r8 * itemp(:)**(-0.2_r8)
982 72960 : call jpl( rate(:,798), m, 0.6_r8, ko, kinf, n )
983 :
984 72960 : end subroutine setrxt
985 :
986 :
987 0 : subroutine setrxt_hrates( rate, temp, m, ncol, kbot )
988 :
989 : use ppgrid, only : pcols, pver
990 :
991 :
992 : use chem_mods, only : rxntot
993 : use mo_jpl, only : jpl
994 :
995 : implicit none
996 :
997 : !-------------------------------------------------------
998 : ! ... dummy arguments
999 : !-------------------------------------------------------
1000 : integer, intent(in) :: ncol
1001 : integer, intent(in) :: kbot
1002 : real(r8), intent(in) :: temp(pcols,pver)
1003 : real(r8), intent(in) :: m(ncol*pver)
1004 : real(r8), intent(inout) :: rate(ncol*pver,max(1,rxntot))
1005 :
1006 : !-------------------------------------------------------
1007 : ! ... local variables
1008 : !-------------------------------------------------------
1009 : integer :: n
1010 : integer :: offset
1011 : integer :: k
1012 0 : real(r8) :: itemp(ncol*kbot)
1013 : real(r8) :: exp_fac(ncol*kbot)
1014 0 : real(r8) :: ko(ncol*kbot)
1015 0 : real(r8) :: kinf(ncol*kbot)
1016 0 : real(r8) :: wrk(ncol*kbot)
1017 :
1018 0 : n = ncol*kbot
1019 :
1020 0 : rate(:n,178) = 6.9e-12_r8
1021 :
1022 0 : do k = 1,kbot
1023 0 : offset = (k-1)*ncol
1024 0 : itemp(offset+1:offset+ncol) = 1._r8 / temp(:ncol,k)
1025 : end do
1026 :
1027 0 : rate(:n,170) = 2.15e-11_r8 * exp( 110._r8 * itemp(:) )
1028 0 : rate(:n,173) = 8e-12_r8 * exp( -2060._r8 * itemp(:) )
1029 0 : rate(:n,182) = 3e-11_r8 * exp( 200._r8 * itemp(:) )
1030 0 : rate(:n,183) = 1e-14_r8 * exp( -490._r8 * itemp(:) )
1031 0 : rate(:n,184) = 1.4e-10_r8 * exp( -470._r8 * itemp(:) )
1032 0 : rate(:n,187) = 4.8e-11_r8 * exp( 250._r8 * itemp(:) )
1033 0 : rate(:n,188) = 1.8e-11_r8 * exp( 180._r8 * itemp(:) )
1034 0 : rate(:n,189) = 1.7e-12_r8 * exp( -940._r8 * itemp(:) )
1035 0 : rate(:n,194) = 2.1e-11_r8 * exp( 100._r8 * itemp(:) )
1036 0 : rate(:n,198) = 3.3e-12_r8 * exp( -3150._r8 * itemp(:) )
1037 0 : rate(:n,199) = 5.1e-12_r8 * exp( 210._r8 * itemp(:) )
1038 0 : rate(:n,207) = 3.44e-12_r8 * exp( 260._r8 * itemp(:) )
1039 0 : rate(:n,208) = 3e-12_r8 * exp( -1500._r8 * itemp(:) )
1040 :
1041 0 : itemp(:) = 300._r8 * itemp(:)
1042 :
1043 0 : ko(:) = 5.3e-32_r8 * itemp(:)**1.8_r8
1044 0 : kinf(:) = 9.5e-11_r8 * itemp(:)**(-0.4_r8)
1045 0 : call jpl( wrk, m, 0.6_r8, ko, kinf, n )
1046 0 : rate(:n,181) = wrk(:)
1047 :
1048 :
1049 :
1050 :
1051 :
1052 :
1053 :
1054 :
1055 :
1056 :
1057 :
1058 :
1059 :
1060 :
1061 :
1062 :
1063 :
1064 :
1065 :
1066 :
1067 :
1068 :
1069 :
1070 :
1071 :
1072 :
1073 :
1074 0 : end subroutine setrxt_hrates
1075 :
1076 : end module mo_setrxt
|