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