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 33692928 : rate(:,144) = 1.2e-10_r8
41 33692928 : rate(:,148) = 1.2e-10_r8
42 33692928 : rate(:,154) = 6.9e-12_r8
43 33692928 : rate(:,155) = 7.2e-11_r8
44 33692928 : rate(:,156) = 1.6e-12_r8
45 33692928 : rate(:,162) = 1.8e-12_r8
46 33692928 : rate(:,166) = 1.8e-12_r8
47 33692928 : rate(:,178) = 3.5e-12_r8
48 33692928 : rate(:,180) = 1.3e-11_r8
49 33692928 : rate(:,181) = 2.2e-11_r8
50 33692928 : rate(:,182) = 5e-11_r8
51 33692928 : rate(:,217) = 1.7e-13_r8
52 33692928 : rate(:,219) = 2.607e-10_r8
53 33692928 : rate(:,220) = 9.75e-11_r8
54 33692928 : rate(:,221) = 2.07e-10_r8
55 33692928 : rate(:,222) = 2.088e-10_r8
56 33692928 : rate(:,223) = 1.17e-10_r8
57 33692928 : rate(:,224) = 4.644e-11_r8
58 33692928 : rate(:,225) = 1.204e-10_r8
59 33692928 : rate(:,226) = 9.9e-11_r8
60 33692928 : rate(:,227) = 3.3e-12_r8
61 33692928 : rate(:,246) = 4.5e-11_r8
62 33692928 : rate(:,247) = 4.62e-10_r8
63 33692928 : rate(:,248) = 1.2e-10_r8
64 33692928 : rate(:,249) = 9e-11_r8
65 33692928 : rate(:,250) = 3e-11_r8
66 33692928 : rate(:,255) = 2.14e-11_r8
67 33692928 : rate(:,256) = 1.9e-10_r8
68 33692928 : rate(:,269) = 2.57e-10_r8
69 33692928 : rate(:,270) = 1.8e-10_r8
70 33692928 : rate(:,271) = 1.794e-10_r8
71 33692928 : rate(:,272) = 1.3e-10_r8
72 33692928 : rate(:,273) = 7.65e-11_r8
73 33692928 : rate(:,286) = 4e-13_r8
74 33692928 : rate(:,290) = 1.31e-10_r8
75 33692928 : rate(:,291) = 3.5e-11_r8
76 33692928 : rate(:,292) = 9e-12_r8
77 33692928 : rate(:,299) = 6.8e-14_r8
78 33692928 : rate(:,300) = 2e-13_r8
79 33692928 : rate(:,315) = 1e-12_r8
80 33692928 : rate(:,319) = 1e-14_r8
81 33692928 : rate(:,320) = 1e-11_r8
82 33692928 : rate(:,321) = 1.15e-11_r8
83 33692928 : rate(:,322) = 4e-14_r8
84 33692928 : rate(:,335) = 3e-12_r8
85 33692928 : rate(:,336) = 6.7e-13_r8
86 33692928 : rate(:,346) = 3.5e-13_r8
87 33692928 : rate(:,347) = 5.4e-11_r8
88 33692928 : rate(:,350) = 2e-12_r8
89 33692928 : rate(:,351) = 1.4e-11_r8
90 33692928 : rate(:,354) = 2.4e-12_r8
91 33692928 : rate(:,365) = 5e-12_r8
92 33692928 : rate(:,375) = 1.6e-12_r8
93 33692928 : rate(:,377) = 6.7e-12_r8
94 33692928 : rate(:,380) = 3.5e-12_r8
95 33692928 : rate(:,383) = 1.3e-11_r8
96 33692928 : rate(:,384) = 1.4e-11_r8
97 33692928 : rate(:,388) = 2.4e-12_r8
98 33692928 : rate(:,389) = 1.4e-11_r8
99 33692928 : rate(:,394) = 2.4e-12_r8
100 33692928 : rate(:,395) = 4e-11_r8
101 33692928 : rate(:,396) = 4e-11_r8
102 33692928 : rate(:,398) = 1.4e-11_r8
103 33692928 : rate(:,402) = 2.4e-12_r8
104 33692928 : rate(:,403) = 4e-11_r8
105 33692928 : rate(:,407) = 7e-11_r8
106 33692928 : rate(:,408) = 1e-10_r8
107 33692928 : rate(:,413) = 2.4e-12_r8
108 33692928 : rate(:,428) = 4.7e-11_r8
109 33692928 : rate(:,441) = 2.1e-12_r8
110 33692928 : rate(:,442) = 2.8e-13_r8
111 33692928 : rate(:,450) = 1.7e-11_r8
112 33692928 : rate(:,456) = 8.4e-11_r8
113 33692928 : rate(:,458) = 1.9e-11_r8
114 33692928 : rate(:,459) = 1.2e-14_r8
115 33692928 : rate(:,460) = 2e-10_r8
116 33692928 : rate(:,467) = 2.4e-12_r8
117 33692928 : rate(:,468) = 2e-11_r8
118 33692928 : rate(:,472) = 2.3e-11_r8
119 33692928 : rate(:,473) = 2e-11_r8
120 33692928 : rate(:,477) = 3.3e-11_r8
121 33692928 : rate(:,478) = 1e-12_r8
122 33692928 : rate(:,479) = 5.7e-11_r8
123 33692928 : rate(:,480) = 3.4e-11_r8
124 33692928 : rate(:,485) = 2.3e-12_r8
125 33692928 : rate(:,487) = 1.2e-11_r8
126 33692928 : rate(:,488) = 5.7e-11_r8
127 33692928 : rate(:,489) = 2.8e-11_r8
128 33692928 : rate(:,490) = 6.6e-11_r8
129 33692928 : rate(:,491) = 1.4e-11_r8
130 33692928 : rate(:,494) = 1.9e-12_r8
131 33692928 : rate(:,506) = 6.34e-08_r8
132 33692928 : rate(:,512) = 1.9e-11_r8
133 33692928 : rate(:,515) = 1.2e-14_r8
134 33692928 : rate(:,516) = 2e-10_r8
135 33692928 : rate(:,527) = 1.34e-11_r8
136 33692928 : rate(:,530) = 1.34e-11_r8
137 33692928 : rate(:,536) = 1.34e-11_r8
138 33692928 : rate(:,537) = 1.34e-11_r8
139 33692928 : rate(:,542) = 1.7e-11_r8
140 33692928 : rate(:,562) = 1.29e-07_r8
141 :
142 2407680 : do n = 1,pver
143 2334720 : offset = (n-1)*ncol
144 36027648 : itemp(offset+1:offset+ncol) = 1._r8 / temp(:ncol,n)
145 : end do
146 :
147 33692928 : rate(:,145) = 1.63e-10_r8 * exp( 60._r8 * itemp(:) )
148 33692928 : rate(:,146) = 2.15e-11_r8 * exp( 110._r8 * itemp(:) )
149 33692928 : rate(:,147) = 3.3e-11_r8 * exp( 55._r8 * itemp(:) )
150 33692928 : rate(:,149) = 8e-12_r8 * exp( -2060._r8 * itemp(:) )
151 33692928 : rate(:,152) = 1.6e-11_r8 * exp( -4570._r8 * itemp(:) )
152 33692928 : exp_fac(:) = exp( -2000._r8 * itemp(:) )
153 33692928 : rate(:,153) = 1.4e-12_r8 * exp_fac(:)
154 33692928 : rate(:,404) = 1.05e-14_r8 * exp_fac(:)
155 33692928 : rate(:,523) = 1.05e-14_r8 * exp_fac(:)
156 33692928 : exp_fac(:) = exp( 200._r8 * itemp(:) )
157 33692928 : rate(:,158) = 3e-11_r8 * exp_fac(:)
158 33692928 : rate(:,244) = 5.5e-12_r8 * exp_fac(:)
159 33692928 : rate(:,283) = 3.8e-12_r8 * exp_fac(:)
160 33692928 : rate(:,304) = 3.8e-12_r8 * exp_fac(:)
161 33692928 : rate(:,331) = 3.8e-12_r8 * exp_fac(:)
162 33692928 : rate(:,339) = 3.8e-12_r8 * exp_fac(:)
163 33692928 : rate(:,343) = 3.8e-12_r8 * exp_fac(:)
164 33692928 : rate(:,359) = 2.3e-11_r8 * exp_fac(:)
165 33692928 : rate(:,369) = 3.8e-12_r8 * exp_fac(:)
166 33692928 : rate(:,379) = 3.8e-12_r8 * exp_fac(:)
167 33692928 : rate(:,406) = 1.52e-11_r8 * exp_fac(:)
168 33692928 : rate(:,414) = 1.52e-12_r8 * exp_fac(:)
169 33692928 : rate(:,420) = 3.8e-12_r8 * exp_fac(:)
170 33692928 : rate(:,423) = 3.8e-12_r8 * exp_fac(:)
171 33692928 : rate(:,427) = 3.8e-12_r8 * exp_fac(:)
172 33692928 : rate(:,443) = 3.8e-12_r8 * exp_fac(:)
173 33692928 : rate(:,447) = 3.8e-12_r8 * exp_fac(:)
174 33692928 : rate(:,453) = 3.8e-12_r8 * exp_fac(:)
175 33692928 : rate(:,457) = 3.8e-12_r8 * exp_fac(:)
176 33692928 : rate(:,159) = 1e-14_r8 * exp( -490._r8 * itemp(:) )
177 33692928 : rate(:,160) = 1.4e-10_r8 * exp( -470._r8 * itemp(:) )
178 33692928 : rate(:,161) = 2.8e-12_r8 * exp( -1800._r8 * itemp(:) )
179 33692928 : exp_fac(:) = exp( 250._r8 * itemp(:) )
180 33692928 : rate(:,163) = 4.8e-11_r8 * exp_fac(:)
181 33692928 : rate(:,242) = 1.7e-11_r8 * exp_fac(:)
182 33692928 : exp_fac(:) = exp( 180._r8 * itemp(:) )
183 33692928 : rate(:,164) = 1.8e-11_r8 * exp_fac(:)
184 33692928 : rate(:,317) = 4.2e-12_r8 * exp_fac(:)
185 33692928 : rate(:,330) = 4.2e-12_r8 * exp_fac(:)
186 33692928 : rate(:,338) = 4.2e-12_r8 * exp_fac(:)
187 33692928 : rate(:,367) = 4.2e-12_r8 * exp_fac(:)
188 33692928 : rate(:,387) = 4.4e-12_r8 * exp_fac(:)
189 33692928 : rate(:,393) = 4.4e-12_r8 * exp_fac(:)
190 33692928 : rate(:,466) = 4.2e-12_r8 * exp_fac(:)
191 33692928 : rate(:,471) = 4.2e-12_r8 * exp_fac(:)
192 33692928 : rate(:,476) = 4.2e-12_r8 * exp_fac(:)
193 33692928 : rate(:,165) = 1.7e-12_r8 * exp( -940._r8 * itemp(:) )
194 33692928 : rate(:,169) = 4.5e-13_r8 * exp( 610._r8 * itemp(:) )
195 33692928 : rate(:,170) = 2.1e-11_r8 * exp( 100._r8 * itemp(:) )
196 33692928 : exp_fac(:) = exp( 220._r8 * itemp(:) )
197 33692928 : rate(:,171) = 2.9e-12_r8 * exp_fac(:)
198 33692928 : rate(:,172) = 1.45e-12_r8 * exp_fac(:)
199 33692928 : rate(:,173) = 1.45e-12_r8 * exp_fac(:)
200 33692928 : rate(:,174) = 3.3e-12_r8 * exp( -3150._r8 * itemp(:) )
201 33692928 : rate(:,175) = 5.1e-12_r8 * exp( 210._r8 * itemp(:) )
202 33692928 : exp_fac(:) = exp( -2450._r8 * itemp(:) )
203 33692928 : rate(:,176) = 1.2e-13_r8 * exp_fac(:)
204 33692928 : rate(:,202) = 3e-11_r8 * exp_fac(:)
205 33692928 : exp_fac(:) = exp( 125._r8 * itemp(:) )
206 33692928 : rate(:,179) = 1.7e-11_r8 * exp_fac(:)
207 33692928 : rate(:,277) = 5.5e-12_r8 * exp_fac(:)
208 33692928 : exp_fac(:) = exp( 260._r8 * itemp(:) )
209 33692928 : rate(:,183) = 3.44e-12_r8 * exp_fac(:)
210 33692928 : rate(:,235) = 2.3e-12_r8 * exp_fac(:)
211 33692928 : rate(:,238) = 8.8e-12_r8 * exp_fac(:)
212 33692928 : exp_fac(:) = exp( -1500._r8 * itemp(:) )
213 33692928 : rate(:,184) = 3e-12_r8 * exp_fac(:)
214 33692928 : rate(:,243) = 5.8e-12_r8 * exp_fac(:)
215 33692928 : exp_fac(:) = exp( 20._r8 * itemp(:) )
216 33692928 : rate(:,186) = 7.26e-11_r8 * exp_fac(:)
217 33692928 : rate(:,187) = 4.64e-11_r8 * exp_fac(:)
218 33692928 : rate(:,194) = 8.1e-11_r8 * exp( -30._r8 * itemp(:) )
219 33692928 : rate(:,195) = 7.1e-12_r8 * exp( -1270._r8 * itemp(:) )
220 33692928 : rate(:,196) = 3.05e-11_r8 * exp( -2270._r8 * itemp(:) )
221 33692928 : rate(:,197) = 1.1e-11_r8 * exp( -980._r8 * itemp(:) )
222 33692928 : exp_fac(:) = exp( 270._r8 * itemp(:) )
223 33692928 : rate(:,198) = 1.4e-11_r8 * exp_fac(:)
224 33692928 : rate(:,212) = 7.4e-12_r8 * exp_fac(:)
225 33692928 : rate(:,313) = 8.1e-12_r8 * exp_fac(:)
226 33692928 : rate(:,199) = 3.6e-11_r8 * exp( -375._r8 * itemp(:) )
227 33692928 : rate(:,200) = 2.3e-11_r8 * exp( -200._r8 * itemp(:) )
228 33692928 : rate(:,201) = 3.3e-12_r8 * exp( -115._r8 * itemp(:) )
229 33692928 : rate(:,203) = 1e-12_r8 * exp( -1590._r8 * itemp(:) )
230 33692928 : rate(:,204) = 3.5e-13_r8 * exp( -1370._r8 * itemp(:) )
231 33692928 : exp_fac(:) = exp( 290._r8 * itemp(:) )
232 33692928 : rate(:,205) = 2.6e-12_r8 * exp_fac(:)
233 33692928 : rate(:,206) = 6.4e-12_r8 * exp_fac(:)
234 33692928 : rate(:,236) = 4.1e-13_r8 * exp_fac(:)
235 33692928 : rate(:,416) = 7.5e-12_r8 * exp_fac(:)
236 33692928 : rate(:,430) = 7.5e-12_r8 * exp_fac(:)
237 33692928 : rate(:,433) = 7.5e-12_r8 * exp_fac(:)
238 33692928 : rate(:,436) = 7.5e-12_r8 * exp_fac(:)
239 33692928 : rate(:,207) = 6.5e-12_r8 * exp( 135._r8 * itemp(:) )
240 33692928 : exp_fac(:) = exp( -840._r8 * itemp(:) )
241 33692928 : rate(:,209) = 3.6e-12_r8 * exp_fac(:)
242 33692928 : rate(:,258) = 2e-12_r8 * exp_fac(:)
243 33692928 : rate(:,210) = 1.2e-12_r8 * exp( -330._r8 * itemp(:) )
244 33692928 : rate(:,211) = 2.8e-11_r8 * exp( 85._r8 * itemp(:) )
245 33692928 : exp_fac(:) = exp( 230._r8 * itemp(:) )
246 33692928 : rate(:,213) = 6e-13_r8 * exp_fac(:)
247 33692928 : rate(:,233) = 1.5e-12_r8 * exp_fac(:)
248 33692928 : rate(:,241) = 1.9e-11_r8 * exp_fac(:)
249 33692928 : rate(:,214) = 1e-11_r8 * exp( -3300._r8 * itemp(:) )
250 33692928 : rate(:,215) = 1.8e-12_r8 * exp( -250._r8 * itemp(:) )
251 33692928 : rate(:,216) = 3.4e-12_r8 * exp( -130._r8 * itemp(:) )
252 33692928 : exp_fac(:) = exp( -500._r8 * itemp(:) )
253 33692928 : rate(:,218) = 3e-12_r8 * exp_fac(:)
254 33692928 : rate(:,252) = 1.4e-10_r8 * exp_fac(:)
255 33692928 : exp_fac(:) = exp( -800._r8 * itemp(:) )
256 33692928 : rate(:,230) = 1.7e-11_r8 * exp_fac(:)
257 33692928 : rate(:,257) = 6.3e-12_r8 * exp_fac(:)
258 33692928 : rate(:,231) = 4.8e-12_r8 * exp( -310._r8 * itemp(:) )
259 33692928 : rate(:,232) = 1.6e-11_r8 * exp( -780._r8 * itemp(:) )
260 33692928 : rate(:,234) = 9.5e-13_r8 * exp( 550._r8 * itemp(:) )
261 33692928 : rate(:,237) = 4.5e-12_r8 * exp( 460._r8 * itemp(:) )
262 33692928 : rate(:,240) = 1.9e-11_r8 * exp( 215._r8 * itemp(:) )
263 33692928 : rate(:,245) = 1.2e-10_r8 * exp( -430._r8 * itemp(:) )
264 33692928 : rate(:,251) = 1.6e-10_r8 * exp( -260._r8 * itemp(:) )
265 33692928 : exp_fac(:) = exp( 0._r8 * itemp(:) )
266 33692928 : rate(:,253) = 1.4e-11_r8 * exp_fac(:)
267 33692928 : rate(:,255) = 2.14e-11_r8 * exp_fac(:)
268 33692928 : rate(:,256) = 1.9e-10_r8 * exp_fac(:)
269 33692928 : rate(:,269) = 2.57e-10_r8 * exp_fac(:)
270 33692928 : rate(:,270) = 1.8e-10_r8 * exp_fac(:)
271 33692928 : rate(:,271) = 1.794e-10_r8 * exp_fac(:)
272 33692928 : rate(:,272) = 1.3e-10_r8 * exp_fac(:)
273 33692928 : rate(:,273) = 7.65e-11_r8 * exp_fac(:)
274 33692928 : rate(:,286) = 4e-13_r8 * exp_fac(:)
275 33692928 : rate(:,290) = 1.31e-10_r8 * exp_fac(:)
276 33692928 : rate(:,291) = 3.5e-11_r8 * exp_fac(:)
277 33692928 : rate(:,292) = 9e-12_r8 * exp_fac(:)
278 33692928 : rate(:,299) = 6.8e-14_r8 * exp_fac(:)
279 33692928 : rate(:,300) = 2e-13_r8 * exp_fac(:)
280 33692928 : rate(:,315) = 1e-12_r8 * exp_fac(:)
281 33692928 : rate(:,319) = 1e-14_r8 * exp_fac(:)
282 33692928 : rate(:,320) = 1e-11_r8 * exp_fac(:)
283 33692928 : rate(:,321) = 1.15e-11_r8 * exp_fac(:)
284 33692928 : rate(:,322) = 4e-14_r8 * exp_fac(:)
285 33692928 : rate(:,335) = 3e-12_r8 * exp_fac(:)
286 33692928 : rate(:,336) = 6.7e-13_r8 * exp_fac(:)
287 33692928 : rate(:,346) = 3.5e-13_r8 * exp_fac(:)
288 33692928 : rate(:,347) = 5.4e-11_r8 * exp_fac(:)
289 33692928 : rate(:,350) = 2e-12_r8 * exp_fac(:)
290 33692928 : rate(:,351) = 1.4e-11_r8 * exp_fac(:)
291 33692928 : rate(:,354) = 2.4e-12_r8 * exp_fac(:)
292 33692928 : rate(:,365) = 5e-12_r8 * exp_fac(:)
293 33692928 : rate(:,375) = 1.6e-12_r8 * exp_fac(:)
294 33692928 : rate(:,377) = 6.7e-12_r8 * exp_fac(:)
295 33692928 : rate(:,380) = 3.5e-12_r8 * exp_fac(:)
296 33692928 : rate(:,383) = 1.3e-11_r8 * exp_fac(:)
297 33692928 : rate(:,384) = 1.4e-11_r8 * exp_fac(:)
298 33692928 : rate(:,388) = 2.4e-12_r8 * exp_fac(:)
299 33692928 : rate(:,389) = 1.4e-11_r8 * exp_fac(:)
300 33692928 : rate(:,394) = 2.4e-12_r8 * exp_fac(:)
301 33692928 : rate(:,395) = 4e-11_r8 * exp_fac(:)
302 33692928 : rate(:,396) = 4e-11_r8 * exp_fac(:)
303 33692928 : rate(:,398) = 1.4e-11_r8 * exp_fac(:)
304 33692928 : rate(:,402) = 2.4e-12_r8 * exp_fac(:)
305 33692928 : rate(:,403) = 4e-11_r8 * exp_fac(:)
306 33692928 : rate(:,407) = 7e-11_r8 * exp_fac(:)
307 33692928 : rate(:,408) = 1e-10_r8 * exp_fac(:)
308 33692928 : rate(:,413) = 2.4e-12_r8 * exp_fac(:)
309 33692928 : rate(:,428) = 4.7e-11_r8 * exp_fac(:)
310 33692928 : rate(:,441) = 2.1e-12_r8 * exp_fac(:)
311 33692928 : rate(:,442) = 2.8e-13_r8 * exp_fac(:)
312 33692928 : rate(:,450) = 1.7e-11_r8 * exp_fac(:)
313 33692928 : rate(:,456) = 8.4e-11_r8 * exp_fac(:)
314 33692928 : rate(:,458) = 1.9e-11_r8 * exp_fac(:)
315 33692928 : rate(:,459) = 1.2e-14_r8 * exp_fac(:)
316 33692928 : rate(:,460) = 2e-10_r8 * exp_fac(:)
317 33692928 : rate(:,467) = 2.4e-12_r8 * exp_fac(:)
318 33692928 : rate(:,468) = 2e-11_r8 * exp_fac(:)
319 33692928 : rate(:,472) = 2.3e-11_r8 * exp_fac(:)
320 33692928 : rate(:,473) = 2e-11_r8 * exp_fac(:)
321 33692928 : rate(:,477) = 3.3e-11_r8 * exp_fac(:)
322 33692928 : rate(:,478) = 1e-12_r8 * exp_fac(:)
323 33692928 : rate(:,479) = 5.7e-11_r8 * exp_fac(:)
324 33692928 : rate(:,480) = 3.4e-11_r8 * exp_fac(:)
325 33692928 : rate(:,485) = 2.3e-12_r8 * exp_fac(:)
326 33692928 : rate(:,487) = 1.2e-11_r8 * exp_fac(:)
327 33692928 : rate(:,488) = 5.7e-11_r8 * exp_fac(:)
328 33692928 : rate(:,489) = 2.8e-11_r8 * exp_fac(:)
329 33692928 : rate(:,490) = 6.6e-11_r8 * exp_fac(:)
330 33692928 : rate(:,491) = 1.4e-11_r8 * exp_fac(:)
331 33692928 : rate(:,494) = 1.9e-12_r8 * exp_fac(:)
332 33692928 : rate(:,506) = 6.34e-08_r8 * exp_fac(:)
333 33692928 : rate(:,512) = 1.9e-11_r8 * exp_fac(:)
334 33692928 : rate(:,515) = 1.2e-14_r8 * exp_fac(:)
335 33692928 : rate(:,516) = 2e-10_r8 * exp_fac(:)
336 33692928 : rate(:,527) = 1.34e-11_r8 * exp_fac(:)
337 33692928 : rate(:,530) = 1.34e-11_r8 * exp_fac(:)
338 33692928 : rate(:,536) = 1.34e-11_r8 * exp_fac(:)
339 33692928 : rate(:,537) = 1.34e-11_r8 * exp_fac(:)
340 33692928 : rate(:,542) = 1.7e-11_r8 * exp_fac(:)
341 33692928 : rate(:,562) = 1.29e-07_r8 * exp_fac(:)
342 33692928 : exp_fac(:) = exp( 400._r8 * itemp(:) )
343 33692928 : rate(:,254) = 6e-12_r8 * exp_fac(:)
344 33692928 : rate(:,352) = 5e-13_r8 * exp_fac(:)
345 33692928 : rate(:,385) = 5e-13_r8 * exp_fac(:)
346 33692928 : rate(:,390) = 5e-13_r8 * exp_fac(:)
347 33692928 : rate(:,399) = 5e-13_r8 * exp_fac(:)
348 33692928 : rate(:,410) = 5e-13_r8 * exp_fac(:)
349 33692928 : rate(:,259) = 1.46e-11_r8 * exp( -1040._r8 * itemp(:) )
350 33692928 : rate(:,260) = 1.42e-12_r8 * exp( -1150._r8 * itemp(:) )
351 33692928 : exp_fac(:) = exp( -1520._r8 * itemp(:) )
352 33692928 : rate(:,261) = 1.64e-12_r8 * exp_fac(:)
353 33692928 : rate(:,371) = 8.5e-16_r8 * exp_fac(:)
354 33692928 : rate(:,262) = 2.03e-11_r8 * exp( -1110._r8 * itemp(:) )
355 33692928 : rate(:,263) = 1.96e-12_r8 * exp( -1200._r8 * itemp(:) )
356 33692928 : rate(:,264) = 4.85e-12_r8 * exp( -850._r8 * itemp(:) )
357 33692928 : rate(:,265) = 9e-13_r8 * exp( -360._r8 * itemp(:) )
358 33692928 : exp_fac(:) = exp( -1600._r8 * itemp(:) )
359 33692928 : rate(:,266) = 1.25e-12_r8 * exp_fac(:)
360 33692928 : rate(:,276) = 3.4e-11_r8 * exp_fac(:)
361 33692928 : rate(:,267) = 1.3e-12_r8 * exp( -1770._r8 * itemp(:) )
362 33692928 : rate(:,268) = 9.2e-13_r8 * exp( -1560._r8 * itemp(:) )
363 33692928 : rate(:,274) = 9.7e-15_r8 * exp( 625._r8 * itemp(:) )
364 33692928 : rate(:,275) = 6e-13_r8 * exp( -2058._r8 * itemp(:) )
365 33692928 : rate(:,278) = 5e-13_r8 * exp( -424._r8 * itemp(:) )
366 33692928 : rate(:,279) = 1.9e-14_r8 * exp( 706._r8 * itemp(:) )
367 33692928 : rate(:,280) = 4.1e-13_r8 * exp( 750._r8 * itemp(:) )
368 33692928 : exp_fac(:) = exp( 300._r8 * itemp(:) )
369 33692928 : rate(:,281) = 2.8e-12_r8 * exp_fac(:)
370 33692928 : rate(:,342) = 2.9e-12_r8 * exp_fac(:)
371 33692928 : rate(:,282) = 2.9e-12_r8 * exp( -345._r8 * itemp(:) )
372 33692928 : rate(:,284) = 2.45e-12_r8 * exp( -1775._r8 * itemp(:) )
373 33692928 : exp_fac(:) = exp( 700._r8 * itemp(:) )
374 33692928 : rate(:,287) = 7.5e-13_r8 * exp_fac(:)
375 33692928 : rate(:,301) = 7.5e-13_r8 * exp_fac(:)
376 33692928 : rate(:,316) = 7.5e-13_r8 * exp_fac(:)
377 33692928 : rate(:,329) = 7.5e-13_r8 * exp_fac(:)
378 33692928 : rate(:,337) = 7.5e-13_r8 * exp_fac(:)
379 33692928 : rate(:,341) = 8.6e-13_r8 * exp_fac(:)
380 33692928 : rate(:,353) = 8e-13_r8 * exp_fac(:)
381 33692928 : rate(:,366) = 7.5e-13_r8 * exp_fac(:)
382 33692928 : rate(:,376) = 7.5e-13_r8 * exp_fac(:)
383 33692928 : rate(:,386) = 8e-13_r8 * exp_fac(:)
384 33692928 : rate(:,391) = 8e-13_r8 * exp_fac(:)
385 33692928 : rate(:,400) = 8e-13_r8 * exp_fac(:)
386 33692928 : rate(:,411) = 8e-13_r8 * exp_fac(:)
387 33692928 : rate(:,418) = 7.5e-13_r8 * exp_fac(:)
388 33692928 : rate(:,422) = 7.5e-13_r8 * exp_fac(:)
389 33692928 : rate(:,425) = 7.5e-13_r8 * exp_fac(:)
390 33692928 : rate(:,438) = 7.5e-13_r8 * exp_fac(:)
391 33692928 : rate(:,445) = 7.5e-13_r8 * exp_fac(:)
392 33692928 : rate(:,451) = 7.5e-13_r8 * exp_fac(:)
393 33692928 : rate(:,454) = 7.5e-13_r8 * exp_fac(:)
394 33692928 : rate(:,465) = 7.5e-13_r8 * exp_fac(:)
395 33692928 : rate(:,470) = 7.5e-13_r8 * exp_fac(:)
396 33692928 : rate(:,475) = 7.5e-13_r8 * exp_fac(:)
397 33692928 : rate(:,518) = 7.5e-13_r8 * exp_fac(:)
398 33692928 : rate(:,525) = 7.5e-13_r8 * exp_fac(:)
399 33692928 : rate(:,528) = 7.5e-13_r8 * exp_fac(:)
400 33692928 : rate(:,539) = 7.5e-13_r8 * exp_fac(:)
401 33692928 : rate(:,543) = 7.5e-13_r8 * exp_fac(:)
402 33692928 : rate(:,288) = 2.4e+12_r8 * exp( -7000._r8 * itemp(:) )
403 33692928 : rate(:,289) = 2.6e-12_r8 * exp( 265._r8 * itemp(:) )
404 33692928 : rate(:,293) = 1.08e-10_r8 * exp( 105._r8 * itemp(:) )
405 33692928 : rate(:,298) = 1.2e-14_r8 * exp( -2630._r8 * itemp(:) )
406 33692928 : exp_fac(:) = exp( 365._r8 * itemp(:) )
407 33692928 : rate(:,302) = 2.6e-12_r8 * exp_fac(:)
408 33692928 : rate(:,419) = 2.6e-12_r8 * exp_fac(:)
409 33692928 : rate(:,424) = 2.6e-12_r8 * exp_fac(:)
410 33692928 : rate(:,426) = 2.6e-12_r8 * exp_fac(:)
411 33692928 : rate(:,439) = 2.6e-12_r8 * exp_fac(:)
412 33692928 : rate(:,446) = 2.6e-12_r8 * exp_fac(:)
413 33692928 : rate(:,452) = 2.6e-12_r8 * exp_fac(:)
414 33692928 : rate(:,455) = 2.6e-12_r8 * exp_fac(:)
415 33692928 : rate(:,519) = 2.6e-12_r8 * exp_fac(:)
416 33692928 : rate(:,526) = 2.6e-12_r8 * exp_fac(:)
417 33692928 : rate(:,529) = 2.6e-12_r8 * exp_fac(:)
418 33692928 : rate(:,540) = 2.6e-12_r8 * exp_fac(:)
419 33692928 : rate(:,544) = 2.6e-12_r8 * exp_fac(:)
420 33692928 : rate(:,303) = 6.9e-12_r8 * exp( -230._r8 * itemp(:) )
421 33692928 : rate(:,305) = 7.2e-11_r8 * exp( -70._r8 * itemp(:) )
422 33692928 : rate(:,306) = 7.66e-12_r8 * exp( -1020._r8 * itemp(:) )
423 33692928 : exp_fac(:) = exp( -1900._r8 * itemp(:) )
424 33692928 : rate(:,307) = 1.4e-12_r8 * exp_fac(:)
425 33692928 : rate(:,327) = 6.5e-15_r8 * exp_fac(:)
426 33692928 : exp_fac(:) = exp( 350._r8 * itemp(:) )
427 33692928 : rate(:,308) = 4.63e-12_r8 * exp_fac(:)
428 33692928 : rate(:,522) = 2.7e-12_r8 * exp_fac(:)
429 33692928 : rate(:,309) = 7.8e-13_r8 * exp( -1050._r8 * itemp(:) )
430 33692928 : exp_fac(:) = exp( 500._r8 * itemp(:) )
431 33692928 : rate(:,310) = 2.9e-12_r8 * exp_fac(:)
432 33692928 : rate(:,311) = 2e-12_r8 * exp_fac(:)
433 33692928 : rate(:,340) = 7.1e-13_r8 * exp_fac(:)
434 33692928 : rate(:,361) = 2e-12_r8 * exp_fac(:)
435 33692928 : rate(:,464) = 2e-12_r8 * exp_fac(:)
436 33692928 : rate(:,469) = 2e-12_r8 * exp_fac(:)
437 33692928 : rate(:,474) = 2e-12_r8 * exp_fac(:)
438 33692928 : exp_fac(:) = exp( 1040._r8 * itemp(:) )
439 33692928 : rate(:,312) = 4.3e-13_r8 * exp_fac(:)
440 33692928 : rate(:,362) = 4.3e-13_r8 * exp_fac(:)
441 33692928 : rate(:,415) = 4.3e-13_r8 * exp_fac(:)
442 33692928 : rate(:,429) = 4.3e-13_r8 * exp_fac(:)
443 33692928 : rate(:,432) = 4.3e-13_r8 * exp_fac(:)
444 33692928 : rate(:,435) = 4.3e-13_r8 * exp_fac(:)
445 33692928 : rate(:,314) = 3.15e-14_r8 * exp( 920._r8 * itemp(:) )
446 33692928 : rate(:,318) = 1.6e+11_r8 * exp( -4150._r8 * itemp(:) )
447 33692928 : rate(:,326) = 4.6e-13_r8 * exp( -1156._r8 * itemp(:) )
448 33692928 : rate(:,328) = 3.75e-13_r8 * exp( -40._r8 * itemp(:) )
449 33692928 : rate(:,332) = 9.19e-12_r8 * exp( -630._r8 * itemp(:) )
450 33692928 : rate(:,333) = 1.4e-12_r8 * exp( -1860._r8 * itemp(:) )
451 33692928 : rate(:,334) = 8.4e-13_r8 * exp( 830._r8 * itemp(:) )
452 33692928 : rate(:,348) = 4.8e-12_r8 * exp( 120._r8 * itemp(:) )
453 33692928 : rate(:,349) = 5.1e-14_r8 * exp( 693._r8 * itemp(:) )
454 33692928 : exp_fac(:) = exp( 360._r8 * itemp(:) )
455 33692928 : rate(:,355) = 2.7e-12_r8 * exp_fac(:)
456 33692928 : rate(:,356) = 1.3e-13_r8 * exp_fac(:)
457 33692928 : rate(:,358) = 9.6e-12_r8 * exp_fac(:)
458 33692928 : rate(:,364) = 5.3e-12_r8 * exp_fac(:)
459 33692928 : rate(:,401) = 2.7e-12_r8 * exp_fac(:)
460 33692928 : rate(:,412) = 2.7e-12_r8 * exp_fac(:)
461 33692928 : rate(:,514) = 2.7e-12_r8 * exp_fac(:)
462 33692928 : rate(:,533) = 2.7e-12_r8 * exp_fac(:)
463 33692928 : rate(:,357) = 1.5e-15_r8 * exp( -2100._r8 * itemp(:) )
464 33692928 : exp_fac(:) = exp( 530._r8 * itemp(:) )
465 33692928 : rate(:,360) = 4.6e-12_r8 * exp_fac(:)
466 33692928 : rate(:,363) = 2.3e-12_r8 * exp_fac(:)
467 33692928 : rate(:,368) = 2.3e-12_r8 * exp( -170._r8 * itemp(:) )
468 33692928 : rate(:,372) = 4.13e-12_r8 * exp( 452._r8 * itemp(:) )
469 33692928 : rate(:,378) = 5.4e-14_r8 * exp( 870._r8 * itemp(:) )
470 33692928 : exp_fac(:) = exp( 175._r8 * itemp(:) )
471 33692928 : rate(:,381) = 1.86e-11_r8 * exp_fac(:)
472 33692928 : rate(:,382) = 1.86e-11_r8 * exp_fac(:)
473 33692928 : rate(:,392) = 1.6e+09_r8 * exp( -8300._r8 * itemp(:) )
474 33692928 : exp_fac(:) = exp( -446._r8 * itemp(:) )
475 33692928 : rate(:,397) = 3.03e-12_r8 * exp_fac(:)
476 33692928 : rate(:,520) = 3.03e-12_r8 * exp_fac(:)
477 33692928 : exp_fac(:) = exp( 410._r8 * itemp(:) )
478 33692928 : rate(:,405) = 2.54e-11_r8 * exp_fac(:)
479 33692928 : rate(:,524) = 2.54e-11_r8 * exp_fac(:)
480 33692928 : rate(:,409) = 1.3e-12_r8 * exp( 640._r8 * itemp(:) )
481 33692928 : exp_fac(:) = exp( -193._r8 * itemp(:) )
482 33692928 : rate(:,417) = 2.3e-12_r8 * exp_fac(:)
483 33692928 : rate(:,517) = 2.3e-12_r8 * exp_fac(:)
484 33692928 : rate(:,421) = 5.9e-12_r8 * exp( 225._r8 * itemp(:) )
485 33692928 : rate(:,440) = 4.7e-13_r8 * exp( 1220._r8 * itemp(:) )
486 33692928 : exp_fac(:) = exp( 352._r8 * itemp(:) )
487 33692928 : rate(:,448) = 1.7e-12_r8 * exp_fac(:)
488 33692928 : rate(:,538) = 1.7e-12_r8 * exp_fac(:)
489 33692928 : exp_fac(:) = exp( 490._r8 * itemp(:) )
490 33692928 : rate(:,461) = 1.2e-12_r8 * exp_fac(:)
491 33692928 : rate(:,531) = 1.2e-12_r8 * exp_fac(:)
492 33692928 : exp_fac(:) = exp( -580._r8 * itemp(:) )
493 33692928 : rate(:,462) = 6.3e-16_r8 * exp_fac(:)
494 33692928 : rate(:,534) = 6.3e-16_r8 * exp_fac(:)
495 33692928 : exp_fac(:) = exp( 440._r8 * itemp(:) )
496 33692928 : rate(:,463) = 1.2e-11_r8 * exp_fac(:)
497 33692928 : rate(:,535) = 1.2e-11_r8 * exp_fac(:)
498 33692928 : rate(:,481) = 1.9e-13_r8 * exp( 520._r8 * itemp(:) )
499 33692928 : rate(:,482) = 1.1e-11_r8 * exp( -280._r8 * itemp(:) )
500 33692928 : rate(:,483) = 2.1e-11_r8 * exp( -2200._r8 * itemp(:) )
501 33692928 : rate(:,484) = 7.2e-14_r8 * exp( -1070._r8 * itemp(:) )
502 33692928 : rate(:,492) = 1.6e-13_r8 * exp( -2280._r8 * itemp(:) )
503 33692928 : rate(:,493) = 3.4e-12_r8 * exp( -1100._r8 * itemp(:) )
504 33692928 : rate(:,495) = 2.6e-11_r8 * exp( 330._r8 * itemp(:) )
505 33692928 : rate(:,498) = 1.7e-12_r8 * exp( -710._r8 * itemp(:) )
506 33692928 : exp_fac(:) = exp( 1300._r8 * itemp(:) )
507 33692928 : rate(:,513) = 2.75e-13_r8 * exp_fac(:)
508 33692928 : rate(:,521) = 2.12e-13_r8 * exp_fac(:)
509 33692928 : rate(:,532) = 2.6e-13_r8 * exp_fac(:)
510 :
511 33692928 : itemp(:) = 300._r8 * itemp(:)
512 :
513 72960 : n = ncol*pver
514 :
515 33692928 : ko(:) = 5.3e-32_r8 * itemp(:)**1.8_r8
516 33692928 : kinf(:) = 9.5e-11_r8 * itemp(:)**(-0.4_r8)
517 72960 : call jpl( rate(:,157), m, 0.6_r8, ko, kinf, n )
518 :
519 33692928 : ko(:) = 6.9e-31_r8 * itemp(:)**1._r8
520 33692928 : kinf(:) = 2.6e-11_r8
521 72960 : call jpl( rate(:,167), m, 0.6_r8, ko, kinf, n )
522 :
523 33692928 : ko(:) = 2.5e-31_r8 * itemp(:)**1.8_r8
524 33692928 : kinf(:) = 2.2e-11_r8 * itemp(:)**0.7_r8
525 72960 : call jpl( rate(:,177), m, 0.6_r8, ko, kinf, n )
526 :
527 33692928 : ko(:) = 9e-32_r8 * itemp(:)**1.5_r8
528 33692928 : kinf(:) = 3e-11_r8
529 72960 : call jpl( rate(:,185), m, 0.6_r8, ko, kinf, n )
530 :
531 33692928 : ko(:) = 1.9e-31_r8 * itemp(:)**3.4_r8
532 33692928 : kinf(:) = 4e-12_r8 * itemp(:)**0.3_r8
533 72960 : call jpl( rate(:,188), m, 0.6_r8, ko, kinf, n )
534 :
535 33692928 : ko(:) = 2.4e-30_r8 * itemp(:)**3._r8
536 33692928 : kinf(:) = 1.6e-12_r8 * itemp(:)**(-0.1_r8)
537 72960 : call jpl( rate(:,189), m, 0.6_r8, ko, kinf, n )
538 :
539 33692928 : ko(:) = 1.8e-30_r8 * itemp(:)**3._r8
540 33692928 : kinf(:) = 2.8e-11_r8
541 72960 : call jpl( rate(:,190), m, 0.6_r8, ko, kinf, n )
542 :
543 33692928 : ko(:) = 1.8e-31_r8 * itemp(:)**3.4_r8
544 33692928 : kinf(:) = 1.5e-11_r8 * itemp(:)**1.9_r8
545 72960 : call jpl( rate(:,208), m, 0.6_r8, ko, kinf, n )
546 :
547 33692928 : ko(:) = 1.9e-32_r8 * itemp(:)**3.6_r8
548 33692928 : kinf(:) = 3.7e-12_r8 * itemp(:)**1.6_r8
549 72960 : call jpl( rate(:,228), m, 0.6_r8, ko, kinf, n )
550 :
551 33692928 : ko(:) = 5.2e-31_r8 * itemp(:)**3.2_r8
552 33692928 : kinf(:) = 6.9e-12_r8 * itemp(:)**2.9_r8
553 72960 : call jpl( rate(:,239), m, 0.6_r8, ko, kinf, n )
554 :
555 33692928 : ko(:) = 6.1e-33_r8 * itemp(:)**1.5_r8
556 33692928 : kinf(:) = 9.8e-15_r8 * itemp(:)**(-4.6_r8)
557 72960 : call jpl( rate(:,285), m, 0.8_r8, ko, kinf, n )
558 :
559 33692928 : ko(:) = 5.2e-30_r8 * itemp(:)**2.4_r8
560 33692928 : kinf(:) = 2.2e-10_r8 * itemp(:)**0.7_r8
561 72960 : call jpl( rate(:,295), m, 0.6_r8, ko, kinf, n )
562 :
563 33692928 : ko(:) = 5.5e-30_r8
564 33692928 : kinf(:) = 8.3e-13_r8 * itemp(:)**(-2._r8)
565 72960 : call jpl( rate(:,296), m, 0.6_r8, ko, kinf, n )
566 :
567 33692928 : ko(:) = 1.6e-29_r8 * itemp(:)**3.3_r8
568 33692928 : kinf(:) = 3.1e-10_r8 * itemp(:)
569 72960 : call jpl( rate(:,297), m, 0.6_r8, ko, kinf, n )
570 :
571 33692928 : ko(:) = 8.6e-29_r8 * itemp(:)**3.1_r8
572 33692928 : kinf(:) = 9e-12_r8 * itemp(:)**0.85_r8
573 72960 : call jpl( rate(:,323), m, 0.48_r8, ko, kinf, n )
574 :
575 33692928 : ko(:) = 7.3e-29_r8 * itemp(:)**4.1_r8
576 33692928 : kinf(:) = 9.5e-12_r8 * itemp(:)**1.6_r8
577 72960 : call jpl( rate(:,324), m, 0.6_r8, ko, kinf, n )
578 :
579 33692928 : ko(:) = 8e-27_r8 * itemp(:)**3.5_r8
580 33692928 : kinf(:) = 3e-11_r8
581 72960 : call jpl( rate(:,344), m, 0.5_r8, ko, kinf, n )
582 :
583 33692928 : ko(:) = 8e-27_r8 * itemp(:)**3.5_r8
584 33692928 : kinf(:) = 3e-11_r8
585 72960 : call jpl( rate(:,370), m, 0.5_r8, ko, kinf, n )
586 :
587 33692928 : ko(:) = 9.7e-29_r8 * itemp(:)**5.6_r8
588 33692928 : kinf(:) = 9.3e-12_r8 * itemp(:)**1.5_r8
589 72960 : call jpl( rate(:,431), m, 0.6_r8, ko, kinf, n )
590 :
591 33692928 : ko(:) = 9.7e-29_r8 * itemp(:)**5.6_r8
592 33692928 : kinf(:) = 9.3e-12_r8 * itemp(:)**1.5_r8
593 72960 : call jpl( rate(:,434), m, 0.6_r8, ko, kinf, n )
594 :
595 33692928 : ko(:) = 9.7e-29_r8 * itemp(:)**5.6_r8
596 33692928 : kinf(:) = 9.3e-12_r8 * itemp(:)**1.5_r8
597 72960 : call jpl( rate(:,437), m, 0.6_r8, ko, kinf, n )
598 :
599 33692928 : ko(:) = 9.7e-29_r8 * itemp(:)**5.6_r8
600 33692928 : kinf(:) = 9.3e-12_r8 * itemp(:)**1.5_r8
601 72960 : call jpl( rate(:,444), m, 0.6_r8, ko, kinf, n )
602 :
603 33692928 : ko(:) = 2.9e-31_r8 * itemp(:)**4.1_r8
604 33692928 : kinf(:) = 1.7e-12_r8 * itemp(:)**(-0.2_r8)
605 72960 : call jpl( rate(:,486), m, 0.6_r8, ko, kinf, n )
606 :
607 72960 : end subroutine setrxt
608 :
609 :
610 0 : subroutine setrxt_hrates( rate, temp, m, ncol, kbot )
611 :
612 : use ppgrid, only : pcols, pver
613 :
614 :
615 : use chem_mods, only : rxntot
616 : use mo_jpl, only : jpl
617 :
618 : implicit none
619 :
620 : !-------------------------------------------------------
621 : ! ... dummy arguments
622 : !-------------------------------------------------------
623 : integer, intent(in) :: ncol
624 : integer, intent(in) :: kbot
625 : real(r8), intent(in) :: temp(pcols,pver)
626 : real(r8), intent(in) :: m(ncol*pver)
627 : real(r8), intent(inout) :: rate(ncol*pver,max(1,rxntot))
628 :
629 : !-------------------------------------------------------
630 : ! ... local variables
631 : !-------------------------------------------------------
632 : integer :: n
633 : integer :: offset
634 : integer :: k
635 0 : real(r8) :: itemp(ncol*kbot)
636 : real(r8) :: exp_fac(ncol*kbot)
637 0 : real(r8) :: ko(ncol*kbot)
638 0 : real(r8) :: kinf(ncol*kbot)
639 0 : real(r8) :: wrk(ncol*kbot)
640 :
641 0 : n = ncol*kbot
642 :
643 0 : rate(:n,154) = 6.9e-12_r8
644 :
645 0 : do k = 1,kbot
646 0 : offset = (k-1)*ncol
647 0 : itemp(offset+1:offset+ncol) = 1._r8 / temp(:ncol,k)
648 : end do
649 :
650 0 : rate(:n,146) = 2.15e-11_r8 * exp( 110._r8 * itemp(:) )
651 0 : rate(:n,149) = 8e-12_r8 * exp( -2060._r8 * itemp(:) )
652 0 : rate(:n,158) = 3e-11_r8 * exp( 200._r8 * itemp(:) )
653 0 : rate(:n,159) = 1e-14_r8 * exp( -490._r8 * itemp(:) )
654 0 : rate(:n,160) = 1.4e-10_r8 * exp( -470._r8 * itemp(:) )
655 0 : rate(:n,163) = 4.8e-11_r8 * exp( 250._r8 * itemp(:) )
656 0 : rate(:n,164) = 1.8e-11_r8 * exp( 180._r8 * itemp(:) )
657 0 : rate(:n,165) = 1.7e-12_r8 * exp( -940._r8 * itemp(:) )
658 0 : rate(:n,170) = 2.1e-11_r8 * exp( 100._r8 * itemp(:) )
659 0 : rate(:n,174) = 3.3e-12_r8 * exp( -3150._r8 * itemp(:) )
660 0 : rate(:n,175) = 5.1e-12_r8 * exp( 210._r8 * itemp(:) )
661 0 : rate(:n,183) = 3.44e-12_r8 * exp( 260._r8 * itemp(:) )
662 0 : rate(:n,184) = 3e-12_r8 * exp( -1500._r8 * itemp(:) )
663 :
664 0 : itemp(:) = 300._r8 * itemp(:)
665 :
666 0 : ko(:) = 5.3e-32_r8 * itemp(:)**1.8_r8
667 0 : kinf(:) = 9.5e-11_r8 * itemp(:)**(-0.4_r8)
668 0 : call jpl( wrk, m, 0.6_r8, ko, kinf, n )
669 0 : rate(:n,157) = wrk(:)
670 :
671 :
672 :
673 :
674 :
675 :
676 :
677 :
678 :
679 :
680 :
681 :
682 :
683 :
684 :
685 :
686 :
687 :
688 :
689 :
690 :
691 :
692 :
693 0 : end subroutine setrxt_hrates
694 :
695 : end module mo_setrxt
|