Line data Source code
1 : module mo_lu_solve
2 : use chem_mods, only: veclen
3 : private
4 : public :: lu_slv
5 : contains
6 9949925 : subroutine lu_slv01( avec_len, lu, b )
7 : use shr_kind_mod, only : r8 => shr_kind_r8
8 : use chem_mods, only : clscnt4, nzcnt
9 : implicit none
10 : !-----------------------------------------------------------------------
11 : ! ... Dummy args
12 : !-----------------------------------------------------------------------
13 : integer, intent(in) :: avec_len
14 : real(r8), intent(in) :: lu(veclen,max(1,nzcnt))
15 : real(r8), intent(inout) :: b(veclen,clscnt4)
16 : !-----------------------------------------------------------------------
17 : ! ... Local variables
18 : !-----------------------------------------------------------------------
19 : integer :: k
20 : !-----------------------------------------------------------------------
21 : ! ... solve L * y = b
22 : !-----------------------------------------------------------------------
23 325174377 : do k = 1,avec_len
24 315224452 : b(k,224) = b(k,224) - lu(k,96) * b(k,54)
25 315224452 : b(k,234) = b(k,234) - lu(k,97) * b(k,54)
26 315224452 : b(k,226) = b(k,226) - lu(k,99) * b(k,55)
27 315224452 : b(k,239) = b(k,239) - lu(k,100) * b(k,55)
28 315224452 : b(k,230) = b(k,230) - lu(k,102) * b(k,56)
29 315224452 : b(k,234) = b(k,234) - lu(k,103) * b(k,56)
30 315224452 : b(k,226) = b(k,226) - lu(k,105) * b(k,57)
31 315224452 : b(k,237) = b(k,237) - lu(k,106) * b(k,57)
32 315224452 : b(k,89) = b(k,89) - lu(k,108) * b(k,58)
33 315224452 : b(k,221) = b(k,221) - lu(k,109) * b(k,58)
34 315224452 : b(k,230) = b(k,230) - lu(k,110) * b(k,58)
35 315224452 : b(k,91) = b(k,91) - lu(k,112) * b(k,59)
36 315224452 : b(k,230) = b(k,230) - lu(k,113) * b(k,59)
37 315224452 : b(k,234) = b(k,234) - lu(k,114) * b(k,59)
38 315224452 : b(k,89) = b(k,89) - lu(k,116) * b(k,60)
39 315224452 : b(k,230) = b(k,230) - lu(k,117) * b(k,60)
40 315224452 : b(k,234) = b(k,234) - lu(k,118) * b(k,60)
41 315224452 : b(k,89) = b(k,89) - lu(k,120) * b(k,61)
42 315224452 : b(k,230) = b(k,230) - lu(k,121) * b(k,61)
43 315224452 : b(k,234) = b(k,234) - lu(k,122) * b(k,61)
44 315224452 : b(k,226) = b(k,226) - lu(k,124) * b(k,62)
45 315224452 : b(k,234) = b(k,234) - lu(k,125) * b(k,62)
46 315224452 : b(k,239) = b(k,239) - lu(k,126) * b(k,62)
47 315224452 : b(k,129) = b(k,129) - lu(k,128) * b(k,63)
48 315224452 : b(k,226) = b(k,226) - lu(k,129) * b(k,63)
49 315224452 : b(k,94) = b(k,94) - lu(k,131) * b(k,64)
50 315224452 : b(k,239) = b(k,239) - lu(k,132) * b(k,64)
51 315224452 : b(k,207) = b(k,207) - lu(k,134) * b(k,65)
52 315224452 : b(k,226) = b(k,226) - lu(k,135) * b(k,65)
53 315224452 : b(k,139) = b(k,139) - lu(k,137) * b(k,66)
54 315224452 : b(k,228) = b(k,228) - lu(k,138) * b(k,66)
55 315224452 : b(k,89) = b(k,89) - lu(k,140) * b(k,67)
56 315224452 : b(k,221) = b(k,221) - lu(k,141) * b(k,67)
57 315224452 : b(k,230) = b(k,230) - lu(k,142) * b(k,67)
58 315224452 : b(k,234) = b(k,234) - lu(k,143) * b(k,67)
59 315224452 : b(k,89) = b(k,89) - lu(k,145) * b(k,68)
60 315224452 : b(k,184) = b(k,184) - lu(k,146) * b(k,68)
61 315224452 : b(k,221) = b(k,221) - lu(k,147) * b(k,68)
62 315224452 : b(k,230) = b(k,230) - lu(k,148) * b(k,68)
63 315224452 : b(k,89) = b(k,89) - lu(k,150) * b(k,69)
64 315224452 : b(k,91) = b(k,91) - lu(k,151) * b(k,69)
65 315224452 : b(k,230) = b(k,230) - lu(k,152) * b(k,69)
66 315224452 : b(k,234) = b(k,234) - lu(k,153) * b(k,69)
67 315224452 : b(k,89) = b(k,89) - lu(k,155) * b(k,70)
68 315224452 : b(k,184) = b(k,184) - lu(k,156) * b(k,70)
69 315224452 : b(k,230) = b(k,230) - lu(k,157) * b(k,70)
70 315224452 : b(k,234) = b(k,234) - lu(k,158) * b(k,70)
71 315224452 : b(k,72) = b(k,72) - lu(k,161) * b(k,71)
72 315224452 : b(k,73) = b(k,73) - lu(k,162) * b(k,71)
73 315224452 : b(k,135) = b(k,135) - lu(k,163) * b(k,71)
74 315224452 : b(k,226) = b(k,226) - lu(k,164) * b(k,71)
75 315224452 : b(k,237) = b(k,237) - lu(k,165) * b(k,71)
76 315224452 : b(k,130) = b(k,130) - lu(k,167) * b(k,72)
77 315224452 : b(k,201) = b(k,201) - lu(k,168) * b(k,72)
78 315224452 : b(k,237) = b(k,237) - lu(k,169) * b(k,72)
79 315224452 : b(k,127) = b(k,127) - lu(k,171) * b(k,73)
80 315224452 : b(k,131) = b(k,131) - lu(k,172) * b(k,73)
81 315224452 : b(k,226) = b(k,226) - lu(k,173) * b(k,73)
82 315224452 : b(k,237) = b(k,237) - lu(k,174) * b(k,73)
83 315224452 : b(k,234) = b(k,234) - lu(k,176) * b(k,74)
84 315224452 : b(k,222) = b(k,222) - lu(k,178) * b(k,75)
85 315224452 : b(k,222) = b(k,222) - lu(k,181) * b(k,76)
86 315224452 : b(k,226) = b(k,226) - lu(k,183) * b(k,77)
87 315224452 : b(k,230) = b(k,230) - lu(k,184) * b(k,77)
88 315224452 : b(k,237) = b(k,237) - lu(k,185) * b(k,77)
89 315224452 : b(k,79) = b(k,79) - lu(k,188) * b(k,78)
90 315224452 : b(k,80) = b(k,80) - lu(k,189) * b(k,78)
91 315224452 : b(k,125) = b(k,125) - lu(k,190) * b(k,78)
92 315224452 : b(k,166) = b(k,166) - lu(k,191) * b(k,78)
93 315224452 : b(k,226) = b(k,226) - lu(k,192) * b(k,78)
94 315224452 : b(k,237) = b(k,237) - lu(k,193) * b(k,78)
95 315224452 : b(k,127) = b(k,127) - lu(k,195) * b(k,79)
96 315224452 : b(k,131) = b(k,131) - lu(k,196) * b(k,79)
97 315224452 : b(k,226) = b(k,226) - lu(k,197) * b(k,79)
98 315224452 : b(k,237) = b(k,237) - lu(k,198) * b(k,79)
99 315224452 : b(k,201) = b(k,201) - lu(k,200) * b(k,80)
100 315224452 : b(k,216) = b(k,216) - lu(k,201) * b(k,80)
101 315224452 : b(k,237) = b(k,237) - lu(k,202) * b(k,80)
102 315224452 : b(k,207) = b(k,207) - lu(k,204) * b(k,81)
103 315224452 : b(k,226) = b(k,226) - lu(k,205) * b(k,81)
104 315224452 : b(k,83) = b(k,83) - lu(k,209) * b(k,82)
105 315224452 : b(k,125) = b(k,125) - lu(k,210) * b(k,82)
106 315224452 : b(k,167) = b(k,167) - lu(k,211) * b(k,82)
107 315224452 : b(k,201) = b(k,201) - lu(k,212) * b(k,82)
108 315224452 : b(k,216) = b(k,216) - lu(k,213) * b(k,82)
109 315224452 : b(k,226) = b(k,226) - lu(k,214) * b(k,82)
110 315224452 : b(k,237) = b(k,237) - lu(k,215) * b(k,82)
111 315224452 : b(k,131) = b(k,131) - lu(k,217) * b(k,83)
112 315224452 : b(k,137) = b(k,137) - lu(k,218) * b(k,83)
113 315224452 : b(k,226) = b(k,226) - lu(k,219) * b(k,83)
114 315224452 : b(k,237) = b(k,237) - lu(k,220) * b(k,83)
115 315224452 : b(k,149) = b(k,149) - lu(k,222) * b(k,84)
116 315224452 : b(k,207) = b(k,207) - lu(k,223) * b(k,84)
117 315224452 : b(k,226) = b(k,226) - lu(k,224) * b(k,84)
118 315224452 : b(k,237) = b(k,237) - lu(k,225) * b(k,84)
119 315224452 : b(k,191) = b(k,191) - lu(k,227) * b(k,85)
120 315224452 : b(k,226) = b(k,226) - lu(k,228) * b(k,85)
121 315224452 : b(k,221) = b(k,221) - lu(k,230) * b(k,86)
122 315224452 : b(k,234) = b(k,234) - lu(k,231) * b(k,86)
123 315224452 : b(k,228) = b(k,228) - lu(k,233) * b(k,87)
124 315224452 : b(k,231) = b(k,231) - lu(k,234) * b(k,87)
125 315224452 : b(k,139) = b(k,139) - lu(k,236) * b(k,88)
126 315224452 : b(k,226) = b(k,226) - lu(k,237) * b(k,88)
127 315224452 : b(k,184) = b(k,184) - lu(k,239) * b(k,89)
128 315224452 : b(k,230) = b(k,230) - lu(k,240) * b(k,89)
129 315224452 : b(k,91) = b(k,91) - lu(k,242) * b(k,90)
130 315224452 : b(k,226) = b(k,226) - lu(k,243) * b(k,90)
131 315224452 : b(k,230) = b(k,230) - lu(k,244) * b(k,90)
132 315224452 : b(k,234) = b(k,234) - lu(k,245) * b(k,90)
133 315224452 : b(k,184) = b(k,184) - lu(k,247) * b(k,91)
134 315224452 : b(k,230) = b(k,230) - lu(k,248) * b(k,91)
135 315224452 : b(k,234) = b(k,234) - lu(k,249) * b(k,91)
136 315224452 : b(k,184) = b(k,184) - lu(k,252) * b(k,92)
137 315224452 : b(k,226) = b(k,226) - lu(k,253) * b(k,92)
138 315224452 : b(k,230) = b(k,230) - lu(k,254) * b(k,92)
139 315224452 : b(k,234) = b(k,234) - lu(k,255) * b(k,92)
140 315224452 : b(k,222) = b(k,222) - lu(k,257) * b(k,93)
141 315224452 : b(k,227) = b(k,227) - lu(k,258) * b(k,93)
142 315224452 : b(k,230) = b(k,230) - lu(k,259) * b(k,93)
143 315224452 : b(k,186) = b(k,186) - lu(k,262) * b(k,94)
144 315224452 : b(k,232) = b(k,232) - lu(k,263) * b(k,94)
145 315224452 : b(k,239) = b(k,239) - lu(k,264) * b(k,94)
146 315224452 : b(k,200) = b(k,200) - lu(k,266) * b(k,95)
147 315224452 : b(k,226) = b(k,226) - lu(k,267) * b(k,95)
148 315224452 : b(k,237) = b(k,237) - lu(k,268) * b(k,95)
149 315224452 : b(k,131) = b(k,131) - lu(k,270) * b(k,96)
150 315224452 : b(k,154) = b(k,154) - lu(k,271) * b(k,96)
151 315224452 : b(k,226) = b(k,226) - lu(k,272) * b(k,96)
152 315224452 : b(k,225) = b(k,225) - lu(k,274) * b(k,97)
153 315224452 : b(k,227) = b(k,227) - lu(k,275) * b(k,97)
154 315224452 : b(k,228) = b(k,228) - lu(k,276) * b(k,97)
155 315224452 : b(k,231) = b(k,231) - lu(k,277) * b(k,97)
156 315224452 : b(k,232) = b(k,232) - lu(k,278) * b(k,97)
157 315224452 : b(k,168) = b(k,168) - lu(k,280) * b(k,98)
158 315224452 : b(k,237) = b(k,237) - lu(k,281) * b(k,98)
159 315224452 : b(k,184) = b(k,184) - lu(k,283) * b(k,99)
160 315224452 : b(k,233) = b(k,233) - lu(k,284) * b(k,99)
161 315224452 : b(k,183) = b(k,183) - lu(k,286) * b(k,100)
162 315224452 : b(k,192) = b(k,192) - lu(k,287) * b(k,100)
163 315224452 : b(k,201) = b(k,201) - lu(k,288) * b(k,100)
164 315224452 : b(k,226) = b(k,226) - lu(k,289) * b(k,100)
165 315224452 : b(k,237) = b(k,237) - lu(k,290) * b(k,100)
166 315224452 : b(k,184) = b(k,184) - lu(k,293) * b(k,101)
167 315224452 : b(k,226) = b(k,226) - lu(k,294) * b(k,101)
168 315224452 : b(k,230) = b(k,230) - lu(k,295) * b(k,101)
169 315224452 : b(k,234) = b(k,234) - lu(k,296) * b(k,101)
170 315224452 : b(k,239) = b(k,239) - lu(k,297) * b(k,101)
171 315224452 : b(k,181) = b(k,181) - lu(k,299) * b(k,102)
172 315224452 : b(k,220) = b(k,220) - lu(k,300) * b(k,102)
173 315224452 : b(k,226) = b(k,226) - lu(k,301) * b(k,102)
174 315224452 : b(k,234) = b(k,234) - lu(k,302) * b(k,102)
175 315224452 : b(k,239) = b(k,239) - lu(k,303) * b(k,102)
176 315224452 : b(k,185) = b(k,185) - lu(k,305) * b(k,103)
177 315224452 : b(k,220) = b(k,220) - lu(k,306) * b(k,103)
178 315224452 : b(k,226) = b(k,226) - lu(k,307) * b(k,103)
179 315224452 : b(k,234) = b(k,234) - lu(k,308) * b(k,103)
180 315224452 : b(k,239) = b(k,239) - lu(k,309) * b(k,103)
181 315224452 : b(k,186) = b(k,186) - lu(k,311) * b(k,104)
182 315224452 : b(k,225) = b(k,225) - lu(k,312) * b(k,104)
183 315224452 : b(k,226) = b(k,226) - lu(k,313) * b(k,104)
184 315224452 : b(k,231) = b(k,231) - lu(k,314) * b(k,104)
185 315224452 : b(k,237) = b(k,237) - lu(k,315) * b(k,104)
186 315224452 : b(k,213) = b(k,213) - lu(k,317) * b(k,105)
187 315224452 : b(k,215) = b(k,215) - lu(k,318) * b(k,105)
188 315224452 : b(k,226) = b(k,226) - lu(k,319) * b(k,105)
189 315224452 : b(k,237) = b(k,237) - lu(k,320) * b(k,105)
190 315224452 : b(k,160) = b(k,160) - lu(k,322) * b(k,106)
191 315224452 : b(k,200) = b(k,200) - lu(k,323) * b(k,106)
192 315224452 : b(k,216) = b(k,216) - lu(k,324) * b(k,106)
193 315224452 : b(k,226) = b(k,226) - lu(k,325) * b(k,106)
194 315224452 : b(k,207) = b(k,207) - lu(k,327) * b(k,107)
195 315224452 : b(k,226) = b(k,226) - lu(k,328) * b(k,107)
196 315224452 : b(k,201) = b(k,201) - lu(k,330) * b(k,108)
197 315224452 : b(k,210) = b(k,210) - lu(k,331) * b(k,108)
198 315224452 : b(k,216) = b(k,216) - lu(k,332) * b(k,108)
199 315224452 : b(k,237) = b(k,237) - lu(k,333) * b(k,108)
200 315224452 : b(k,186) = b(k,186) - lu(k,335) * b(k,109)
201 315224452 : b(k,219) = b(k,219) - lu(k,336) * b(k,109)
202 315224452 : b(k,224) = b(k,224) - lu(k,337) * b(k,109)
203 315224452 : b(k,232) = b(k,232) - lu(k,338) * b(k,109)
204 315224452 : b(k,127) = b(k,127) - lu(k,340) * b(k,110)
205 315224452 : b(k,192) = b(k,192) - lu(k,341) * b(k,110)
206 315224452 : b(k,226) = b(k,226) - lu(k,342) * b(k,110)
207 315224452 : b(k,237) = b(k,237) - lu(k,343) * b(k,110)
208 315224452 : b(k,125) = b(k,125) - lu(k,346) * b(k,111)
209 315224452 : b(k,139) = b(k,139) - lu(k,347) * b(k,111)
210 315224452 : b(k,226) = b(k,226) - lu(k,348) * b(k,111)
211 315224452 : b(k,237) = b(k,237) - lu(k,349) * b(k,111)
212 315224452 : b(k,181) = b(k,181) - lu(k,351) * b(k,112)
213 315224452 : b(k,200) = b(k,200) - lu(k,352) * b(k,112)
214 315224452 : b(k,226) = b(k,226) - lu(k,353) * b(k,112)
215 315224452 : b(k,237) = b(k,237) - lu(k,354) * b(k,112)
216 315224452 : b(k,146) = b(k,146) - lu(k,356) * b(k,113)
217 315224452 : b(k,195) = b(k,195) - lu(k,357) * b(k,113)
218 315224452 : b(k,200) = b(k,200) - lu(k,358) * b(k,113)
219 315224452 : b(k,225) = b(k,225) - lu(k,359) * b(k,113)
220 315224452 : b(k,226) = b(k,226) - lu(k,360) * b(k,113)
221 315224452 : b(k,228) = b(k,228) - lu(k,361) * b(k,113)
222 315224452 : b(k,229) = b(k,229) - lu(k,362) * b(k,113)
223 315224452 : b(k,163) = b(k,163) - lu(k,364) * b(k,114)
224 315224452 : b(k,186) = b(k,186) - lu(k,365) * b(k,114)
225 315224452 : b(k,201) = b(k,201) - lu(k,366) * b(k,114)
226 315224452 : b(k,219) = b(k,219) - lu(k,367) * b(k,114)
227 315224452 : b(k,226) = b(k,226) - lu(k,368) * b(k,114)
228 315224452 : b(k,232) = b(k,232) - lu(k,369) * b(k,114)
229 325174377 : b(k,233) = b(k,233) - lu(k,370) * b(k,114)
230 : end do
231 9949925 : end subroutine lu_slv01
232 9949925 : subroutine lu_slv02( avec_len, lu, b )
233 : use shr_kind_mod, only : r8 => shr_kind_r8
234 : use chem_mods, only : clscnt4, nzcnt
235 : implicit none
236 : !-----------------------------------------------------------------------
237 : ! ... Dummy args
238 : !-----------------------------------------------------------------------
239 : integer, intent(in) :: avec_len
240 : real(r8), intent(in) :: lu(veclen,max(1,nzcnt))
241 : real(r8), intent(inout) :: b(veclen,clscnt4)
242 : !-----------------------------------------------------------------------
243 : ! ... Local variables
244 : !-----------------------------------------------------------------------
245 : integer :: k
246 : !-----------------------------------------------------------------------
247 : ! ... solve L * y = b
248 : !-----------------------------------------------------------------------
249 325174377 : do k = 1,avec_len
250 315224452 : b(k,201) = b(k,201) - lu(k,372) * b(k,115)
251 315224452 : b(k,220) = b(k,220) - lu(k,373) * b(k,115)
252 315224452 : b(k,226) = b(k,226) - lu(k,374) * b(k,115)
253 315224452 : b(k,234) = b(k,234) - lu(k,375) * b(k,115)
254 315224452 : b(k,237) = b(k,237) - lu(k,376) * b(k,115)
255 315224452 : b(k,238) = b(k,238) - lu(k,377) * b(k,115)
256 315224452 : b(k,239) = b(k,239) - lu(k,378) * b(k,115)
257 315224452 : b(k,192) = b(k,192) - lu(k,380) * b(k,116)
258 315224452 : b(k,201) = b(k,201) - lu(k,381) * b(k,116)
259 315224452 : b(k,210) = b(k,210) - lu(k,382) * b(k,116)
260 315224452 : b(k,216) = b(k,216) - lu(k,383) * b(k,116)
261 315224452 : b(k,237) = b(k,237) - lu(k,384) * b(k,116)
262 315224452 : b(k,220) = b(k,220) - lu(k,386) * b(k,117)
263 315224452 : b(k,221) = b(k,221) - lu(k,387) * b(k,117)
264 315224452 : b(k,226) = b(k,226) - lu(k,388) * b(k,117)
265 315224452 : b(k,230) = b(k,230) - lu(k,389) * b(k,117)
266 315224452 : b(k,234) = b(k,234) - lu(k,390) * b(k,117)
267 315224452 : b(k,194) = b(k,194) - lu(k,392) * b(k,118)
268 315224452 : b(k,197) = b(k,197) - lu(k,393) * b(k,118)
269 315224452 : b(k,226) = b(k,226) - lu(k,394) * b(k,118)
270 315224452 : b(k,228) = b(k,228) - lu(k,395) * b(k,118)
271 315224452 : b(k,231) = b(k,231) - lu(k,396) * b(k,118)
272 315224452 : b(k,206) = b(k,206) - lu(k,398) * b(k,119)
273 315224452 : b(k,216) = b(k,216) - lu(k,399) * b(k,119)
274 315224452 : b(k,226) = b(k,226) - lu(k,400) * b(k,119)
275 315224452 : b(k,229) = b(k,229) - lu(k,401) * b(k,119)
276 315224452 : b(k,239) = b(k,239) - lu(k,402) * b(k,119)
277 315224452 : b(k,162) = b(k,162) - lu(k,404) * b(k,120)
278 315224452 : b(k,180) = b(k,180) - lu(k,405) * b(k,120)
279 315224452 : b(k,226) = b(k,226) - lu(k,406) * b(k,120)
280 315224452 : b(k,228) = b(k,228) - lu(k,407) * b(k,120)
281 315224452 : b(k,237) = b(k,237) - lu(k,408) * b(k,120)
282 315224452 : b(k,185) = b(k,185) - lu(k,410) * b(k,121)
283 315224452 : b(k,195) = b(k,195) - lu(k,411) * b(k,121)
284 315224452 : b(k,226) = b(k,226) - lu(k,412) * b(k,121)
285 315224452 : b(k,237) = b(k,237) - lu(k,413) * b(k,121)
286 315224452 : b(k,239) = b(k,239) - lu(k,414) * b(k,121)
287 315224452 : b(k,130) = b(k,130) - lu(k,416) * b(k,122)
288 315224452 : b(k,135) = b(k,135) - lu(k,417) * b(k,122)
289 315224452 : b(k,192) = b(k,192) - lu(k,418) * b(k,122)
290 315224452 : b(k,226) = b(k,226) - lu(k,419) * b(k,122)
291 315224452 : b(k,237) = b(k,237) - lu(k,420) * b(k,122)
292 315224452 : b(k,137) = b(k,137) - lu(k,422) * b(k,123)
293 315224452 : b(k,192) = b(k,192) - lu(k,423) * b(k,123)
294 315224452 : b(k,210) = b(k,210) - lu(k,424) * b(k,123)
295 315224452 : b(k,226) = b(k,226) - lu(k,425) * b(k,123)
296 315224452 : b(k,237) = b(k,237) - lu(k,426) * b(k,123)
297 315224452 : b(k,180) = b(k,180) - lu(k,428) * b(k,124)
298 315224452 : b(k,202) = b(k,202) - lu(k,429) * b(k,124)
299 315224452 : b(k,226) = b(k,226) - lu(k,430) * b(k,124)
300 315224452 : b(k,231) = b(k,231) - lu(k,431) * b(k,124)
301 315224452 : b(k,237) = b(k,237) - lu(k,432) * b(k,124)
302 315224452 : b(k,139) = b(k,139) - lu(k,436) * b(k,125)
303 315224452 : b(k,226) = b(k,226) - lu(k,437) * b(k,125)
304 315224452 : b(k,227) = b(k,227) - lu(k,438) * b(k,125)
305 315224452 : b(k,228) = b(k,228) - lu(k,439) * b(k,125)
306 315224452 : b(k,237) = b(k,237) - lu(k,440) * b(k,125)
307 315224452 : b(k,183) = b(k,183) - lu(k,442) * b(k,126)
308 315224452 : b(k,227) = b(k,227) - lu(k,443) * b(k,126)
309 315224452 : b(k,228) = b(k,228) - lu(k,444) * b(k,126)
310 315224452 : b(k,229) = b(k,229) - lu(k,445) * b(k,126)
311 315224452 : b(k,237) = b(k,237) - lu(k,446) * b(k,126)
312 315224452 : b(k,192) = b(k,192) - lu(k,449) * b(k,127)
313 315224452 : b(k,226) = b(k,226) - lu(k,450) * b(k,127)
314 315224452 : b(k,227) = b(k,227) - lu(k,451) * b(k,127)
315 315224452 : b(k,228) = b(k,228) - lu(k,452) * b(k,127)
316 315224452 : b(k,237) = b(k,237) - lu(k,453) * b(k,127)
317 315224452 : b(k,226) = b(k,226) - lu(k,455) * b(k,128)
318 315224452 : b(k,229) = b(k,229) - lu(k,456) * b(k,128)
319 315224452 : b(k,233) = b(k,233) - lu(k,457) * b(k,128)
320 315224452 : b(k,238) = b(k,238) - lu(k,458) * b(k,128)
321 315224452 : b(k,239) = b(k,239) - lu(k,459) * b(k,128)
322 315224452 : b(k,202) = b(k,202) - lu(k,461) * b(k,129)
323 315224452 : b(k,222) = b(k,222) - lu(k,462) * b(k,129)
324 315224452 : b(k,229) = b(k,229) - lu(k,463) * b(k,129)
325 315224452 : b(k,237) = b(k,237) - lu(k,464) * b(k,129)
326 315224452 : b(k,165) = b(k,165) - lu(k,466) * b(k,130)
327 315224452 : b(k,237) = b(k,237) - lu(k,467) * b(k,130)
328 315224452 : b(k,154) = b(k,154) - lu(k,469) * b(k,131)
329 315224452 : b(k,228) = b(k,228) - lu(k,470) * b(k,131)
330 315224452 : b(k,235) = b(k,235) - lu(k,471) * b(k,131)
331 315224452 : b(k,222) = b(k,222) - lu(k,473) * b(k,132)
332 315224452 : b(k,225) = b(k,225) - lu(k,474) * b(k,132)
333 315224452 : b(k,226) = b(k,226) - lu(k,475) * b(k,132)
334 315224452 : b(k,228) = b(k,228) - lu(k,476) * b(k,132)
335 315224452 : b(k,237) = b(k,237) - lu(k,477) * b(k,132)
336 315224452 : b(k,239) = b(k,239) - lu(k,478) * b(k,132)
337 315224452 : b(k,220) = b(k,220) - lu(k,480) * b(k,133)
338 315224452 : b(k,221) = b(k,221) - lu(k,481) * b(k,133)
339 315224452 : b(k,226) = b(k,226) - lu(k,482) * b(k,133)
340 315224452 : b(k,230) = b(k,230) - lu(k,483) * b(k,133)
341 315224452 : b(k,234) = b(k,234) - lu(k,484) * b(k,133)
342 315224452 : b(k,239) = b(k,239) - lu(k,485) * b(k,133)
343 315224452 : b(k,193) = b(k,193) - lu(k,487) * b(k,134)
344 315224452 : b(k,194) = b(k,194) - lu(k,488) * b(k,134)
345 315224452 : b(k,195) = b(k,195) - lu(k,489) * b(k,134)
346 315224452 : b(k,226) = b(k,226) - lu(k,490) * b(k,134)
347 315224452 : b(k,229) = b(k,229) - lu(k,491) * b(k,134)
348 315224452 : b(k,237) = b(k,237) - lu(k,492) * b(k,134)
349 315224452 : b(k,165) = b(k,165) - lu(k,496) * b(k,135)
350 315224452 : b(k,192) = b(k,192) - lu(k,497) * b(k,135)
351 315224452 : b(k,226) = b(k,226) - lu(k,498) * b(k,135)
352 315224452 : b(k,227) = b(k,227) - lu(k,499) * b(k,135)
353 315224452 : b(k,228) = b(k,228) - lu(k,500) * b(k,135)
354 315224452 : b(k,237) = b(k,237) - lu(k,501) * b(k,135)
355 315224452 : b(k,171) = b(k,171) - lu(k,503) * b(k,136)
356 315224452 : b(k,176) = b(k,176) - lu(k,504) * b(k,136)
357 315224452 : b(k,177) = b(k,177) - lu(k,505) * b(k,136)
358 315224452 : b(k,189) = b(k,189) - lu(k,506) * b(k,136)
359 315224452 : b(k,222) = b(k,222) - lu(k,507) * b(k,136)
360 315224452 : b(k,232) = b(k,232) - lu(k,508) * b(k,136)
361 315224452 : b(k,192) = b(k,192) - lu(k,511) * b(k,137)
362 315224452 : b(k,210) = b(k,210) - lu(k,512) * b(k,137)
363 315224452 : b(k,226) = b(k,226) - lu(k,513) * b(k,137)
364 315224452 : b(k,227) = b(k,227) - lu(k,514) * b(k,137)
365 315224452 : b(k,228) = b(k,228) - lu(k,515) * b(k,137)
366 315224452 : b(k,237) = b(k,237) - lu(k,516) * b(k,137)
367 315224452 : b(k,160) = b(k,160) - lu(k,518) * b(k,138)
368 315224452 : b(k,181) = b(k,181) - lu(k,519) * b(k,138)
369 315224452 : b(k,216) = b(k,216) - lu(k,520) * b(k,138)
370 315224452 : b(k,226) = b(k,226) - lu(k,521) * b(k,138)
371 315224452 : b(k,154) = b(k,154) - lu(k,524) * b(k,139)
372 315224452 : b(k,226) = b(k,226) - lu(k,525) * b(k,139)
373 315224452 : b(k,227) = b(k,227) - lu(k,526) * b(k,139)
374 315224452 : b(k,228) = b(k,228) - lu(k,527) * b(k,139)
375 315224452 : b(k,237) = b(k,237) - lu(k,528) * b(k,139)
376 315224452 : b(k,189) = b(k,189) - lu(k,530) * b(k,140)
377 315224452 : b(k,222) = b(k,222) - lu(k,531) * b(k,140)
378 315224452 : b(k,227) = b(k,227) - lu(k,532) * b(k,140)
379 315224452 : b(k,230) = b(k,230) - lu(k,533) * b(k,140)
380 315224452 : b(k,216) = b(k,216) - lu(k,535) * b(k,141)
381 315224452 : b(k,217) = b(k,217) - lu(k,536) * b(k,141)
382 315224452 : b(k,225) = b(k,225) - lu(k,537) * b(k,141)
383 315224452 : b(k,226) = b(k,226) - lu(k,538) * b(k,141)
384 315224452 : b(k,228) = b(k,228) - lu(k,539) * b(k,141)
385 315224452 : b(k,229) = b(k,229) - lu(k,540) * b(k,141)
386 315224452 : b(k,238) = b(k,238) - lu(k,541) * b(k,141)
387 315224452 : b(k,172) = b(k,172) - lu(k,543) * b(k,142)
388 315224452 : b(k,200) = b(k,200) - lu(k,544) * b(k,142)
389 315224452 : b(k,205) = b(k,205) - lu(k,545) * b(k,142)
390 315224452 : b(k,226) = b(k,226) - lu(k,546) * b(k,142)
391 315224452 : b(k,229) = b(k,229) - lu(k,547) * b(k,142)
392 315224452 : b(k,237) = b(k,237) - lu(k,548) * b(k,142)
393 315224452 : b(k,239) = b(k,239) - lu(k,549) * b(k,142)
394 315224452 : b(k,170) = b(k,170) - lu(k,551) * b(k,143)
395 315224452 : b(k,183) = b(k,183) - lu(k,552) * b(k,143)
396 315224452 : b(k,201) = b(k,201) - lu(k,553) * b(k,143)
397 315224452 : b(k,226) = b(k,226) - lu(k,554) * b(k,143)
398 315224452 : b(k,229) = b(k,229) - lu(k,555) * b(k,143)
399 315224452 : b(k,235) = b(k,235) - lu(k,556) * b(k,143)
400 315224452 : b(k,237) = b(k,237) - lu(k,557) * b(k,143)
401 315224452 : b(k,180) = b(k,180) - lu(k,559) * b(k,144)
402 315224452 : b(k,202) = b(k,202) - lu(k,560) * b(k,144)
403 315224452 : b(k,204) = b(k,204) - lu(k,561) * b(k,144)
404 315224452 : b(k,205) = b(k,205) - lu(k,562) * b(k,144)
405 315224452 : b(k,226) = b(k,226) - lu(k,563) * b(k,144)
406 315224452 : b(k,231) = b(k,231) - lu(k,564) * b(k,144)
407 315224452 : b(k,237) = b(k,237) - lu(k,565) * b(k,144)
408 315224452 : b(k,169) = b(k,169) - lu(k,567) * b(k,145)
409 315224452 : b(k,221) = b(k,221) - lu(k,568) * b(k,145)
410 315224452 : b(k,223) = b(k,223) - lu(k,569) * b(k,145)
411 315224452 : b(k,225) = b(k,225) - lu(k,570) * b(k,145)
412 315224452 : b(k,228) = b(k,228) - lu(k,571) * b(k,145)
413 315224452 : b(k,231) = b(k,231) - lu(k,572) * b(k,145)
414 315224452 : b(k,232) = b(k,232) - lu(k,573) * b(k,145)
415 315224452 : b(k,195) = b(k,195) - lu(k,575) * b(k,146)
416 315224452 : b(k,200) = b(k,200) - lu(k,576) * b(k,146)
417 315224452 : b(k,204) = b(k,204) - lu(k,577) * b(k,146)
418 315224452 : b(k,227) = b(k,227) - lu(k,578) * b(k,146)
419 315224452 : b(k,228) = b(k,228) - lu(k,579) * b(k,146)
420 315224452 : b(k,229) = b(k,229) - lu(k,580) * b(k,146)
421 315224452 : b(k,237) = b(k,237) - lu(k,581) * b(k,146)
422 315224452 : b(k,220) = b(k,220) - lu(k,583) * b(k,147)
423 315224452 : b(k,226) = b(k,226) - lu(k,584) * b(k,147)
424 315224452 : b(k,232) = b(k,232) - lu(k,585) * b(k,147)
425 315224452 : b(k,234) = b(k,234) - lu(k,586) * b(k,147)
426 315224452 : b(k,237) = b(k,237) - lu(k,587) * b(k,147)
427 315224452 : b(k,239) = b(k,239) - lu(k,588) * b(k,147)
428 315224452 : b(k,202) = b(k,202) - lu(k,590) * b(k,148)
429 315224452 : b(k,204) = b(k,204) - lu(k,591) * b(k,148)
430 315224452 : b(k,205) = b(k,205) - lu(k,592) * b(k,148)
431 315224452 : b(k,226) = b(k,226) - lu(k,593) * b(k,148)
432 315224452 : b(k,228) = b(k,228) - lu(k,594) * b(k,148)
433 315224452 : b(k,229) = b(k,229) - lu(k,595) * b(k,148)
434 315224452 : b(k,231) = b(k,231) - lu(k,596) * b(k,148)
435 315224452 : b(k,237) = b(k,237) - lu(k,597) * b(k,148)
436 315224452 : b(k,182) = b(k,182) - lu(k,599) * b(k,149)
437 315224452 : b(k,201) = b(k,201) - lu(k,600) * b(k,149)
438 315224452 : b(k,237) = b(k,237) - lu(k,601) * b(k,149)
439 315224452 : b(k,220) = b(k,220) - lu(k,603) * b(k,150)
440 315224452 : b(k,221) = b(k,221) - lu(k,604) * b(k,150)
441 315224452 : b(k,226) = b(k,226) - lu(k,605) * b(k,150)
442 315224452 : b(k,230) = b(k,230) - lu(k,606) * b(k,150)
443 315224452 : b(k,234) = b(k,234) - lu(k,607) * b(k,150)
444 315224452 : b(k,237) = b(k,237) - lu(k,608) * b(k,150)
445 315224452 : b(k,238) = b(k,238) - lu(k,609) * b(k,150)
446 315224452 : b(k,239) = b(k,239) - lu(k,610) * b(k,150)
447 315224452 : b(k,205) = b(k,205) - lu(k,612) * b(k,151)
448 315224452 : b(k,215) = b(k,215) - lu(k,613) * b(k,151)
449 315224452 : b(k,217) = b(k,217) - lu(k,614) * b(k,151)
450 315224452 : b(k,225) = b(k,225) - lu(k,615) * b(k,151)
451 315224452 : b(k,226) = b(k,226) - lu(k,616) * b(k,151)
452 315224452 : b(k,228) = b(k,228) - lu(k,617) * b(k,151)
453 315224452 : b(k,229) = b(k,229) - lu(k,618) * b(k,151)
454 325174377 : b(k,237) = b(k,237) - lu(k,619) * b(k,151)
455 : end do
456 9949925 : end subroutine lu_slv02
457 9949925 : subroutine lu_slv03( avec_len, lu, b )
458 : use shr_kind_mod, only : r8 => shr_kind_r8
459 : use chem_mods, only : clscnt4, nzcnt
460 : implicit none
461 : !-----------------------------------------------------------------------
462 : ! ... Dummy args
463 : !-----------------------------------------------------------------------
464 : integer, intent(in) :: avec_len
465 : real(r8), intent(in) :: lu(veclen,max(1,nzcnt))
466 : real(r8), intent(inout) :: b(veclen,clscnt4)
467 : !-----------------------------------------------------------------------
468 : ! ... Local variables
469 : !-----------------------------------------------------------------------
470 : integer :: k
471 : !-----------------------------------------------------------------------
472 : ! ... solve L * y = b
473 : !-----------------------------------------------------------------------
474 325174377 : do k = 1,avec_len
475 315224452 : b(k,216) = b(k,216) - lu(k,621) * b(k,152)
476 315224452 : b(k,217) = b(k,217) - lu(k,622) * b(k,152)
477 315224452 : b(k,226) = b(k,226) - lu(k,623) * b(k,152)
478 315224452 : b(k,229) = b(k,229) - lu(k,624) * b(k,152)
479 315224452 : b(k,238) = b(k,238) - lu(k,625) * b(k,152)
480 315224452 : b(k,239) = b(k,239) - lu(k,626) * b(k,152)
481 315224452 : b(k,194) = b(k,194) - lu(k,628) * b(k,153)
482 315224452 : b(k,226) = b(k,226) - lu(k,629) * b(k,153)
483 315224452 : b(k,228) = b(k,228) - lu(k,630) * b(k,153)
484 315224452 : b(k,231) = b(k,231) - lu(k,631) * b(k,153)
485 315224452 : b(k,237) = b(k,237) - lu(k,632) * b(k,153)
486 315224452 : b(k,226) = b(k,226) - lu(k,636) * b(k,154)
487 315224452 : b(k,227) = b(k,227) - lu(k,637) * b(k,154)
488 315224452 : b(k,228) = b(k,228) - lu(k,638) * b(k,154)
489 315224452 : b(k,235) = b(k,235) - lu(k,639) * b(k,154)
490 315224452 : b(k,237) = b(k,237) - lu(k,640) * b(k,154)
491 315224452 : b(k,160) = b(k,160) - lu(k,643) * b(k,155)
492 315224452 : b(k,181) = b(k,181) - lu(k,644) * b(k,155)
493 315224452 : b(k,195) = b(k,195) - lu(k,645) * b(k,155)
494 315224452 : b(k,200) = b(k,200) - lu(k,646) * b(k,155)
495 315224452 : b(k,216) = b(k,216) - lu(k,647) * b(k,155)
496 315224452 : b(k,226) = b(k,226) - lu(k,648) * b(k,155)
497 315224452 : b(k,228) = b(k,228) - lu(k,649) * b(k,155)
498 315224452 : b(k,229) = b(k,229) - lu(k,650) * b(k,155)
499 315224452 : b(k,237) = b(k,237) - lu(k,651) * b(k,155)
500 315224452 : b(k,195) = b(k,195) - lu(k,653) * b(k,156)
501 315224452 : b(k,198) = b(k,198) - lu(k,654) * b(k,156)
502 315224452 : b(k,201) = b(k,201) - lu(k,655) * b(k,156)
503 315224452 : b(k,202) = b(k,202) - lu(k,656) * b(k,156)
504 315224452 : b(k,203) = b(k,203) - lu(k,657) * b(k,156)
505 315224452 : b(k,217) = b(k,217) - lu(k,658) * b(k,156)
506 315224452 : b(k,226) = b(k,226) - lu(k,659) * b(k,156)
507 315224452 : b(k,229) = b(k,229) - lu(k,660) * b(k,156)
508 315224452 : b(k,237) = b(k,237) - lu(k,661) * b(k,156)
509 315224452 : b(k,165) = b(k,165) - lu(k,666) * b(k,157)
510 315224452 : b(k,166) = b(k,166) - lu(k,667) * b(k,157)
511 315224452 : b(k,168) = b(k,168) - lu(k,668) * b(k,157)
512 315224452 : b(k,182) = b(k,182) - lu(k,669) * b(k,157)
513 315224452 : b(k,192) = b(k,192) - lu(k,670) * b(k,157)
514 315224452 : b(k,201) = b(k,201) - lu(k,671) * b(k,157)
515 315224452 : b(k,210) = b(k,210) - lu(k,672) * b(k,157)
516 315224452 : b(k,226) = b(k,226) - lu(k,673) * b(k,157)
517 315224452 : b(k,237) = b(k,237) - lu(k,674) * b(k,157)
518 315224452 : b(k,171) = b(k,171) - lu(k,677) * b(k,158)
519 315224452 : b(k,176) = b(k,176) - lu(k,678) * b(k,158)
520 315224452 : b(k,177) = b(k,177) - lu(k,679) * b(k,158)
521 315224452 : b(k,178) = b(k,178) - lu(k,680) * b(k,158)
522 315224452 : b(k,189) = b(k,189) - lu(k,681) * b(k,158)
523 315224452 : b(k,222) = b(k,222) - lu(k,682) * b(k,158)
524 315224452 : b(k,227) = b(k,227) - lu(k,683) * b(k,158)
525 315224452 : b(k,230) = b(k,230) - lu(k,684) * b(k,158)
526 315224452 : b(k,232) = b(k,232) - lu(k,685) * b(k,158)
527 315224452 : b(k,160) = b(k,160) - lu(k,688) * b(k,159)
528 315224452 : b(k,181) = b(k,181) - lu(k,689) * b(k,159)
529 315224452 : b(k,191) = b(k,191) - lu(k,690) * b(k,159)
530 315224452 : b(k,195) = b(k,195) - lu(k,691) * b(k,159)
531 315224452 : b(k,200) = b(k,200) - lu(k,692) * b(k,159)
532 315224452 : b(k,216) = b(k,216) - lu(k,693) * b(k,159)
533 315224452 : b(k,226) = b(k,226) - lu(k,694) * b(k,159)
534 315224452 : b(k,229) = b(k,229) - lu(k,695) * b(k,159)
535 315224452 : b(k,237) = b(k,237) - lu(k,696) * b(k,159)
536 315224452 : b(k,200) = b(k,200) - lu(k,699) * b(k,160)
537 315224452 : b(k,216) = b(k,216) - lu(k,700) * b(k,160)
538 315224452 : b(k,226) = b(k,226) - lu(k,701) * b(k,160)
539 315224452 : b(k,227) = b(k,227) - lu(k,702) * b(k,160)
540 315224452 : b(k,228) = b(k,228) - lu(k,703) * b(k,160)
541 315224452 : b(k,237) = b(k,237) - lu(k,704) * b(k,160)
542 315224452 : b(k,217) = b(k,217) - lu(k,706) * b(k,161)
543 315224452 : b(k,226) = b(k,226) - lu(k,707) * b(k,161)
544 315224452 : b(k,238) = b(k,238) - lu(k,708) * b(k,161)
545 315224452 : b(k,239) = b(k,239) - lu(k,709) * b(k,161)
546 315224452 : b(k,207) = b(k,207) - lu(k,712) * b(k,162)
547 315224452 : b(k,209) = b(k,209) - lu(k,713) * b(k,162)
548 315224452 : b(k,214) = b(k,214) - lu(k,714) * b(k,162)
549 315224452 : b(k,226) = b(k,226) - lu(k,715) * b(k,162)
550 315224452 : b(k,229) = b(k,229) - lu(k,716) * b(k,162)
551 315224452 : b(k,237) = b(k,237) - lu(k,717) * b(k,162)
552 315224452 : b(k,219) = b(k,219) - lu(k,719) * b(k,163)
553 315224452 : b(k,222) = b(k,222) - lu(k,720) * b(k,163)
554 315224452 : b(k,226) = b(k,226) - lu(k,721) * b(k,163)
555 315224452 : b(k,232) = b(k,232) - lu(k,722) * b(k,163)
556 315224452 : b(k,233) = b(k,233) - lu(k,723) * b(k,163)
557 315224452 : b(k,235) = b(k,235) - lu(k,724) * b(k,163)
558 315224452 : b(k,165) = b(k,165) - lu(k,730) * b(k,164)
559 315224452 : b(k,167) = b(k,167) - lu(k,731) * b(k,164)
560 315224452 : b(k,168) = b(k,168) - lu(k,732) * b(k,164)
561 315224452 : b(k,182) = b(k,182) - lu(k,733) * b(k,164)
562 315224452 : b(k,192) = b(k,192) - lu(k,734) * b(k,164)
563 315224452 : b(k,201) = b(k,201) - lu(k,735) * b(k,164)
564 315224452 : b(k,210) = b(k,210) - lu(k,736) * b(k,164)
565 315224452 : b(k,216) = b(k,216) - lu(k,737) * b(k,164)
566 315224452 : b(k,226) = b(k,226) - lu(k,738) * b(k,164)
567 315224452 : b(k,237) = b(k,237) - lu(k,739) * b(k,164)
568 315224452 : b(k,192) = b(k,192) - lu(k,741) * b(k,165)
569 315224452 : b(k,201) = b(k,201) - lu(k,742) * b(k,165)
570 315224452 : b(k,227) = b(k,227) - lu(k,743) * b(k,165)
571 315224452 : b(k,228) = b(k,228) - lu(k,744) * b(k,165)
572 315224452 : b(k,237) = b(k,237) - lu(k,745) * b(k,165)
573 315224452 : b(k,168) = b(k,168) - lu(k,752) * b(k,166)
574 315224452 : b(k,182) = b(k,182) - lu(k,753) * b(k,166)
575 315224452 : b(k,192) = b(k,192) - lu(k,754) * b(k,166)
576 315224452 : b(k,201) = b(k,201) - lu(k,755) * b(k,166)
577 315224452 : b(k,210) = b(k,210) - lu(k,756) * b(k,166)
578 315224452 : b(k,226) = b(k,226) - lu(k,757) * b(k,166)
579 315224452 : b(k,227) = b(k,227) - lu(k,758) * b(k,166)
580 315224452 : b(k,228) = b(k,228) - lu(k,759) * b(k,166)
581 315224452 : b(k,237) = b(k,237) - lu(k,760) * b(k,166)
582 315224452 : b(k,168) = b(k,168) - lu(k,768) * b(k,167)
583 315224452 : b(k,182) = b(k,182) - lu(k,769) * b(k,167)
584 315224452 : b(k,192) = b(k,192) - lu(k,770) * b(k,167)
585 315224452 : b(k,201) = b(k,201) - lu(k,771) * b(k,167)
586 315224452 : b(k,210) = b(k,210) - lu(k,772) * b(k,167)
587 315224452 : b(k,216) = b(k,216) - lu(k,773) * b(k,167)
588 315224452 : b(k,226) = b(k,226) - lu(k,774) * b(k,167)
589 315224452 : b(k,227) = b(k,227) - lu(k,775) * b(k,167)
590 315224452 : b(k,228) = b(k,228) - lu(k,776) * b(k,167)
591 315224452 : b(k,237) = b(k,237) - lu(k,777) * b(k,167)
592 315224452 : b(k,201) = b(k,201) - lu(k,779) * b(k,168)
593 315224452 : b(k,210) = b(k,210) - lu(k,780) * b(k,168)
594 315224452 : b(k,226) = b(k,226) - lu(k,781) * b(k,168)
595 315224452 : b(k,227) = b(k,227) - lu(k,782) * b(k,168)
596 315224452 : b(k,228) = b(k,228) - lu(k,783) * b(k,168)
597 315224452 : b(k,237) = b(k,237) - lu(k,784) * b(k,168)
598 315224452 : b(k,238) = b(k,238) - lu(k,785) * b(k,168)
599 315224452 : b(k,220) = b(k,220) - lu(k,788) * b(k,169)
600 315224452 : b(k,221) = b(k,221) - lu(k,789) * b(k,169)
601 315224452 : b(k,223) = b(k,223) - lu(k,790) * b(k,169)
602 315224452 : b(k,226) = b(k,226) - lu(k,791) * b(k,169)
603 315224452 : b(k,232) = b(k,232) - lu(k,792) * b(k,169)
604 315224452 : b(k,234) = b(k,234) - lu(k,793) * b(k,169)
605 315224452 : b(k,239) = b(k,239) - lu(k,794) * b(k,169)
606 315224452 : b(k,202) = b(k,202) - lu(k,798) * b(k,170)
607 315224452 : b(k,222) = b(k,222) - lu(k,799) * b(k,170)
608 315224452 : b(k,226) = b(k,226) - lu(k,800) * b(k,170)
609 315224452 : b(k,227) = b(k,227) - lu(k,801) * b(k,170)
610 315224452 : b(k,228) = b(k,228) - lu(k,802) * b(k,170)
611 315224452 : b(k,229) = b(k,229) - lu(k,803) * b(k,170)
612 315224452 : b(k,237) = b(k,237) - lu(k,804) * b(k,170)
613 315224452 : b(k,176) = b(k,176) - lu(k,806) * b(k,171)
614 315224452 : b(k,177) = b(k,177) - lu(k,807) * b(k,171)
615 315224452 : b(k,189) = b(k,189) - lu(k,808) * b(k,171)
616 315224452 : b(k,201) = b(k,201) - lu(k,809) * b(k,171)
617 315224452 : b(k,217) = b(k,217) - lu(k,810) * b(k,171)
618 315224452 : b(k,222) = b(k,222) - lu(k,811) * b(k,171)
619 315224452 : b(k,232) = b(k,232) - lu(k,812) * b(k,171)
620 315224452 : b(k,200) = b(k,200) - lu(k,815) * b(k,172)
621 315224452 : b(k,205) = b(k,205) - lu(k,816) * b(k,172)
622 315224452 : b(k,222) = b(k,222) - lu(k,817) * b(k,172)
623 315224452 : b(k,226) = b(k,226) - lu(k,818) * b(k,172)
624 315224452 : b(k,227) = b(k,227) - lu(k,819) * b(k,172)
625 315224452 : b(k,228) = b(k,228) - lu(k,820) * b(k,172)
626 315224452 : b(k,229) = b(k,229) - lu(k,821) * b(k,172)
627 315224452 : b(k,237) = b(k,237) - lu(k,822) * b(k,172)
628 315224452 : b(k,239) = b(k,239) - lu(k,823) * b(k,172)
629 315224452 : b(k,226) = b(k,226) - lu(k,825) * b(k,173)
630 315224452 : b(k,229) = b(k,229) - lu(k,826) * b(k,173)
631 315224452 : b(k,237) = b(k,237) - lu(k,827) * b(k,173)
632 315224452 : b(k,221) = b(k,221) - lu(k,829) * b(k,174)
633 315224452 : b(k,223) = b(k,223) - lu(k,830) * b(k,174)
634 315224452 : b(k,226) = b(k,226) - lu(k,831) * b(k,174)
635 315224452 : b(k,230) = b(k,230) - lu(k,832) * b(k,174)
636 315224452 : b(k,232) = b(k,232) - lu(k,833) * b(k,174)
637 315224452 : b(k,233) = b(k,233) - lu(k,834) * b(k,174)
638 315224452 : b(k,239) = b(k,239) - lu(k,835) * b(k,174)
639 315224452 : b(k,220) = b(k,220) - lu(k,838) * b(k,175)
640 315224452 : b(k,224) = b(k,224) - lu(k,839) * b(k,175)
641 315224452 : b(k,226) = b(k,226) - lu(k,840) * b(k,175)
642 315224452 : b(k,232) = b(k,232) - lu(k,841) * b(k,175)
643 315224452 : b(k,234) = b(k,234) - lu(k,842) * b(k,175)
644 315224452 : b(k,239) = b(k,239) - lu(k,843) * b(k,175)
645 315224452 : b(k,178) = b(k,178) - lu(k,846) * b(k,176)
646 315224452 : b(k,189) = b(k,189) - lu(k,847) * b(k,176)
647 315224452 : b(k,222) = b(k,222) - lu(k,848) * b(k,176)
648 315224452 : b(k,227) = b(k,227) - lu(k,849) * b(k,176)
649 315224452 : b(k,230) = b(k,230) - lu(k,850) * b(k,176)
650 315224452 : b(k,232) = b(k,232) - lu(k,851) * b(k,176)
651 315224452 : b(k,178) = b(k,178) - lu(k,854) * b(k,177)
652 315224452 : b(k,189) = b(k,189) - lu(k,855) * b(k,177)
653 315224452 : b(k,222) = b(k,222) - lu(k,856) * b(k,177)
654 315224452 : b(k,227) = b(k,227) - lu(k,857) * b(k,177)
655 315224452 : b(k,230) = b(k,230) - lu(k,858) * b(k,177)
656 315224452 : b(k,232) = b(k,232) - lu(k,859) * b(k,177)
657 315224452 : b(k,189) = b(k,189) - lu(k,866) * b(k,178)
658 315224452 : b(k,201) = b(k,201) - lu(k,867) * b(k,178)
659 315224452 : b(k,217) = b(k,217) - lu(k,868) * b(k,178)
660 315224452 : b(k,222) = b(k,222) - lu(k,869) * b(k,178)
661 315224452 : b(k,227) = b(k,227) - lu(k,870) * b(k,178)
662 315224452 : b(k,230) = b(k,230) - lu(k,871) * b(k,178)
663 315224452 : b(k,232) = b(k,232) - lu(k,872) * b(k,178)
664 315224452 : b(k,183) = b(k,183) - lu(k,877) * b(k,179)
665 315224452 : b(k,199) = b(k,199) - lu(k,878) * b(k,179)
666 315224452 : b(k,201) = b(k,201) - lu(k,879) * b(k,179)
667 315224452 : b(k,208) = b(k,208) - lu(k,880) * b(k,179)
668 315224452 : b(k,209) = b(k,209) - lu(k,881) * b(k,179)
669 315224452 : b(k,211) = b(k,211) - lu(k,882) * b(k,179)
670 315224452 : b(k,212) = b(k,212) - lu(k,883) * b(k,179)
671 315224452 : b(k,214) = b(k,214) - lu(k,884) * b(k,179)
672 315224452 : b(k,216) = b(k,216) - lu(k,885) * b(k,179)
673 315224452 : b(k,225) = b(k,225) - lu(k,886) * b(k,179)
674 315224452 : b(k,226) = b(k,226) - lu(k,887) * b(k,179)
675 315224452 : b(k,229) = b(k,229) - lu(k,888) * b(k,179)
676 315224452 : b(k,235) = b(k,235) - lu(k,889) * b(k,179)
677 315224452 : b(k,237) = b(k,237) - lu(k,890) * b(k,179)
678 325174377 : b(k,238) = b(k,238) - lu(k,891) * b(k,179)
679 : end do
680 9949925 : end subroutine lu_slv03
681 9949925 : subroutine lu_slv04( avec_len, lu, b )
682 : use shr_kind_mod, only : r8 => shr_kind_r8
683 : use chem_mods, only : clscnt4, nzcnt
684 : implicit none
685 : !-----------------------------------------------------------------------
686 : ! ... Dummy args
687 : !-----------------------------------------------------------------------
688 : integer, intent(in) :: avec_len
689 : real(r8), intent(in) :: lu(veclen,max(1,nzcnt))
690 : real(r8), intent(inout) :: b(veclen,clscnt4)
691 : !-----------------------------------------------------------------------
692 : ! ... Local variables
693 : !-----------------------------------------------------------------------
694 : integer :: k
695 : !-----------------------------------------------------------------------
696 : ! ... solve L * y = b
697 : !-----------------------------------------------------------------------
698 325174377 : do k = 1,avec_len
699 315224452 : b(k,210) = b(k,210) - lu(k,893) * b(k,180)
700 315224452 : b(k,216) = b(k,216) - lu(k,894) * b(k,180)
701 315224452 : b(k,226) = b(k,226) - lu(k,895) * b(k,180)
702 315224452 : b(k,228) = b(k,228) - lu(k,896) * b(k,180)
703 315224452 : b(k,229) = b(k,229) - lu(k,897) * b(k,180)
704 315224452 : b(k,200) = b(k,200) - lu(k,902) * b(k,181)
705 315224452 : b(k,222) = b(k,222) - lu(k,903) * b(k,181)
706 315224452 : b(k,226) = b(k,226) - lu(k,904) * b(k,181)
707 315224452 : b(k,227) = b(k,227) - lu(k,905) * b(k,181)
708 315224452 : b(k,228) = b(k,228) - lu(k,906) * b(k,181)
709 315224452 : b(k,229) = b(k,229) - lu(k,907) * b(k,181)
710 315224452 : b(k,237) = b(k,237) - lu(k,908) * b(k,181)
711 315224452 : b(k,238) = b(k,238) - lu(k,909) * b(k,181)
712 315224452 : b(k,192) = b(k,192) - lu(k,911) * b(k,182)
713 315224452 : b(k,201) = b(k,201) - lu(k,912) * b(k,182)
714 315224452 : b(k,210) = b(k,210) - lu(k,913) * b(k,182)
715 315224452 : b(k,226) = b(k,226) - lu(k,914) * b(k,182)
716 315224452 : b(k,227) = b(k,227) - lu(k,915) * b(k,182)
717 315224452 : b(k,228) = b(k,228) - lu(k,916) * b(k,182)
718 315224452 : b(k,237) = b(k,237) - lu(k,917) * b(k,182)
719 315224452 : b(k,238) = b(k,238) - lu(k,918) * b(k,182)
720 315224452 : b(k,217) = b(k,217) - lu(k,920) * b(k,183)
721 315224452 : b(k,226) = b(k,226) - lu(k,921) * b(k,183)
722 315224452 : b(k,237) = b(k,237) - lu(k,922) * b(k,183)
723 315224452 : b(k,239) = b(k,239) - lu(k,923) * b(k,183)
724 315224452 : b(k,218) = b(k,218) - lu(k,926) * b(k,184)
725 315224452 : b(k,225) = b(k,225) - lu(k,927) * b(k,184)
726 315224452 : b(k,226) = b(k,226) - lu(k,928) * b(k,184)
727 315224452 : b(k,231) = b(k,231) - lu(k,929) * b(k,184)
728 315224452 : b(k,233) = b(k,233) - lu(k,930) * b(k,184)
729 315224452 : b(k,236) = b(k,236) - lu(k,931) * b(k,184)
730 315224452 : b(k,238) = b(k,238) - lu(k,932) * b(k,184)
731 315224452 : b(k,239) = b(k,239) - lu(k,933) * b(k,184)
732 315224452 : b(k,195) = b(k,195) - lu(k,936) * b(k,185)
733 315224452 : b(k,200) = b(k,200) - lu(k,937) * b(k,185)
734 315224452 : b(k,222) = b(k,222) - lu(k,938) * b(k,185)
735 315224452 : b(k,226) = b(k,226) - lu(k,939) * b(k,185)
736 315224452 : b(k,227) = b(k,227) - lu(k,940) * b(k,185)
737 315224452 : b(k,228) = b(k,228) - lu(k,941) * b(k,185)
738 315224452 : b(k,229) = b(k,229) - lu(k,942) * b(k,185)
739 315224452 : b(k,237) = b(k,237) - lu(k,943) * b(k,185)
740 315224452 : b(k,238) = b(k,238) - lu(k,944) * b(k,185)
741 315224452 : b(k,239) = b(k,239) - lu(k,945) * b(k,185)
742 315224452 : b(k,219) = b(k,219) - lu(k,948) * b(k,186)
743 315224452 : b(k,226) = b(k,226) - lu(k,949) * b(k,186)
744 315224452 : b(k,232) = b(k,232) - lu(k,950) * b(k,186)
745 315224452 : b(k,237) = b(k,237) - lu(k,951) * b(k,186)
746 315224452 : b(k,239) = b(k,239) - lu(k,952) * b(k,186)
747 315224452 : b(k,220) = b(k,220) - lu(k,956) * b(k,187)
748 315224452 : b(k,224) = b(k,224) - lu(k,957) * b(k,187)
749 315224452 : b(k,225) = b(k,225) - lu(k,958) * b(k,187)
750 315224452 : b(k,226) = b(k,226) - lu(k,959) * b(k,187)
751 315224452 : b(k,228) = b(k,228) - lu(k,960) * b(k,187)
752 315224452 : b(k,231) = b(k,231) - lu(k,961) * b(k,187)
753 315224452 : b(k,232) = b(k,232) - lu(k,962) * b(k,187)
754 315224452 : b(k,234) = b(k,234) - lu(k,963) * b(k,187)
755 315224452 : b(k,239) = b(k,239) - lu(k,964) * b(k,187)
756 315224452 : b(k,191) = b(k,191) - lu(k,975) * b(k,188)
757 315224452 : b(k,192) = b(k,192) - lu(k,976) * b(k,188)
758 315224452 : b(k,193) = b(k,193) - lu(k,977) * b(k,188)
759 315224452 : b(k,194) = b(k,194) - lu(k,978) * b(k,188)
760 315224452 : b(k,195) = b(k,195) - lu(k,979) * b(k,188)
761 315224452 : b(k,197) = b(k,197) - lu(k,980) * b(k,188)
762 315224452 : b(k,198) = b(k,198) - lu(k,981) * b(k,188)
763 315224452 : b(k,201) = b(k,201) - lu(k,982) * b(k,188)
764 315224452 : b(k,206) = b(k,206) - lu(k,983) * b(k,188)
765 315224452 : b(k,210) = b(k,210) - lu(k,984) * b(k,188)
766 315224452 : b(k,216) = b(k,216) - lu(k,985) * b(k,188)
767 315224452 : b(k,217) = b(k,217) - lu(k,986) * b(k,188)
768 315224452 : b(k,225) = b(k,225) - lu(k,987) * b(k,188)
769 315224452 : b(k,226) = b(k,226) - lu(k,988) * b(k,188)
770 315224452 : b(k,229) = b(k,229) - lu(k,989) * b(k,188)
771 315224452 : b(k,235) = b(k,235) - lu(k,990) * b(k,188)
772 315224452 : b(k,237) = b(k,237) - lu(k,991) * b(k,188)
773 315224452 : b(k,239) = b(k,239) - lu(k,992) * b(k,188)
774 315224452 : b(k,201) = b(k,201) - lu(k,1000) * b(k,189)
775 315224452 : b(k,217) = b(k,217) - lu(k,1001) * b(k,189)
776 315224452 : b(k,222) = b(k,222) - lu(k,1002) * b(k,189)
777 315224452 : b(k,226) = b(k,226) - lu(k,1003) * b(k,189)
778 315224452 : b(k,227) = b(k,227) - lu(k,1004) * b(k,189)
779 315224452 : b(k,228) = b(k,228) - lu(k,1005) * b(k,189)
780 315224452 : b(k,230) = b(k,230) - lu(k,1006) * b(k,189)
781 315224452 : b(k,232) = b(k,232) - lu(k,1007) * b(k,189)
782 315224452 : b(k,233) = b(k,233) - lu(k,1008) * b(k,189)
783 315224452 : b(k,191) = b(k,191) - lu(k,1019) * b(k,190)
784 315224452 : b(k,192) = b(k,192) - lu(k,1020) * b(k,190)
785 315224452 : b(k,193) = b(k,193) - lu(k,1021) * b(k,190)
786 315224452 : b(k,194) = b(k,194) - lu(k,1022) * b(k,190)
787 315224452 : b(k,195) = b(k,195) - lu(k,1023) * b(k,190)
788 315224452 : b(k,197) = b(k,197) - lu(k,1024) * b(k,190)
789 315224452 : b(k,198) = b(k,198) - lu(k,1025) * b(k,190)
790 315224452 : b(k,201) = b(k,201) - lu(k,1026) * b(k,190)
791 315224452 : b(k,206) = b(k,206) - lu(k,1027) * b(k,190)
792 315224452 : b(k,210) = b(k,210) - lu(k,1028) * b(k,190)
793 315224452 : b(k,216) = b(k,216) - lu(k,1029) * b(k,190)
794 315224452 : b(k,217) = b(k,217) - lu(k,1030) * b(k,190)
795 315224452 : b(k,225) = b(k,225) - lu(k,1031) * b(k,190)
796 315224452 : b(k,226) = b(k,226) - lu(k,1032) * b(k,190)
797 315224452 : b(k,229) = b(k,229) - lu(k,1033) * b(k,190)
798 315224452 : b(k,235) = b(k,235) - lu(k,1034) * b(k,190)
799 315224452 : b(k,237) = b(k,237) - lu(k,1035) * b(k,190)
800 315224452 : b(k,239) = b(k,239) - lu(k,1036) * b(k,190)
801 315224452 : b(k,195) = b(k,195) - lu(k,1043) * b(k,191)
802 315224452 : b(k,200) = b(k,200) - lu(k,1044) * b(k,191)
803 315224452 : b(k,216) = b(k,216) - lu(k,1045) * b(k,191)
804 315224452 : b(k,222) = b(k,222) - lu(k,1046) * b(k,191)
805 315224452 : b(k,226) = b(k,226) - lu(k,1047) * b(k,191)
806 315224452 : b(k,227) = b(k,227) - lu(k,1048) * b(k,191)
807 315224452 : b(k,228) = b(k,228) - lu(k,1049) * b(k,191)
808 315224452 : b(k,229) = b(k,229) - lu(k,1050) * b(k,191)
809 315224452 : b(k,237) = b(k,237) - lu(k,1051) * b(k,191)
810 315224452 : b(k,238) = b(k,238) - lu(k,1052) * b(k,191)
811 315224452 : b(k,201) = b(k,201) - lu(k,1055) * b(k,192)
812 315224452 : b(k,217) = b(k,217) - lu(k,1056) * b(k,192)
813 315224452 : b(k,226) = b(k,226) - lu(k,1057) * b(k,192)
814 315224452 : b(k,237) = b(k,237) - lu(k,1058) * b(k,192)
815 315224452 : b(k,194) = b(k,194) - lu(k,1063) * b(k,193)
816 315224452 : b(k,195) = b(k,195) - lu(k,1064) * b(k,193)
817 315224452 : b(k,226) = b(k,226) - lu(k,1065) * b(k,193)
818 315224452 : b(k,227) = b(k,227) - lu(k,1066) * b(k,193)
819 315224452 : b(k,228) = b(k,228) - lu(k,1067) * b(k,193)
820 315224452 : b(k,229) = b(k,229) - lu(k,1068) * b(k,193)
821 315224452 : b(k,231) = b(k,231) - lu(k,1069) * b(k,193)
822 315224452 : b(k,237) = b(k,237) - lu(k,1070) * b(k,193)
823 315224452 : b(k,238) = b(k,238) - lu(k,1071) * b(k,193)
824 315224452 : b(k,197) = b(k,197) - lu(k,1073) * b(k,194)
825 315224452 : b(k,198) = b(k,198) - lu(k,1074) * b(k,194)
826 315224452 : b(k,201) = b(k,201) - lu(k,1075) * b(k,194)
827 315224452 : b(k,203) = b(k,203) - lu(k,1076) * b(k,194)
828 315224452 : b(k,225) = b(k,225) - lu(k,1077) * b(k,194)
829 315224452 : b(k,226) = b(k,226) - lu(k,1078) * b(k,194)
830 315224452 : b(k,237) = b(k,237) - lu(k,1079) * b(k,194)
831 315224452 : b(k,206) = b(k,206) - lu(k,1081) * b(k,195)
832 315224452 : b(k,216) = b(k,216) - lu(k,1082) * b(k,195)
833 315224452 : b(k,226) = b(k,226) - lu(k,1083) * b(k,195)
834 315224452 : b(k,238) = b(k,238) - lu(k,1084) * b(k,195)
835 315224452 : b(k,239) = b(k,239) - lu(k,1085) * b(k,195)
836 315224452 : b(k,201) = b(k,201) - lu(k,1091) * b(k,196)
837 315224452 : b(k,210) = b(k,210) - lu(k,1092) * b(k,196)
838 315224452 : b(k,216) = b(k,216) - lu(k,1093) * b(k,196)
839 315224452 : b(k,217) = b(k,217) - lu(k,1094) * b(k,196)
840 315224452 : b(k,226) = b(k,226) - lu(k,1095) * b(k,196)
841 315224452 : b(k,227) = b(k,227) - lu(k,1096) * b(k,196)
842 315224452 : b(k,228) = b(k,228) - lu(k,1097) * b(k,196)
843 315224452 : b(k,229) = b(k,229) - lu(k,1098) * b(k,196)
844 315224452 : b(k,231) = b(k,231) - lu(k,1099) * b(k,196)
845 315224452 : b(k,237) = b(k,237) - lu(k,1100) * b(k,196)
846 315224452 : b(k,238) = b(k,238) - lu(k,1101) * b(k,196)
847 315224452 : b(k,198) = b(k,198) - lu(k,1107) * b(k,197)
848 315224452 : b(k,201) = b(k,201) - lu(k,1108) * b(k,197)
849 315224452 : b(k,203) = b(k,203) - lu(k,1109) * b(k,197)
850 315224452 : b(k,225) = b(k,225) - lu(k,1110) * b(k,197)
851 315224452 : b(k,226) = b(k,226) - lu(k,1111) * b(k,197)
852 315224452 : b(k,227) = b(k,227) - lu(k,1112) * b(k,197)
853 315224452 : b(k,228) = b(k,228) - lu(k,1113) * b(k,197)
854 315224452 : b(k,229) = b(k,229) - lu(k,1114) * b(k,197)
855 315224452 : b(k,231) = b(k,231) - lu(k,1115) * b(k,197)
856 315224452 : b(k,237) = b(k,237) - lu(k,1116) * b(k,197)
857 315224452 : b(k,238) = b(k,238) - lu(k,1117) * b(k,197)
858 315224452 : b(k,201) = b(k,201) - lu(k,1120) * b(k,198)
859 315224452 : b(k,206) = b(k,206) - lu(k,1121) * b(k,198)
860 315224452 : b(k,216) = b(k,216) - lu(k,1122) * b(k,198)
861 315224452 : b(k,217) = b(k,217) - lu(k,1123) * b(k,198)
862 315224452 : b(k,226) = b(k,226) - lu(k,1124) * b(k,198)
863 315224452 : b(k,229) = b(k,229) - lu(k,1125) * b(k,198)
864 315224452 : b(k,237) = b(k,237) - lu(k,1126) * b(k,198)
865 315224452 : b(k,238) = b(k,238) - lu(k,1127) * b(k,198)
866 315224452 : b(k,239) = b(k,239) - lu(k,1128) * b(k,198)
867 315224452 : b(k,200) = b(k,200) - lu(k,1134) * b(k,199)
868 315224452 : b(k,201) = b(k,201) - lu(k,1135) * b(k,199)
869 315224452 : b(k,205) = b(k,205) - lu(k,1136) * b(k,199)
870 315224452 : b(k,210) = b(k,210) - lu(k,1137) * b(k,199)
871 315224452 : b(k,216) = b(k,216) - lu(k,1138) * b(k,199)
872 315224452 : b(k,217) = b(k,217) - lu(k,1139) * b(k,199)
873 315224452 : b(k,222) = b(k,222) - lu(k,1140) * b(k,199)
874 315224452 : b(k,225) = b(k,225) - lu(k,1141) * b(k,199)
875 315224452 : b(k,226) = b(k,226) - lu(k,1142) * b(k,199)
876 315224452 : b(k,227) = b(k,227) - lu(k,1143) * b(k,199)
877 315224452 : b(k,228) = b(k,228) - lu(k,1144) * b(k,199)
878 315224452 : b(k,229) = b(k,229) - lu(k,1145) * b(k,199)
879 315224452 : b(k,235) = b(k,235) - lu(k,1146) * b(k,199)
880 315224452 : b(k,236) = b(k,236) - lu(k,1147) * b(k,199)
881 315224452 : b(k,237) = b(k,237) - lu(k,1148) * b(k,199)
882 315224452 : b(k,238) = b(k,238) - lu(k,1149) * b(k,199)
883 315224452 : b(k,239) = b(k,239) - lu(k,1150) * b(k,199)
884 315224452 : b(k,201) = b(k,201) - lu(k,1152) * b(k,200)
885 315224452 : b(k,216) = b(k,216) - lu(k,1153) * b(k,200)
886 315224452 : b(k,225) = b(k,225) - lu(k,1154) * b(k,200)
887 315224452 : b(k,226) = b(k,226) - lu(k,1155) * b(k,200)
888 315224452 : b(k,231) = b(k,231) - lu(k,1156) * b(k,200)
889 315224452 : b(k,237) = b(k,237) - lu(k,1157) * b(k,200)
890 315224452 : b(k,238) = b(k,238) - lu(k,1158) * b(k,200)
891 315224452 : b(k,239) = b(k,239) - lu(k,1159) * b(k,200)
892 315224452 : b(k,217) = b(k,217) - lu(k,1161) * b(k,201)
893 315224452 : b(k,226) = b(k,226) - lu(k,1162) * b(k,201)
894 315224452 : b(k,237) = b(k,237) - lu(k,1163) * b(k,201)
895 315224452 : b(k,217) = b(k,217) - lu(k,1167) * b(k,202)
896 315224452 : b(k,226) = b(k,226) - lu(k,1168) * b(k,202)
897 315224452 : b(k,229) = b(k,229) - lu(k,1169) * b(k,202)
898 315224452 : b(k,237) = b(k,237) - lu(k,1170) * b(k,202)
899 315224452 : b(k,206) = b(k,206) - lu(k,1179) * b(k,203)
900 315224452 : b(k,216) = b(k,216) - lu(k,1180) * b(k,203)
901 315224452 : b(k,217) = b(k,217) - lu(k,1181) * b(k,203)
902 315224452 : b(k,226) = b(k,226) - lu(k,1182) * b(k,203)
903 315224452 : b(k,227) = b(k,227) - lu(k,1183) * b(k,203)
904 315224452 : b(k,228) = b(k,228) - lu(k,1184) * b(k,203)
905 315224452 : b(k,229) = b(k,229) - lu(k,1185) * b(k,203)
906 315224452 : b(k,231) = b(k,231) - lu(k,1186) * b(k,203)
907 315224452 : b(k,237) = b(k,237) - lu(k,1187) * b(k,203)
908 315224452 : b(k,238) = b(k,238) - lu(k,1188) * b(k,203)
909 325174377 : b(k,239) = b(k,239) - lu(k,1189) * b(k,203)
910 : end do
911 9949925 : end subroutine lu_slv04
912 9949925 : subroutine lu_slv05( avec_len, lu, b )
913 : use shr_kind_mod, only : r8 => shr_kind_r8
914 : use chem_mods, only : clscnt4, nzcnt
915 : implicit none
916 : !-----------------------------------------------------------------------
917 : ! ... Dummy args
918 : !-----------------------------------------------------------------------
919 : integer, intent(in) :: avec_len
920 : real(r8), intent(in) :: lu(veclen,max(1,nzcnt))
921 : real(r8), intent(inout) :: b(veclen,clscnt4)
922 : !-----------------------------------------------------------------------
923 : ! ... Local variables
924 : !-----------------------------------------------------------------------
925 : integer :: k
926 : !-----------------------------------------------------------------------
927 : ! ... solve L * y = b
928 : !-----------------------------------------------------------------------
929 325174377 : do k = 1,avec_len
930 315224452 : b(k,205) = b(k,205) - lu(k,1196) * b(k,204)
931 315224452 : b(k,206) = b(k,206) - lu(k,1197) * b(k,204)
932 315224452 : b(k,216) = b(k,216) - lu(k,1198) * b(k,204)
933 315224452 : b(k,217) = b(k,217) - lu(k,1199) * b(k,204)
934 315224452 : b(k,225) = b(k,225) - lu(k,1200) * b(k,204)
935 315224452 : b(k,226) = b(k,226) - lu(k,1201) * b(k,204)
936 315224452 : b(k,228) = b(k,228) - lu(k,1202) * b(k,204)
937 315224452 : b(k,229) = b(k,229) - lu(k,1203) * b(k,204)
938 315224452 : b(k,231) = b(k,231) - lu(k,1204) * b(k,204)
939 315224452 : b(k,237) = b(k,237) - lu(k,1205) * b(k,204)
940 315224452 : b(k,238) = b(k,238) - lu(k,1206) * b(k,204)
941 315224452 : b(k,239) = b(k,239) - lu(k,1207) * b(k,204)
942 315224452 : b(k,210) = b(k,210) - lu(k,1209) * b(k,205)
943 315224452 : b(k,216) = b(k,216) - lu(k,1210) * b(k,205)
944 315224452 : b(k,226) = b(k,226) - lu(k,1211) * b(k,205)
945 315224452 : b(k,229) = b(k,229) - lu(k,1212) * b(k,205)
946 315224452 : b(k,237) = b(k,237) - lu(k,1213) * b(k,205)
947 315224452 : b(k,210) = b(k,210) - lu(k,1218) * b(k,206)
948 315224452 : b(k,216) = b(k,216) - lu(k,1219) * b(k,206)
949 315224452 : b(k,226) = b(k,226) - lu(k,1220) * b(k,206)
950 315224452 : b(k,227) = b(k,227) - lu(k,1221) * b(k,206)
951 315224452 : b(k,228) = b(k,228) - lu(k,1222) * b(k,206)
952 315224452 : b(k,229) = b(k,229) - lu(k,1223) * b(k,206)
953 315224452 : b(k,237) = b(k,237) - lu(k,1224) * b(k,206)
954 315224452 : b(k,238) = b(k,238) - lu(k,1225) * b(k,206)
955 315224452 : b(k,239) = b(k,239) - lu(k,1226) * b(k,206)
956 315224452 : b(k,210) = b(k,210) - lu(k,1234) * b(k,207)
957 315224452 : b(k,216) = b(k,216) - lu(k,1235) * b(k,207)
958 315224452 : b(k,217) = b(k,217) - lu(k,1236) * b(k,207)
959 315224452 : b(k,225) = b(k,225) - lu(k,1237) * b(k,207)
960 315224452 : b(k,226) = b(k,226) - lu(k,1238) * b(k,207)
961 315224452 : b(k,227) = b(k,227) - lu(k,1239) * b(k,207)
962 315224452 : b(k,228) = b(k,228) - lu(k,1240) * b(k,207)
963 315224452 : b(k,229) = b(k,229) - lu(k,1241) * b(k,207)
964 315224452 : b(k,237) = b(k,237) - lu(k,1242) * b(k,207)
965 315224452 : b(k,238) = b(k,238) - lu(k,1243) * b(k,207)
966 315224452 : b(k,209) = b(k,209) - lu(k,1254) * b(k,208)
967 315224452 : b(k,210) = b(k,210) - lu(k,1255) * b(k,208)
968 315224452 : b(k,214) = b(k,214) - lu(k,1256) * b(k,208)
969 315224452 : b(k,216) = b(k,216) - lu(k,1257) * b(k,208)
970 315224452 : b(k,217) = b(k,217) - lu(k,1258) * b(k,208)
971 315224452 : b(k,225) = b(k,225) - lu(k,1259) * b(k,208)
972 315224452 : b(k,226) = b(k,226) - lu(k,1260) * b(k,208)
973 315224452 : b(k,227) = b(k,227) - lu(k,1261) * b(k,208)
974 315224452 : b(k,228) = b(k,228) - lu(k,1262) * b(k,208)
975 315224452 : b(k,229) = b(k,229) - lu(k,1263) * b(k,208)
976 315224452 : b(k,231) = b(k,231) - lu(k,1264) * b(k,208)
977 315224452 : b(k,237) = b(k,237) - lu(k,1265) * b(k,208)
978 315224452 : b(k,238) = b(k,238) - lu(k,1266) * b(k,208)
979 315224452 : b(k,210) = b(k,210) - lu(k,1270) * b(k,209)
980 315224452 : b(k,213) = b(k,213) - lu(k,1271) * b(k,209)
981 315224452 : b(k,215) = b(k,215) - lu(k,1272) * b(k,209)
982 315224452 : b(k,216) = b(k,216) - lu(k,1273) * b(k,209)
983 315224452 : b(k,217) = b(k,217) - lu(k,1274) * b(k,209)
984 315224452 : b(k,226) = b(k,226) - lu(k,1275) * b(k,209)
985 315224452 : b(k,229) = b(k,229) - lu(k,1276) * b(k,209)
986 315224452 : b(k,235) = b(k,235) - lu(k,1277) * b(k,209)
987 315224452 : b(k,237) = b(k,237) - lu(k,1278) * b(k,209)
988 315224452 : b(k,239) = b(k,239) - lu(k,1279) * b(k,209)
989 315224452 : b(k,216) = b(k,216) - lu(k,1282) * b(k,210)
990 315224452 : b(k,217) = b(k,217) - lu(k,1283) * b(k,210)
991 315224452 : b(k,225) = b(k,225) - lu(k,1284) * b(k,210)
992 315224452 : b(k,226) = b(k,226) - lu(k,1285) * b(k,210)
993 315224452 : b(k,231) = b(k,231) - lu(k,1286) * b(k,210)
994 315224452 : b(k,237) = b(k,237) - lu(k,1287) * b(k,210)
995 315224452 : b(k,239) = b(k,239) - lu(k,1288) * b(k,210)
996 315224452 : b(k,213) = b(k,213) - lu(k,1300) * b(k,211)
997 315224452 : b(k,214) = b(k,214) - lu(k,1301) * b(k,211)
998 315224452 : b(k,215) = b(k,215) - lu(k,1302) * b(k,211)
999 315224452 : b(k,216) = b(k,216) - lu(k,1303) * b(k,211)
1000 315224452 : b(k,217) = b(k,217) - lu(k,1304) * b(k,211)
1001 315224452 : b(k,225) = b(k,225) - lu(k,1305) * b(k,211)
1002 315224452 : b(k,226) = b(k,226) - lu(k,1306) * b(k,211)
1003 315224452 : b(k,227) = b(k,227) - lu(k,1307) * b(k,211)
1004 315224452 : b(k,228) = b(k,228) - lu(k,1308) * b(k,211)
1005 315224452 : b(k,229) = b(k,229) - lu(k,1309) * b(k,211)
1006 315224452 : b(k,231) = b(k,231) - lu(k,1310) * b(k,211)
1007 315224452 : b(k,235) = b(k,235) - lu(k,1311) * b(k,211)
1008 315224452 : b(k,237) = b(k,237) - lu(k,1312) * b(k,211)
1009 315224452 : b(k,238) = b(k,238) - lu(k,1313) * b(k,211)
1010 315224452 : b(k,239) = b(k,239) - lu(k,1314) * b(k,211)
1011 315224452 : b(k,213) = b(k,213) - lu(k,1333) * b(k,212)
1012 315224452 : b(k,214) = b(k,214) - lu(k,1334) * b(k,212)
1013 315224452 : b(k,215) = b(k,215) - lu(k,1335) * b(k,212)
1014 315224452 : b(k,216) = b(k,216) - lu(k,1336) * b(k,212)
1015 315224452 : b(k,217) = b(k,217) - lu(k,1337) * b(k,212)
1016 315224452 : b(k,225) = b(k,225) - lu(k,1338) * b(k,212)
1017 315224452 : b(k,226) = b(k,226) - lu(k,1339) * b(k,212)
1018 315224452 : b(k,227) = b(k,227) - lu(k,1340) * b(k,212)
1019 315224452 : b(k,228) = b(k,228) - lu(k,1341) * b(k,212)
1020 315224452 : b(k,229) = b(k,229) - lu(k,1342) * b(k,212)
1021 315224452 : b(k,231) = b(k,231) - lu(k,1343) * b(k,212)
1022 315224452 : b(k,235) = b(k,235) - lu(k,1344) * b(k,212)
1023 315224452 : b(k,237) = b(k,237) - lu(k,1345) * b(k,212)
1024 315224452 : b(k,238) = b(k,238) - lu(k,1346) * b(k,212)
1025 315224452 : b(k,239) = b(k,239) - lu(k,1347) * b(k,212)
1026 315224452 : b(k,215) = b(k,215) - lu(k,1357) * b(k,213)
1027 315224452 : b(k,216) = b(k,216) - lu(k,1358) * b(k,213)
1028 315224452 : b(k,217) = b(k,217) - lu(k,1359) * b(k,213)
1029 315224452 : b(k,225) = b(k,225) - lu(k,1360) * b(k,213)
1030 315224452 : b(k,226) = b(k,226) - lu(k,1361) * b(k,213)
1031 315224452 : b(k,227) = b(k,227) - lu(k,1362) * b(k,213)
1032 315224452 : b(k,228) = b(k,228) - lu(k,1363) * b(k,213)
1033 315224452 : b(k,229) = b(k,229) - lu(k,1364) * b(k,213)
1034 315224452 : b(k,231) = b(k,231) - lu(k,1365) * b(k,213)
1035 315224452 : b(k,237) = b(k,237) - lu(k,1366) * b(k,213)
1036 315224452 : b(k,238) = b(k,238) - lu(k,1367) * b(k,213)
1037 315224452 : b(k,239) = b(k,239) - lu(k,1368) * b(k,213)
1038 315224452 : b(k,215) = b(k,215) - lu(k,1377) * b(k,214)
1039 315224452 : b(k,216) = b(k,216) - lu(k,1378) * b(k,214)
1040 315224452 : b(k,217) = b(k,217) - lu(k,1379) * b(k,214)
1041 315224452 : b(k,222) = b(k,222) - lu(k,1380) * b(k,214)
1042 315224452 : b(k,225) = b(k,225) - lu(k,1381) * b(k,214)
1043 315224452 : b(k,226) = b(k,226) - lu(k,1382) * b(k,214)
1044 315224452 : b(k,227) = b(k,227) - lu(k,1383) * b(k,214)
1045 315224452 : b(k,228) = b(k,228) - lu(k,1384) * b(k,214)
1046 315224452 : b(k,229) = b(k,229) - lu(k,1385) * b(k,214)
1047 315224452 : b(k,231) = b(k,231) - lu(k,1386) * b(k,214)
1048 315224452 : b(k,235) = b(k,235) - lu(k,1387) * b(k,214)
1049 315224452 : b(k,236) = b(k,236) - lu(k,1388) * b(k,214)
1050 315224452 : b(k,237) = b(k,237) - lu(k,1389) * b(k,214)
1051 315224452 : b(k,238) = b(k,238) - lu(k,1390) * b(k,214)
1052 315224452 : b(k,239) = b(k,239) - lu(k,1391) * b(k,214)
1053 315224452 : b(k,216) = b(k,216) - lu(k,1398) * b(k,215)
1054 315224452 : b(k,217) = b(k,217) - lu(k,1399) * b(k,215)
1055 315224452 : b(k,225) = b(k,225) - lu(k,1400) * b(k,215)
1056 315224452 : b(k,226) = b(k,226) - lu(k,1401) * b(k,215)
1057 315224452 : b(k,227) = b(k,227) - lu(k,1402) * b(k,215)
1058 315224452 : b(k,228) = b(k,228) - lu(k,1403) * b(k,215)
1059 315224452 : b(k,229) = b(k,229) - lu(k,1404) * b(k,215)
1060 315224452 : b(k,231) = b(k,231) - lu(k,1405) * b(k,215)
1061 315224452 : b(k,235) = b(k,235) - lu(k,1406) * b(k,215)
1062 315224452 : b(k,237) = b(k,237) - lu(k,1407) * b(k,215)
1063 315224452 : b(k,238) = b(k,238) - lu(k,1408) * b(k,215)
1064 315224452 : b(k,239) = b(k,239) - lu(k,1409) * b(k,215)
1065 315224452 : b(k,217) = b(k,217) - lu(k,1429) * b(k,216)
1066 315224452 : b(k,222) = b(k,222) - lu(k,1430) * b(k,216)
1067 315224452 : b(k,225) = b(k,225) - lu(k,1431) * b(k,216)
1068 315224452 : b(k,226) = b(k,226) - lu(k,1432) * b(k,216)
1069 315224452 : b(k,227) = b(k,227) - lu(k,1433) * b(k,216)
1070 315224452 : b(k,228) = b(k,228) - lu(k,1434) * b(k,216)
1071 315224452 : b(k,229) = b(k,229) - lu(k,1435) * b(k,216)
1072 315224452 : b(k,231) = b(k,231) - lu(k,1436) * b(k,216)
1073 315224452 : b(k,235) = b(k,235) - lu(k,1437) * b(k,216)
1074 315224452 : b(k,236) = b(k,236) - lu(k,1438) * b(k,216)
1075 315224452 : b(k,237) = b(k,237) - lu(k,1439) * b(k,216)
1076 315224452 : b(k,238) = b(k,238) - lu(k,1440) * b(k,216)
1077 315224452 : b(k,239) = b(k,239) - lu(k,1441) * b(k,216)
1078 315224452 : b(k,222) = b(k,222) - lu(k,1451) * b(k,217)
1079 315224452 : b(k,226) = b(k,226) - lu(k,1452) * b(k,217)
1080 315224452 : b(k,227) = b(k,227) - lu(k,1453) * b(k,217)
1081 315224452 : b(k,228) = b(k,228) - lu(k,1454) * b(k,217)
1082 315224452 : b(k,230) = b(k,230) - lu(k,1455) * b(k,217)
1083 315224452 : b(k,232) = b(k,232) - lu(k,1456) * b(k,217)
1084 315224452 : b(k,233) = b(k,233) - lu(k,1457) * b(k,217)
1085 315224452 : b(k,237) = b(k,237) - lu(k,1458) * b(k,217)
1086 315224452 : b(k,220) = b(k,220) - lu(k,1462) * b(k,218)
1087 315224452 : b(k,225) = b(k,225) - lu(k,1463) * b(k,218)
1088 315224452 : b(k,226) = b(k,226) - lu(k,1464) * b(k,218)
1089 315224452 : b(k,230) = b(k,230) - lu(k,1465) * b(k,218)
1090 315224452 : b(k,231) = b(k,231) - lu(k,1466) * b(k,218)
1091 315224452 : b(k,232) = b(k,232) - lu(k,1467) * b(k,218)
1092 315224452 : b(k,233) = b(k,233) - lu(k,1468) * b(k,218)
1093 315224452 : b(k,234) = b(k,234) - lu(k,1469) * b(k,218)
1094 315224452 : b(k,236) = b(k,236) - lu(k,1470) * b(k,218)
1095 315224452 : b(k,238) = b(k,238) - lu(k,1471) * b(k,218)
1096 315224452 : b(k,239) = b(k,239) - lu(k,1472) * b(k,218)
1097 315224452 : b(k,221) = b(k,221) - lu(k,1477) * b(k,219)
1098 315224452 : b(k,222) = b(k,222) - lu(k,1478) * b(k,219)
1099 315224452 : b(k,223) = b(k,223) - lu(k,1479) * b(k,219)
1100 315224452 : b(k,224) = b(k,224) - lu(k,1480) * b(k,219)
1101 315224452 : b(k,226) = b(k,226) - lu(k,1481) * b(k,219)
1102 315224452 : b(k,227) = b(k,227) - lu(k,1482) * b(k,219)
1103 315224452 : b(k,228) = b(k,228) - lu(k,1483) * b(k,219)
1104 315224452 : b(k,232) = b(k,232) - lu(k,1484) * b(k,219)
1105 315224452 : b(k,233) = b(k,233) - lu(k,1485) * b(k,219)
1106 315224452 : b(k,234) = b(k,234) - lu(k,1486) * b(k,219)
1107 315224452 : b(k,235) = b(k,235) - lu(k,1487) * b(k,219)
1108 315224452 : b(k,237) = b(k,237) - lu(k,1488) * b(k,219)
1109 315224452 : b(k,239) = b(k,239) - lu(k,1489) * b(k,219)
1110 315224452 : b(k,221) = b(k,221) - lu(k,1496) * b(k,220)
1111 315224452 : b(k,223) = b(k,223) - lu(k,1497) * b(k,220)
1112 315224452 : b(k,224) = b(k,224) - lu(k,1498) * b(k,220)
1113 315224452 : b(k,225) = b(k,225) - lu(k,1499) * b(k,220)
1114 315224452 : b(k,226) = b(k,226) - lu(k,1500) * b(k,220)
1115 315224452 : b(k,228) = b(k,228) - lu(k,1501) * b(k,220)
1116 315224452 : b(k,230) = b(k,230) - lu(k,1502) * b(k,220)
1117 315224452 : b(k,231) = b(k,231) - lu(k,1503) * b(k,220)
1118 315224452 : b(k,232) = b(k,232) - lu(k,1504) * b(k,220)
1119 315224452 : b(k,233) = b(k,233) - lu(k,1505) * b(k,220)
1120 315224452 : b(k,234) = b(k,234) - lu(k,1506) * b(k,220)
1121 315224452 : b(k,239) = b(k,239) - lu(k,1507) * b(k,220)
1122 315224452 : b(k,222) = b(k,222) - lu(k,1512) * b(k,221)
1123 315224452 : b(k,223) = b(k,223) - lu(k,1513) * b(k,221)
1124 315224452 : b(k,226) = b(k,226) - lu(k,1514) * b(k,221)
1125 315224452 : b(k,227) = b(k,227) - lu(k,1515) * b(k,221)
1126 315224452 : b(k,228) = b(k,228) - lu(k,1516) * b(k,221)
1127 315224452 : b(k,229) = b(k,229) - lu(k,1517) * b(k,221)
1128 315224452 : b(k,230) = b(k,230) - lu(k,1518) * b(k,221)
1129 315224452 : b(k,232) = b(k,232) - lu(k,1519) * b(k,221)
1130 315224452 : b(k,233) = b(k,233) - lu(k,1520) * b(k,221)
1131 315224452 : b(k,235) = b(k,235) - lu(k,1521) * b(k,221)
1132 315224452 : b(k,237) = b(k,237) - lu(k,1522) * b(k,221)
1133 325174377 : b(k,239) = b(k,239) - lu(k,1523) * b(k,221)
1134 : end do
1135 9949925 : end subroutine lu_slv05
1136 9949925 : subroutine lu_slv06( avec_len, lu, b )
1137 : use shr_kind_mod, only : r8 => shr_kind_r8
1138 : use chem_mods, only : clscnt4, nzcnt
1139 : implicit none
1140 : !-----------------------------------------------------------------------
1141 : ! ... Dummy args
1142 : !-----------------------------------------------------------------------
1143 : integer, intent(in) :: avec_len
1144 : real(r8), intent(in) :: lu(veclen,max(1,nzcnt))
1145 : real(r8), intent(inout) :: b(veclen,clscnt4)
1146 : !-----------------------------------------------------------------------
1147 : ! ... Local variables
1148 : !-----------------------------------------------------------------------
1149 : integer :: k
1150 : !-----------------------------------------------------------------------
1151 : ! ... solve L * y = b
1152 : !-----------------------------------------------------------------------
1153 325174377 : do k = 1,avec_len
1154 315224452 : b(k,223) = b(k,223) - lu(k,1543) * b(k,222)
1155 315224452 : b(k,224) = b(k,224) - lu(k,1544) * b(k,222)
1156 315224452 : b(k,226) = b(k,226) - lu(k,1545) * b(k,222)
1157 315224452 : b(k,227) = b(k,227) - lu(k,1546) * b(k,222)
1158 315224452 : b(k,228) = b(k,228) - lu(k,1547) * b(k,222)
1159 315224452 : b(k,229) = b(k,229) - lu(k,1548) * b(k,222)
1160 315224452 : b(k,230) = b(k,230) - lu(k,1549) * b(k,222)
1161 315224452 : b(k,232) = b(k,232) - lu(k,1550) * b(k,222)
1162 315224452 : b(k,233) = b(k,233) - lu(k,1551) * b(k,222)
1163 315224452 : b(k,234) = b(k,234) - lu(k,1552) * b(k,222)
1164 315224452 : b(k,235) = b(k,235) - lu(k,1553) * b(k,222)
1165 315224452 : b(k,237) = b(k,237) - lu(k,1554) * b(k,222)
1166 315224452 : b(k,239) = b(k,239) - lu(k,1555) * b(k,222)
1167 315224452 : b(k,224) = b(k,224) - lu(k,1566) * b(k,223)
1168 315224452 : b(k,225) = b(k,225) - lu(k,1567) * b(k,223)
1169 315224452 : b(k,226) = b(k,226) - lu(k,1568) * b(k,223)
1170 315224452 : b(k,227) = b(k,227) - lu(k,1569) * b(k,223)
1171 315224452 : b(k,228) = b(k,228) - lu(k,1570) * b(k,223)
1172 315224452 : b(k,229) = b(k,229) - lu(k,1571) * b(k,223)
1173 315224452 : b(k,230) = b(k,230) - lu(k,1572) * b(k,223)
1174 315224452 : b(k,231) = b(k,231) - lu(k,1573) * b(k,223)
1175 315224452 : b(k,232) = b(k,232) - lu(k,1574) * b(k,223)
1176 315224452 : b(k,233) = b(k,233) - lu(k,1575) * b(k,223)
1177 315224452 : b(k,234) = b(k,234) - lu(k,1576) * b(k,223)
1178 315224452 : b(k,235) = b(k,235) - lu(k,1577) * b(k,223)
1179 315224452 : b(k,237) = b(k,237) - lu(k,1578) * b(k,223)
1180 315224452 : b(k,239) = b(k,239) - lu(k,1579) * b(k,223)
1181 315224452 : b(k,225) = b(k,225) - lu(k,1593) * b(k,224)
1182 315224452 : b(k,226) = b(k,226) - lu(k,1594) * b(k,224)
1183 315224452 : b(k,227) = b(k,227) - lu(k,1595) * b(k,224)
1184 315224452 : b(k,228) = b(k,228) - lu(k,1596) * b(k,224)
1185 315224452 : b(k,229) = b(k,229) - lu(k,1597) * b(k,224)
1186 315224452 : b(k,230) = b(k,230) - lu(k,1598) * b(k,224)
1187 315224452 : b(k,231) = b(k,231) - lu(k,1599) * b(k,224)
1188 315224452 : b(k,232) = b(k,232) - lu(k,1600) * b(k,224)
1189 315224452 : b(k,233) = b(k,233) - lu(k,1601) * b(k,224)
1190 315224452 : b(k,234) = b(k,234) - lu(k,1602) * b(k,224)
1191 315224452 : b(k,235) = b(k,235) - lu(k,1603) * b(k,224)
1192 315224452 : b(k,237) = b(k,237) - lu(k,1604) * b(k,224)
1193 315224452 : b(k,238) = b(k,238) - lu(k,1605) * b(k,224)
1194 315224452 : b(k,239) = b(k,239) - lu(k,1606) * b(k,224)
1195 315224452 : b(k,226) = b(k,226) - lu(k,1653) * b(k,225)
1196 315224452 : b(k,227) = b(k,227) - lu(k,1654) * b(k,225)
1197 315224452 : b(k,228) = b(k,228) - lu(k,1655) * b(k,225)
1198 315224452 : b(k,229) = b(k,229) - lu(k,1656) * b(k,225)
1199 315224452 : b(k,230) = b(k,230) - lu(k,1657) * b(k,225)
1200 315224452 : b(k,231) = b(k,231) - lu(k,1658) * b(k,225)
1201 315224452 : b(k,232) = b(k,232) - lu(k,1659) * b(k,225)
1202 315224452 : b(k,233) = b(k,233) - lu(k,1660) * b(k,225)
1203 315224452 : b(k,234) = b(k,234) - lu(k,1661) * b(k,225)
1204 315224452 : b(k,235) = b(k,235) - lu(k,1662) * b(k,225)
1205 315224452 : b(k,236) = b(k,236) - lu(k,1663) * b(k,225)
1206 315224452 : b(k,237) = b(k,237) - lu(k,1664) * b(k,225)
1207 315224452 : b(k,238) = b(k,238) - lu(k,1665) * b(k,225)
1208 315224452 : b(k,239) = b(k,239) - lu(k,1666) * b(k,225)
1209 315224452 : b(k,227) = b(k,227) - lu(k,1819) * b(k,226)
1210 315224452 : b(k,228) = b(k,228) - lu(k,1820) * b(k,226)
1211 315224452 : b(k,229) = b(k,229) - lu(k,1821) * b(k,226)
1212 315224452 : b(k,230) = b(k,230) - lu(k,1822) * b(k,226)
1213 315224452 : b(k,231) = b(k,231) - lu(k,1823) * b(k,226)
1214 315224452 : b(k,232) = b(k,232) - lu(k,1824) * b(k,226)
1215 315224452 : b(k,233) = b(k,233) - lu(k,1825) * b(k,226)
1216 315224452 : b(k,234) = b(k,234) - lu(k,1826) * b(k,226)
1217 315224452 : b(k,235) = b(k,235) - lu(k,1827) * b(k,226)
1218 315224452 : b(k,236) = b(k,236) - lu(k,1828) * b(k,226)
1219 315224452 : b(k,237) = b(k,237) - lu(k,1829) * b(k,226)
1220 315224452 : b(k,238) = b(k,238) - lu(k,1830) * b(k,226)
1221 315224452 : b(k,239) = b(k,239) - lu(k,1831) * b(k,226)
1222 315224452 : b(k,228) = b(k,228) - lu(k,1915) * b(k,227)
1223 315224452 : b(k,229) = b(k,229) - lu(k,1916) * b(k,227)
1224 315224452 : b(k,230) = b(k,230) - lu(k,1917) * b(k,227)
1225 315224452 : b(k,231) = b(k,231) - lu(k,1918) * b(k,227)
1226 315224452 : b(k,232) = b(k,232) - lu(k,1919) * b(k,227)
1227 315224452 : b(k,233) = b(k,233) - lu(k,1920) * b(k,227)
1228 315224452 : b(k,234) = b(k,234) - lu(k,1921) * b(k,227)
1229 315224452 : b(k,235) = b(k,235) - lu(k,1922) * b(k,227)
1230 315224452 : b(k,236) = b(k,236) - lu(k,1923) * b(k,227)
1231 315224452 : b(k,237) = b(k,237) - lu(k,1924) * b(k,227)
1232 315224452 : b(k,238) = b(k,238) - lu(k,1925) * b(k,227)
1233 315224452 : b(k,239) = b(k,239) - lu(k,1926) * b(k,227)
1234 315224452 : b(k,229) = b(k,229) - lu(k,1961) * b(k,228)
1235 315224452 : b(k,230) = b(k,230) - lu(k,1962) * b(k,228)
1236 315224452 : b(k,231) = b(k,231) - lu(k,1963) * b(k,228)
1237 315224452 : b(k,232) = b(k,232) - lu(k,1964) * b(k,228)
1238 315224452 : b(k,233) = b(k,233) - lu(k,1965) * b(k,228)
1239 315224452 : b(k,234) = b(k,234) - lu(k,1966) * b(k,228)
1240 315224452 : b(k,235) = b(k,235) - lu(k,1967) * b(k,228)
1241 315224452 : b(k,236) = b(k,236) - lu(k,1968) * b(k,228)
1242 315224452 : b(k,237) = b(k,237) - lu(k,1969) * b(k,228)
1243 315224452 : b(k,238) = b(k,238) - lu(k,1970) * b(k,228)
1244 315224452 : b(k,239) = b(k,239) - lu(k,1971) * b(k,228)
1245 315224452 : b(k,230) = b(k,230) - lu(k,1988) * b(k,229)
1246 315224452 : b(k,231) = b(k,231) - lu(k,1989) * b(k,229)
1247 315224452 : b(k,232) = b(k,232) - lu(k,1990) * b(k,229)
1248 315224452 : b(k,233) = b(k,233) - lu(k,1991) * b(k,229)
1249 315224452 : b(k,234) = b(k,234) - lu(k,1992) * b(k,229)
1250 315224452 : b(k,235) = b(k,235) - lu(k,1993) * b(k,229)
1251 315224452 : b(k,236) = b(k,236) - lu(k,1994) * b(k,229)
1252 315224452 : b(k,237) = b(k,237) - lu(k,1995) * b(k,229)
1253 315224452 : b(k,238) = b(k,238) - lu(k,1996) * b(k,229)
1254 315224452 : b(k,239) = b(k,239) - lu(k,1997) * b(k,229)
1255 315224452 : b(k,231) = b(k,231) - lu(k,2032) * b(k,230)
1256 315224452 : b(k,232) = b(k,232) - lu(k,2033) * b(k,230)
1257 315224452 : b(k,233) = b(k,233) - lu(k,2034) * b(k,230)
1258 315224452 : b(k,234) = b(k,234) - lu(k,2035) * b(k,230)
1259 315224452 : b(k,235) = b(k,235) - lu(k,2036) * b(k,230)
1260 315224452 : b(k,236) = b(k,236) - lu(k,2037) * b(k,230)
1261 315224452 : b(k,237) = b(k,237) - lu(k,2038) * b(k,230)
1262 315224452 : b(k,238) = b(k,238) - lu(k,2039) * b(k,230)
1263 315224452 : b(k,239) = b(k,239) - lu(k,2040) * b(k,230)
1264 315224452 : b(k,232) = b(k,232) - lu(k,2056) * b(k,231)
1265 315224452 : b(k,233) = b(k,233) - lu(k,2057) * b(k,231)
1266 315224452 : b(k,234) = b(k,234) - lu(k,2058) * b(k,231)
1267 315224452 : b(k,235) = b(k,235) - lu(k,2059) * b(k,231)
1268 315224452 : b(k,236) = b(k,236) - lu(k,2060) * b(k,231)
1269 315224452 : b(k,237) = b(k,237) - lu(k,2061) * b(k,231)
1270 315224452 : b(k,238) = b(k,238) - lu(k,2062) * b(k,231)
1271 315224452 : b(k,239) = b(k,239) - lu(k,2063) * b(k,231)
1272 315224452 : b(k,233) = b(k,233) - lu(k,2100) * b(k,232)
1273 315224452 : b(k,234) = b(k,234) - lu(k,2101) * b(k,232)
1274 315224452 : b(k,235) = b(k,235) - lu(k,2102) * b(k,232)
1275 315224452 : b(k,236) = b(k,236) - lu(k,2103) * b(k,232)
1276 315224452 : b(k,237) = b(k,237) - lu(k,2104) * b(k,232)
1277 315224452 : b(k,238) = b(k,238) - lu(k,2105) * b(k,232)
1278 315224452 : b(k,239) = b(k,239) - lu(k,2106) * b(k,232)
1279 315224452 : b(k,234) = b(k,234) - lu(k,2122) * b(k,233)
1280 315224452 : b(k,235) = b(k,235) - lu(k,2123) * b(k,233)
1281 315224452 : b(k,236) = b(k,236) - lu(k,2124) * b(k,233)
1282 315224452 : b(k,237) = b(k,237) - lu(k,2125) * b(k,233)
1283 315224452 : b(k,238) = b(k,238) - lu(k,2126) * b(k,233)
1284 315224452 : b(k,239) = b(k,239) - lu(k,2127) * b(k,233)
1285 315224452 : b(k,235) = b(k,235) - lu(k,2169) * b(k,234)
1286 315224452 : b(k,236) = b(k,236) - lu(k,2170) * b(k,234)
1287 315224452 : b(k,237) = b(k,237) - lu(k,2171) * b(k,234)
1288 315224452 : b(k,238) = b(k,238) - lu(k,2172) * b(k,234)
1289 315224452 : b(k,239) = b(k,239) - lu(k,2173) * b(k,234)
1290 315224452 : b(k,236) = b(k,236) - lu(k,2234) * b(k,235)
1291 315224452 : b(k,237) = b(k,237) - lu(k,2235) * b(k,235)
1292 315224452 : b(k,238) = b(k,238) - lu(k,2236) * b(k,235)
1293 315224452 : b(k,239) = b(k,239) - lu(k,2237) * b(k,235)
1294 315224452 : b(k,237) = b(k,237) - lu(k,2260) * b(k,236)
1295 315224452 : b(k,238) = b(k,238) - lu(k,2261) * b(k,236)
1296 315224452 : b(k,239) = b(k,239) - lu(k,2262) * b(k,236)
1297 315224452 : b(k,238) = b(k,238) - lu(k,2369) * b(k,237)
1298 315224452 : b(k,239) = b(k,239) - lu(k,2370) * b(k,237)
1299 325174377 : b(k,239) = b(k,239) - lu(k,2422) * b(k,238)
1300 : end do
1301 9949925 : end subroutine lu_slv06
1302 9949925 : subroutine lu_slv07( avec_len, lu, b )
1303 : use shr_kind_mod, only : r8 => shr_kind_r8
1304 : use chem_mods, only : clscnt4, nzcnt
1305 : implicit none
1306 : !-----------------------------------------------------------------------
1307 : ! ... Dummy args
1308 : !-----------------------------------------------------------------------
1309 : integer, intent(in) :: avec_len
1310 : real(r8), intent(in) :: lu(veclen,max(1,nzcnt))
1311 : real(r8), intent(inout) :: b(veclen,clscnt4)
1312 : !-----------------------------------------------------------------------
1313 : ! ... Local variables
1314 : !-----------------------------------------------------------------------
1315 : integer :: k
1316 : !-----------------------------------------------------------------------
1317 : ! ... solve L * y = b
1318 : !-----------------------------------------------------------------------
1319 325174377 : do k = 1,avec_len
1320 : !-----------------------------------------------------------------------
1321 : ! ... Solve U * x = y
1322 : !-----------------------------------------------------------------------
1323 315224452 : b(k,239) = b(k,239) * lu(k,2449)
1324 315224452 : b(k,238) = b(k,238) - lu(k,2448) * b(k,239)
1325 315224452 : b(k,237) = b(k,237) - lu(k,2447) * b(k,239)
1326 315224452 : b(k,236) = b(k,236) - lu(k,2446) * b(k,239)
1327 315224452 : b(k,235) = b(k,235) - lu(k,2445) * b(k,239)
1328 315224452 : b(k,234) = b(k,234) - lu(k,2444) * b(k,239)
1329 315224452 : b(k,233) = b(k,233) - lu(k,2443) * b(k,239)
1330 315224452 : b(k,232) = b(k,232) - lu(k,2442) * b(k,239)
1331 315224452 : b(k,231) = b(k,231) - lu(k,2441) * b(k,239)
1332 315224452 : b(k,230) = b(k,230) - lu(k,2440) * b(k,239)
1333 315224452 : b(k,229) = b(k,229) - lu(k,2439) * b(k,239)
1334 315224452 : b(k,228) = b(k,228) - lu(k,2438) * b(k,239)
1335 315224452 : b(k,227) = b(k,227) - lu(k,2437) * b(k,239)
1336 315224452 : b(k,226) = b(k,226) - lu(k,2436) * b(k,239)
1337 315224452 : b(k,225) = b(k,225) - lu(k,2435) * b(k,239)
1338 315224452 : b(k,224) = b(k,224) - lu(k,2434) * b(k,239)
1339 315224452 : b(k,223) = b(k,223) - lu(k,2433) * b(k,239)
1340 315224452 : b(k,222) = b(k,222) - lu(k,2432) * b(k,239)
1341 315224452 : b(k,221) = b(k,221) - lu(k,2431) * b(k,239)
1342 315224452 : b(k,220) = b(k,220) - lu(k,2430) * b(k,239)
1343 315224452 : b(k,219) = b(k,219) - lu(k,2429) * b(k,239)
1344 315224452 : b(k,218) = b(k,218) - lu(k,2428) * b(k,239)
1345 315224452 : b(k,186) = b(k,186) - lu(k,2427) * b(k,239)
1346 315224452 : b(k,184) = b(k,184) - lu(k,2426) * b(k,239)
1347 315224452 : b(k,99) = b(k,99) - lu(k,2425) * b(k,239)
1348 315224452 : b(k,94) = b(k,94) - lu(k,2424) * b(k,239)
1349 315224452 : b(k,64) = b(k,64) - lu(k,2423) * b(k,239)
1350 315224452 : b(k,238) = b(k,238) * lu(k,2421)
1351 315224452 : b(k,237) = b(k,237) - lu(k,2420) * b(k,238)
1352 315224452 : b(k,236) = b(k,236) - lu(k,2419) * b(k,238)
1353 315224452 : b(k,235) = b(k,235) - lu(k,2418) * b(k,238)
1354 315224452 : b(k,234) = b(k,234) - lu(k,2417) * b(k,238)
1355 315224452 : b(k,233) = b(k,233) - lu(k,2416) * b(k,238)
1356 315224452 : b(k,232) = b(k,232) - lu(k,2415) * b(k,238)
1357 315224452 : b(k,231) = b(k,231) - lu(k,2414) * b(k,238)
1358 315224452 : b(k,230) = b(k,230) - lu(k,2413) * b(k,238)
1359 315224452 : b(k,229) = b(k,229) - lu(k,2412) * b(k,238)
1360 315224452 : b(k,228) = b(k,228) - lu(k,2411) * b(k,238)
1361 315224452 : b(k,227) = b(k,227) - lu(k,2410) * b(k,238)
1362 315224452 : b(k,226) = b(k,226) - lu(k,2409) * b(k,238)
1363 315224452 : b(k,225) = b(k,225) - lu(k,2408) * b(k,238)
1364 315224452 : b(k,224) = b(k,224) - lu(k,2407) * b(k,238)
1365 315224452 : b(k,223) = b(k,223) - lu(k,2406) * b(k,238)
1366 315224452 : b(k,222) = b(k,222) - lu(k,2405) * b(k,238)
1367 315224452 : b(k,217) = b(k,217) - lu(k,2404) * b(k,238)
1368 315224452 : b(k,216) = b(k,216) - lu(k,2403) * b(k,238)
1369 315224452 : b(k,215) = b(k,215) - lu(k,2402) * b(k,238)
1370 315224452 : b(k,214) = b(k,214) - lu(k,2401) * b(k,238)
1371 315224452 : b(k,213) = b(k,213) - lu(k,2400) * b(k,238)
1372 315224452 : b(k,212) = b(k,212) - lu(k,2399) * b(k,238)
1373 315224452 : b(k,211) = b(k,211) - lu(k,2398) * b(k,238)
1374 315224452 : b(k,210) = b(k,210) - lu(k,2397) * b(k,238)
1375 315224452 : b(k,209) = b(k,209) - lu(k,2396) * b(k,238)
1376 315224452 : b(k,208) = b(k,208) - lu(k,2395) * b(k,238)
1377 315224452 : b(k,207) = b(k,207) - lu(k,2394) * b(k,238)
1378 315224452 : b(k,206) = b(k,206) - lu(k,2393) * b(k,238)
1379 315224452 : b(k,205) = b(k,205) - lu(k,2392) * b(k,238)
1380 315224452 : b(k,203) = b(k,203) - lu(k,2391) * b(k,238)
1381 315224452 : b(k,202) = b(k,202) - lu(k,2390) * b(k,238)
1382 315224452 : b(k,201) = b(k,201) - lu(k,2389) * b(k,238)
1383 315224452 : b(k,200) = b(k,200) - lu(k,2388) * b(k,238)
1384 315224452 : b(k,198) = b(k,198) - lu(k,2387) * b(k,238)
1385 315224452 : b(k,197) = b(k,197) - lu(k,2386) * b(k,238)
1386 315224452 : b(k,196) = b(k,196) - lu(k,2385) * b(k,238)
1387 315224452 : b(k,195) = b(k,195) - lu(k,2384) * b(k,238)
1388 315224452 : b(k,194) = b(k,194) - lu(k,2383) * b(k,238)
1389 315224452 : b(k,193) = b(k,193) - lu(k,2382) * b(k,238)
1390 315224452 : b(k,192) = b(k,192) - lu(k,2381) * b(k,238)
1391 315224452 : b(k,185) = b(k,185) - lu(k,2380) * b(k,238)
1392 315224452 : b(k,181) = b(k,181) - lu(k,2379) * b(k,238)
1393 315224452 : b(k,180) = b(k,180) - lu(k,2378) * b(k,238)
1394 315224452 : b(k,173) = b(k,173) - lu(k,2377) * b(k,238)
1395 315224452 : b(k,161) = b(k,161) - lu(k,2376) * b(k,238)
1396 315224452 : b(k,153) = b(k,153) - lu(k,2375) * b(k,238)
1397 315224452 : b(k,128) = b(k,128) - lu(k,2374) * b(k,238)
1398 315224452 : b(k,124) = b(k,124) - lu(k,2373) * b(k,238)
1399 315224452 : b(k,107) = b(k,107) - lu(k,2372) * b(k,238)
1400 315224452 : b(k,95) = b(k,95) - lu(k,2371) * b(k,238)
1401 315224452 : b(k,237) = b(k,237) * lu(k,2368)
1402 315224452 : b(k,236) = b(k,236) - lu(k,2367) * b(k,237)
1403 315224452 : b(k,235) = b(k,235) - lu(k,2366) * b(k,237)
1404 315224452 : b(k,234) = b(k,234) - lu(k,2365) * b(k,237)
1405 315224452 : b(k,233) = b(k,233) - lu(k,2364) * b(k,237)
1406 315224452 : b(k,232) = b(k,232) - lu(k,2363) * b(k,237)
1407 315224452 : b(k,231) = b(k,231) - lu(k,2362) * b(k,237)
1408 315224452 : b(k,230) = b(k,230) - lu(k,2361) * b(k,237)
1409 315224452 : b(k,229) = b(k,229) - lu(k,2360) * b(k,237)
1410 315224452 : b(k,228) = b(k,228) - lu(k,2359) * b(k,237)
1411 315224452 : b(k,227) = b(k,227) - lu(k,2358) * b(k,237)
1412 315224452 : b(k,226) = b(k,226) - lu(k,2357) * b(k,237)
1413 315224452 : b(k,225) = b(k,225) - lu(k,2356) * b(k,237)
1414 315224452 : b(k,224) = b(k,224) - lu(k,2355) * b(k,237)
1415 315224452 : b(k,223) = b(k,223) - lu(k,2354) * b(k,237)
1416 315224452 : b(k,222) = b(k,222) - lu(k,2353) * b(k,237)
1417 315224452 : b(k,221) = b(k,221) - lu(k,2352) * b(k,237)
1418 315224452 : b(k,220) = b(k,220) - lu(k,2351) * b(k,237)
1419 315224452 : b(k,218) = b(k,218) - lu(k,2350) * b(k,237)
1420 315224452 : b(k,217) = b(k,217) - lu(k,2349) * b(k,237)
1421 315224452 : b(k,216) = b(k,216) - lu(k,2348) * b(k,237)
1422 315224452 : b(k,215) = b(k,215) - lu(k,2347) * b(k,237)
1423 315224452 : b(k,214) = b(k,214) - lu(k,2346) * b(k,237)
1424 315224452 : b(k,213) = b(k,213) - lu(k,2345) * b(k,237)
1425 315224452 : b(k,212) = b(k,212) - lu(k,2344) * b(k,237)
1426 315224452 : b(k,211) = b(k,211) - lu(k,2343) * b(k,237)
1427 315224452 : b(k,210) = b(k,210) - lu(k,2342) * b(k,237)
1428 315224452 : b(k,209) = b(k,209) - lu(k,2341) * b(k,237)
1429 315224452 : b(k,208) = b(k,208) - lu(k,2340) * b(k,237)
1430 315224452 : b(k,207) = b(k,207) - lu(k,2339) * b(k,237)
1431 315224452 : b(k,206) = b(k,206) - lu(k,2338) * b(k,237)
1432 315224452 : b(k,205) = b(k,205) - lu(k,2337) * b(k,237)
1433 315224452 : b(k,203) = b(k,203) - lu(k,2336) * b(k,237)
1434 315224452 : b(k,202) = b(k,202) - lu(k,2335) * b(k,237)
1435 315224452 : b(k,201) = b(k,201) - lu(k,2334) * b(k,237)
1436 315224452 : b(k,200) = b(k,200) - lu(k,2333) * b(k,237)
1437 315224452 : b(k,198) = b(k,198) - lu(k,2332) * b(k,237)
1438 315224452 : b(k,197) = b(k,197) - lu(k,2331) * b(k,237)
1439 315224452 : b(k,195) = b(k,195) - lu(k,2330) * b(k,237)
1440 315224452 : b(k,194) = b(k,194) - lu(k,2329) * b(k,237)
1441 315224452 : b(k,193) = b(k,193) - lu(k,2328) * b(k,237)
1442 315224452 : b(k,192) = b(k,192) - lu(k,2327) * b(k,237)
1443 315224452 : b(k,191) = b(k,191) - lu(k,2326) * b(k,237)
1444 315224452 : b(k,185) = b(k,185) - lu(k,2325) * b(k,237)
1445 315224452 : b(k,183) = b(k,183) - lu(k,2324) * b(k,237)
1446 315224452 : b(k,182) = b(k,182) - lu(k,2323) * b(k,237)
1447 315224452 : b(k,181) = b(k,181) - lu(k,2322) * b(k,237)
1448 315224452 : b(k,180) = b(k,180) - lu(k,2321) * b(k,237)
1449 315224452 : b(k,175) = b(k,175) - lu(k,2320) * b(k,237)
1450 315224452 : b(k,174) = b(k,174) - lu(k,2319) * b(k,237)
1451 315224452 : b(k,172) = b(k,172) - lu(k,2318) * b(k,237)
1452 315224452 : b(k,170) = b(k,170) - lu(k,2317) * b(k,237)
1453 315224452 : b(k,169) = b(k,169) - lu(k,2316) * b(k,237)
1454 315224452 : b(k,168) = b(k,168) - lu(k,2315) * b(k,237)
1455 315224452 : b(k,167) = b(k,167) - lu(k,2314) * b(k,237)
1456 315224452 : b(k,166) = b(k,166) - lu(k,2313) * b(k,237)
1457 315224452 : b(k,165) = b(k,165) - lu(k,2312) * b(k,237)
1458 315224452 : b(k,164) = b(k,164) - lu(k,2311) * b(k,237)
1459 315224452 : b(k,162) = b(k,162) - lu(k,2310) * b(k,237)
1460 315224452 : b(k,161) = b(k,161) - lu(k,2309) * b(k,237)
1461 315224452 : b(k,160) = b(k,160) - lu(k,2308) * b(k,237)
1462 315224452 : b(k,159) = b(k,159) - lu(k,2307) * b(k,237)
1463 315224452 : b(k,157) = b(k,157) - lu(k,2306) * b(k,237)
1464 315224452 : b(k,156) = b(k,156) - lu(k,2305) * b(k,237)
1465 315224452 : b(k,154) = b(k,154) - lu(k,2304) * b(k,237)
1466 315224452 : b(k,152) = b(k,152) - lu(k,2303) * b(k,237)
1467 315224452 : b(k,147) = b(k,147) - lu(k,2302) * b(k,237)
1468 315224452 : b(k,142) = b(k,142) - lu(k,2301) * b(k,237)
1469 315224452 : b(k,139) = b(k,139) - lu(k,2300) * b(k,237)
1470 315224452 : b(k,137) = b(k,137) - lu(k,2299) * b(k,237)
1471 315224452 : b(k,135) = b(k,135) - lu(k,2298) * b(k,237)
1472 315224452 : b(k,134) = b(k,134) - lu(k,2297) * b(k,237)
1473 315224452 : b(k,132) = b(k,132) - lu(k,2296) * b(k,237)
1474 315224452 : b(k,131) = b(k,131) - lu(k,2295) * b(k,237)
1475 315224452 : b(k,130) = b(k,130) - lu(k,2294) * b(k,237)
1476 315224452 : b(k,129) = b(k,129) - lu(k,2293) * b(k,237)
1477 315224452 : b(k,128) = b(k,128) - lu(k,2292) * b(k,237)
1478 315224452 : b(k,127) = b(k,127) - lu(k,2291) * b(k,237)
1479 315224452 : b(k,126) = b(k,126) - lu(k,2290) * b(k,237)
1480 315224452 : b(k,125) = b(k,125) - lu(k,2289) * b(k,237)
1481 315224452 : b(k,123) = b(k,123) - lu(k,2288) * b(k,237)
1482 315224452 : b(k,122) = b(k,122) - lu(k,2287) * b(k,237)
1483 315224452 : b(k,121) = b(k,121) - lu(k,2286) * b(k,237)
1484 315224452 : b(k,120) = b(k,120) - lu(k,2285) * b(k,237)
1485 315224452 : b(k,119) = b(k,119) - lu(k,2284) * b(k,237)
1486 315224452 : b(k,118) = b(k,118) - lu(k,2283) * b(k,237)
1487 315224452 : b(k,112) = b(k,112) - lu(k,2282) * b(k,237)
1488 315224452 : b(k,111) = b(k,111) - lu(k,2281) * b(k,237)
1489 315224452 : b(k,110) = b(k,110) - lu(k,2280) * b(k,237)
1490 315224452 : b(k,106) = b(k,106) - lu(k,2279) * b(k,237)
1491 315224452 : b(k,105) = b(k,105) - lu(k,2278) * b(k,237)
1492 315224452 : b(k,96) = b(k,96) - lu(k,2277) * b(k,237)
1493 315224452 : b(k,81) = b(k,81) - lu(k,2276) * b(k,237)
1494 315224452 : b(k,63) = b(k,63) - lu(k,2275) * b(k,237)
1495 315224452 : b(k,53) = b(k,53) - lu(k,2274) * b(k,237)
1496 315224452 : b(k,52) = b(k,52) - lu(k,2273) * b(k,237)
1497 315224452 : b(k,51) = b(k,51) - lu(k,2272) * b(k,237)
1498 315224452 : b(k,49) = b(k,49) - lu(k,2271) * b(k,237)
1499 315224452 : b(k,48) = b(k,48) - lu(k,2270) * b(k,237)
1500 315224452 : b(k,47) = b(k,47) - lu(k,2269) * b(k,237)
1501 315224452 : b(k,46) = b(k,46) - lu(k,2268) * b(k,237)
1502 315224452 : b(k,43) = b(k,43) - lu(k,2267) * b(k,237)
1503 315224452 : b(k,42) = b(k,42) - lu(k,2266) * b(k,237)
1504 315224452 : b(k,41) = b(k,41) - lu(k,2265) * b(k,237)
1505 315224452 : b(k,40) = b(k,40) - lu(k,2264) * b(k,237)
1506 315224452 : b(k,39) = b(k,39) - lu(k,2263) * b(k,237)
1507 315224452 : b(k,236) = b(k,236) * lu(k,2259)
1508 315224452 : b(k,235) = b(k,235) - lu(k,2258) * b(k,236)
1509 315224452 : b(k,234) = b(k,234) - lu(k,2257) * b(k,236)
1510 315224452 : b(k,233) = b(k,233) - lu(k,2256) * b(k,236)
1511 315224452 : b(k,232) = b(k,232) - lu(k,2255) * b(k,236)
1512 315224452 : b(k,231) = b(k,231) - lu(k,2254) * b(k,236)
1513 315224452 : b(k,230) = b(k,230) - lu(k,2253) * b(k,236)
1514 315224452 : b(k,229) = b(k,229) - lu(k,2252) * b(k,236)
1515 315224452 : b(k,228) = b(k,228) - lu(k,2251) * b(k,236)
1516 315224452 : b(k,227) = b(k,227) - lu(k,2250) * b(k,236)
1517 315224452 : b(k,226) = b(k,226) - lu(k,2249) * b(k,236)
1518 315224452 : b(k,225) = b(k,225) - lu(k,2248) * b(k,236)
1519 315224452 : b(k,224) = b(k,224) - lu(k,2247) * b(k,236)
1520 315224452 : b(k,223) = b(k,223) - lu(k,2246) * b(k,236)
1521 315224452 : b(k,222) = b(k,222) - lu(k,2245) * b(k,236)
1522 315224452 : b(k,221) = b(k,221) - lu(k,2244) * b(k,236)
1523 315224452 : b(k,220) = b(k,220) - lu(k,2243) * b(k,236)
1524 315224452 : b(k,218) = b(k,218) - lu(k,2242) * b(k,236)
1525 315224452 : b(k,217) = b(k,217) - lu(k,2241) * b(k,236)
1526 315224452 : b(k,201) = b(k,201) - lu(k,2240) * b(k,236)
1527 315224452 : b(k,184) = b(k,184) - lu(k,2239) * b(k,236)
1528 325174377 : b(k,99) = b(k,99) - lu(k,2238) * b(k,236)
1529 : end do
1530 9949925 : end subroutine lu_slv07
1531 9949925 : subroutine lu_slv08( avec_len, lu, b )
1532 : use shr_kind_mod, only : r8 => shr_kind_r8
1533 : use chem_mods, only : clscnt4, nzcnt
1534 : implicit none
1535 : !-----------------------------------------------------------------------
1536 : ! ... Dummy args
1537 : !-----------------------------------------------------------------------
1538 : integer, intent(in) :: avec_len
1539 : real(r8), intent(in) :: lu(veclen,max(1,nzcnt))
1540 : real(r8), intent(inout) :: b(veclen,clscnt4)
1541 : !-----------------------------------------------------------------------
1542 : ! ... Local variables
1543 : !-----------------------------------------------------------------------
1544 : integer :: k
1545 : !-----------------------------------------------------------------------
1546 : ! ... solve L * y = b
1547 : !-----------------------------------------------------------------------
1548 325174377 : do k = 1,avec_len
1549 315224452 : b(k,235) = b(k,235) * lu(k,2233)
1550 315224452 : b(k,234) = b(k,234) - lu(k,2232) * b(k,235)
1551 315224452 : b(k,233) = b(k,233) - lu(k,2231) * b(k,235)
1552 315224452 : b(k,232) = b(k,232) - lu(k,2230) * b(k,235)
1553 315224452 : b(k,231) = b(k,231) - lu(k,2229) * b(k,235)
1554 315224452 : b(k,230) = b(k,230) - lu(k,2228) * b(k,235)
1555 315224452 : b(k,229) = b(k,229) - lu(k,2227) * b(k,235)
1556 315224452 : b(k,228) = b(k,228) - lu(k,2226) * b(k,235)
1557 315224452 : b(k,227) = b(k,227) - lu(k,2225) * b(k,235)
1558 315224452 : b(k,226) = b(k,226) - lu(k,2224) * b(k,235)
1559 315224452 : b(k,225) = b(k,225) - lu(k,2223) * b(k,235)
1560 315224452 : b(k,224) = b(k,224) - lu(k,2222) * b(k,235)
1561 315224452 : b(k,223) = b(k,223) - lu(k,2221) * b(k,235)
1562 315224452 : b(k,222) = b(k,222) - lu(k,2220) * b(k,235)
1563 315224452 : b(k,221) = b(k,221) - lu(k,2219) * b(k,235)
1564 315224452 : b(k,219) = b(k,219) - lu(k,2218) * b(k,235)
1565 315224452 : b(k,217) = b(k,217) - lu(k,2217) * b(k,235)
1566 315224452 : b(k,216) = b(k,216) - lu(k,2216) * b(k,235)
1567 315224452 : b(k,215) = b(k,215) - lu(k,2215) * b(k,235)
1568 315224452 : b(k,214) = b(k,214) - lu(k,2214) * b(k,235)
1569 315224452 : b(k,213) = b(k,213) - lu(k,2213) * b(k,235)
1570 315224452 : b(k,212) = b(k,212) - lu(k,2212) * b(k,235)
1571 315224452 : b(k,211) = b(k,211) - lu(k,2211) * b(k,235)
1572 315224452 : b(k,210) = b(k,210) - lu(k,2210) * b(k,235)
1573 315224452 : b(k,209) = b(k,209) - lu(k,2209) * b(k,235)
1574 315224452 : b(k,208) = b(k,208) - lu(k,2208) * b(k,235)
1575 315224452 : b(k,206) = b(k,206) - lu(k,2207) * b(k,235)
1576 315224452 : b(k,205) = b(k,205) - lu(k,2206) * b(k,235)
1577 315224452 : b(k,203) = b(k,203) - lu(k,2205) * b(k,235)
1578 315224452 : b(k,202) = b(k,202) - lu(k,2204) * b(k,235)
1579 315224452 : b(k,201) = b(k,201) - lu(k,2203) * b(k,235)
1580 315224452 : b(k,200) = b(k,200) - lu(k,2202) * b(k,235)
1581 315224452 : b(k,199) = b(k,199) - lu(k,2201) * b(k,235)
1582 315224452 : b(k,198) = b(k,198) - lu(k,2200) * b(k,235)
1583 315224452 : b(k,197) = b(k,197) - lu(k,2199) * b(k,235)
1584 315224452 : b(k,195) = b(k,195) - lu(k,2198) * b(k,235)
1585 315224452 : b(k,194) = b(k,194) - lu(k,2197) * b(k,235)
1586 315224452 : b(k,193) = b(k,193) - lu(k,2196) * b(k,235)
1587 315224452 : b(k,192) = b(k,192) - lu(k,2195) * b(k,235)
1588 315224452 : b(k,191) = b(k,191) - lu(k,2194) * b(k,235)
1589 315224452 : b(k,190) = b(k,190) - lu(k,2193) * b(k,235)
1590 315224452 : b(k,188) = b(k,188) - lu(k,2192) * b(k,235)
1591 315224452 : b(k,186) = b(k,186) - lu(k,2191) * b(k,235)
1592 315224452 : b(k,183) = b(k,183) - lu(k,2190) * b(k,235)
1593 315224452 : b(k,179) = b(k,179) - lu(k,2189) * b(k,235)
1594 315224452 : b(k,170) = b(k,170) - lu(k,2188) * b(k,235)
1595 315224452 : b(k,163) = b(k,163) - lu(k,2187) * b(k,235)
1596 315224452 : b(k,161) = b(k,161) - lu(k,2186) * b(k,235)
1597 315224452 : b(k,154) = b(k,154) - lu(k,2185) * b(k,235)
1598 315224452 : b(k,143) = b(k,143) - lu(k,2184) * b(k,235)
1599 315224452 : b(k,131) = b(k,131) - lu(k,2183) * b(k,235)
1600 315224452 : b(k,116) = b(k,116) - lu(k,2182) * b(k,235)
1601 315224452 : b(k,85) = b(k,85) - lu(k,2181) * b(k,235)
1602 315224452 : b(k,76) = b(k,76) - lu(k,2180) * b(k,235)
1603 315224452 : b(k,75) = b(k,75) - lu(k,2179) * b(k,235)
1604 315224452 : b(k,43) = b(k,43) - lu(k,2178) * b(k,235)
1605 315224452 : b(k,42) = b(k,42) - lu(k,2177) * b(k,235)
1606 315224452 : b(k,41) = b(k,41) - lu(k,2176) * b(k,235)
1607 315224452 : b(k,40) = b(k,40) - lu(k,2175) * b(k,235)
1608 315224452 : b(k,39) = b(k,39) - lu(k,2174) * b(k,235)
1609 315224452 : b(k,234) = b(k,234) * lu(k,2168)
1610 315224452 : b(k,233) = b(k,233) - lu(k,2167) * b(k,234)
1611 315224452 : b(k,232) = b(k,232) - lu(k,2166) * b(k,234)
1612 315224452 : b(k,231) = b(k,231) - lu(k,2165) * b(k,234)
1613 315224452 : b(k,230) = b(k,230) - lu(k,2164) * b(k,234)
1614 315224452 : b(k,229) = b(k,229) - lu(k,2163) * b(k,234)
1615 315224452 : b(k,228) = b(k,228) - lu(k,2162) * b(k,234)
1616 315224452 : b(k,227) = b(k,227) - lu(k,2161) * b(k,234)
1617 315224452 : b(k,226) = b(k,226) - lu(k,2160) * b(k,234)
1618 315224452 : b(k,225) = b(k,225) - lu(k,2159) * b(k,234)
1619 315224452 : b(k,224) = b(k,224) - lu(k,2158) * b(k,234)
1620 315224452 : b(k,223) = b(k,223) - lu(k,2157) * b(k,234)
1621 315224452 : b(k,222) = b(k,222) - lu(k,2156) * b(k,234)
1622 315224452 : b(k,221) = b(k,221) - lu(k,2155) * b(k,234)
1623 315224452 : b(k,220) = b(k,220) - lu(k,2154) * b(k,234)
1624 315224452 : b(k,218) = b(k,218) - lu(k,2153) * b(k,234)
1625 315224452 : b(k,217) = b(k,217) - lu(k,2152) * b(k,234)
1626 315224452 : b(k,216) = b(k,216) - lu(k,2151) * b(k,234)
1627 315224452 : b(k,210) = b(k,210) - lu(k,2150) * b(k,234)
1628 315224452 : b(k,206) = b(k,206) - lu(k,2149) * b(k,234)
1629 315224452 : b(k,202) = b(k,202) - lu(k,2148) * b(k,234)
1630 315224452 : b(k,201) = b(k,201) - lu(k,2147) * b(k,234)
1631 315224452 : b(k,200) = b(k,200) - lu(k,2146) * b(k,234)
1632 315224452 : b(k,195) = b(k,195) - lu(k,2145) * b(k,234)
1633 315224452 : b(k,192) = b(k,192) - lu(k,2144) * b(k,234)
1634 315224452 : b(k,187) = b(k,187) - lu(k,2143) * b(k,234)
1635 315224452 : b(k,185) = b(k,185) - lu(k,2142) * b(k,234)
1636 315224452 : b(k,183) = b(k,183) - lu(k,2141) * b(k,234)
1637 315224452 : b(k,181) = b(k,181) - lu(k,2140) * b(k,234)
1638 315224452 : b(k,175) = b(k,175) - lu(k,2139) * b(k,234)
1639 315224452 : b(k,170) = b(k,170) - lu(k,2138) * b(k,234)
1640 315224452 : b(k,150) = b(k,150) - lu(k,2137) * b(k,234)
1641 315224452 : b(k,147) = b(k,147) - lu(k,2136) * b(k,234)
1642 315224452 : b(k,143) = b(k,143) - lu(k,2135) * b(k,234)
1643 315224452 : b(k,133) = b(k,133) - lu(k,2134) * b(k,234)
1644 315224452 : b(k,117) = b(k,117) - lu(k,2133) * b(k,234)
1645 315224452 : b(k,115) = b(k,115) - lu(k,2132) * b(k,234)
1646 315224452 : b(k,103) = b(k,103) - lu(k,2131) * b(k,234)
1647 315224452 : b(k,102) = b(k,102) - lu(k,2130) * b(k,234)
1648 315224452 : b(k,100) = b(k,100) - lu(k,2129) * b(k,234)
1649 315224452 : b(k,74) = b(k,74) - lu(k,2128) * b(k,234)
1650 315224452 : b(k,233) = b(k,233) * lu(k,2121)
1651 315224452 : b(k,232) = b(k,232) - lu(k,2120) * b(k,233)
1652 315224452 : b(k,231) = b(k,231) - lu(k,2119) * b(k,233)
1653 315224452 : b(k,230) = b(k,230) - lu(k,2118) * b(k,233)
1654 315224452 : b(k,229) = b(k,229) - lu(k,2117) * b(k,233)
1655 315224452 : b(k,228) = b(k,228) - lu(k,2116) * b(k,233)
1656 315224452 : b(k,227) = b(k,227) - lu(k,2115) * b(k,233)
1657 315224452 : b(k,226) = b(k,226) - lu(k,2114) * b(k,233)
1658 315224452 : b(k,225) = b(k,225) - lu(k,2113) * b(k,233)
1659 315224452 : b(k,224) = b(k,224) - lu(k,2112) * b(k,233)
1660 315224452 : b(k,223) = b(k,223) - lu(k,2111) * b(k,233)
1661 315224452 : b(k,222) = b(k,222) - lu(k,2110) * b(k,233)
1662 315224452 : b(k,221) = b(k,221) - lu(k,2109) * b(k,233)
1663 315224452 : b(k,220) = b(k,220) - lu(k,2108) * b(k,233)
1664 315224452 : b(k,218) = b(k,218) - lu(k,2107) * b(k,233)
1665 315224452 : b(k,232) = b(k,232) * lu(k,2099)
1666 315224452 : b(k,231) = b(k,231) - lu(k,2098) * b(k,232)
1667 315224452 : b(k,230) = b(k,230) - lu(k,2097) * b(k,232)
1668 315224452 : b(k,229) = b(k,229) - lu(k,2096) * b(k,232)
1669 315224452 : b(k,228) = b(k,228) - lu(k,2095) * b(k,232)
1670 315224452 : b(k,227) = b(k,227) - lu(k,2094) * b(k,232)
1671 315224452 : b(k,226) = b(k,226) - lu(k,2093) * b(k,232)
1672 315224452 : b(k,225) = b(k,225) - lu(k,2092) * b(k,232)
1673 315224452 : b(k,224) = b(k,224) - lu(k,2091) * b(k,232)
1674 315224452 : b(k,223) = b(k,223) - lu(k,2090) * b(k,232)
1675 315224452 : b(k,222) = b(k,222) - lu(k,2089) * b(k,232)
1676 315224452 : b(k,221) = b(k,221) - lu(k,2088) * b(k,232)
1677 315224452 : b(k,220) = b(k,220) - lu(k,2087) * b(k,232)
1678 315224452 : b(k,219) = b(k,219) - lu(k,2086) * b(k,232)
1679 315224452 : b(k,218) = b(k,218) - lu(k,2085) * b(k,232)
1680 315224452 : b(k,217) = b(k,217) - lu(k,2084) * b(k,232)
1681 315224452 : b(k,201) = b(k,201) - lu(k,2083) * b(k,232)
1682 315224452 : b(k,189) = b(k,189) - lu(k,2082) * b(k,232)
1683 315224452 : b(k,187) = b(k,187) - lu(k,2081) * b(k,232)
1684 315224452 : b(k,186) = b(k,186) - lu(k,2080) * b(k,232)
1685 315224452 : b(k,178) = b(k,178) - lu(k,2079) * b(k,232)
1686 315224452 : b(k,177) = b(k,177) - lu(k,2078) * b(k,232)
1687 315224452 : b(k,176) = b(k,176) - lu(k,2077) * b(k,232)
1688 315224452 : b(k,175) = b(k,175) - lu(k,2076) * b(k,232)
1689 315224452 : b(k,174) = b(k,174) - lu(k,2075) * b(k,232)
1690 315224452 : b(k,171) = b(k,171) - lu(k,2074) * b(k,232)
1691 315224452 : b(k,169) = b(k,169) - lu(k,2073) * b(k,232)
1692 315224452 : b(k,163) = b(k,163) - lu(k,2072) * b(k,232)
1693 315224452 : b(k,158) = b(k,158) - lu(k,2071) * b(k,232)
1694 315224452 : b(k,147) = b(k,147) - lu(k,2070) * b(k,232)
1695 315224452 : b(k,145) = b(k,145) - lu(k,2069) * b(k,232)
1696 315224452 : b(k,140) = b(k,140) - lu(k,2068) * b(k,232)
1697 315224452 : b(k,136) = b(k,136) - lu(k,2067) * b(k,232)
1698 315224452 : b(k,114) = b(k,114) - lu(k,2066) * b(k,232)
1699 315224452 : b(k,76) = b(k,76) - lu(k,2065) * b(k,232)
1700 315224452 : b(k,75) = b(k,75) - lu(k,2064) * b(k,232)
1701 315224452 : b(k,231) = b(k,231) * lu(k,2055)
1702 315224452 : b(k,230) = b(k,230) - lu(k,2054) * b(k,231)
1703 315224452 : b(k,229) = b(k,229) - lu(k,2053) * b(k,231)
1704 315224452 : b(k,228) = b(k,228) - lu(k,2052) * b(k,231)
1705 315224452 : b(k,227) = b(k,227) - lu(k,2051) * b(k,231)
1706 315224452 : b(k,226) = b(k,226) - lu(k,2050) * b(k,231)
1707 315224452 : b(k,225) = b(k,225) - lu(k,2049) * b(k,231)
1708 315224452 : b(k,224) = b(k,224) - lu(k,2048) * b(k,231)
1709 315224452 : b(k,223) = b(k,223) - lu(k,2047) * b(k,231)
1710 315224452 : b(k,222) = b(k,222) - lu(k,2046) * b(k,231)
1711 315224452 : b(k,221) = b(k,221) - lu(k,2045) * b(k,231)
1712 315224452 : b(k,220) = b(k,220) - lu(k,2044) * b(k,231)
1713 315224452 : b(k,218) = b(k,218) - lu(k,2043) * b(k,231)
1714 315224452 : b(k,184) = b(k,184) - lu(k,2042) * b(k,231)
1715 315224452 : b(k,99) = b(k,99) - lu(k,2041) * b(k,231)
1716 315224452 : b(k,230) = b(k,230) * lu(k,2031)
1717 315224452 : b(k,229) = b(k,229) - lu(k,2030) * b(k,230)
1718 315224452 : b(k,228) = b(k,228) - lu(k,2029) * b(k,230)
1719 315224452 : b(k,227) = b(k,227) - lu(k,2028) * b(k,230)
1720 315224452 : b(k,226) = b(k,226) - lu(k,2027) * b(k,230)
1721 315224452 : b(k,225) = b(k,225) - lu(k,2026) * b(k,230)
1722 315224452 : b(k,224) = b(k,224) - lu(k,2025) * b(k,230)
1723 315224452 : b(k,223) = b(k,223) - lu(k,2024) * b(k,230)
1724 315224452 : b(k,222) = b(k,222) - lu(k,2023) * b(k,230)
1725 315224452 : b(k,221) = b(k,221) - lu(k,2022) * b(k,230)
1726 315224452 : b(k,220) = b(k,220) - lu(k,2021) * b(k,230)
1727 315224452 : b(k,218) = b(k,218) - lu(k,2020) * b(k,230)
1728 315224452 : b(k,184) = b(k,184) - lu(k,2019) * b(k,230)
1729 315224452 : b(k,174) = b(k,174) - lu(k,2018) * b(k,230)
1730 315224452 : b(k,150) = b(k,150) - lu(k,2017) * b(k,230)
1731 315224452 : b(k,133) = b(k,133) - lu(k,2016) * b(k,230)
1732 315224452 : b(k,117) = b(k,117) - lu(k,2015) * b(k,230)
1733 315224452 : b(k,101) = b(k,101) - lu(k,2014) * b(k,230)
1734 315224452 : b(k,93) = b(k,93) - lu(k,2013) * b(k,230)
1735 315224452 : b(k,92) = b(k,92) - lu(k,2012) * b(k,230)
1736 315224452 : b(k,91) = b(k,91) - lu(k,2011) * b(k,230)
1737 315224452 : b(k,90) = b(k,90) - lu(k,2010) * b(k,230)
1738 315224452 : b(k,89) = b(k,89) - lu(k,2009) * b(k,230)
1739 315224452 : b(k,77) = b(k,77) - lu(k,2008) * b(k,230)
1740 315224452 : b(k,76) = b(k,76) - lu(k,2007) * b(k,230)
1741 315224452 : b(k,70) = b(k,70) - lu(k,2006) * b(k,230)
1742 315224452 : b(k,69) = b(k,69) - lu(k,2005) * b(k,230)
1743 315224452 : b(k,68) = b(k,68) - lu(k,2004) * b(k,230)
1744 315224452 : b(k,67) = b(k,67) - lu(k,2003) * b(k,230)
1745 315224452 : b(k,61) = b(k,61) - lu(k,2002) * b(k,230)
1746 315224452 : b(k,60) = b(k,60) - lu(k,2001) * b(k,230)
1747 315224452 : b(k,59) = b(k,59) - lu(k,2000) * b(k,230)
1748 315224452 : b(k,58) = b(k,58) - lu(k,1999) * b(k,230)
1749 325174377 : b(k,56) = b(k,56) - lu(k,1998) * b(k,230)
1750 : end do
1751 9949925 : end subroutine lu_slv08
1752 9949925 : subroutine lu_slv09( avec_len, lu, b )
1753 : use shr_kind_mod, only : r8 => shr_kind_r8
1754 : use chem_mods, only : clscnt4, nzcnt
1755 : implicit none
1756 : !-----------------------------------------------------------------------
1757 : ! ... Dummy args
1758 : !-----------------------------------------------------------------------
1759 : integer, intent(in) :: avec_len
1760 : real(r8), intent(in) :: lu(veclen,max(1,nzcnt))
1761 : real(r8), intent(inout) :: b(veclen,clscnt4)
1762 : !-----------------------------------------------------------------------
1763 : ! ... Local variables
1764 : !-----------------------------------------------------------------------
1765 : integer :: k
1766 : !-----------------------------------------------------------------------
1767 : ! ... solve L * y = b
1768 : !-----------------------------------------------------------------------
1769 325174377 : do k = 1,avec_len
1770 315224452 : b(k,229) = b(k,229) * lu(k,1987)
1771 315224452 : b(k,228) = b(k,228) - lu(k,1986) * b(k,229)
1772 315224452 : b(k,227) = b(k,227) - lu(k,1985) * b(k,229)
1773 315224452 : b(k,226) = b(k,226) - lu(k,1984) * b(k,229)
1774 315224452 : b(k,225) = b(k,225) - lu(k,1983) * b(k,229)
1775 315224452 : b(k,224) = b(k,224) - lu(k,1982) * b(k,229)
1776 315224452 : b(k,223) = b(k,223) - lu(k,1981) * b(k,229)
1777 315224452 : b(k,222) = b(k,222) - lu(k,1980) * b(k,229)
1778 315224452 : b(k,221) = b(k,221) - lu(k,1979) * b(k,229)
1779 315224452 : b(k,220) = b(k,220) - lu(k,1978) * b(k,229)
1780 315224452 : b(k,218) = b(k,218) - lu(k,1977) * b(k,229)
1781 315224452 : b(k,217) = b(k,217) - lu(k,1976) * b(k,229)
1782 315224452 : b(k,201) = b(k,201) - lu(k,1975) * b(k,229)
1783 315224452 : b(k,183) = b(k,183) - lu(k,1974) * b(k,229)
1784 315224452 : b(k,174) = b(k,174) - lu(k,1973) * b(k,229)
1785 315224452 : b(k,126) = b(k,126) - lu(k,1972) * b(k,229)
1786 315224452 : b(k,228) = b(k,228) * lu(k,1960)
1787 315224452 : b(k,227) = b(k,227) - lu(k,1959) * b(k,228)
1788 315224452 : b(k,226) = b(k,226) - lu(k,1958) * b(k,228)
1789 315224452 : b(k,225) = b(k,225) - lu(k,1957) * b(k,228)
1790 315224452 : b(k,224) = b(k,224) - lu(k,1956) * b(k,228)
1791 315224452 : b(k,223) = b(k,223) - lu(k,1955) * b(k,228)
1792 315224452 : b(k,222) = b(k,222) - lu(k,1954) * b(k,228)
1793 315224452 : b(k,221) = b(k,221) - lu(k,1953) * b(k,228)
1794 315224452 : b(k,220) = b(k,220) - lu(k,1952) * b(k,228)
1795 315224452 : b(k,219) = b(k,219) - lu(k,1951) * b(k,228)
1796 315224452 : b(k,217) = b(k,217) - lu(k,1950) * b(k,228)
1797 315224452 : b(k,216) = b(k,216) - lu(k,1949) * b(k,228)
1798 315224452 : b(k,215) = b(k,215) - lu(k,1948) * b(k,228)
1799 315224452 : b(k,210) = b(k,210) - lu(k,1947) * b(k,228)
1800 315224452 : b(k,205) = b(k,205) - lu(k,1946) * b(k,228)
1801 315224452 : b(k,201) = b(k,201) - lu(k,1945) * b(k,228)
1802 315224452 : b(k,192) = b(k,192) - lu(k,1944) * b(k,228)
1803 315224452 : b(k,189) = b(k,189) - lu(k,1943) * b(k,228)
1804 315224452 : b(k,187) = b(k,187) - lu(k,1942) * b(k,228)
1805 315224452 : b(k,186) = b(k,186) - lu(k,1941) * b(k,228)
1806 315224452 : b(k,182) = b(k,182) - lu(k,1940) * b(k,228)
1807 315224452 : b(k,169) = b(k,169) - lu(k,1939) * b(k,228)
1808 315224452 : b(k,168) = b(k,168) - lu(k,1938) * b(k,228)
1809 315224452 : b(k,165) = b(k,165) - lu(k,1937) * b(k,228)
1810 315224452 : b(k,154) = b(k,154) - lu(k,1936) * b(k,228)
1811 315224452 : b(k,151) = b(k,151) - lu(k,1935) * b(k,228)
1812 315224452 : b(k,145) = b(k,145) - lu(k,1934) * b(k,228)
1813 315224452 : b(k,141) = b(k,141) - lu(k,1933) * b(k,228)
1814 315224452 : b(k,139) = b(k,139) - lu(k,1932) * b(k,228)
1815 315224452 : b(k,132) = b(k,132) - lu(k,1931) * b(k,228)
1816 315224452 : b(k,131) = b(k,131) - lu(k,1930) * b(k,228)
1817 315224452 : b(k,97) = b(k,97) - lu(k,1929) * b(k,228)
1818 315224452 : b(k,93) = b(k,93) - lu(k,1928) * b(k,228)
1819 315224452 : b(k,66) = b(k,66) - lu(k,1927) * b(k,228)
1820 315224452 : b(k,227) = b(k,227) * lu(k,1914)
1821 315224452 : b(k,226) = b(k,226) - lu(k,1913) * b(k,227)
1822 315224452 : b(k,225) = b(k,225) - lu(k,1912) * b(k,227)
1823 315224452 : b(k,224) = b(k,224) - lu(k,1911) * b(k,227)
1824 315224452 : b(k,223) = b(k,223) - lu(k,1910) * b(k,227)
1825 315224452 : b(k,222) = b(k,222) - lu(k,1909) * b(k,227)
1826 315224452 : b(k,221) = b(k,221) - lu(k,1908) * b(k,227)
1827 315224452 : b(k,217) = b(k,217) - lu(k,1907) * b(k,227)
1828 315224452 : b(k,216) = b(k,216) - lu(k,1906) * b(k,227)
1829 315224452 : b(k,215) = b(k,215) - lu(k,1905) * b(k,227)
1830 315224452 : b(k,214) = b(k,214) - lu(k,1904) * b(k,227)
1831 315224452 : b(k,213) = b(k,213) - lu(k,1903) * b(k,227)
1832 315224452 : b(k,212) = b(k,212) - lu(k,1902) * b(k,227)
1833 315224452 : b(k,211) = b(k,211) - lu(k,1901) * b(k,227)
1834 315224452 : b(k,210) = b(k,210) - lu(k,1900) * b(k,227)
1835 315224452 : b(k,209) = b(k,209) - lu(k,1899) * b(k,227)
1836 315224452 : b(k,208) = b(k,208) - lu(k,1898) * b(k,227)
1837 315224452 : b(k,207) = b(k,207) - lu(k,1897) * b(k,227)
1838 315224452 : b(k,206) = b(k,206) - lu(k,1896) * b(k,227)
1839 315224452 : b(k,205) = b(k,205) - lu(k,1895) * b(k,227)
1840 315224452 : b(k,204) = b(k,204) - lu(k,1894) * b(k,227)
1841 315224452 : b(k,203) = b(k,203) - lu(k,1893) * b(k,227)
1842 315224452 : b(k,202) = b(k,202) - lu(k,1892) * b(k,227)
1843 315224452 : b(k,201) = b(k,201) - lu(k,1891) * b(k,227)
1844 315224452 : b(k,200) = b(k,200) - lu(k,1890) * b(k,227)
1845 315224452 : b(k,198) = b(k,198) - lu(k,1889) * b(k,227)
1846 315224452 : b(k,197) = b(k,197) - lu(k,1888) * b(k,227)
1847 315224452 : b(k,196) = b(k,196) - lu(k,1887) * b(k,227)
1848 315224452 : b(k,195) = b(k,195) - lu(k,1886) * b(k,227)
1849 315224452 : b(k,194) = b(k,194) - lu(k,1885) * b(k,227)
1850 315224452 : b(k,193) = b(k,193) - lu(k,1884) * b(k,227)
1851 315224452 : b(k,192) = b(k,192) - lu(k,1883) * b(k,227)
1852 315224452 : b(k,191) = b(k,191) - lu(k,1882) * b(k,227)
1853 315224452 : b(k,189) = b(k,189) - lu(k,1881) * b(k,227)
1854 315224452 : b(k,185) = b(k,185) - lu(k,1880) * b(k,227)
1855 315224452 : b(k,183) = b(k,183) - lu(k,1879) * b(k,227)
1856 315224452 : b(k,182) = b(k,182) - lu(k,1878) * b(k,227)
1857 315224452 : b(k,181) = b(k,181) - lu(k,1877) * b(k,227)
1858 315224452 : b(k,180) = b(k,180) - lu(k,1876) * b(k,227)
1859 315224452 : b(k,178) = b(k,178) - lu(k,1875) * b(k,227)
1860 315224452 : b(k,177) = b(k,177) - lu(k,1874) * b(k,227)
1861 315224452 : b(k,176) = b(k,176) - lu(k,1873) * b(k,227)
1862 315224452 : b(k,172) = b(k,172) - lu(k,1872) * b(k,227)
1863 315224452 : b(k,170) = b(k,170) - lu(k,1871) * b(k,227)
1864 315224452 : b(k,168) = b(k,168) - lu(k,1870) * b(k,227)
1865 315224452 : b(k,167) = b(k,167) - lu(k,1869) * b(k,227)
1866 315224452 : b(k,166) = b(k,166) - lu(k,1868) * b(k,227)
1867 315224452 : b(k,165) = b(k,165) - lu(k,1867) * b(k,227)
1868 315224452 : b(k,160) = b(k,160) - lu(k,1866) * b(k,227)
1869 315224452 : b(k,155) = b(k,155) - lu(k,1865) * b(k,227)
1870 315224452 : b(k,154) = b(k,154) - lu(k,1864) * b(k,227)
1871 315224452 : b(k,153) = b(k,153) - lu(k,1863) * b(k,227)
1872 315224452 : b(k,149) = b(k,149) - lu(k,1862) * b(k,227)
1873 315224452 : b(k,148) = b(k,148) - lu(k,1861) * b(k,227)
1874 315224452 : b(k,146) = b(k,146) - lu(k,1860) * b(k,227)
1875 315224452 : b(k,144) = b(k,144) - lu(k,1859) * b(k,227)
1876 315224452 : b(k,139) = b(k,139) - lu(k,1858) * b(k,227)
1877 315224452 : b(k,138) = b(k,138) - lu(k,1857) * b(k,227)
1878 315224452 : b(k,137) = b(k,137) - lu(k,1856) * b(k,227)
1879 315224452 : b(k,135) = b(k,135) - lu(k,1855) * b(k,227)
1880 315224452 : b(k,131) = b(k,131) - lu(k,1854) * b(k,227)
1881 315224452 : b(k,130) = b(k,130) - lu(k,1853) * b(k,227)
1882 315224452 : b(k,129) = b(k,129) - lu(k,1852) * b(k,227)
1883 315224452 : b(k,127) = b(k,127) - lu(k,1851) * b(k,227)
1884 315224452 : b(k,126) = b(k,126) - lu(k,1850) * b(k,227)
1885 315224452 : b(k,125) = b(k,125) - lu(k,1849) * b(k,227)
1886 315224452 : b(k,108) = b(k,108) - lu(k,1848) * b(k,227)
1887 315224452 : b(k,107) = b(k,107) - lu(k,1847) * b(k,227)
1888 315224452 : b(k,98) = b(k,98) - lu(k,1846) * b(k,227)
1889 315224452 : b(k,88) = b(k,88) - lu(k,1845) * b(k,227)
1890 315224452 : b(k,87) = b(k,87) - lu(k,1844) * b(k,227)
1891 315224452 : b(k,53) = b(k,53) - lu(k,1843) * b(k,227)
1892 315224452 : b(k,52) = b(k,52) - lu(k,1842) * b(k,227)
1893 315224452 : b(k,51) = b(k,51) - lu(k,1841) * b(k,227)
1894 315224452 : b(k,49) = b(k,49) - lu(k,1840) * b(k,227)
1895 315224452 : b(k,48) = b(k,48) - lu(k,1839) * b(k,227)
1896 315224452 : b(k,47) = b(k,47) - lu(k,1838) * b(k,227)
1897 315224452 : b(k,46) = b(k,46) - lu(k,1837) * b(k,227)
1898 315224452 : b(k,43) = b(k,43) - lu(k,1836) * b(k,227)
1899 315224452 : b(k,42) = b(k,42) - lu(k,1835) * b(k,227)
1900 315224452 : b(k,41) = b(k,41) - lu(k,1834) * b(k,227)
1901 315224452 : b(k,40) = b(k,40) - lu(k,1833) * b(k,227)
1902 315224452 : b(k,39) = b(k,39) - lu(k,1832) * b(k,227)
1903 315224452 : b(k,226) = b(k,226) * lu(k,1818)
1904 315224452 : b(k,225) = b(k,225) - lu(k,1817) * b(k,226)
1905 315224452 : b(k,224) = b(k,224) - lu(k,1816) * b(k,226)
1906 315224452 : b(k,223) = b(k,223) - lu(k,1815) * b(k,226)
1907 315224452 : b(k,222) = b(k,222) - lu(k,1814) * b(k,226)
1908 315224452 : b(k,221) = b(k,221) - lu(k,1813) * b(k,226)
1909 315224452 : b(k,220) = b(k,220) - lu(k,1812) * b(k,226)
1910 315224452 : b(k,219) = b(k,219) - lu(k,1811) * b(k,226)
1911 315224452 : b(k,218) = b(k,218) - lu(k,1810) * b(k,226)
1912 315224452 : b(k,217) = b(k,217) - lu(k,1809) * b(k,226)
1913 315224452 : b(k,216) = b(k,216) - lu(k,1808) * b(k,226)
1914 315224452 : b(k,215) = b(k,215) - lu(k,1807) * b(k,226)
1915 315224452 : b(k,214) = b(k,214) - lu(k,1806) * b(k,226)
1916 315224452 : b(k,213) = b(k,213) - lu(k,1805) * b(k,226)
1917 315224452 : b(k,212) = b(k,212) - lu(k,1804) * b(k,226)
1918 315224452 : b(k,211) = b(k,211) - lu(k,1803) * b(k,226)
1919 315224452 : b(k,210) = b(k,210) - lu(k,1802) * b(k,226)
1920 315224452 : b(k,209) = b(k,209) - lu(k,1801) * b(k,226)
1921 315224452 : b(k,208) = b(k,208) - lu(k,1800) * b(k,226)
1922 315224452 : b(k,207) = b(k,207) - lu(k,1799) * b(k,226)
1923 315224452 : b(k,206) = b(k,206) - lu(k,1798) * b(k,226)
1924 315224452 : b(k,205) = b(k,205) - lu(k,1797) * b(k,226)
1925 315224452 : b(k,204) = b(k,204) - lu(k,1796) * b(k,226)
1926 315224452 : b(k,203) = b(k,203) - lu(k,1795) * b(k,226)
1927 315224452 : b(k,202) = b(k,202) - lu(k,1794) * b(k,226)
1928 315224452 : b(k,201) = b(k,201) - lu(k,1793) * b(k,226)
1929 315224452 : b(k,200) = b(k,200) - lu(k,1792) * b(k,226)
1930 315224452 : b(k,199) = b(k,199) - lu(k,1791) * b(k,226)
1931 315224452 : b(k,198) = b(k,198) - lu(k,1790) * b(k,226)
1932 315224452 : b(k,197) = b(k,197) - lu(k,1789) * b(k,226)
1933 315224452 : b(k,196) = b(k,196) - lu(k,1788) * b(k,226)
1934 315224452 : b(k,195) = b(k,195) - lu(k,1787) * b(k,226)
1935 315224452 : b(k,194) = b(k,194) - lu(k,1786) * b(k,226)
1936 315224452 : b(k,193) = b(k,193) - lu(k,1785) * b(k,226)
1937 315224452 : b(k,192) = b(k,192) - lu(k,1784) * b(k,226)
1938 315224452 : b(k,191) = b(k,191) - lu(k,1783) * b(k,226)
1939 315224452 : b(k,190) = b(k,190) - lu(k,1782) * b(k,226)
1940 315224452 : b(k,189) = b(k,189) - lu(k,1781) * b(k,226)
1941 315224452 : b(k,188) = b(k,188) - lu(k,1780) * b(k,226)
1942 315224452 : b(k,187) = b(k,187) - lu(k,1779) * b(k,226)
1943 315224452 : b(k,186) = b(k,186) - lu(k,1778) * b(k,226)
1944 315224452 : b(k,185) = b(k,185) - lu(k,1777) * b(k,226)
1945 315224452 : b(k,184) = b(k,184) - lu(k,1776) * b(k,226)
1946 315224452 : b(k,183) = b(k,183) - lu(k,1775) * b(k,226)
1947 315224452 : b(k,182) = b(k,182) - lu(k,1774) * b(k,226)
1948 315224452 : b(k,181) = b(k,181) - lu(k,1773) * b(k,226)
1949 315224452 : b(k,180) = b(k,180) - lu(k,1772) * b(k,226)
1950 315224452 : b(k,179) = b(k,179) - lu(k,1771) * b(k,226)
1951 315224452 : b(k,175) = b(k,175) - lu(k,1770) * b(k,226)
1952 315224452 : b(k,174) = b(k,174) - lu(k,1769) * b(k,226)
1953 315224452 : b(k,173) = b(k,173) - lu(k,1768) * b(k,226)
1954 315224452 : b(k,172) = b(k,172) - lu(k,1767) * b(k,226)
1955 315224452 : b(k,170) = b(k,170) - lu(k,1766) * b(k,226)
1956 315224452 : b(k,168) = b(k,168) - lu(k,1765) * b(k,226)
1957 315224452 : b(k,167) = b(k,167) - lu(k,1764) * b(k,226)
1958 315224452 : b(k,166) = b(k,166) - lu(k,1763) * b(k,226)
1959 315224452 : b(k,165) = b(k,165) - lu(k,1762) * b(k,226)
1960 315224452 : b(k,164) = b(k,164) - lu(k,1761) * b(k,226)
1961 315224452 : b(k,163) = b(k,163) - lu(k,1760) * b(k,226)
1962 315224452 : b(k,162) = b(k,162) - lu(k,1759) * b(k,226)
1963 315224452 : b(k,161) = b(k,161) - lu(k,1758) * b(k,226)
1964 315224452 : b(k,160) = b(k,160) - lu(k,1757) * b(k,226)
1965 315224452 : b(k,159) = b(k,159) - lu(k,1756) * b(k,226)
1966 315224452 : b(k,157) = b(k,157) - lu(k,1755) * b(k,226)
1967 315224452 : b(k,156) = b(k,156) - lu(k,1754) * b(k,226)
1968 315224452 : b(k,155) = b(k,155) - lu(k,1753) * b(k,226)
1969 315224452 : b(k,154) = b(k,154) - lu(k,1752) * b(k,226)
1970 315224452 : b(k,153) = b(k,153) - lu(k,1751) * b(k,226)
1971 315224452 : b(k,152) = b(k,152) - lu(k,1750) * b(k,226)
1972 315224452 : b(k,151) = b(k,151) - lu(k,1749) * b(k,226)
1973 315224452 : b(k,150) = b(k,150) - lu(k,1748) * b(k,226)
1974 315224452 : b(k,149) = b(k,149) - lu(k,1747) * b(k,226)
1975 315224452 : b(k,148) = b(k,148) - lu(k,1746) * b(k,226)
1976 315224452 : b(k,147) = b(k,147) - lu(k,1745) * b(k,226)
1977 315224452 : b(k,146) = b(k,146) - lu(k,1744) * b(k,226)
1978 315224452 : b(k,144) = b(k,144) - lu(k,1743) * b(k,226)
1979 315224452 : b(k,143) = b(k,143) - lu(k,1742) * b(k,226)
1980 315224452 : b(k,142) = b(k,142) - lu(k,1741) * b(k,226)
1981 315224452 : b(k,141) = b(k,141) - lu(k,1740) * b(k,226)
1982 315224452 : b(k,139) = b(k,139) - lu(k,1739) * b(k,226)
1983 315224452 : b(k,138) = b(k,138) - lu(k,1738) * b(k,226)
1984 315224452 : b(k,137) = b(k,137) - lu(k,1737) * b(k,226)
1985 315224452 : b(k,135) = b(k,135) - lu(k,1736) * b(k,226)
1986 315224452 : b(k,134) = b(k,134) - lu(k,1735) * b(k,226)
1987 315224452 : b(k,133) = b(k,133) - lu(k,1734) * b(k,226)
1988 315224452 : b(k,132) = b(k,132) - lu(k,1733) * b(k,226)
1989 315224452 : b(k,131) = b(k,131) - lu(k,1732) * b(k,226)
1990 315224452 : b(k,130) = b(k,130) - lu(k,1731) * b(k,226)
1991 315224452 : b(k,128) = b(k,128) - lu(k,1730) * b(k,226)
1992 315224452 : b(k,127) = b(k,127) - lu(k,1729) * b(k,226)
1993 315224452 : b(k,125) = b(k,125) - lu(k,1728) * b(k,226)
1994 315224452 : b(k,124) = b(k,124) - lu(k,1727) * b(k,226)
1995 315224452 : b(k,123) = b(k,123) - lu(k,1726) * b(k,226)
1996 315224452 : b(k,122) = b(k,122) - lu(k,1725) * b(k,226)
1997 315224452 : b(k,121) = b(k,121) - lu(k,1724) * b(k,226)
1998 315224452 : b(k,120) = b(k,120) - lu(k,1723) * b(k,226)
1999 315224452 : b(k,119) = b(k,119) - lu(k,1722) * b(k,226)
2000 315224452 : b(k,118) = b(k,118) - lu(k,1721) * b(k,226)
2001 315224452 : b(k,117) = b(k,117) - lu(k,1720) * b(k,226)
2002 315224452 : b(k,115) = b(k,115) - lu(k,1719) * b(k,226)
2003 315224452 : b(k,114) = b(k,114) - lu(k,1718) * b(k,226)
2004 315224452 : b(k,113) = b(k,113) - lu(k,1717) * b(k,226)
2005 315224452 : b(k,112) = b(k,112) - lu(k,1716) * b(k,226)
2006 315224452 : b(k,111) = b(k,111) - lu(k,1715) * b(k,226)
2007 315224452 : b(k,110) = b(k,110) - lu(k,1714) * b(k,226)
2008 315224452 : b(k,107) = b(k,107) - lu(k,1713) * b(k,226)
2009 315224452 : b(k,106) = b(k,106) - lu(k,1712) * b(k,226)
2010 315224452 : b(k,105) = b(k,105) - lu(k,1711) * b(k,226)
2011 315224452 : b(k,104) = b(k,104) - lu(k,1710) * b(k,226)
2012 315224452 : b(k,103) = b(k,103) - lu(k,1709) * b(k,226)
2013 315224452 : b(k,102) = b(k,102) - lu(k,1708) * b(k,226)
2014 315224452 : b(k,101) = b(k,101) - lu(k,1707) * b(k,226)
2015 315224452 : b(k,100) = b(k,100) - lu(k,1706) * b(k,226)
2016 315224452 : b(k,96) = b(k,96) - lu(k,1705) * b(k,226)
2017 315224452 : b(k,95) = b(k,95) - lu(k,1704) * b(k,226)
2018 315224452 : b(k,94) = b(k,94) - lu(k,1703) * b(k,226)
2019 315224452 : b(k,92) = b(k,92) - lu(k,1702) * b(k,226)
2020 315224452 : b(k,91) = b(k,91) - lu(k,1701) * b(k,226)
2021 315224452 : b(k,90) = b(k,90) - lu(k,1700) * b(k,226)
2022 315224452 : b(k,89) = b(k,89) - lu(k,1699) * b(k,226)
2023 315224452 : b(k,88) = b(k,88) - lu(k,1698) * b(k,226)
2024 315224452 : b(k,87) = b(k,87) - lu(k,1697) * b(k,226)
2025 315224452 : b(k,85) = b(k,85) - lu(k,1696) * b(k,226)
2026 315224452 : b(k,84) = b(k,84) - lu(k,1695) * b(k,226)
2027 315224452 : b(k,83) = b(k,83) - lu(k,1694) * b(k,226)
2028 315224452 : b(k,82) = b(k,82) - lu(k,1693) * b(k,226)
2029 315224452 : b(k,81) = b(k,81) - lu(k,1692) * b(k,226)
2030 315224452 : b(k,80) = b(k,80) - lu(k,1691) * b(k,226)
2031 315224452 : b(k,79) = b(k,79) - lu(k,1690) * b(k,226)
2032 315224452 : b(k,78) = b(k,78) - lu(k,1689) * b(k,226)
2033 315224452 : b(k,77) = b(k,77) - lu(k,1688) * b(k,226)
2034 315224452 : b(k,73) = b(k,73) - lu(k,1687) * b(k,226)
2035 315224452 : b(k,72) = b(k,72) - lu(k,1686) * b(k,226)
2036 315224452 : b(k,71) = b(k,71) - lu(k,1685) * b(k,226)
2037 315224452 : b(k,65) = b(k,65) - lu(k,1684) * b(k,226)
2038 315224452 : b(k,62) = b(k,62) - lu(k,1683) * b(k,226)
2039 315224452 : b(k,57) = b(k,57) - lu(k,1682) * b(k,226)
2040 315224452 : b(k,55) = b(k,55) - lu(k,1681) * b(k,226)
2041 315224452 : b(k,53) = b(k,53) - lu(k,1680) * b(k,226)
2042 315224452 : b(k,52) = b(k,52) - lu(k,1679) * b(k,226)
2043 315224452 : b(k,51) = b(k,51) - lu(k,1678) * b(k,226)
2044 315224452 : b(k,50) = b(k,50) - lu(k,1677) * b(k,226)
2045 315224452 : b(k,49) = b(k,49) - lu(k,1676) * b(k,226)
2046 315224452 : b(k,48) = b(k,48) - lu(k,1675) * b(k,226)
2047 315224452 : b(k,47) = b(k,47) - lu(k,1674) * b(k,226)
2048 315224452 : b(k,46) = b(k,46) - lu(k,1673) * b(k,226)
2049 315224452 : b(k,45) = b(k,45) - lu(k,1672) * b(k,226)
2050 315224452 : b(k,43) = b(k,43) - lu(k,1671) * b(k,226)
2051 315224452 : b(k,42) = b(k,42) - lu(k,1670) * b(k,226)
2052 315224452 : b(k,41) = b(k,41) - lu(k,1669) * b(k,226)
2053 315224452 : b(k,40) = b(k,40) - lu(k,1668) * b(k,226)
2054 325174377 : b(k,39) = b(k,39) - lu(k,1667) * b(k,226)
2055 : end do
2056 9949925 : end subroutine lu_slv09
2057 9949925 : subroutine lu_slv10( avec_len, lu, b )
2058 : use shr_kind_mod, only : r8 => shr_kind_r8
2059 : use chem_mods, only : clscnt4, nzcnt
2060 : implicit none
2061 : !-----------------------------------------------------------------------
2062 : ! ... Dummy args
2063 : !-----------------------------------------------------------------------
2064 : integer, intent(in) :: avec_len
2065 : real(r8), intent(in) :: lu(veclen,max(1,nzcnt))
2066 : real(r8), intent(inout) :: b(veclen,clscnt4)
2067 : !-----------------------------------------------------------------------
2068 : ! ... Local variables
2069 : !-----------------------------------------------------------------------
2070 : integer :: k
2071 : !-----------------------------------------------------------------------
2072 : ! ... solve L * y = b
2073 : !-----------------------------------------------------------------------
2074 325174377 : do k = 1,avec_len
2075 315224452 : b(k,225) = b(k,225) * lu(k,1652)
2076 315224452 : b(k,224) = b(k,224) - lu(k,1651) * b(k,225)
2077 315224452 : b(k,223) = b(k,223) - lu(k,1650) * b(k,225)
2078 315224452 : b(k,222) = b(k,222) - lu(k,1649) * b(k,225)
2079 315224452 : b(k,221) = b(k,221) - lu(k,1648) * b(k,225)
2080 315224452 : b(k,219) = b(k,219) - lu(k,1647) * b(k,225)
2081 315224452 : b(k,217) = b(k,217) - lu(k,1646) * b(k,225)
2082 315224452 : b(k,216) = b(k,216) - lu(k,1645) * b(k,225)
2083 315224452 : b(k,215) = b(k,215) - lu(k,1644) * b(k,225)
2084 315224452 : b(k,214) = b(k,214) - lu(k,1643) * b(k,225)
2085 315224452 : b(k,213) = b(k,213) - lu(k,1642) * b(k,225)
2086 315224452 : b(k,212) = b(k,212) - lu(k,1641) * b(k,225)
2087 315224452 : b(k,211) = b(k,211) - lu(k,1640) * b(k,225)
2088 315224452 : b(k,210) = b(k,210) - lu(k,1639) * b(k,225)
2089 315224452 : b(k,209) = b(k,209) - lu(k,1638) * b(k,225)
2090 315224452 : b(k,208) = b(k,208) - lu(k,1637) * b(k,225)
2091 315224452 : b(k,207) = b(k,207) - lu(k,1636) * b(k,225)
2092 315224452 : b(k,206) = b(k,206) - lu(k,1635) * b(k,225)
2093 315224452 : b(k,205) = b(k,205) - lu(k,1634) * b(k,225)
2094 315224452 : b(k,204) = b(k,204) - lu(k,1633) * b(k,225)
2095 315224452 : b(k,203) = b(k,203) - lu(k,1632) * b(k,225)
2096 315224452 : b(k,202) = b(k,202) - lu(k,1631) * b(k,225)
2097 315224452 : b(k,201) = b(k,201) - lu(k,1630) * b(k,225)
2098 315224452 : b(k,200) = b(k,200) - lu(k,1629) * b(k,225)
2099 315224452 : b(k,199) = b(k,199) - lu(k,1628) * b(k,225)
2100 315224452 : b(k,198) = b(k,198) - lu(k,1627) * b(k,225)
2101 315224452 : b(k,197) = b(k,197) - lu(k,1626) * b(k,225)
2102 315224452 : b(k,196) = b(k,196) - lu(k,1625) * b(k,225)
2103 315224452 : b(k,195) = b(k,195) - lu(k,1624) * b(k,225)
2104 315224452 : b(k,194) = b(k,194) - lu(k,1623) * b(k,225)
2105 315224452 : b(k,193) = b(k,193) - lu(k,1622) * b(k,225)
2106 315224452 : b(k,192) = b(k,192) - lu(k,1621) * b(k,225)
2107 315224452 : b(k,191) = b(k,191) - lu(k,1620) * b(k,225)
2108 315224452 : b(k,190) = b(k,190) - lu(k,1619) * b(k,225)
2109 315224452 : b(k,188) = b(k,188) - lu(k,1618) * b(k,225)
2110 315224452 : b(k,186) = b(k,186) - lu(k,1617) * b(k,225)
2111 315224452 : b(k,183) = b(k,183) - lu(k,1616) * b(k,225)
2112 315224452 : b(k,180) = b(k,180) - lu(k,1615) * b(k,225)
2113 315224452 : b(k,179) = b(k,179) - lu(k,1614) * b(k,225)
2114 315224452 : b(k,146) = b(k,146) - lu(k,1613) * b(k,225)
2115 315224452 : b(k,113) = b(k,113) - lu(k,1612) * b(k,225)
2116 315224452 : b(k,107) = b(k,107) - lu(k,1611) * b(k,225)
2117 315224452 : b(k,104) = b(k,104) - lu(k,1610) * b(k,225)
2118 315224452 : b(k,97) = b(k,97) - lu(k,1609) * b(k,225)
2119 315224452 : b(k,43) = b(k,43) - lu(k,1608) * b(k,225)
2120 315224452 : b(k,42) = b(k,42) - lu(k,1607) * b(k,225)
2121 315224452 : b(k,224) = b(k,224) * lu(k,1592)
2122 315224452 : b(k,223) = b(k,223) - lu(k,1591) * b(k,224)
2123 315224452 : b(k,222) = b(k,222) - lu(k,1590) * b(k,224)
2124 315224452 : b(k,221) = b(k,221) - lu(k,1589) * b(k,224)
2125 315224452 : b(k,220) = b(k,220) - lu(k,1588) * b(k,224)
2126 315224452 : b(k,219) = b(k,219) - lu(k,1587) * b(k,224)
2127 315224452 : b(k,187) = b(k,187) - lu(k,1586) * b(k,224)
2128 315224452 : b(k,186) = b(k,186) - lu(k,1585) * b(k,224)
2129 315224452 : b(k,175) = b(k,175) - lu(k,1584) * b(k,224)
2130 315224452 : b(k,109) = b(k,109) - lu(k,1583) * b(k,224)
2131 315224452 : b(k,86) = b(k,86) - lu(k,1582) * b(k,224)
2132 315224452 : b(k,74) = b(k,74) - lu(k,1581) * b(k,224)
2133 315224452 : b(k,54) = b(k,54) - lu(k,1580) * b(k,224)
2134 315224452 : b(k,223) = b(k,223) * lu(k,1565)
2135 315224452 : b(k,222) = b(k,222) - lu(k,1564) * b(k,223)
2136 315224452 : b(k,221) = b(k,221) - lu(k,1563) * b(k,223)
2137 315224452 : b(k,220) = b(k,220) - lu(k,1562) * b(k,223)
2138 315224452 : b(k,219) = b(k,219) - lu(k,1561) * b(k,223)
2139 315224452 : b(k,186) = b(k,186) - lu(k,1560) * b(k,223)
2140 315224452 : b(k,169) = b(k,169) - lu(k,1559) * b(k,223)
2141 315224452 : b(k,145) = b(k,145) - lu(k,1558) * b(k,223)
2142 315224452 : b(k,109) = b(k,109) - lu(k,1557) * b(k,223)
2143 315224452 : b(k,86) = b(k,86) - lu(k,1556) * b(k,223)
2144 315224452 : b(k,222) = b(k,222) * lu(k,1542)
2145 315224452 : b(k,221) = b(k,221) - lu(k,1541) * b(k,222)
2146 315224452 : b(k,219) = b(k,219) - lu(k,1540) * b(k,222)
2147 315224452 : b(k,217) = b(k,217) - lu(k,1539) * b(k,222)
2148 315224452 : b(k,202) = b(k,202) - lu(k,1538) * b(k,222)
2149 315224452 : b(k,201) = b(k,201) - lu(k,1537) * b(k,222)
2150 315224452 : b(k,189) = b(k,189) - lu(k,1536) * b(k,222)
2151 315224452 : b(k,186) = b(k,186) - lu(k,1535) * b(k,222)
2152 315224452 : b(k,178) = b(k,178) - lu(k,1534) * b(k,222)
2153 315224452 : b(k,177) = b(k,177) - lu(k,1533) * b(k,222)
2154 315224452 : b(k,176) = b(k,176) - lu(k,1532) * b(k,222)
2155 315224452 : b(k,171) = b(k,171) - lu(k,1531) * b(k,222)
2156 315224452 : b(k,163) = b(k,163) - lu(k,1530) * b(k,222)
2157 315224452 : b(k,158) = b(k,158) - lu(k,1529) * b(k,222)
2158 315224452 : b(k,140) = b(k,140) - lu(k,1528) * b(k,222)
2159 315224452 : b(k,136) = b(k,136) - lu(k,1527) * b(k,222)
2160 315224452 : b(k,129) = b(k,129) - lu(k,1526) * b(k,222)
2161 315224452 : b(k,76) = b(k,76) - lu(k,1525) * b(k,222)
2162 315224452 : b(k,75) = b(k,75) - lu(k,1524) * b(k,222)
2163 315224452 : b(k,221) = b(k,221) * lu(k,1511)
2164 315224452 : b(k,217) = b(k,217) - lu(k,1510) * b(k,221)
2165 315224452 : b(k,201) = b(k,201) - lu(k,1509) * b(k,221)
2166 315224452 : b(k,174) = b(k,174) - lu(k,1508) * b(k,221)
2167 315224452 : b(k,220) = b(k,220) * lu(k,1495)
2168 315224452 : b(k,187) = b(k,187) - lu(k,1494) * b(k,220)
2169 315224452 : b(k,175) = b(k,175) - lu(k,1493) * b(k,220)
2170 315224452 : b(k,169) = b(k,169) - lu(k,1492) * b(k,220)
2171 315224452 : b(k,86) = b(k,86) - lu(k,1491) * b(k,220)
2172 315224452 : b(k,74) = b(k,74) - lu(k,1490) * b(k,220)
2173 315224452 : b(k,219) = b(k,219) * lu(k,1476)
2174 315224452 : b(k,186) = b(k,186) - lu(k,1475) * b(k,219)
2175 315224452 : b(k,163) = b(k,163) - lu(k,1474) * b(k,219)
2176 315224452 : b(k,109) = b(k,109) - lu(k,1473) * b(k,219)
2177 315224452 : b(k,218) = b(k,218) * lu(k,1461)
2178 315224452 : b(k,184) = b(k,184) - lu(k,1460) * b(k,218)
2179 315224452 : b(k,99) = b(k,99) - lu(k,1459) * b(k,218)
2180 315224452 : b(k,217) = b(k,217) * lu(k,1450)
2181 315224452 : b(k,201) = b(k,201) - lu(k,1449) * b(k,217)
2182 315224452 : b(k,189) = b(k,189) - lu(k,1448) * b(k,217)
2183 315224452 : b(k,178) = b(k,178) - lu(k,1447) * b(k,217)
2184 315224452 : b(k,177) = b(k,177) - lu(k,1446) * b(k,217)
2185 315224452 : b(k,176) = b(k,176) - lu(k,1445) * b(k,217)
2186 315224452 : b(k,171) = b(k,171) - lu(k,1444) * b(k,217)
2187 315224452 : b(k,76) = b(k,76) - lu(k,1443) * b(k,217)
2188 315224452 : b(k,75) = b(k,75) - lu(k,1442) * b(k,217)
2189 315224452 : b(k,216) = b(k,216) * lu(k,1428)
2190 315224452 : b(k,215) = b(k,215) - lu(k,1427) * b(k,216)
2191 315224452 : b(k,214) = b(k,214) - lu(k,1426) * b(k,216)
2192 315224452 : b(k,213) = b(k,213) - lu(k,1425) * b(k,216)
2193 315224452 : b(k,212) = b(k,212) - lu(k,1424) * b(k,216)
2194 315224452 : b(k,211) = b(k,211) - lu(k,1423) * b(k,216)
2195 315224452 : b(k,210) = b(k,210) - lu(k,1422) * b(k,216)
2196 315224452 : b(k,209) = b(k,209) - lu(k,1421) * b(k,216)
2197 315224452 : b(k,208) = b(k,208) - lu(k,1420) * b(k,216)
2198 315224452 : b(k,207) = b(k,207) - lu(k,1419) * b(k,216)
2199 315224452 : b(k,205) = b(k,205) - lu(k,1418) * b(k,216)
2200 315224452 : b(k,202) = b(k,202) - lu(k,1417) * b(k,216)
2201 315224452 : b(k,201) = b(k,201) - lu(k,1416) * b(k,216)
2202 315224452 : b(k,196) = b(k,196) - lu(k,1415) * b(k,216)
2203 315224452 : b(k,192) = b(k,192) - lu(k,1414) * b(k,216)
2204 315224452 : b(k,161) = b(k,161) - lu(k,1413) * b(k,216)
2205 315224452 : b(k,152) = b(k,152) - lu(k,1412) * b(k,216)
2206 315224452 : b(k,141) = b(k,141) - lu(k,1411) * b(k,216)
2207 315224452 : b(k,107) = b(k,107) - lu(k,1410) * b(k,216)
2208 315224452 : b(k,215) = b(k,215) * lu(k,1397)
2209 315224452 : b(k,210) = b(k,210) - lu(k,1396) * b(k,215)
2210 315224452 : b(k,205) = b(k,205) - lu(k,1395) * b(k,215)
2211 315224452 : b(k,161) = b(k,161) - lu(k,1394) * b(k,215)
2212 315224452 : b(k,152) = b(k,152) - lu(k,1393) * b(k,215)
2213 315224452 : b(k,151) = b(k,151) - lu(k,1392) * b(k,215)
2214 315224452 : b(k,214) = b(k,214) * lu(k,1376)
2215 315224452 : b(k,213) = b(k,213) - lu(k,1375) * b(k,214)
2216 315224452 : b(k,210) = b(k,210) - lu(k,1374) * b(k,214)
2217 315224452 : b(k,205) = b(k,205) - lu(k,1373) * b(k,214)
2218 315224452 : b(k,201) = b(k,201) - lu(k,1372) * b(k,214)
2219 315224452 : b(k,200) = b(k,200) - lu(k,1371) * b(k,214)
2220 315224452 : b(k,199) = b(k,199) - lu(k,1370) * b(k,214)
2221 315224452 : b(k,183) = b(k,183) - lu(k,1369) * b(k,214)
2222 315224452 : b(k,213) = b(k,213) * lu(k,1356)
2223 315224452 : b(k,210) = b(k,210) - lu(k,1355) * b(k,213)
2224 315224452 : b(k,206) = b(k,206) - lu(k,1354) * b(k,213)
2225 315224452 : b(k,205) = b(k,205) - lu(k,1353) * b(k,213)
2226 315224452 : b(k,204) = b(k,204) - lu(k,1352) * b(k,213)
2227 315224452 : b(k,202) = b(k,202) - lu(k,1351) * b(k,213)
2228 315224452 : b(k,201) = b(k,201) - lu(k,1350) * b(k,213)
2229 315224452 : b(k,173) = b(k,173) - lu(k,1349) * b(k,213)
2230 315224452 : b(k,105) = b(k,105) - lu(k,1348) * b(k,213)
2231 315224452 : b(k,212) = b(k,212) * lu(k,1332)
2232 315224452 : b(k,210) = b(k,210) - lu(k,1331) * b(k,212)
2233 315224452 : b(k,209) = b(k,209) - lu(k,1330) * b(k,212)
2234 315224452 : b(k,207) = b(k,207) - lu(k,1329) * b(k,212)
2235 315224452 : b(k,206) = b(k,206) - lu(k,1328) * b(k,212)
2236 315224452 : b(k,205) = b(k,205) - lu(k,1327) * b(k,212)
2237 315224452 : b(k,204) = b(k,204) - lu(k,1326) * b(k,212)
2238 315224452 : b(k,202) = b(k,202) - lu(k,1325) * b(k,212)
2239 315224452 : b(k,201) = b(k,201) - lu(k,1324) * b(k,212)
2240 315224452 : b(k,192) = b(k,192) - lu(k,1323) * b(k,212)
2241 315224452 : b(k,182) = b(k,182) - lu(k,1322) * b(k,212)
2242 315224452 : b(k,180) = b(k,180) - lu(k,1321) * b(k,212)
2243 315224452 : b(k,173) = b(k,173) - lu(k,1320) * b(k,212)
2244 315224452 : b(k,162) = b(k,162) - lu(k,1319) * b(k,212)
2245 315224452 : b(k,149) = b(k,149) - lu(k,1318) * b(k,212)
2246 315224452 : b(k,144) = b(k,144) - lu(k,1317) * b(k,212)
2247 315224452 : b(k,107) = b(k,107) - lu(k,1316) * b(k,212)
2248 315224452 : b(k,84) = b(k,84) - lu(k,1315) * b(k,212)
2249 315224452 : b(k,211) = b(k,211) * lu(k,1299)
2250 315224452 : b(k,210) = b(k,210) - lu(k,1298) * b(k,211)
2251 315224452 : b(k,209) = b(k,209) - lu(k,1297) * b(k,211)
2252 315224452 : b(k,207) = b(k,207) - lu(k,1296) * b(k,211)
2253 315224452 : b(k,206) = b(k,206) - lu(k,1295) * b(k,211)
2254 315224452 : b(k,205) = b(k,205) - lu(k,1294) * b(k,211)
2255 315224452 : b(k,204) = b(k,204) - lu(k,1293) * b(k,211)
2256 315224452 : b(k,202) = b(k,202) - lu(k,1292) * b(k,211)
2257 315224452 : b(k,173) = b(k,173) - lu(k,1291) * b(k,211)
2258 315224452 : b(k,162) = b(k,162) - lu(k,1290) * b(k,211)
2259 315224452 : b(k,148) = b(k,148) - lu(k,1289) * b(k,211)
2260 315224452 : b(k,210) = b(k,210) * lu(k,1281)
2261 315224452 : b(k,201) = b(k,201) - lu(k,1280) * b(k,210)
2262 315224452 : b(k,209) = b(k,209) * lu(k,1269)
2263 315224452 : b(k,201) = b(k,201) - lu(k,1268) * b(k,209)
2264 315224452 : b(k,183) = b(k,183) - lu(k,1267) * b(k,209)
2265 315224452 : b(k,208) = b(k,208) * lu(k,1253)
2266 315224452 : b(k,207) = b(k,207) - lu(k,1252) * b(k,208)
2267 315224452 : b(k,202) = b(k,202) - lu(k,1251) * b(k,208)
2268 315224452 : b(k,201) = b(k,201) - lu(k,1250) * b(k,208)
2269 315224452 : b(k,196) = b(k,196) - lu(k,1249) * b(k,208)
2270 315224452 : b(k,180) = b(k,180) - lu(k,1248) * b(k,208)
2271 315224452 : b(k,173) = b(k,173) - lu(k,1247) * b(k,208)
2272 315224452 : b(k,162) = b(k,162) - lu(k,1246) * b(k,208)
2273 315224452 : b(k,124) = b(k,124) - lu(k,1245) * b(k,208)
2274 315224452 : b(k,120) = b(k,120) - lu(k,1244) * b(k,208)
2275 315224452 : b(k,207) = b(k,207) * lu(k,1233)
2276 315224452 : b(k,205) = b(k,205) - lu(k,1232) * b(k,207)
2277 315224452 : b(k,202) = b(k,202) - lu(k,1231) * b(k,207)
2278 315224452 : b(k,201) = b(k,201) - lu(k,1230) * b(k,207)
2279 315224452 : b(k,192) = b(k,192) - lu(k,1229) * b(k,207)
2280 315224452 : b(k,173) = b(k,173) - lu(k,1228) * b(k,207)
2281 325174377 : b(k,81) = b(k,81) - lu(k,1227) * b(k,207)
2282 : end do
2283 9949925 : end subroutine lu_slv10
2284 9949925 : subroutine lu_slv11( avec_len, lu, b )
2285 : use shr_kind_mod, only : r8 => shr_kind_r8
2286 : use chem_mods, only : clscnt4, nzcnt
2287 : implicit none
2288 : !-----------------------------------------------------------------------
2289 : ! ... Dummy args
2290 : !-----------------------------------------------------------------------
2291 : integer, intent(in) :: avec_len
2292 : real(r8), intent(in) :: lu(veclen,max(1,nzcnt))
2293 : real(r8), intent(inout) :: b(veclen,clscnt4)
2294 : !-----------------------------------------------------------------------
2295 : ! ... Local variables
2296 : !-----------------------------------------------------------------------
2297 : integer :: k
2298 : !-----------------------------------------------------------------------
2299 : ! ... solve L * y = b
2300 : !-----------------------------------------------------------------------
2301 325174377 : do k = 1,avec_len
2302 315224452 : b(k,206) = b(k,206) * lu(k,1217)
2303 315224452 : b(k,205) = b(k,205) - lu(k,1216) * b(k,206)
2304 315224452 : b(k,173) = b(k,173) - lu(k,1215) * b(k,206)
2305 315224452 : b(k,119) = b(k,119) - lu(k,1214) * b(k,206)
2306 315224452 : b(k,205) = b(k,205) * lu(k,1208)
2307 315224452 : b(k,204) = b(k,204) * lu(k,1195)
2308 315224452 : b(k,202) = b(k,202) - lu(k,1194) * b(k,204)
2309 315224452 : b(k,201) = b(k,201) - lu(k,1193) * b(k,204)
2310 315224452 : b(k,200) = b(k,200) - lu(k,1192) * b(k,204)
2311 315224452 : b(k,195) = b(k,195) - lu(k,1191) * b(k,204)
2312 315224452 : b(k,87) = b(k,87) - lu(k,1190) * b(k,204)
2313 315224452 : b(k,203) = b(k,203) * lu(k,1178)
2314 315224452 : b(k,202) = b(k,202) - lu(k,1177) * b(k,203)
2315 315224452 : b(k,201) = b(k,201) - lu(k,1176) * b(k,203)
2316 315224452 : b(k,198) = b(k,198) - lu(k,1175) * b(k,203)
2317 315224452 : b(k,195) = b(k,195) - lu(k,1174) * b(k,203)
2318 315224452 : b(k,173) = b(k,173) - lu(k,1173) * b(k,203)
2319 315224452 : b(k,156) = b(k,156) - lu(k,1172) * b(k,203)
2320 315224452 : b(k,87) = b(k,87) - lu(k,1171) * b(k,203)
2321 315224452 : b(k,202) = b(k,202) * lu(k,1166)
2322 315224452 : b(k,201) = b(k,201) - lu(k,1165) * b(k,202)
2323 315224452 : b(k,192) = b(k,192) - lu(k,1164) * b(k,202)
2324 315224452 : b(k,201) = b(k,201) * lu(k,1160)
2325 315224452 : b(k,200) = b(k,200) * lu(k,1151)
2326 315224452 : b(k,199) = b(k,199) * lu(k,1133)
2327 315224452 : b(k,183) = b(k,183) - lu(k,1132) * b(k,199)
2328 315224452 : b(k,180) = b(k,180) - lu(k,1131) * b(k,199)
2329 315224452 : b(k,172) = b(k,172) - lu(k,1130) * b(k,199)
2330 315224452 : b(k,161) = b(k,161) - lu(k,1129) * b(k,199)
2331 315224452 : b(k,198) = b(k,198) * lu(k,1119)
2332 315224452 : b(k,195) = b(k,195) - lu(k,1118) * b(k,198)
2333 315224452 : b(k,197) = b(k,197) * lu(k,1106)
2334 315224452 : b(k,194) = b(k,194) - lu(k,1105) * b(k,197)
2335 315224452 : b(k,173) = b(k,173) - lu(k,1104) * b(k,197)
2336 315224452 : b(k,153) = b(k,153) - lu(k,1103) * b(k,197)
2337 315224452 : b(k,118) = b(k,118) - lu(k,1102) * b(k,197)
2338 315224452 : b(k,196) = b(k,196) * lu(k,1090)
2339 315224452 : b(k,192) = b(k,192) - lu(k,1089) * b(k,196)
2340 315224452 : b(k,182) = b(k,182) - lu(k,1088) * b(k,196)
2341 315224452 : b(k,180) = b(k,180) - lu(k,1087) * b(k,196)
2342 315224452 : b(k,149) = b(k,149) - lu(k,1086) * b(k,196)
2343 315224452 : b(k,195) = b(k,195) * lu(k,1080)
2344 315224452 : b(k,194) = b(k,194) * lu(k,1072)
2345 315224452 : b(k,193) = b(k,193) * lu(k,1062)
2346 315224452 : b(k,173) = b(k,173) - lu(k,1061) * b(k,193)
2347 315224452 : b(k,153) = b(k,153) - lu(k,1060) * b(k,193)
2348 315224452 : b(k,134) = b(k,134) - lu(k,1059) * b(k,193)
2349 315224452 : b(k,192) = b(k,192) * lu(k,1054)
2350 315224452 : b(k,39) = b(k,39) - lu(k,1053) * b(k,192)
2351 315224452 : b(k,191) = b(k,191) * lu(k,1042)
2352 315224452 : b(k,181) = b(k,181) - lu(k,1041) * b(k,191)
2353 315224452 : b(k,160) = b(k,160) - lu(k,1040) * b(k,191)
2354 315224452 : b(k,159) = b(k,159) - lu(k,1039) * b(k,191)
2355 315224452 : b(k,155) = b(k,155) - lu(k,1038) * b(k,191)
2356 315224452 : b(k,138) = b(k,138) - lu(k,1037) * b(k,191)
2357 315224452 : b(k,190) = b(k,190) * lu(k,1018)
2358 315224452 : b(k,183) = b(k,183) - lu(k,1017) * b(k,190)
2359 315224452 : b(k,116) = b(k,116) - lu(k,1016) * b(k,190)
2360 315224452 : b(k,85) = b(k,85) - lu(k,1015) * b(k,190)
2361 315224452 : b(k,51) = b(k,51) - lu(k,1014) * b(k,190)
2362 315224452 : b(k,43) = b(k,43) - lu(k,1013) * b(k,190)
2363 315224452 : b(k,42) = b(k,42) - lu(k,1012) * b(k,190)
2364 315224452 : b(k,41) = b(k,41) - lu(k,1011) * b(k,190)
2365 315224452 : b(k,40) = b(k,40) - lu(k,1010) * b(k,190)
2366 315224452 : b(k,39) = b(k,39) - lu(k,1009) * b(k,190)
2367 315224452 : b(k,189) = b(k,189) * lu(k,999)
2368 315224452 : b(k,178) = b(k,178) - lu(k,998) * b(k,189)
2369 315224452 : b(k,177) = b(k,177) - lu(k,997) * b(k,189)
2370 315224452 : b(k,176) = b(k,176) - lu(k,996) * b(k,189)
2371 315224452 : b(k,171) = b(k,171) - lu(k,995) * b(k,189)
2372 315224452 : b(k,136) = b(k,136) - lu(k,994) * b(k,189)
2373 315224452 : b(k,93) = b(k,93) - lu(k,993) * b(k,189)
2374 315224452 : b(k,188) = b(k,188) * lu(k,974)
2375 315224452 : b(k,183) = b(k,183) - lu(k,973) * b(k,188)
2376 315224452 : b(k,116) = b(k,116) - lu(k,972) * b(k,188)
2377 315224452 : b(k,85) = b(k,85) - lu(k,971) * b(k,188)
2378 315224452 : b(k,46) = b(k,46) - lu(k,970) * b(k,188)
2379 315224452 : b(k,43) = b(k,43) - lu(k,969) * b(k,188)
2380 315224452 : b(k,42) = b(k,42) - lu(k,968) * b(k,188)
2381 315224452 : b(k,41) = b(k,41) - lu(k,967) * b(k,188)
2382 315224452 : b(k,40) = b(k,40) - lu(k,966) * b(k,188)
2383 315224452 : b(k,39) = b(k,39) - lu(k,965) * b(k,188)
2384 315224452 : b(k,187) = b(k,187) * lu(k,955)
2385 315224452 : b(k,175) = b(k,175) - lu(k,954) * b(k,187)
2386 315224452 : b(k,74) = b(k,74) - lu(k,953) * b(k,187)
2387 315224452 : b(k,186) = b(k,186) * lu(k,947)
2388 315224452 : b(k,94) = b(k,94) - lu(k,946) * b(k,186)
2389 315224452 : b(k,185) = b(k,185) * lu(k,935)
2390 315224452 : b(k,121) = b(k,121) - lu(k,934) * b(k,185)
2391 315224452 : b(k,184) = b(k,184) * lu(k,925)
2392 315224452 : b(k,99) = b(k,99) - lu(k,924) * b(k,184)
2393 315224452 : b(k,183) = b(k,183) * lu(k,919)
2394 315224452 : b(k,182) = b(k,182) * lu(k,910)
2395 315224452 : b(k,181) = b(k,181) * lu(k,901)
2396 315224452 : b(k,173) = b(k,173) - lu(k,900) * b(k,181)
2397 315224452 : b(k,112) = b(k,112) - lu(k,899) * b(k,181)
2398 315224452 : b(k,95) = b(k,95) - lu(k,898) * b(k,181)
2399 315224452 : b(k,180) = b(k,180) * lu(k,892)
2400 315224452 : b(k,179) = b(k,179) * lu(k,876)
2401 315224452 : b(k,48) = b(k,48) - lu(k,875) * b(k,179)
2402 315224452 : b(k,43) = b(k,43) - lu(k,874) * b(k,179)
2403 315224452 : b(k,42) = b(k,42) - lu(k,873) * b(k,179)
2404 315224452 : b(k,178) = b(k,178) * lu(k,865)
2405 315224452 : b(k,177) = b(k,177) - lu(k,864) * b(k,178)
2406 315224452 : b(k,176) = b(k,176) - lu(k,863) * b(k,178)
2407 315224452 : b(k,171) = b(k,171) - lu(k,862) * b(k,178)
2408 315224452 : b(k,158) = b(k,158) - lu(k,861) * b(k,178)
2409 315224452 : b(k,140) = b(k,140) - lu(k,860) * b(k,178)
2410 315224452 : b(k,177) = b(k,177) * lu(k,853)
2411 315224452 : b(k,176) = b(k,176) - lu(k,852) * b(k,177)
2412 315224452 : b(k,176) = b(k,176) * lu(k,845)
2413 315224452 : b(k,140) = b(k,140) - lu(k,844) * b(k,176)
2414 315224452 : b(k,175) = b(k,175) * lu(k,837)
2415 315224452 : b(k,74) = b(k,74) - lu(k,836) * b(k,175)
2416 315224452 : b(k,174) = b(k,174) * lu(k,828)
2417 315224452 : b(k,173) = b(k,173) * lu(k,824)
2418 315224452 : b(k,172) = b(k,172) * lu(k,814)
2419 315224452 : b(k,142) = b(k,142) - lu(k,813) * b(k,172)
2420 315224452 : b(k,171) = b(k,171) * lu(k,805)
2421 315224452 : b(k,170) = b(k,170) * lu(k,797)
2422 315224452 : b(k,129) = b(k,129) - lu(k,796) * b(k,170)
2423 315224452 : b(k,63) = b(k,63) - lu(k,795) * b(k,170)
2424 315224452 : b(k,169) = b(k,169) * lu(k,787)
2425 315224452 : b(k,86) = b(k,86) - lu(k,786) * b(k,169)
2426 315224452 : b(k,168) = b(k,168) * lu(k,778)
2427 315224452 : b(k,167) = b(k,167) * lu(k,767)
2428 315224452 : b(k,165) = b(k,165) - lu(k,766) * b(k,167)
2429 315224452 : b(k,164) = b(k,164) - lu(k,765) * b(k,167)
2430 315224452 : b(k,149) = b(k,149) - lu(k,764) * b(k,167)
2431 315224452 : b(k,130) = b(k,130) - lu(k,763) * b(k,167)
2432 315224452 : b(k,108) = b(k,108) - lu(k,762) * b(k,167)
2433 315224452 : b(k,98) = b(k,98) - lu(k,761) * b(k,167)
2434 315224452 : b(k,166) = b(k,166) * lu(k,751)
2435 315224452 : b(k,165) = b(k,165) - lu(k,750) * b(k,166)
2436 315224452 : b(k,157) = b(k,157) - lu(k,749) * b(k,166)
2437 315224452 : b(k,149) = b(k,149) - lu(k,748) * b(k,166)
2438 315224452 : b(k,130) = b(k,130) - lu(k,747) * b(k,166)
2439 315224452 : b(k,98) = b(k,98) - lu(k,746) * b(k,166)
2440 315224452 : b(k,165) = b(k,165) * lu(k,740)
2441 315224452 : b(k,164) = b(k,164) * lu(k,729)
2442 315224452 : b(k,149) = b(k,149) - lu(k,728) * b(k,164)
2443 315224452 : b(k,130) = b(k,130) - lu(k,727) * b(k,164)
2444 315224452 : b(k,108) = b(k,108) - lu(k,726) * b(k,164)
2445 315224452 : b(k,98) = b(k,98) - lu(k,725) * b(k,164)
2446 315224452 : b(k,163) = b(k,163) * lu(k,718)
2447 315224452 : b(k,162) = b(k,162) * lu(k,711)
2448 315224452 : b(k,65) = b(k,65) - lu(k,710) * b(k,162)
2449 315224452 : b(k,161) = b(k,161) * lu(k,705)
2450 315224452 : b(k,160) = b(k,160) * lu(k,698)
2451 315224452 : b(k,106) = b(k,106) - lu(k,697) * b(k,160)
2452 315224452 : b(k,159) = b(k,159) * lu(k,687)
2453 315224452 : b(k,138) = b(k,138) - lu(k,686) * b(k,159)
2454 315224452 : b(k,158) = b(k,158) * lu(k,676)
2455 315224452 : b(k,140) = b(k,140) - lu(k,675) * b(k,158)
2456 315224452 : b(k,157) = b(k,157) * lu(k,665)
2457 315224452 : b(k,149) = b(k,149) - lu(k,664) * b(k,157)
2458 315224452 : b(k,130) = b(k,130) - lu(k,663) * b(k,157)
2459 315224452 : b(k,98) = b(k,98) - lu(k,662) * b(k,157)
2460 315224452 : b(k,156) = b(k,156) * lu(k,652)
2461 315224452 : b(k,155) = b(k,155) * lu(k,642)
2462 315224452 : b(k,138) = b(k,138) - lu(k,641) * b(k,155)
2463 315224452 : b(k,154) = b(k,154) * lu(k,635)
2464 315224452 : b(k,131) = b(k,131) - lu(k,634) * b(k,154)
2465 315224452 : b(k,96) = b(k,96) - lu(k,633) * b(k,154)
2466 315224452 : b(k,153) = b(k,153) * lu(k,627)
2467 315224452 : b(k,152) = b(k,152) * lu(k,620)
2468 315224452 : b(k,151) = b(k,151) * lu(k,611)
2469 315224452 : b(k,150) = b(k,150) * lu(k,602)
2470 315224452 : b(k,149) = b(k,149) * lu(k,598)
2471 315224452 : b(k,148) = b(k,148) * lu(k,589)
2472 315224452 : b(k,147) = b(k,147) * lu(k,582)
2473 315224452 : b(k,146) = b(k,146) * lu(k,574)
2474 315224452 : b(k,145) = b(k,145) * lu(k,566)
2475 315224452 : b(k,144) = b(k,144) * lu(k,558)
2476 315224452 : b(k,143) = b(k,143) * lu(k,550)
2477 315224452 : b(k,142) = b(k,142) * lu(k,542)
2478 315224452 : b(k,141) = b(k,141) * lu(k,534)
2479 315224452 : b(k,140) = b(k,140) * lu(k,529)
2480 315224452 : b(k,139) = b(k,139) * lu(k,523)
2481 315224452 : b(k,66) = b(k,66) - lu(k,522) * b(k,139)
2482 315224452 : b(k,138) = b(k,138) * lu(k,517)
2483 315224452 : b(k,137) = b(k,137) * lu(k,510)
2484 315224452 : b(k,123) = b(k,123) - lu(k,509) * b(k,137)
2485 315224452 : b(k,136) = b(k,136) * lu(k,502)
2486 315224452 : b(k,135) = b(k,135) * lu(k,495)
2487 315224452 : b(k,130) = b(k,130) - lu(k,494) * b(k,135)
2488 315224452 : b(k,122) = b(k,122) - lu(k,493) * b(k,135)
2489 315224452 : b(k,134) = b(k,134) * lu(k,486)
2490 315224452 : b(k,133) = b(k,133) * lu(k,479)
2491 315224452 : b(k,132) = b(k,132) * lu(k,472)
2492 315224452 : b(k,131) = b(k,131) * lu(k,468)
2493 315224452 : b(k,130) = b(k,130) * lu(k,465)
2494 315224452 : b(k,129) = b(k,129) * lu(k,460)
2495 315224452 : b(k,128) = b(k,128) * lu(k,454)
2496 315224452 : b(k,127) = b(k,127) * lu(k,448)
2497 315224452 : b(k,110) = b(k,110) - lu(k,447) * b(k,127)
2498 315224452 : b(k,126) = b(k,126) * lu(k,441)
2499 315224452 : b(k,125) = b(k,125) * lu(k,435)
2500 315224452 : b(k,111) = b(k,111) - lu(k,434) * b(k,125)
2501 315224452 : b(k,88) = b(k,88) - lu(k,433) * b(k,125)
2502 325174377 : b(k,124) = b(k,124) * lu(k,427)
2503 : end do
2504 9949925 : end subroutine lu_slv11
2505 9949925 : subroutine lu_slv12( avec_len, lu, b )
2506 : use shr_kind_mod, only : r8 => shr_kind_r8
2507 : use chem_mods, only : clscnt4, nzcnt
2508 : implicit none
2509 : !-----------------------------------------------------------------------
2510 : ! ... Dummy args
2511 : !-----------------------------------------------------------------------
2512 : integer, intent(in) :: avec_len
2513 : real(r8), intent(in) :: lu(veclen,max(1,nzcnt))
2514 : real(r8), intent(inout) :: b(veclen,clscnt4)
2515 : !-----------------------------------------------------------------------
2516 : ! ... Local variables
2517 : !-----------------------------------------------------------------------
2518 : integer :: k
2519 : !-----------------------------------------------------------------------
2520 : ! ... solve L * y = b
2521 : !-----------------------------------------------------------------------
2522 325174377 : do k = 1,avec_len
2523 315224452 : b(k,123) = b(k,123) * lu(k,421)
2524 315224452 : b(k,122) = b(k,122) * lu(k,415)
2525 315224452 : b(k,121) = b(k,121) * lu(k,409)
2526 315224452 : b(k,120) = b(k,120) * lu(k,403)
2527 315224452 : b(k,119) = b(k,119) * lu(k,397)
2528 315224452 : b(k,118) = b(k,118) * lu(k,391)
2529 315224452 : b(k,117) = b(k,117) * lu(k,385)
2530 315224452 : b(k,116) = b(k,116) * lu(k,379)
2531 315224452 : b(k,115) = b(k,115) * lu(k,371)
2532 315224452 : b(k,114) = b(k,114) * lu(k,363)
2533 315224452 : b(k,113) = b(k,113) * lu(k,355)
2534 315224452 : b(k,112) = b(k,112) * lu(k,350)
2535 315224452 : b(k,111) = b(k,111) * lu(k,345)
2536 315224452 : b(k,88) = b(k,88) - lu(k,344) * b(k,111)
2537 315224452 : b(k,110) = b(k,110) * lu(k,339)
2538 315224452 : b(k,109) = b(k,109) * lu(k,334)
2539 315224452 : b(k,108) = b(k,108) * lu(k,329)
2540 315224452 : b(k,107) = b(k,107) * lu(k,326)
2541 315224452 : b(k,106) = b(k,106) * lu(k,321)
2542 315224452 : b(k,105) = b(k,105) * lu(k,316)
2543 315224452 : b(k,104) = b(k,104) * lu(k,310)
2544 315224452 : b(k,103) = b(k,103) * lu(k,304)
2545 315224452 : b(k,102) = b(k,102) * lu(k,298)
2546 315224452 : b(k,101) = b(k,101) * lu(k,292)
2547 315224452 : b(k,89) = b(k,89) - lu(k,291) * b(k,101)
2548 315224452 : b(k,100) = b(k,100) * lu(k,285)
2549 315224452 : b(k,99) = b(k,99) * lu(k,282)
2550 315224452 : b(k,98) = b(k,98) * lu(k,279)
2551 315224452 : b(k,97) = b(k,97) * lu(k,273)
2552 315224452 : b(k,96) = b(k,96) * lu(k,269)
2553 315224452 : b(k,95) = b(k,95) * lu(k,265)
2554 315224452 : b(k,94) = b(k,94) * lu(k,261)
2555 315224452 : b(k,64) = b(k,64) - lu(k,260) * b(k,94)
2556 315224452 : b(k,93) = b(k,93) * lu(k,256)
2557 315224452 : b(k,92) = b(k,92) * lu(k,251)
2558 315224452 : b(k,89) = b(k,89) - lu(k,250) * b(k,92)
2559 315224452 : b(k,91) = b(k,91) * lu(k,246)
2560 315224452 : b(k,90) = b(k,90) * lu(k,241)
2561 315224452 : b(k,89) = b(k,89) * lu(k,238)
2562 315224452 : b(k,88) = b(k,88) * lu(k,235)
2563 315224452 : b(k,87) = b(k,87) * lu(k,232)
2564 315224452 : b(k,86) = b(k,86) * lu(k,229)
2565 315224452 : b(k,85) = b(k,85) * lu(k,226)
2566 315224452 : b(k,84) = b(k,84) * lu(k,221)
2567 315224452 : b(k,83) = b(k,83) * lu(k,216)
2568 315224452 : b(k,82) = b(k,82) * lu(k,208)
2569 315224452 : b(k,80) = b(k,80) - lu(k,207) * b(k,82)
2570 315224452 : b(k,53) = b(k,53) - lu(k,206) * b(k,82)
2571 315224452 : b(k,81) = b(k,81) * lu(k,203)
2572 315224452 : b(k,80) = b(k,80) * lu(k,199)
2573 315224452 : b(k,79) = b(k,79) * lu(k,194)
2574 315224452 : b(k,78) = b(k,78) * lu(k,187)
2575 315224452 : b(k,52) = b(k,52) - lu(k,186) * b(k,78)
2576 315224452 : b(k,77) = b(k,77) * lu(k,182)
2577 315224452 : b(k,76) = b(k,76) * lu(k,180)
2578 315224452 : b(k,75) = b(k,75) - lu(k,179) * b(k,76)
2579 315224452 : b(k,75) = b(k,75) * lu(k,177)
2580 315224452 : b(k,74) = b(k,74) * lu(k,175)
2581 315224452 : b(k,73) = b(k,73) * lu(k,170)
2582 315224452 : b(k,72) = b(k,72) * lu(k,166)
2583 315224452 : b(k,71) = b(k,71) * lu(k,160)
2584 315224452 : b(k,47) = b(k,47) - lu(k,159) * b(k,71)
2585 315224452 : b(k,70) = b(k,70) * lu(k,154)
2586 315224452 : b(k,69) = b(k,69) * lu(k,149)
2587 315224452 : b(k,68) = b(k,68) * lu(k,144)
2588 315224452 : b(k,67) = b(k,67) * lu(k,139)
2589 315224452 : b(k,66) = b(k,66) * lu(k,136)
2590 315224452 : b(k,65) = b(k,65) * lu(k,133)
2591 315224452 : b(k,64) = b(k,64) * lu(k,130)
2592 315224452 : b(k,63) = b(k,63) * lu(k,127)
2593 315224452 : b(k,62) = b(k,62) * lu(k,123)
2594 315224452 : b(k,61) = b(k,61) * lu(k,119)
2595 315224452 : b(k,60) = b(k,60) * lu(k,115)
2596 315224452 : b(k,59) = b(k,59) * lu(k,111)
2597 315224452 : b(k,58) = b(k,58) * lu(k,107)
2598 315224452 : b(k,57) = b(k,57) * lu(k,104)
2599 315224452 : b(k,56) = b(k,56) * lu(k,101)
2600 315224452 : b(k,55) = b(k,55) * lu(k,98)
2601 315224452 : b(k,54) = b(k,54) * lu(k,95)
2602 315224452 : b(k,53) = b(k,53) * lu(k,94)
2603 315224452 : b(k,43) = b(k,43) - lu(k,93) * b(k,53)
2604 315224452 : b(k,42) = b(k,42) - lu(k,92) * b(k,53)
2605 315224452 : b(k,41) = b(k,41) - lu(k,91) * b(k,53)
2606 315224452 : b(k,40) = b(k,40) - lu(k,90) * b(k,53)
2607 315224452 : b(k,39) = b(k,39) - lu(k,89) * b(k,53)
2608 315224452 : b(k,52) = b(k,52) * lu(k,88)
2609 315224452 : b(k,43) = b(k,43) - lu(k,87) * b(k,52)
2610 315224452 : b(k,42) = b(k,42) - lu(k,86) * b(k,52)
2611 315224452 : b(k,41) = b(k,41) - lu(k,85) * b(k,52)
2612 315224452 : b(k,40) = b(k,40) - lu(k,84) * b(k,52)
2613 315224452 : b(k,39) = b(k,39) - lu(k,83) * b(k,52)
2614 315224452 : b(k,51) = b(k,51) * lu(k,82)
2615 315224452 : b(k,43) = b(k,43) - lu(k,81) * b(k,51)
2616 315224452 : b(k,42) = b(k,42) - lu(k,80) * b(k,51)
2617 315224452 : b(k,41) = b(k,41) - lu(k,79) * b(k,51)
2618 315224452 : b(k,40) = b(k,40) - lu(k,78) * b(k,51)
2619 315224452 : b(k,39) = b(k,39) - lu(k,77) * b(k,51)
2620 315224452 : b(k,50) = b(k,50) * lu(k,76)
2621 315224452 : b(k,49) = b(k,49) - lu(k,75) * b(k,50)
2622 315224452 : b(k,49) = b(k,49) * lu(k,74)
2623 315224452 : b(k,43) = b(k,43) - lu(k,73) * b(k,49)
2624 315224452 : b(k,42) = b(k,42) - lu(k,72) * b(k,49)
2625 315224452 : b(k,41) = b(k,41) - lu(k,71) * b(k,49)
2626 315224452 : b(k,40) = b(k,40) - lu(k,70) * b(k,49)
2627 315224452 : b(k,39) = b(k,39) - lu(k,69) * b(k,49)
2628 315224452 : b(k,48) = b(k,48) * lu(k,68)
2629 315224452 : b(k,43) = b(k,43) - lu(k,67) * b(k,48)
2630 315224452 : b(k,42) = b(k,42) - lu(k,66) * b(k,48)
2631 315224452 : b(k,41) = b(k,41) - lu(k,65) * b(k,48)
2632 315224452 : b(k,40) = b(k,40) - lu(k,64) * b(k,48)
2633 315224452 : b(k,39) = b(k,39) - lu(k,63) * b(k,48)
2634 315224452 : b(k,47) = b(k,47) * lu(k,62)
2635 315224452 : b(k,43) = b(k,43) - lu(k,61) * b(k,47)
2636 315224452 : b(k,42) = b(k,42) - lu(k,60) * b(k,47)
2637 315224452 : b(k,41) = b(k,41) - lu(k,59) * b(k,47)
2638 315224452 : b(k,40) = b(k,40) - lu(k,58) * b(k,47)
2639 315224452 : b(k,39) = b(k,39) - lu(k,57) * b(k,47)
2640 315224452 : b(k,46) = b(k,46) * lu(k,56)
2641 315224452 : b(k,43) = b(k,43) - lu(k,55) * b(k,46)
2642 315224452 : b(k,42) = b(k,42) - lu(k,54) * b(k,46)
2643 315224452 : b(k,41) = b(k,41) - lu(k,53) * b(k,46)
2644 315224452 : b(k,40) = b(k,40) - lu(k,52) * b(k,46)
2645 315224452 : b(k,39) = b(k,39) - lu(k,51) * b(k,46)
2646 315224452 : b(k,45) = b(k,45) * lu(k,50)
2647 315224452 : b(k,43) = b(k,43) - lu(k,49) * b(k,45)
2648 315224452 : b(k,42) = b(k,42) - lu(k,48) * b(k,45)
2649 315224452 : b(k,41) = b(k,41) - lu(k,47) * b(k,45)
2650 315224452 : b(k,40) = b(k,40) - lu(k,46) * b(k,45)
2651 315224452 : b(k,39) = b(k,39) - lu(k,45) * b(k,45)
2652 315224452 : b(k,44) = b(k,44) * lu(k,44)
2653 315224452 : b(k,43) = b(k,43) * lu(k,43)
2654 315224452 : b(k,42) = b(k,42) * lu(k,42)
2655 315224452 : b(k,41) = b(k,41) * lu(k,41)
2656 315224452 : b(k,40) = b(k,40) * lu(k,40)
2657 315224452 : b(k,39) = b(k,39) * lu(k,39)
2658 315224452 : b(k,38) = b(k,38) * lu(k,38)
2659 315224452 : b(k,37) = b(k,37) * lu(k,37)
2660 315224452 : b(k,36) = b(k,36) * lu(k,36)
2661 315224452 : b(k,35) = b(k,35) * lu(k,35)
2662 315224452 : b(k,34) = b(k,34) * lu(k,34)
2663 315224452 : b(k,33) = b(k,33) * lu(k,33)
2664 315224452 : b(k,32) = b(k,32) * lu(k,32)
2665 315224452 : b(k,31) = b(k,31) * lu(k,31)
2666 315224452 : b(k,30) = b(k,30) * lu(k,30)
2667 315224452 : b(k,29) = b(k,29) * lu(k,29)
2668 315224452 : b(k,28) = b(k,28) * lu(k,28)
2669 315224452 : b(k,27) = b(k,27) * lu(k,27)
2670 315224452 : b(k,26) = b(k,26) * lu(k,26)
2671 315224452 : b(k,25) = b(k,25) * lu(k,25)
2672 315224452 : b(k,24) = b(k,24) * lu(k,24)
2673 315224452 : b(k,23) = b(k,23) * lu(k,23)
2674 315224452 : b(k,22) = b(k,22) * lu(k,22)
2675 315224452 : b(k,21) = b(k,21) * lu(k,21)
2676 315224452 : b(k,20) = b(k,20) * lu(k,20)
2677 315224452 : b(k,19) = b(k,19) * lu(k,19)
2678 315224452 : b(k,18) = b(k,18) * lu(k,18)
2679 315224452 : b(k,17) = b(k,17) * lu(k,17)
2680 315224452 : b(k,16) = b(k,16) * lu(k,16)
2681 315224452 : b(k,15) = b(k,15) * lu(k,15)
2682 315224452 : b(k,14) = b(k,14) * lu(k,14)
2683 315224452 : b(k,13) = b(k,13) * lu(k,13)
2684 315224452 : b(k,12) = b(k,12) * lu(k,12)
2685 315224452 : b(k,11) = b(k,11) * lu(k,11)
2686 315224452 : b(k,10) = b(k,10) * lu(k,10)
2687 315224452 : b(k,9) = b(k,9) * lu(k,9)
2688 315224452 : b(k,8) = b(k,8) * lu(k,8)
2689 315224452 : b(k,7) = b(k,7) * lu(k,7)
2690 315224452 : b(k,6) = b(k,6) * lu(k,6)
2691 315224452 : b(k,5) = b(k,5) * lu(k,5)
2692 315224452 : b(k,4) = b(k,4) * lu(k,4)
2693 315224452 : b(k,3) = b(k,3) * lu(k,3)
2694 315224452 : b(k,2) = b(k,2) * lu(k,2)
2695 325174377 : b(k,1) = b(k,1) * lu(k,1)
2696 : end do
2697 9949925 : end subroutine lu_slv12
2698 9949925 : subroutine lu_slv( avec_len, lu, b )
2699 : use shr_kind_mod, only : r8 => shr_kind_r8
2700 : use chem_mods, only : clscnt4, nzcnt
2701 : implicit none
2702 : !-----------------------------------------------------------------------
2703 : ! ... Dummy args
2704 : !-----------------------------------------------------------------------
2705 : integer, intent(in) :: avec_len
2706 : real(r8), intent(in) :: lu(veclen,max(1,nzcnt))
2707 : real(r8), intent(inout) :: b(veclen,clscnt4)
2708 9949925 : call lu_slv01( avec_len, lu, b )
2709 9949925 : call lu_slv02( avec_len, lu, b )
2710 9949925 : call lu_slv03( avec_len, lu, b )
2711 9949925 : call lu_slv04( avec_len, lu, b )
2712 9949925 : call lu_slv05( avec_len, lu, b )
2713 9949925 : call lu_slv06( avec_len, lu, b )
2714 9949925 : call lu_slv07( avec_len, lu, b )
2715 9949925 : call lu_slv08( avec_len, lu, b )
2716 9949925 : call lu_slv09( avec_len, lu, b )
2717 9949925 : call lu_slv10( avec_len, lu, b )
2718 9949925 : call lu_slv11( avec_len, lu, b )
2719 9949925 : call lu_slv12( avec_len, lu, b )
2720 9949925 : end subroutine lu_slv
2721 : end module mo_lu_solve
|