Line data Source code
1 : module mo_adjrxt
2 : private
3 : public :: adjrxt
4 : contains
5 153600 : subroutine adjrxt( rate, inv, m, ncol, nlev )
6 : use shr_kind_mod, only : r8 => shr_kind_r8
7 : use chem_mods, only : nfs, rxntot
8 : implicit none
9 : !--------------------------------------------------------------------
10 : ! ... dummy arguments
11 : !--------------------------------------------------------------------
12 : integer, intent(in) :: ncol, nlev
13 : real(r8), intent(in) :: inv(ncol,nlev,nfs)
14 : real(r8), intent(in) :: m(ncol,nlev)
15 : real(r8), intent(inout) :: rate(ncol,nlev,rxntot)
16 165734400 : rate(:,:, 96) = rate(:,:, 96) * inv(:,:, 2)
17 165734400 : rate(:,:, 100) = rate(:,:, 100) * inv(:,:, 2)
18 165734400 : rate(:,:, 104) = rate(:,:, 104) * inv(:,:, 2)
19 165734400 : rate(:,:, 109) = rate(:,:, 109) * inv(:,:, 1)
20 165734400 : rate(:,:, 110) = rate(:,:, 110) * inv(:,:, 1)
21 165734400 : rate(:,:, 116) = rate(:,:, 116) * inv(:,:, 1)
22 165734400 : rate(:,:, 126) = rate(:,:, 126) * inv(:,:, 1)
23 165734400 : rate(:,:, 138) = rate(:,:, 138) * inv(:,:, 1)
24 165734400 : rate(:,:, 146) = rate(:,:, 146) * inv(:,:, 1)
25 165734400 : rate(:,:, 149) = rate(:,:, 149) * inv(:,:, 1)
26 165734400 : rate(:,:, 150) = rate(:,:, 150) * inv(:,:, 1)
27 165734400 : rate(:,:, 151) = rate(:,:, 151) * inv(:,:, 1)
28 165734400 : rate(:,:, 153) = rate(:,:, 153) * inv(:,:, 1)
29 165734400 : rate(:,:, 154) = rate(:,:, 154) * inv(:,:, 1)
30 165734400 : rate(:,:, 169) = rate(:,:, 169) * inv(:,:, 1)
31 165734400 : rate(:,:, 189) = rate(:,:, 189) * inv(:,:, 1)
32 165734400 : rate(:,:, 190) = rate(:,:, 190) * inv(:,:, 1)
33 165734400 : rate(:,:, 200) = rate(:,:, 200) * inv(:,:, 1)
34 165734400 : rate(:,:, 242) = rate(:,:, 242) * inv(:,:, 1)
35 165734400 : rate(:,:, 298) = rate(:,:, 298) * inv(:,:, 2)
36 165734400 : rate(:,:, 301) = rate(:,:, 301) * inv(:,:, 2)
37 165734400 : rate(:,:, 305) = rate(:,:, 305) * inv(:,:, 2)
38 165734400 : rate(:,:, 310) = rate(:,:, 310) * inv(:,:, 2)
39 165734400 : rate(:,:, 311) = rate(:,:, 311) * inv(:,:, 2)
40 165734400 : rate(:,:, 94) = rate(:,:, 94) * m(:,:)
41 165734400 : rate(:,:, 95) = rate(:,:, 95) * m(:,:)
42 165734400 : rate(:,:, 97) = rate(:,:, 97) * m(:,:)
43 165734400 : rate(:,:, 98) = rate(:,:, 98) * m(:,:)
44 165734400 : rate(:,:, 99) = rate(:,:, 99) * m(:,:)
45 165734400 : rate(:,:, 101) = rate(:,:, 101) * m(:,:)
46 165734400 : rate(:,:, 102) = rate(:,:, 102) * m(:,:)
47 165734400 : rate(:,:, 103) = rate(:,:, 103) * m(:,:)
48 165734400 : rate(:,:, 105) = rate(:,:, 105) * m(:,:)
49 165734400 : rate(:,:, 106) = rate(:,:, 106) * m(:,:)
50 165734400 : rate(:,:, 107) = rate(:,:, 107) * m(:,:)
51 165734400 : rate(:,:, 108) = rate(:,:, 108) * m(:,:)
52 165734400 : rate(:,:, 109) = rate(:,:, 109) * m(:,:)
53 165734400 : rate(:,:, 110) = rate(:,:, 110) * m(:,:)
54 165734400 : rate(:,:, 111) = rate(:,:, 111) * m(:,:)
55 165734400 : rate(:,:, 112) = rate(:,:, 112) * m(:,:)
56 165734400 : rate(:,:, 113) = rate(:,:, 113) * m(:,:)
57 165734400 : rate(:,:, 114) = rate(:,:, 114) * m(:,:)
58 165734400 : rate(:,:, 115) = rate(:,:, 115) * m(:,:)
59 165734400 : rate(:,:, 116) = rate(:,:, 116) * m(:,:)
60 165734400 : rate(:,:, 117) = rate(:,:, 117) * m(:,:)
61 165734400 : rate(:,:, 118) = rate(:,:, 118) * m(:,:)
62 165734400 : rate(:,:, 119) = rate(:,:, 119) * m(:,:)
63 165734400 : rate(:,:, 120) = rate(:,:, 120) * m(:,:)
64 165734400 : rate(:,:, 121) = rate(:,:, 121) * m(:,:)
65 165734400 : rate(:,:, 122) = rate(:,:, 122) * m(:,:)
66 165734400 : rate(:,:, 123) = rate(:,:, 123) * m(:,:)
67 165734400 : rate(:,:, 124) = rate(:,:, 124) * m(:,:)
68 165734400 : rate(:,:, 125) = rate(:,:, 125) * m(:,:)
69 165734400 : rate(:,:, 126) = rate(:,:, 126) * m(:,:)
70 165734400 : rate(:,:, 127) = rate(:,:, 127) * m(:,:)
71 165734400 : rate(:,:, 128) = rate(:,:, 128) * m(:,:)
72 165734400 : rate(:,:, 129) = rate(:,:, 129) * m(:,:)
73 165734400 : rate(:,:, 130) = rate(:,:, 130) * m(:,:)
74 165734400 : rate(:,:, 131) = rate(:,:, 131) * m(:,:)
75 165734400 : rate(:,:, 132) = rate(:,:, 132) * m(:,:)
76 165734400 : rate(:,:, 133) = rate(:,:, 133) * m(:,:)
77 165734400 : rate(:,:, 134) = rate(:,:, 134) * m(:,:)
78 165734400 : rate(:,:, 135) = rate(:,:, 135) * m(:,:)
79 165734400 : rate(:,:, 136) = rate(:,:, 136) * m(:,:)
80 165734400 : rate(:,:, 137) = rate(:,:, 137) * m(:,:)
81 165734400 : rate(:,:, 138) = rate(:,:, 138) * m(:,:)
82 165734400 : rate(:,:, 139) = rate(:,:, 139) * m(:,:)
83 165734400 : rate(:,:, 140) = rate(:,:, 140) * m(:,:)
84 165734400 : rate(:,:, 141) = rate(:,:, 141) * m(:,:)
85 165734400 : rate(:,:, 142) = rate(:,:, 142) * m(:,:)
86 165734400 : rate(:,:, 143) = rate(:,:, 143) * m(:,:)
87 165734400 : rate(:,:, 144) = rate(:,:, 144) * m(:,:)
88 165734400 : rate(:,:, 145) = rate(:,:, 145) * m(:,:)
89 165734400 : rate(:,:, 146) = rate(:,:, 146) * m(:,:)
90 165734400 : rate(:,:, 147) = rate(:,:, 147) * m(:,:)
91 165734400 : rate(:,:, 148) = rate(:,:, 148) * m(:,:)
92 165734400 : rate(:,:, 149) = rate(:,:, 149) * m(:,:)
93 165734400 : rate(:,:, 150) = rate(:,:, 150) * m(:,:)
94 165734400 : rate(:,:, 151) = rate(:,:, 151) * m(:,:)
95 165734400 : rate(:,:, 152) = rate(:,:, 152) * m(:,:)
96 165734400 : rate(:,:, 155) = rate(:,:, 155) * m(:,:)
97 165734400 : rate(:,:, 156) = rate(:,:, 156) * m(:,:)
98 165734400 : rate(:,:, 157) = rate(:,:, 157) * m(:,:)
99 165734400 : rate(:,:, 158) = rate(:,:, 158) * m(:,:)
100 165734400 : rate(:,:, 159) = rate(:,:, 159) * m(:,:)
101 165734400 : rate(:,:, 160) = rate(:,:, 160) * m(:,:)
102 165734400 : rate(:,:, 161) = rate(:,:, 161) * m(:,:)
103 165734400 : rate(:,:, 162) = rate(:,:, 162) * m(:,:)
104 165734400 : rate(:,:, 163) = rate(:,:, 163) * m(:,:)
105 165734400 : rate(:,:, 164) = rate(:,:, 164) * m(:,:)
106 165734400 : rate(:,:, 165) = rate(:,:, 165) * m(:,:)
107 165734400 : rate(:,:, 166) = rate(:,:, 166) * m(:,:)
108 165734400 : rate(:,:, 167) = rate(:,:, 167) * m(:,:)
109 165734400 : rate(:,:, 168) = rate(:,:, 168) * m(:,:)
110 165734400 : rate(:,:, 169) = rate(:,:, 169) * m(:,:)
111 165734400 : rate(:,:, 170) = rate(:,:, 170) * m(:,:)
112 165734400 : rate(:,:, 171) = rate(:,:, 171) * m(:,:)
113 165734400 : rate(:,:, 172) = rate(:,:, 172) * m(:,:)
114 165734400 : rate(:,:, 173) = rate(:,:, 173) * m(:,:)
115 165734400 : rate(:,:, 174) = rate(:,:, 174) * m(:,:)
116 165734400 : rate(:,:, 175) = rate(:,:, 175) * m(:,:)
117 165734400 : rate(:,:, 176) = rate(:,:, 176) * m(:,:)
118 165734400 : rate(:,:, 177) = rate(:,:, 177) * m(:,:)
119 165734400 : rate(:,:, 178) = rate(:,:, 178) * m(:,:)
120 165734400 : rate(:,:, 179) = rate(:,:, 179) * m(:,:)
121 165734400 : rate(:,:, 180) = rate(:,:, 180) * m(:,:)
122 165734400 : rate(:,:, 181) = rate(:,:, 181) * m(:,:)
123 165734400 : rate(:,:, 182) = rate(:,:, 182) * m(:,:)
124 165734400 : rate(:,:, 183) = rate(:,:, 183) * m(:,:)
125 165734400 : rate(:,:, 184) = rate(:,:, 184) * m(:,:)
126 165734400 : rate(:,:, 185) = rate(:,:, 185) * m(:,:)
127 165734400 : rate(:,:, 186) = rate(:,:, 186) * m(:,:)
128 165734400 : rate(:,:, 187) = rate(:,:, 187) * m(:,:)
129 165734400 : rate(:,:, 188) = rate(:,:, 188) * m(:,:)
130 165734400 : rate(:,:, 189) = rate(:,:, 189) * m(:,:)
131 165734400 : rate(:,:, 191) = rate(:,:, 191) * m(:,:)
132 165734400 : rate(:,:, 192) = rate(:,:, 192) * m(:,:)
133 165734400 : rate(:,:, 193) = rate(:,:, 193) * m(:,:)
134 165734400 : rate(:,:, 194) = rate(:,:, 194) * m(:,:)
135 165734400 : rate(:,:, 195) = rate(:,:, 195) * m(:,:)
136 165734400 : rate(:,:, 196) = rate(:,:, 196) * m(:,:)
137 165734400 : rate(:,:, 197) = rate(:,:, 197) * m(:,:)
138 165734400 : rate(:,:, 198) = rate(:,:, 198) * m(:,:)
139 165734400 : rate(:,:, 199) = rate(:,:, 199) * m(:,:)
140 165734400 : rate(:,:, 200) = rate(:,:, 200) * m(:,:)
141 165734400 : rate(:,:, 201) = rate(:,:, 201) * m(:,:)
142 165734400 : rate(:,:, 202) = rate(:,:, 202) * m(:,:)
143 165734400 : rate(:,:, 203) = rate(:,:, 203) * m(:,:)
144 165734400 : rate(:,:, 204) = rate(:,:, 204) * m(:,:)
145 165734400 : rate(:,:, 205) = rate(:,:, 205) * m(:,:)
146 165734400 : rate(:,:, 206) = rate(:,:, 206) * m(:,:)
147 165734400 : rate(:,:, 207) = rate(:,:, 207) * m(:,:)
148 165734400 : rate(:,:, 208) = rate(:,:, 208) * m(:,:)
149 165734400 : rate(:,:, 209) = rate(:,:, 209) * m(:,:)
150 165734400 : rate(:,:, 210) = rate(:,:, 210) * m(:,:)
151 165734400 : rate(:,:, 211) = rate(:,:, 211) * m(:,:)
152 165734400 : rate(:,:, 212) = rate(:,:, 212) * m(:,:)
153 165734400 : rate(:,:, 213) = rate(:,:, 213) * m(:,:)
154 165734400 : rate(:,:, 214) = rate(:,:, 214) * m(:,:)
155 165734400 : rate(:,:, 215) = rate(:,:, 215) * m(:,:)
156 165734400 : rate(:,:, 216) = rate(:,:, 216) * m(:,:)
157 165734400 : rate(:,:, 217) = rate(:,:, 217) * m(:,:)
158 165734400 : rate(:,:, 218) = rate(:,:, 218) * m(:,:)
159 165734400 : rate(:,:, 219) = rate(:,:, 219) * m(:,:)
160 165734400 : rate(:,:, 220) = rate(:,:, 220) * m(:,:)
161 165734400 : rate(:,:, 221) = rate(:,:, 221) * m(:,:)
162 165734400 : rate(:,:, 222) = rate(:,:, 222) * m(:,:)
163 165734400 : rate(:,:, 223) = rate(:,:, 223) * m(:,:)
164 165734400 : rate(:,:, 224) = rate(:,:, 224) * m(:,:)
165 165734400 : rate(:,:, 225) = rate(:,:, 225) * m(:,:)
166 165734400 : rate(:,:, 226) = rate(:,:, 226) * m(:,:)
167 165734400 : rate(:,:, 227) = rate(:,:, 227) * m(:,:)
168 165734400 : rate(:,:, 228) = rate(:,:, 228) * m(:,:)
169 165734400 : rate(:,:, 229) = rate(:,:, 229) * m(:,:)
170 165734400 : rate(:,:, 230) = rate(:,:, 230) * m(:,:)
171 165734400 : rate(:,:, 231) = rate(:,:, 231) * m(:,:)
172 165734400 : rate(:,:, 232) = rate(:,:, 232) * m(:,:)
173 165734400 : rate(:,:, 233) = rate(:,:, 233) * m(:,:)
174 165734400 : rate(:,:, 234) = rate(:,:, 234) * m(:,:)
175 165734400 : rate(:,:, 235) = rate(:,:, 235) * m(:,:)
176 165734400 : rate(:,:, 236) = rate(:,:, 236) * m(:,:)
177 165734400 : rate(:,:, 237) = rate(:,:, 237) * m(:,:)
178 165734400 : rate(:,:, 238) = rate(:,:, 238) * m(:,:)
179 165734400 : rate(:,:, 239) = rate(:,:, 239) * m(:,:)
180 165734400 : rate(:,:, 240) = rate(:,:, 240) * m(:,:)
181 165734400 : rate(:,:, 241) = rate(:,:, 241) * m(:,:)
182 165734400 : rate(:,:, 242) = rate(:,:, 242) * m(:,:)
183 165734400 : rate(:,:, 243) = rate(:,:, 243) * m(:,:)
184 165734400 : rate(:,:, 244) = rate(:,:, 244) * m(:,:)
185 165734400 : rate(:,:, 245) = rate(:,:, 245) * m(:,:)
186 165734400 : rate(:,:, 246) = rate(:,:, 246) * m(:,:)
187 165734400 : rate(:,:, 247) = rate(:,:, 247) * m(:,:)
188 165734400 : rate(:,:, 248) = rate(:,:, 248) * m(:,:)
189 165734400 : rate(:,:, 249) = rate(:,:, 249) * m(:,:)
190 165734400 : rate(:,:, 250) = rate(:,:, 250) * m(:,:)
191 165734400 : rate(:,:, 251) = rate(:,:, 251) * m(:,:)
192 165734400 : rate(:,:, 252) = rate(:,:, 252) * m(:,:)
193 165734400 : rate(:,:, 253) = rate(:,:, 253) * m(:,:)
194 165734400 : rate(:,:, 254) = rate(:,:, 254) * m(:,:)
195 165734400 : rate(:,:, 255) = rate(:,:, 255) * m(:,:)
196 165734400 : rate(:,:, 256) = rate(:,:, 256) * m(:,:)
197 165734400 : rate(:,:, 257) = rate(:,:, 257) * m(:,:)
198 165734400 : rate(:,:, 258) = rate(:,:, 258) * m(:,:)
199 165734400 : rate(:,:, 259) = rate(:,:, 259) * m(:,:)
200 165734400 : rate(:,:, 260) = rate(:,:, 260) * m(:,:)
201 165734400 : rate(:,:, 261) = rate(:,:, 261) * m(:,:)
202 165734400 : rate(:,:, 262) = rate(:,:, 262) * m(:,:)
203 165734400 : rate(:,:, 263) = rate(:,:, 263) * m(:,:)
204 165734400 : rate(:,:, 269) = rate(:,:, 269) * m(:,:)
205 165734400 : rate(:,:, 274) = rate(:,:, 274) * m(:,:)
206 165734400 : rate(:,:, 275) = rate(:,:, 275) * m(:,:)
207 165734400 : rate(:,:, 276) = rate(:,:, 276) * m(:,:)
208 165734400 : rate(:,:, 279) = rate(:,:, 279) * m(:,:)
209 165734400 : rate(:,:, 280) = rate(:,:, 280) * m(:,:)
210 165734400 : rate(:,:, 281) = rate(:,:, 281) * m(:,:)
211 165734400 : rate(:,:, 284) = rate(:,:, 284) * m(:,:)
212 165734400 : rate(:,:, 288) = rate(:,:, 288) * m(:,:)
213 165734400 : rate(:,:, 289) = rate(:,:, 289) * m(:,:)
214 165734400 : rate(:,:, 290) = rate(:,:, 290) * m(:,:)
215 165734400 : rate(:,:, 291) = rate(:,:, 291) * m(:,:)
216 165734400 : rate(:,:, 292) = rate(:,:, 292) * m(:,:)
217 165734400 : rate(:,:, 293) = rate(:,:, 293) * m(:,:)
218 165734400 : rate(:,:, 294) = rate(:,:, 294) * m(:,:)
219 165734400 : rate(:,:, 295) = rate(:,:, 295) * m(:,:)
220 165734400 : rate(:,:, 296) = rate(:,:, 296) * m(:,:)
221 165734400 : rate(:,:, 297) = rate(:,:, 297) * m(:,:)
222 165734400 : rate(:,:, 299) = rate(:,:, 299) * m(:,:)
223 165734400 : rate(:,:, 300) = rate(:,:, 300) * m(:,:)
224 165734400 : rate(:,:, 302) = rate(:,:, 302) * m(:,:)
225 165734400 : rate(:,:, 303) = rate(:,:, 303) * m(:,:)
226 165734400 : rate(:,:, 304) = rate(:,:, 304) * m(:,:)
227 165734400 : rate(:,:, 306) = rate(:,:, 306) * m(:,:)
228 165734400 : rate(:,:, 307) = rate(:,:, 307) * m(:,:)
229 165734400 : rate(:,:, 308) = rate(:,:, 308) * m(:,:)
230 165734400 : rate(:,:, 309) = rate(:,:, 309) * m(:,:)
231 165734400 : rate(:,:, 312) = rate(:,:, 312) * m(:,:)
232 153600 : end subroutine adjrxt
233 : end module mo_adjrxt
|
