27         SUBROUTINE lu_slv01_vec(ncol,nb,nz,lu, b)  argument
35             REAL(KIND=r8), intent(in)    :: lu(ncol,nz)  local
47          b(:,125) = b(:,125) - lu(:,18) * b(:,17)
48          b(:,131) = b(:,131) - lu(:,19) * b(:,17)
50          b(:,124) = b(:,124) - lu(:,21) * b(:,18)
51          b(:,126) = b(:,126) - lu(:,22) * b(:,18)
55          b(:,125) = b(:,125) - lu(:,18) * b(:,17)
56          b(:,131) = b(:,131) - lu(:,19) * b(:,17)
58          b(i,124) = b(i,124) - lu(i,21) * b(i,18)
59          b(i,126) = b(i,126) - lu(i,22) * b(i,18)
61          b(i,79) = b(i,79) - lu(i,24) * b(i,19)
62          b(i,131) = b(i,131) - lu(i,25) * b(i,19)
63          b(i,41) = b(i,41) - lu(i,27) * b(i,20)
64          b(i,131) = b(i,131) - lu(i,28) * b(i,20)
65          b(i,96) = b(i,96) - lu(i,30) * b(i,21)
66          b(i,131) = b(i,131) - lu(i,31) * b(i,21)
67          b(i,134) = b(i,134) - lu(i,32) * b(i,21)
68          b(i,23) = b(i,23) - lu(i,34) * b(i,22)
69          b(i,65) = b(i,65) - lu(i,35) * b(i,22)
70          b(i,125) = b(i,125) - lu(i,36) * b(i,22)
71          b(i,131) = b(i,131) - lu(i,37) * b(i,22)
72          b(i,31) = b(i,31) - lu(i,39) * b(i,23)
73          b(i,131) = b(i,131) - lu(i,40) * b(i,23)
74          b(i,56) = b(i,56) - lu(i,42) * b(i,24)
75          b(i,131) = b(i,131) - lu(i,43) * b(i,24)
76          b(i,88) = b(i,88) - lu(i,45) * b(i,25)
77          b(i,122) = b(i,122) - lu(i,46) * b(i,25)
78          b(i,36) = b(i,36) - lu(i,48) * b(i,26)
79          b(i,134) = b(i,134) - lu(i,49) * b(i,26)
80          b(i,120) = b(i,120) - lu(i,51) * b(i,27)
81          b(i,120) = b(i,120) - lu(i,54) * b(i,28)
82          b(i,126) = b(i,126) - lu(i,56) * b(i,29)
83          b(i,122) = b(i,122) - lu(i,58) * b(i,30)
84          b(i,125) = b(i,125) - lu(i,59) * b(i,30)
85          b(i,131) = b(i,131) - lu(i,60) * b(i,30)
86          b(i,66) = b(i,66) - lu(i,62) * b(i,31)
87          b(i,125) = b(i,125) - lu(i,63) * b(i,31)
88          b(i,130) = b(i,130) - lu(i,64) * b(i,31)
89          b(i,88) = b(i,88) - lu(i,66) * b(i,32)
90          b(i,122) = b(i,122) - lu(i,67) * b(i,32)
91          b(i,126) = b(i,126) - lu(i,68) * b(i,32)
92          b(i,118) = b(i,118) - lu(i,70) * b(i,33)
93          b(i,126) = b(i,126) - lu(i,71) * b(i,33)
94          b(i,88) = b(i,88) - lu(i,73) * b(i,34)
95          b(i,127) = b(i,127) - lu(i,74) * b(i,34)
96          b(i,104) = b(i,104) - lu(i,76) * b(i,35)
97          b(i,125) = b(i,125) - lu(i,77) * b(i,35)
98          b(i,131) = b(i,131) - lu(i,78) * b(i,35)
99          b(i,99) = b(i,99) - lu(i,81) * b(i,36)
100          b(i,121) = b(i,121) - lu(i,82) * b(i,36)
101          b(i,134) = b(i,134) - lu(i,83) * b(i,36)
102          b(i,91) = b(i,91) - lu(i,85) * b(i,37)
103          b(i,117) = b(i,117) - lu(i,86) * b(i,37)
104          b(i,126) = b(i,126) - lu(i,87) * b(i,37)
105          b(i,131) = b(i,131) - lu(i,88) * b(i,37)
106          b(i,134) = b(i,134) - lu(i,89) * b(i,37)
107          b(i,64) = b(i,64) - lu(i,91) * b(i,38)
108          b(i,81) = b(i,81) - lu(i,92) * b(i,38)
109          b(i,103) = b(i,103) - lu(i,93) * b(i,38)
110          b(i,125) = b(i,125) - lu(i,94) * b(i,38)
111          b(i,131) = b(i,131) - lu(i,95) * b(i,38)
112          b(i,99) = b(i,99) - lu(i,97) * b(i,39)
113          b(i,125) = b(i,125) - lu(i,98) * b(i,39)
114          b(i,131) = b(i,131) - lu(i,99) * b(i,39)
115          b(i,132) = b(i,132) - lu(i,100) * b(i,39)
116          b(i,133) = b(i,133) - lu(i,101) * b(i,39)
117          b(i,121) = b(i,121) - lu(i,103) * b(i,40)
118          b(i,129) = b(i,129) - lu(i,104) * b(i,40)
119          b(i,130) = b(i,130) - lu(i,105) * b(i,40)
120          b(i,132) = b(i,132) - lu(i,106) * b(i,40)
121          b(i,133) = b(i,133) - lu(i,107) * b(i,40)
122          b(i,80) = b(i,80) - lu(i,109) * b(i,41)
123          b(i,104) = b(i,104) - lu(i,110) * b(i,41)
124          b(i,125) = b(i,125) - lu(i,111) * b(i,41)
125          b(i,129) = b(i,129) - lu(i,112) * b(i,41)
126          b(i,130) = b(i,130) - lu(i,113) * b(i,41)
127          b(i,135) = b(i,135) - lu(i,114) * b(i,41)
128          b(i,77) = b(i,77) - lu(i,116) * b(i,42)
129          b(i,104) = b(i,104) - lu(i,117) * b(i,42)
130          b(i,115) = b(i,115) - lu(i,118) * b(i,42)
131          b(i,131) = b(i,131) - lu(i,119) * b(i,42)
132          b(i,112) = b(i,112) - lu(i,121) * b(i,43)
133          b(i,114) = b(i,114) - lu(i,122) * b(i,43)
134          b(i,125) = b(i,125) - lu(i,123) * b(i,43)
135          b(i,131) = b(i,131) - lu(i,124) * b(i,43)
136          b(i,91) = b(i,91) - lu(i,126) * b(i,44)
137          b(i,104) = b(i,104) - lu(i,127) * b(i,44)
138          b(i,125) = b(i,125) - lu(i,128) * b(i,44)
139          b(i,131) = b(i,131) - lu(i,129) * b(i,44)
140          b(i,110) = b(i,110) - lu(i,131) * b(i,45)
141          b(i,131) = b(i,131) - lu(i,132) * b(i,45)
142          b(i,134) = b(i,134) - lu(i,133) * b(i,45)
143          b(i,99) = b(i,99) - lu(i,135) * b(i,46)
144          b(i,116) = b(i,116) - lu(i,136) * b(i,46)
145          b(i,121) = b(i,121) - lu(i,137) * b(i,46)
146          b(i,124) = b(i,124) - lu(i,138) * b(i,46)
147          b(i,110) = b(i,110) - lu(i,140) * b(i,47)
148          b(i,131) = b(i,131) - lu(i,141) * b(i,47)
149          b(i,82) = b(i,82) - lu(i,143) * b(i,48)
150          b(i,99) = b(i,99) - lu(i,144) * b(i,48)
151          b(i,103) = b(i,103) - lu(i,145) * b(i,48)
152          b(i,116) = b(i,116) - lu(i,146) * b(i,48)
153          b(i,121) = b(i,121) - lu(i,147) * b(i,48)
154          b(i,127) = b(i,127) - lu(i,148) * b(i,48)
155          b(i,131) = b(i,131) - lu(i,149) * b(i,48)
156          b(i,109) = b(i,109) - lu(i,151) * b(i,49)
157          b(i,130) = b(i,130) - lu(i,152) * b(i,49)
158          b(i,131) = b(i,131) - lu(i,153) * b(i,49)
159          b(i,119) = b(i,119) - lu(i,155) * b(i,50)
160          b(i,127) = b(i,127) - lu(i,156) * b(i,50)
161          b(i,131) = b(i,131) - lu(i,157) * b(i,50)
162          b(i,134) = b(i,134) - lu(i,158) * b(i,50)
163          b(i,135) = b(i,135) - lu(i,159) * b(i,50)
164          b(i,65) = b(i,65) - lu(i,161) * b(i,51)
165          b(i,66) = b(i,66) - lu(i,162) * b(i,51)
166          b(i,81) = b(i,81) - lu(i,163) * b(i,51)
167          b(i,109) = b(i,109) - lu(i,164) * b(i,51)
168          b(i,131) = b(i,131) - lu(i,165) * b(i,51)
169          b(i,80) = b(i,80) - lu(i,167) * b(i,52)
170          b(i,96) = b(i,96) - lu(i,168) * b(i,52)
171          b(i,125) = b(i,125) - lu(i,169) * b(i,52)
172          b(i,131) = b(i,131) - lu(i,170) * b(i,52)
173          b(i,134) = b(i,134) - lu(i,171) * b(i,52)
174          b(i,106) = b(i,106) - lu(i,173) * b(i,53)
175          b(i,115) = b(i,115) - lu(i,174) * b(i,53)
176          b(i,131) = b(i,131) - lu(i,175) * b(i,53)
177          b(i,134) = b(i,134) - lu(i,176) * b(i,53)
178          b(i,135) = b(i,135) - lu(i,177) * b(i,53)
179          b(i,64) = b(i,64) - lu(i,179) * b(i,54)
180          b(i,125) = b(i,125) - lu(i,180) * b(i,54)
181          b(i,129) = b(i,129) - lu(i,181) * b(i,54)
182          b(i,130) = b(i,130) - lu(i,182) * b(i,54)
183          b(i,135) = b(i,135) - lu(i,183) * b(i,54)
184          b(i,77) = b(i,77) - lu(i,185) * b(i,55)
185          b(i,91) = b(i,91) - lu(i,186) * b(i,55)
186          b(i,115) = b(i,115) - lu(i,187) * b(i,55)
187          b(i,131) = b(i,131) - lu(i,188) * b(i,55)
188          b(i,95) = b(i,95) - lu(i,190) * b(i,56)
189          b(i,120) = b(i,120) - lu(i,191) * b(i,56)
190          b(i,125) = b(i,125) - lu(i,192) * b(i,56)
191          b(i,135) = b(i,135) - lu(i,193) * b(i,56)
192          b(i,115) = b(i,115) - lu(i,195) * b(i,57)
193          b(i,119) = b(i,119) - lu(i,196) * b(i,57)
194          b(i,130) = b(i,130) - lu(i,197) * b(i,57)
195          b(i,131) = b(i,131) - lu(i,198) * b(i,57)
196          b(i,132) = b(i,132) - lu(i,199) * b(i,57)
197          b(i,135) = b(i,135) - lu(i,200) * b(i,57)
198          b(i,72) = b(i,72) - lu(i,202) * b(i,58)
199          b(i,85) = b(i,85) - lu(i,203) * b(i,58)
200          b(i,86) = b(i,86) - lu(i,204) * b(i,58)
201          b(i,92) = b(i,92) - lu(i,205) * b(i,58)
202          b(i,120) = b(i,120) - lu(i,206) * b(i,58)
203          b(i,121) = b(i,121) - lu(i,207) * b(i,58)
204          b(i,80) = b(i,80) - lu(i,209) * b(i,59)
205          b(i,98) = b(i,98) - lu(i,210) * b(i,59)
206          b(i,107) = b(i,107) - lu(i,211) * b(i,59)
207          b(i,113) = b(i,113) - lu(i,212) * b(i,59)
208          b(i,125) = b(i,125) - lu(i,213) * b(i,59)
209          b(i,131) = b(i,131) - lu(i,214) * b(i,59)
210          b(i,120) = b(i,120) - lu(i,216) * b(i,60)
211          b(i,125) = b(i,125) - lu(i,217) * b(i,60)
212          b(i,130) = b(i,130) - lu(i,218) * b(i,60)
213          b(i,131) = b(i,131) - lu(i,219) * b(i,60)
214          b(i,132) = b(i,132) - lu(i,220) * b(i,60)
215          b(i,134) = b(i,134) - lu(i,221) * b(i,60)
216          b(i,92) = b(i,92) - lu(i,223) * b(i,61)
217          b(i,120) = b(i,120) - lu(i,224) * b(i,61)
218          b(i,122) = b(i,122) - lu(i,225) * b(i,61)
219          b(i,129) = b(i,129) - lu(i,226) * b(i,61)
220          b(i,115) = b(i,115) - lu(i,228) * b(i,62)
221          b(i,119) = b(i,119) - lu(i,229) * b(i,62)
222          b(i,131) = b(i,131) - lu(i,230) * b(i,62)
223          b(i,134) = b(i,134) - lu(i,231) * b(i,62)
224          b(i,135) = b(i,135) - lu(i,232) * b(i,62)
225          b(i,64) = b(i,64) - lu(i,234) * b(i,63)
226          b(i,83) = b(i,83) - lu(i,235) * b(i,63)
227          b(i,103) = b(i,103) - lu(i,236) * b(i,63)
228          b(i,123) = b(i,123) - lu(i,237) * b(i,63)
229          b(i,125) = b(i,125) - lu(i,238) * b(i,63)
230          b(i,131) = b(i,131) - lu(i,239) * b(i,63)
231          b(i,135) = b(i,135) - lu(i,240) * b(i,63)
232          b(i,125) = b(i,125) - lu(i,242) * b(i,64)
233          b(i,131) = b(i,131) - lu(i,243) * b(i,64)
234          b(i,134) = b(i,134) - lu(i,244) * b(i,64)
235          b(i,66) = b(i,66) - lu(i,247) * b(i,65)
236          b(i,81) = b(i,81) - lu(i,248) * b(i,65)
237          b(i,109) = b(i,109) - lu(i,249) * b(i,65)
238          b(i,125) = b(i,125) - lu(i,250) * b(i,65)
239          b(i,129) = b(i,129) - lu(i,251) * b(i,65)
240          b(i,130) = b(i,130) - lu(i,252) * b(i,65)
241          b(i,131) = b(i,131) - lu(i,253) * b(i,65)
242          b(i,81) = b(i,81) - lu(i,255) * b(i,66)
243          b(i,103) = b(i,103) - lu(i,256) * b(i,66)
244          b(i,109) = b(i,109) - lu(i,257) * b(i,66)
245          b(i,115) = b(i,115) - lu(i,258) * b(i,66)
246          b(i,125) = b(i,125) - lu(i,259) * b(i,66)
247          b(i,89) = b(i,89) - lu(i,261) * b(i,67)
248          b(i,104) = b(i,104) - lu(i,262) * b(i,67)
249          b(i,105) = b(i,105) - lu(i,263) * b(i,67)
250          b(i,125) = b(i,125) - lu(i,264) * b(i,67)
251          b(i,131) = b(i,131) - lu(i,265) * b(i,67)
252          b(i,134) = b(i,134) - lu(i,266) * b(i,67)
253          b(i,135) = b(i,135) - lu(i,267) * b(i,67)
254          b(i,125) = b(i,125) - lu(i,269) * b(i,68)
255          b(i,131) = b(i,131) - lu(i,270) * b(i,68)
256          b(i,135) = b(i,135) - lu(i,271) * b(i,68)
257          b(i,107) = b(i,107) - lu(i,273) * b(i,69)
258          b(i,110) = b(i,110) - lu(i,274) * b(i,69)
259          b(i,111) = b(i,111) - lu(i,275) * b(i,69)
260          b(i,113) = b(i,113) - lu(i,276) * b(i,69)
261          b(i,125) = b(i,125) - lu(i,277) * b(i,69)
262          b(i,131) = b(i,131) - lu(i,278) * b(i,69)
263          b(i,135) = b(i,135) - lu(i,279) * b(i,69)
272         SUBROUTINE lu_slv02_vec(ncol,nb,nz,lu, b)  argument
279             REAL(KIND=r8), intent(in) :: lu(ncol,nz)  local
291          b(i,84) =   b(i,84) - lu(i,281) * b(i,70)
292          b(i,118) = b(i,118) - lu(i,282) * b(i,70)
293          b(i,121) = b(i,121) - lu(i,283) * b(i,70)
294          b(i,128) = b(i,128) - lu(i,284) * b(i,70)
295          b(i,130) = b(i,130) - lu(i,285) * b(i,70)
296          b(i,132) = b(i,132) - lu(i,286) * b(i,70)
297          b(i,133) = b(i,133) - lu(i,287) * b(i,70)
301          b(i,105) = b(i,105) - lu(i,289) * b(i,71)
302          b(i,114) = b(i,114) - lu(i,290) * b(i,71)
303          b(i,125) = b(i,125) - lu(i,291) * b(i,71)
304          b(i,130) = b(i,130) - lu(i,292) * b(i,71)
305          b(i,131) = b(i,131) - lu(i,293) * b(i,71)
306          b(i,132) = b(i,132) - lu(i,294) * b(i,71)
307          b(i,135) = b(i,135) - lu(i,295) * b(i,71)
308          b(i,85) = b(i,85) - lu(i,297) * b(i,72)
309          b(i,86) = b(i,86) - lu(i,298) * b(i,72)
310          b(i,92) = b(i,92) - lu(i,299) * b(i,72)
311          b(i,103) = b(i,103) - lu(i,300) * b(i,72)
312          b(i,120) = b(i,120) - lu(i,301) * b(i,72)
313          b(i,121) = b(i,121) - lu(i,302) * b(i,72)
314          b(i,98) = b(i,98) - lu(i,304) * b(i,73)
315          b(i,107) = b(i,107) - lu(i,305) * b(i,73)
316          b(i,113) = b(i,113) - lu(i,306) * b(i,73)
317          b(i,123) = b(i,123) - lu(i,307) * b(i,73)
318          b(i,125) = b(i,125) - lu(i,308) * b(i,73)
319          b(i,130) = b(i,130) - lu(i,309) * b(i,73)
320          b(i,131) = b(i,131) - lu(i,310) * b(i,73)
321          b(i,132) = b(i,132) - lu(i,311) * b(i,73)
322          b(i,117) = b(i,117) - lu(i,313) * b(i,74)
323          b(i,121) = b(i,121) - lu(i,314) * b(i,74)
324          b(i,125) = b(i,125) - lu(i,315) * b(i,74)
325          b(i,126) = b(i,126) - lu(i,316) * b(i,74)
326          b(i,131) = b(i,131) - lu(i,317) * b(i,74)
327          b(i,134) = b(i,134) - lu(i,318) * b(i,74)
328          b(i,119) = b(i,119) - lu(i,320) * b(i,75)
329          b(i,131) = b(i,131) - lu(i,321) * b(i,75)
330          b(i,134) = b(i,134) - lu(i,322) * b(i,75)
331          b(i,77) = b(i,77) - lu(i,325) * b(i,76)
332          b(i,79) = b(i,79) - lu(i,326) * b(i,76)
333          b(i,80) = b(i,80) - lu(i,327) * b(i,76)
334          b(i,91) = b(i,91) - lu(i,328) * b(i,76)
335          b(i,104) = b(i,104) - lu(i,329) * b(i,76)
336          b(i,115) = b(i,115) - lu(i,330) * b(i,76)
337          b(i,125) = b(i,125) - lu(i,331) * b(i,76)
338          b(i,131) = b(i,131) - lu(i,332) * b(i,76)
339          b(i,135) = b(i,135) - lu(i,333) * b(i,76)
340          b(i,104) = b(i,104) - lu(i,336) * b(i,77)
341          b(i,115) = b(i,115) - lu(i,337) * b(i,77)
342          b(i,125) = b(i,125) - lu(i,338) * b(i,77)
343          b(i,129) = b(i,129) - lu(i,339) * b(i,77)
344          b(i,130) = b(i,130) - lu(i,340) * b(i,77)
345          b(i,131) = b(i,131) - lu(i,341) * b(i,77)
346          b(i,85) = b(i,85) - lu(i,345) * b(i,78)
347          b(i,86) = b(i,86) - lu(i,346) * b(i,78)
348          b(i,87) = b(i,87) - lu(i,347) * b(i,78)
349          b(i,92) = b(i,92) - lu(i,348) * b(i,78)
350          b(i,103) = b(i,103) - lu(i,349) * b(i,78)
351          b(i,120) = b(i,120) - lu(i,350) * b(i,78)
352          b(i,121) = b(i,121) - lu(i,351) * b(i,78)
353          b(i,122) = b(i,122) - lu(i,352) * b(i,78)
354          b(i,129) = b(i,129) - lu(i,353) * b(i,78)
355          b(i,80) = b(i,80) - lu(i,359) * b(i,79)
356          b(i,91) = b(i,91) - lu(i,360) * b(i,79)
357          b(i,104) = b(i,104) - lu(i,361) * b(i,79)
358          b(i,109) = b(i,109) - lu(i,362) * b(i,79)
359          b(i,115) = b(i,115) - lu(i,363) * b(i,79)
360          b(i,125) = b(i,125) - lu(i,364) * b(i,79)
361          b(i,129) = b(i,129) - lu(i,365) * b(i,79)
362          b(i,130) = b(i,130) - lu(i,366) * b(i,79)
363          b(i,131) = b(i,131) - lu(i,367) * b(i,79)
364          b(i,135) = b(i,135) - lu(i,368) * b(i,79)
365          b(i,106) = b(i,106) - lu(i,370) * b(i,80)
366          b(i,115) = b(i,115) - lu(i,371) * b(i,80)
367          b(i,119) = b(i,119) - lu(i,372) * b(i,80)
368          b(i,131) = b(i,131) - lu(i,373) * b(i,80)
369          b(i,134) = b(i,134) - lu(i,374) * b(i,80)
370          b(i,103) = b(i,103) - lu(i,376) * b(i,81)
371          b(i,125) = b(i,125) - lu(i,377) * b(i,81)
372          b(i,131) = b(i,131) - lu(i,378) * b(i,81)
373          b(i,116) = b(i,116) - lu(i,380) * b(i,82)
374          b(i,120) = b(i,120) - lu(i,381) * b(i,82)
375          b(i,121) = b(i,121) - lu(i,382) * b(i,82)
376          b(i,123) = b(i,123) - lu(i,383) * b(i,82)
377          b(i,127) = b(i,127) - lu(i,384) * b(i,82)
378          b(i,131) = b(i,131) - lu(i,385) * b(i,82)
379          b(i,95) = b(i,95) - lu(i,389) * b(i,83)
380          b(i,120) = b(i,120) - lu(i,390) * b(i,83)
381          b(i,125) = b(i,125) - lu(i,391) * b(i,83)
382          b(i,129) = b(i,129) - lu(i,392) * b(i,83)
383          b(i,130) = b(i,130) - lu(i,393) * b(i,83)
384          b(i,131) = b(i,131) - lu(i,394) * b(i,83)
385          b(i,135) = b(i,135) - lu(i,395) * b(i,83)
386          b(i,117) = b(i,117) - lu(i,398) * b(i,84)
387          b(i,118) = b(i,118) - lu(i,399) * b(i,84)
388          b(i,121) = b(i,121) - lu(i,400) * b(i,84)
389          b(i,126) = b(i,126) - lu(i,401) * b(i,84)
390          b(i,128) = b(i,128) - lu(i,402) * b(i,84)
391          b(i,131) = b(i,131) - lu(i,403) * b(i,84)
392          b(i,134) = b(i,134) - lu(i,404) * b(i,84)
393          b(i,86) = b(i,86) - lu(i,406) * b(i,85)
394          b(i,87) = b(i,87) - lu(i,407) * b(i,85)
395          b(i,92) = b(i,92) - lu(i,408) * b(i,85)
396          b(i,120) = b(i,120) - lu(i,409) * b(i,85)
397          b(i,121) = b(i,121) - lu(i,410) * b(i,85)
398          b(i,122) = b(i,122) - lu(i,411) * b(i,85)
399          b(i,129) = b(i,129) - lu(i,412) * b(i,85)
400          b(i,87) = b(i,87) - lu(i,415) * b(i,86)
401          b(i,92) = b(i,92) - lu(i,416) * b(i,86)
402          b(i,120) = b(i,120) - lu(i,417) * b(i,86)
403          b(i,121) = b(i,121) - lu(i,418) * b(i,86)
404          b(i,122) = b(i,122) - lu(i,419) * b(i,86)
405          b(i,129) = b(i,129) - lu(i,420) * b(i,86)
406          b(i,92) = b(i,92) - lu(i,426) * b(i,87)
407          b(i,103) = b(i,103) - lu(i,427) * b(i,87)
408          b(i,120) = b(i,120) - lu(i,428) * b(i,87)
409          b(i,121) = b(i,121) - lu(i,429) * b(i,87)
410          b(i,122) = b(i,122) - lu(i,430) * b(i,87)
411          b(i,129) = b(i,129) - lu(i,431) * b(i,87)
412          b(i,108) = b(i,108) - lu(i,434) * b(i,88)
413          b(i,119) = b(i,119) - lu(i,435) * b(i,88)
414          b(i,127) = b(i,127) - lu(i,436) * b(i,88)
415          b(i,131) = b(i,131) - lu(i,437) * b(i,88)
416          b(i,132) = b(i,132) - lu(i,438) * b(i,88)
417          b(i,133) = b(i,133) - lu(i,439) * b(i,88)
418          b(i,134) = b(i,134) - lu(i,440) * b(i,88)
419          b(i,104) = b(i,104) - lu(i,443) * b(i,89)
420          b(i,105) = b(i,105) - lu(i,444) * b(i,89)
421          b(i,120) = b(i,120) - lu(i,445) * b(i,89)
422          b(i,125) = b(i,125) - lu(i,446) * b(i,89)
423          b(i,129) = b(i,129) - lu(i,447) * b(i,89)
424          b(i,130) = b(i,130) - lu(i,448) * b(i,89)
425          b(i,131) = b(i,131) - lu(i,449) * b(i,89)
426          b(i,134) = b(i,134) - lu(i,450) * b(i,89)
427          b(i,135) = b(i,135) - lu(i,451) * b(i,89)
428          b(i,118) = b(i,118) - lu(i,453) * b(i,90)
429          b(i,121) = b(i,121) - lu(i,454) * b(i,90)
430          b(i,122) = b(i,122) - lu(i,455) * b(i,90)
431          b(i,127) = b(i,127) - lu(i,456) * b(i,90)
432          b(i,131) = b(i,131) - lu(i,457) * b(i,90)
433          b(i,134) = b(i,134) - lu(i,458) * b(i,90)
434          b(i,104) = b(i,104) - lu(i,463) * b(i,91)
435          b(i,119) = b(i,119) - lu(i,464) * b(i,91)
436          b(i,120) = b(i,120) - lu(i,465) * b(i,91)
437          b(i,125) = b(i,125) - lu(i,466) * b(i,91)
438          b(i,129) = b(i,129) - lu(i,467) * b(i,91)
439          b(i,130) = b(i,130) - lu(i,468) * b(i,91)
440          b(i,131) = b(i,131) - lu(i,469) * b(i,91)
441          b(i,135) = b(i,135) - lu(i,470) * b(i,91)
442          b(i,103) = b(i,103) - lu(i,477) * b(i,92)
443          b(i,120) = b(i,120) - lu(i,478) * b(i,92)
444          b(i,121) = b(i,121) - lu(i,479) * b(i,92)
445          b(i,122) = b(i,122) - lu(i,480) * b(i,92)
446          b(i,127) = b(i,127) - lu(i,481) * b(i,92)
447          b(i,129) = b(i,129) - lu(i,482) * b(i,92)
448          b(i,130) = b(i,130) - lu(i,483) * b(i,92)
449          b(i,131) = b(i,131) - lu(i,484) * b(i,92)
450          b(i,117) = b(i,117) - lu(i,487) * b(i,93)
451          b(i,121) = b(i,121) - lu(i,488) * b(i,93)
452          b(i,124) = b(i,124) - lu(i,489) * b(i,93)
453          b(i,126) = b(i,126) - lu(i,490) * b(i,93)
454          b(i,131) = b(i,131) - lu(i,491) * b(i,93)
455          b(i,134) = b(i,134) - lu(i,492) * b(i,93)
456          b(i,101) = b(i,101) - lu(i,495) * b(i,94)
457          b(i,102) = b(i,102) - lu(i,496) * b(i,94)
458          b(i,103) = b(i,103) - lu(i,497) * b(i,94)
459          b(i,107) = b(i,107) - lu(i,498) * b(i,94)
460          b(i,111) = b(i,111) - lu(i,499) * b(i,94)
461          b(i,113) = b(i,113) - lu(i,500) * b(i,94)
462          b(i,114) = b(i,114) - lu(i,501) * b(i,94)
463          b(i,119) = b(i,119) - lu(i,502) * b(i,94)
464          b(i,123) = b(i,123) - lu(i,503) * b(i,94)
465          b(i,125) = b(i,125) - lu(i,504) * b(i,94)
466          b(i,131) = b(i,131) - lu(i,505) * b(i,94)
467          b(i,132) = b(i,132) - lu(i,506) * b(i,94)
468          b(i,134) = b(i,134) - lu(i,507) * b(i,94)
469          b(i,135) = b(i,135) - lu(i,508) * b(i,94)
470          b(i,103) = b(i,103) - lu(i,511) * b(i,95)
471          b(i,125) = b(i,125) - lu(i,512) * b(i,95)
472          b(i,131) = b(i,131) - lu(i,513) * b(i,95)
473          b(i,135) = b(i,135) - lu(i,514) * b(i,95)
474          b(i,104) = b(i,104) - lu(i,518) * b(i,96)
475          b(i,106) = b(i,106) - lu(i,519) * b(i,96)
476          b(i,115) = b(i,115) - lu(i,520) * b(i,96)
477          b(i,119) = b(i,119) - lu(i,521) * b(i,96)
478          b(i,120) = b(i,120) - lu(i,522) * b(i,96)
479          b(i,125) = b(i,125) - lu(i,523) * b(i,96)
480          b(i,129) = b(i,129) - lu(i,524) * b(i,96)
481          b(i,130) = b(i,130) - lu(i,525) * b(i,96)
482          b(i,131) = b(i,131) - lu(i,526) * b(i,96)
483          b(i,134) = b(i,134) - lu(i,527) * b(i,96)
484          b(i,135) = b(i,135) - lu(i,528) * b(i,96)
485          b(i,103) = b(i,103) - lu(i,531) * b(i,97)
486          b(i,110) = b(i,110) - lu(i,532) * b(i,97)
487          b(i,125) = b(i,125) - lu(i,533) * b(i,97)
488          b(i,130) = b(i,130) - lu(i,534) * b(i,97)
489          b(i,131) = b(i,131) - lu(i,535) * b(i,97)
490          b(i,132) = b(i,132) - lu(i,536) * b(i,97)
491          b(i,135) = b(i,135) - lu(i,537) * b(i,97)
492          b(i,106) = b(i,106) - lu(i,541) * b(i,98)
493          b(i,107) = b(i,107) - lu(i,542) * b(i,98)
494          b(i,113) = b(i,113) - lu(i,543) * b(i,98)
495          b(i,115) = b(i,115) - lu(i,544) * b(i,98)
496          b(i,119) = b(i,119) - lu(i,545) * b(i,98)
497          b(i,125) = b(i,125) - lu(i,546) * b(i,98)
498          b(i,129) = b(i,129) - lu(i,547) * b(i,98)
499          b(i,130) = b(i,130) - lu(i,548) * b(i,98)
500          b(i,131) = b(i,131) - lu(i,549) * b(i,98)
501          b(i,134) = b(i,134) - lu(i,550) * b(i,98)
510         SUBROUTINE lu_slv03_vec(ncol,nb,nz,lu, b)  argument
517             REAL(KIND=r8), intent(in) :: lu(ncol,nz)  local
529          b(i,116) = b(i,116) - lu(i,553) * b(i,99)
530          b(i,121) = b(i,121) - lu(i,554) * b(i,99)
531          b(i,125) = b(i,125) - lu(i,555) * b(i,99)
532          b(i,131) = b(i,131) - lu(i,556) * b(i,99)
533          b(i,134) = b(i,134) - lu(i,557) * b(i,99)
534          b(i,117) = b(i,117) - lu(i,561) * b(i,100)
535          b(i,121) = b(i,121) - lu(i,562) * b(i,100)
536          b(i,124) = b(i,124) - lu(i,563) * b(i,100)
537          b(i,126) = b(i,126) - lu(i,564) * b(i,100)
538          b(i,130) = b(i,130) - lu(i,565) * b(i,100)
539          b(i,131) = b(i,131) - lu(i,566) * b(i,100)
540          b(i,132) = b(i,132) - lu(i,567) * b(i,100)
541          b(i,133) = b(i,133) - lu(i,568) * b(i,100)
542          b(i,134) = b(i,134) - lu(i,569) * b(i,100)
543          b(i,103) = b(i,103) - lu(i,573) * b(i,101)
544          b(i,107) = b(i,107) - lu(i,574) * b(i,101)
545          b(i,110) = b(i,110) - lu(i,575) * b(i,101)
546          b(i,113) = b(i,113) - lu(i,576) * b(i,101)
547          b(i,125) = b(i,125) - lu(i,577) * b(i,101)
548          b(i,129) = b(i,129) - lu(i,578) * b(i,101)
549          b(i,130) = b(i,130) - lu(i,579) * b(i,101)
550          b(i,131) = b(i,131) - lu(i,580) * b(i,101)
551          b(i,132) = b(i,132) - lu(i,581) * b(i,101)
552          b(i,134) = b(i,134) - lu(i,582) * b(i,101)
553          b(i,135) = b(i,135) - lu(i,583) * b(i,101)
554          b(i,103) = b(i,103) - lu(i,588) * b(i,102)
555          b(i,104) = b(i,104) - lu(i,589) * b(i,102)
556          b(i,105) = b(i,105) - lu(i,590) * b(i,102)
557          b(i,109) = b(i,109) - lu(i,591) * b(i,102)
558          b(i,119) = b(i,119) - lu(i,592) * b(i,102)
559          b(i,120) = b(i,120) - lu(i,593) * b(i,102)
560          b(i,123) = b(i,123) - lu(i,594) * b(i,102)
561          b(i,125) = b(i,125) - lu(i,595) * b(i,102)
562          b(i,129) = b(i,129) - lu(i,596) * b(i,102)
563          b(i,130) = b(i,130) - lu(i,597) * b(i,102)
564          b(i,131) = b(i,131) - lu(i,598) * b(i,102)
565          b(i,132) = b(i,132) - lu(i,599) * b(i,102)
566          b(i,134) = b(i,134) - lu(i,600) * b(i,102)
567          b(i,135) = b(i,135) - lu(i,601) * b(i,102)
568          b(i,125) = b(i,125) - lu(i,603) * b(i,103)
569          b(i,127) = b(i,127) - lu(i,604) * b(i,103)
570          b(i,131) = b(i,131) - lu(i,605) * b(i,103)
571          b(i,115) = b(i,115) - lu(i,608) * b(i,104)
572          b(i,119) = b(i,119) - lu(i,609) * b(i,104)
573          b(i,125) = b(i,125) - lu(i,610) * b(i,104)
574          b(i,127) = b(i,127) - lu(i,611) * b(i,104)
575          b(i,131) = b(i,131) - lu(i,612) * b(i,104)
576          b(i,132) = b(i,132) - lu(i,613) * b(i,104)
577          b(i,133) = b(i,133) - lu(i,614) * b(i,104)
578          b(i,134) = b(i,134) - lu(i,615) * b(i,104)
579          b(i,109) = b(i,109) - lu(i,617) * b(i,105)
580          b(i,115) = b(i,115) - lu(i,618) * b(i,105)
581          b(i,125) = b(i,125) - lu(i,619) * b(i,105)
582          b(i,131) = b(i,131) - lu(i,620) * b(i,105)
583          b(i,135) = b(i,135) - lu(i,621) * b(i,105)
584          b(i,109) = b(i,109) - lu(i,626) * b(i,106)
585          b(i,115) = b(i,115) - lu(i,627) * b(i,106)
586          b(i,119) = b(i,119) - lu(i,628) * b(i,106)
587          b(i,120) = b(i,120) - lu(i,629) * b(i,106)
588          b(i,125) = b(i,125) - lu(i,630) * b(i,106)
589          b(i,129) = b(i,129) - lu(i,631) * b(i,106)
590          b(i,130) = b(i,130) - lu(i,632) * b(i,106)
591          b(i,131) = b(i,131) - lu(i,633) * b(i,106)
592          b(i,134) = b(i,134) - lu(i,634) * b(i,106)
593          b(i,135) = b(i,135) - lu(i,635) * b(i,106)
594          b(i,109) = b(i,109) - lu(i,638) * b(i,107)
595          b(i,112) = b(i,112) - lu(i,639) * b(i,107)
596          b(i,114) = b(i,114) - lu(i,640) * b(i,107)
597          b(i,115) = b(i,115) - lu(i,641) * b(i,107)
598          b(i,123) = b(i,123) - lu(i,642) * b(i,107)
599          b(i,125) = b(i,125) - lu(i,643) * b(i,107)
600          b(i,127) = b(i,127) - lu(i,644) * b(i,107)
601          b(i,131) = b(i,131) - lu(i,645) * b(i,107)
602          b(i,134) = b(i,134) - lu(i,646) * b(i,107)
603          b(i,135) = b(i,135) - lu(i,647) * b(i,107)
604          b(i,117) = b(i,117) - lu(i,651) * b(i,108)
605          b(i,119) = b(i,119) - lu(i,652) * b(i,108)
606          b(i,121) = b(i,121) - lu(i,653) * b(i,108)
607          b(i,122) = b(i,122) - lu(i,654) * b(i,108)
608          b(i,126) = b(i,126) - lu(i,655) * b(i,108)
609          b(i,127) = b(i,127) - lu(i,656) * b(i,108)
610          b(i,131) = b(i,131) - lu(i,657) * b(i,108)
611          b(i,132) = b(i,132) - lu(i,658) * b(i,108)
612          b(i,133) = b(i,133) - lu(i,659) * b(i,108)
613          b(i,134) = b(i,134) - lu(i,660) * b(i,108)
614          b(i,115) = b(i,115) - lu(i,663) * b(i,109)
615          b(i,125) = b(i,125) - lu(i,664) * b(i,109)
616          b(i,127) = b(i,127) - lu(i,665) * b(i,109)
617          b(i,131) = b(i,131) - lu(i,666) * b(i,109)
618          b(i,132) = b(i,132) - lu(i,667) * b(i,109)
619          b(i,133) = b(i,133) - lu(i,668) * b(i,109)
620          b(i,134) = b(i,134) - lu(i,669) * b(i,109)
621          b(i,115) = b(i,115) - lu(i,678) * b(i,110)
622          b(i,119) = b(i,119) - lu(i,679) * b(i,110)
623          b(i,125) = b(i,125) - lu(i,680) * b(i,110)
624          b(i,127) = b(i,127) - lu(i,681) * b(i,110)
625          b(i,129) = b(i,129) - lu(i,682) * b(i,110)
626          b(i,130) = b(i,130) - lu(i,683) * b(i,110)
627          b(i,131) = b(i,131) - lu(i,684) * b(i,110)
628          b(i,132) = b(i,132) - lu(i,685) * b(i,110)
629          b(i,133) = b(i,133) - lu(i,686) * b(i,110)
630          b(i,134) = b(i,134) - lu(i,687) * b(i,110)
631          b(i,135) = b(i,135) - lu(i,688) * b(i,110)
632          b(i,112) = b(i,112) - lu(i,698) * b(i,111)
633          b(i,113) = b(i,113) - lu(i,699) * b(i,111)
634          b(i,114) = b(i,114) - lu(i,700) * b(i,111)
635          b(i,115) = b(i,115) - lu(i,701) * b(i,111)
636          b(i,119) = b(i,119) - lu(i,702) * b(i,111)
637          b(i,123) = b(i,123) - lu(i,703) * b(i,111)
638          b(i,125) = b(i,125) - lu(i,704) * b(i,111)
639          b(i,127) = b(i,127) - lu(i,705) * b(i,111)
640          b(i,129) = b(i,129) - lu(i,706) * b(i,111)
641          b(i,130) = b(i,130) - lu(i,707) * b(i,111)
642          b(i,131) = b(i,131) - lu(i,708) * b(i,111)
643          b(i,132) = b(i,132) - lu(i,709) * b(i,111)
644          b(i,133) = b(i,133) - lu(i,710) * b(i,111)
645          b(i,134) = b(i,134) - lu(i,711) * b(i,111)
646          b(i,135) = b(i,135) - lu(i,712) * b(i,111)
647          b(i,114) = b(i,114) - lu(i,722) * b(i,112)
648          b(i,115) = b(i,115) - lu(i,723) * b(i,112)
649          b(i,119) = b(i,119) - lu(i,724) * b(i,112)
650          b(i,125) = b(i,125) - lu(i,725) * b(i,112)
651          b(i,127) = b(i,127) - lu(i,726) * b(i,112)
652          b(i,129) = b(i,129) - lu(i,727) * b(i,112)
653          b(i,130) = b(i,130) - lu(i,728) * b(i,112)
654          b(i,131) = b(i,131) - lu(i,729) * b(i,112)
655          b(i,132) = b(i,132) - lu(i,730) * b(i,112)
656          b(i,133) = b(i,133) - lu(i,731) * b(i,112)
657          b(i,134) = b(i,134) - lu(i,732) * b(i,112)
658          b(i,135) = b(i,135) - lu(i,733) * b(i,112)
659          b(i,114) = b(i,114) - lu(i,741) * b(i,113)
660          b(i,115) = b(i,115) - lu(i,742) * b(i,113)
661          b(i,119) = b(i,119) - lu(i,743) * b(i,113)
662          b(i,120) = b(i,120) - lu(i,744) * b(i,113)
663          b(i,123) = b(i,123) - lu(i,745) * b(i,113)
664          b(i,125) = b(i,125) - lu(i,746) * b(i,113)
665          b(i,127) = b(i,127) - lu(i,747) * b(i,113)
666          b(i,129) = b(i,129) - lu(i,748) * b(i,113)
667          b(i,130) = b(i,130) - lu(i,749) * b(i,113)
668          b(i,131) = b(i,131) - lu(i,750) * b(i,113)
669          b(i,132) = b(i,132) - lu(i,751) * b(i,113)
670          b(i,133) = b(i,133) - lu(i,752) * b(i,113)
671          b(i,134) = b(i,134) - lu(i,753) * b(i,113)
672          b(i,135) = b(i,135) - lu(i,754) * b(i,113)
673          b(i,115) = b(i,115) - lu(i,761) * b(i,114)
674          b(i,119) = b(i,119) - lu(i,762) * b(i,114)
675          b(i,120) = b(i,120) - lu(i,763) * b(i,114)
676          b(i,123) = b(i,123) - lu(i,764) * b(i,114)
677          b(i,125) = b(i,125) - lu(i,765) * b(i,114)
678          b(i,127) = b(i,127) - lu(i,766) * b(i,114)
679          b(i,129) = b(i,129) - lu(i,767) * b(i,114)
680          b(i,130) = b(i,130) - lu(i,768) * b(i,114)
681          b(i,131) = b(i,131) - lu(i,769) * b(i,114)
682          b(i,132) = b(i,132) - lu(i,770) * b(i,114)
683          b(i,133) = b(i,133) - lu(i,771) * b(i,114)
684          b(i,134) = b(i,134) - lu(i,772) * b(i,114)
685          b(i,135) = b(i,135) - lu(i,773) * b(i,114)
686          b(i,119) = b(i,119) - lu(i,790) * b(i,115)
687          b(i,120) = b(i,120) - lu(i,791) * b(i,115)
688          b(i,123) = b(i,123) - lu(i,792) * b(i,115)
689          b(i,125) = b(i,125) - lu(i,793) * b(i,115)
690          b(i,127) = b(i,127) - lu(i,794) * b(i,115)
691          b(i,129) = b(i,129) - lu(i,795) * b(i,115)
692          b(i,130) = b(i,130) - lu(i,796) * b(i,115)
693          b(i,131) = b(i,131) - lu(i,797) * b(i,115)
694          b(i,132) = b(i,132) - lu(i,798) * b(i,115)
695          b(i,133) = b(i,133) - lu(i,799) * b(i,115)
696          b(i,134) = b(i,134) - lu(i,800) * b(i,115)
697          b(i,135) = b(i,135) - lu(i,801) * b(i,115)
698          b(i,118) = b(i,118) - lu(i,806) * b(i,116)
699          b(i,120) = b(i,120) - lu(i,807) * b(i,116)
700          b(i,121) = b(i,121) - lu(i,808) * b(i,116)
701          b(i,123) = b(i,123) - lu(i,809) * b(i,116)
702          b(i,124) = b(i,124) - lu(i,810) * b(i,116)
703          b(i,125) = b(i,125) - lu(i,811) * b(i,116)
704          b(i,126) = b(i,126) - lu(i,812) * b(i,116)
705          b(i,127) = b(i,127) - lu(i,813) * b(i,116)
706          b(i,128) = b(i,128) - lu(i,814) * b(i,116)
707          b(i,129) = b(i,129) - lu(i,815) * b(i,116)
708          b(i,130) = b(i,130) - lu(i,816) * b(i,116)
709          b(i,131) = b(i,131) - lu(i,817) * b(i,116)
710          b(i,134) = b(i,134) - lu(i,818) * b(i,116)
711          b(i,118) = b(i,118) - lu(i,825) * b(i,117)
712          b(i,121) = b(i,121) - lu(i,826) * b(i,117)
713          b(i,122) = b(i,122) - lu(i,827) * b(i,117)
714          b(i,124) = b(i,124) - lu(i,828) * b(i,117)
715          b(i,126) = b(i,126) - lu(i,829) * b(i,117)
716          b(i,127) = b(i,127) - lu(i,830) * b(i,117)
717          b(i,128) = b(i,128) - lu(i,831) * b(i,117)
718          b(i,130) = b(i,130) - lu(i,832) * b(i,117)
719          b(i,131) = b(i,131) - lu(i,833) * b(i,117)
720          b(i,132) = b(i,132) - lu(i,834) * b(i,117)
721          b(i,133) = b(i,133) - lu(i,835) * b(i,117)
722          b(i,134) = b(i,134) - lu(i,836) * b(i,117)
723          b(i,120) = b(i,120) - lu(i,840) * b(i,118)
724          b(i,121) = b(i,121) - lu(i,841) * b(i,118)
725          b(i,122) = b(i,122) - lu(i,842) * b(i,118)
726          b(i,123) = b(i,123) - lu(i,843) * b(i,118)
727          b(i,125) = b(i,125) - lu(i,844) * b(i,118)
728          b(i,127) = b(i,127) - lu(i,845) * b(i,118)
729          b(i,128) = b(i,128) - lu(i,846) * b(i,118)
730          b(i,131) = b(i,131) - lu(i,847) * b(i,118)
731          b(i,134) = b(i,134) - lu(i,848) * b(i,118)
732          b(i,135) = b(i,135) - lu(i,849) * b(i,118)
741         SUBROUTINE lu_slv04_vec(ncol,nb,nz,lu, b)  argument
748             REAL(KIND=r8), intent(in) :: lu(ncol,nz)  local
760          b(i,120) = b(i,120) - lu(i,873) * b(i,119)
761          b(i,123) = b(i,123) - lu(i,874) * b(i,119)
762          b(i,124) = b(i,124) - lu(i,875) * b(i,119)
763          b(i,125) = b(i,125) - lu(i,876) * b(i,119)
764          b(i,126) = b(i,126) - lu(i,877) * b(i,119)
765          b(i,127) = b(i,127) - lu(i,878) * b(i,119)
766          b(i,129) = b(i,129) - lu(i,879) * b(i,119)
767          b(i,130) = b(i,130) - lu(i,880) * b(i,119)
768          b(i,131) = b(i,131) - lu(i,881) * b(i,119)
769          b(i,132) = b(i,132) - lu(i,882) * b(i,119)
770          b(i,133) = b(i,133) - lu(i,883) * b(i,119)
771          b(i,134) = b(i,134) - lu(i,884) * b(i,119)
772          b(i,135) = b(i,135) - lu(i,885) * b(i,119)
773          b(i,121) = b(i,121) - lu(i,904) * b(i,120)
774          b(i,122) = b(i,122) - lu(i,905) * b(i,120)
775          b(i,123) = b(i,123) - lu(i,906) * b(i,120)
776          b(i,124) = b(i,124) - lu(i,907) * b(i,120)
777          b(i,125) = b(i,125) - lu(i,908) * b(i,120)
778          b(i,126) = b(i,126) - lu(i,909) * b(i,120)
779          b(i,127) = b(i,127) - lu(i,910) * b(i,120)
780          b(i,128) = b(i,128) - lu(i,911) * b(i,120)
781          b(i,129) = b(i,129) - lu(i,912) * b(i,120)
782          b(i,130) = b(i,130) - lu(i,913) * b(i,120)
783          b(i,131) = b(i,131) - lu(i,914) * b(i,120)
784          b(i,134) = b(i,134) - lu(i,915) * b(i,120)
785          b(i,135) = b(i,135) - lu(i,916) * b(i,120)
786          b(i,122) = b(i,122) - lu(i,944) * b(i,121)
787          b(i,123) = b(i,123) - lu(i,945) * b(i,121)
788          b(i,124) = b(i,124) - lu(i,946) * b(i,121)
789          b(i,125) = b(i,125) - lu(i,947) * b(i,121)
790          b(i,126) = b(i,126) - lu(i,948) * b(i,121)
791          b(i,127) = b(i,127) - lu(i,949) * b(i,121)
792          b(i,128) = b(i,128) - lu(i,950) * b(i,121)
793          b(i,129) = b(i,129) - lu(i,951) * b(i,121)
794          b(i,130) = b(i,130) - lu(i,952) * b(i,121)
795          b(i,131) = b(i,131) - lu(i,953) * b(i,121)
796          b(i,132) = b(i,132) - lu(i,954) * b(i,121)
797          b(i,133) = b(i,133) - lu(i,955) * b(i,121)
798          b(i,134) = b(i,134) - lu(i,956) * b(i,121)
799          b(i,135) = b(i,135) - lu(i,957) * b(i,121)
800          b(i,123) = b(i,123) - lu(i,971) * b(i,122)
801          b(i,124) = b(i,124) - lu(i,972) * b(i,122)
802          b(i,125) = b(i,125) - lu(i,973) * b(i,122)
803          b(i,126) = b(i,126) - lu(i,974) * b(i,122)
804          b(i,127) = b(i,127) - lu(i,975) * b(i,122)
805          b(i,128) = b(i,128) - lu(i,976) * b(i,122)
806          b(i,129) = b(i,129) - lu(i,977) * b(i,122)
807          b(i,130) = b(i,130) - lu(i,978) * b(i,122)
808          b(i,131) = b(i,131) - lu(i,979) * b(i,122)
809          b(i,132) = b(i,132) - lu(i,980) * b(i,122)
810          b(i,133) = b(i,133) - lu(i,981) * b(i,122)
811          b(i,134) = b(i,134) - lu(i,982) * b(i,122)
812          b(i,135) = b(i,135) - lu(i,983) * b(i,122)
813          b(i,124) = b(i,124) - lu(i,1017) * b(i,123)
814          b(i,125) = b(i,125) - lu(i,1018) * b(i,123)
815          b(i,126) = b(i,126) - lu(i,1019) * b(i,123)
816          b(i,127) = b(i,127) - lu(i,1020) * b(i,123)
817          b(i,128) = b(i,128) - lu(i,1021) * b(i,123)
818          b(i,129) = b(i,129) - lu(i,1022) * b(i,123)
819          b(i,130) = b(i,130) - lu(i,1023) * b(i,123)
820          b(i,131) = b(i,131) - lu(i,1024) * b(i,123)
821          b(i,132) = b(i,132) - lu(i,1025) * b(i,123)
822          b(i,133) = b(i,133) - lu(i,1026) * b(i,123)
823          b(i,134) = b(i,134) - lu(i,1027) * b(i,123)
824          b(i,135) = b(i,135) - lu(i,1028) * b(i,123)
825          b(i,125) = b(i,125) - lu(i,1045) * b(i,124)
826          b(i,126) = b(i,126) - lu(i,1046) * b(i,124)
827          b(i,127) = b(i,127) - lu(i,1047) * b(i,124)
828          b(i,128) = b(i,128) - lu(i,1048) * b(i,124)
829          b(i,129) = b(i,129) - lu(i,1049) * b(i,124)
830          b(i,130) = b(i,130) - lu(i,1050) * b(i,124)
831          b(i,131) = b(i,131) - lu(i,1051) * b(i,124)
832          b(i,132) = b(i,132) - lu(i,1052) * b(i,124)
833          b(i,133) = b(i,133) - lu(i,1053) * b(i,124)
834          b(i,134) = b(i,134) - lu(i,1054) * b(i,124)
835          b(i,135) = b(i,135) - lu(i,1055) * b(i,124)
836          b(i,126) = b(i,126) - lu(i,1115) * b(i,125)
837          b(i,127) = b(i,127) - lu(i,1116) * b(i,125)
838          b(i,128) = b(i,128) - lu(i,1117) * b(i,125)
839          b(i,129) = b(i,129) - lu(i,1118) * b(i,125)
840          b(i,130) = b(i,130) - lu(i,1119) * b(i,125)
841          b(i,131) = b(i,131) - lu(i,1120) * b(i,125)
842          b(i,132) = b(i,132) - lu(i,1121) * b(i,125)
843          b(i,133) = b(i,133) - lu(i,1122) * b(i,125)
844          b(i,134) = b(i,134) - lu(i,1123) * b(i,125)
845          b(i,135) = b(i,135) - lu(i,1124) * b(i,125)
846          b(i,127) = b(i,127) - lu(i,1151) * b(i,126)
847          b(i,128) = b(i,128) - lu(i,1152) * b(i,126)
848          b(i,129) = b(i,129) - lu(i,1153) * b(i,126)
849          b(i,130) = b(i,130) - lu(i,1154) * b(i,126)
850          b(i,131) = b(i,131) - lu(i,1155) * b(i,126)
851          b(i,132) = b(i,132) - lu(i,1156) * b(i,126)
852          b(i,133) = b(i,133) - lu(i,1157) * b(i,126)
853          b(i,134) = b(i,134) - lu(i,1158) * b(i,126)
854          b(i,135) = b(i,135) - lu(i,1159) * b(i,126)
855          b(i,128) = b(i,128) - lu(i,1172) * b(i,127)
856          b(i,129) = b(i,129) - lu(i,1173) * b(i,127)
857          b(i,130) = b(i,130) - lu(i,1174) * b(i,127)
858          b(i,131) = b(i,131) - lu(i,1175) * b(i,127)
859          b(i,132) = b(i,132) - lu(i,1176) * b(i,127)
860          b(i,133) = b(i,133) - lu(i,1177) * b(i,127)
861          b(i,134) = b(i,134) - lu(i,1178) * b(i,127)
862          b(i,135) = b(i,135) - lu(i,1179) * b(i,127)
863          b(i,129) = b(i,129) - lu(i,1197) * b(i,128)
864          b(i,130) = b(i,130) - lu(i,1198) * b(i,128)
865          b(i,131) = b(i,131) - lu(i,1199) * b(i,128)
866          b(i,132) = b(i,132) - lu(i,1200) * b(i,128)
867          b(i,133) = b(i,133) - lu(i,1201) * b(i,128)
868          b(i,134) = b(i,134) - lu(i,1202) * b(i,128)
869          b(i,135) = b(i,135) - lu(i,1203) * b(i,128)
870          b(i,130) = b(i,130) - lu(i,1253) * b(i,129)
871          b(i,131) = b(i,131) - lu(i,1254) * b(i,129)
872          b(i,132) = b(i,132) - lu(i,1255) * b(i,129)
873          b(i,133) = b(i,133) - lu(i,1256) * b(i,129)
874          b(i,134) = b(i,134) - lu(i,1257) * b(i,129)
875          b(i,135) = b(i,135) - lu(i,1258) * b(i,129)
876          b(i,131) = b(i,131) - lu(i,1291) * b(i,130)
877          b(i,132) = b(i,132) - lu(i,1292) * b(i,130)
878          b(i,133) = b(i,133) - lu(i,1293) * b(i,130)
879          b(i,134) = b(i,134) - lu(i,1294) * b(i,130)
880          b(i,135) = b(i,135) - lu(i,1295) * b(i,130)
881          b(i,132) = b(i,132) - lu(i,1390) * b(i,131)
882          b(i,133) = b(i,133) - lu(i,1391) * b(i,131)
883          b(i,134) = b(i,134) - lu(i,1392) * b(i,131)
884          b(i,135) = b(i,135) - lu(i,1393) * b(i,131)
885          b(i,133) = b(i,133) - lu(i,1435) * b(i,132)
886          b(i,134) = b(i,134) - lu(i,1436) * b(i,132)
887          b(i,135) = b(i,135) - lu(i,1437) * b(i,132)
888          b(i,134) = b(i,134) - lu(i,1458) * b(i,133)
889          b(i,135) = b(i,135) - lu(i,1459) * b(i,133)
890          b(i,135) = b(i,135) - lu(i,1485) * b(i,134)
899         SUBROUTINE lu_slv05_vec(ncol,nb,nz,lu, b)  argument
906             REAL(KIND=r8), intent(in) :: lu(ncol,nz)  local
922          b(i,135) = b(i,135) * lu(i,1509)
923          b(i,134) = b(i,134) - lu(i,1508) * b(i,135)
924          b(i,133) = b(i,133) - lu(i,1507) * b(i,135)
925          b(i,132) = b(i,132) - lu(i,1506) * b(i,135)
926          b(i,131) = b(i,131) - lu(i,1505) * b(i,135)
927          b(i,130) = b(i,130) - lu(i,1504) * b(i,135)
928          b(i,129) = b(i,129) - lu(i,1503) * b(i,135)
929          b(i,128) = b(i,128) - lu(i,1502) * b(i,135)
930          b(i,127) = b(i,127) - lu(i,1501) * b(i,135)
931          b(i,126) = b(i,126) - lu(i,1500) * b(i,135)
932          b(i,125) = b(i,125) - lu(i,1499) * b(i,135)
933          b(i,124) = b(i,124) - lu(i,1498) * b(i,135)
934          b(i,123) = b(i,123) - lu(i,1497) * b(i,135)
935          b(i,122) = b(i,122) - lu(i,1496) * b(i,135)
936          b(i,121) = b(i,121) - lu(i,1495) * b(i,135)
937          b(i,120) = b(i,120) - lu(i,1494) * b(i,135)
938          b(i,119) = b(i,119) - lu(i,1493) * b(i,135)
939          b(i,118) = b(i,118) - lu(i,1492) * b(i,135)
940          b(i,117) = b(i,117) - lu(i,1491) * b(i,135)
941          b(i,108) = b(i,108) - lu(i,1490) * b(i,135)
942          b(i,103) = b(i,103) - lu(i,1489) * b(i,135)
943          b(i,90) = b(i,90) - lu(i,1488) * b(i,135)
944          b(i,64) = b(i,64) - lu(i,1487) * b(i,135)
945          b(i,54) = b(i,54) - lu(i,1486) * b(i,135)
946          b(i,134) = b(i,134) * lu(i,1484)
947          b(i,133) = b(i,133) - lu(i,1483) * b(i,134)
948          b(i,132) = b(i,132) - lu(i,1482) * b(i,134)
949          b(i,131) = b(i,131) - lu(i,1481) * b(i,134)
950          b(i,130) = b(i,130) - lu(i,1480) * b(i,134)
951          b(i,129) = b(i,129) - lu(i,1479) * b(i,134)
952          b(i,128) = b(i,128) - lu(i,1478) * b(i,134)
953          b(i,127) = b(i,127) - lu(i,1477) * b(i,134)
954          b(i,126) = b(i,126) - lu(i,1476) * b(i,134)
955          b(i,125) = b(i,125) - lu(i,1475) * b(i,134)
956          b(i,124) = b(i,124) - lu(i,1474) * b(i,134)
957          b(i,123) = b(i,123) - lu(i,1473) * b(i,134)
958          b(i,122) = b(i,122) - lu(i,1472) * b(i,134)
959          b(i,121) = b(i,121) - lu(i,1471) * b(i,134)
960          b(i,120) = b(i,120) - lu(i,1470) * b(i,134)
961          b(i,119) = b(i,119) - lu(i,1469) * b(i,134)
962          b(i,118) = b(i,118) - lu(i,1468) * b(i,134)
963          b(i,117) = b(i,117) - lu(i,1467) * b(i,134)
964          b(i,116) = b(i,116) - lu(i,1466) * b(i,134)
965          b(i,108) = b(i,108) - lu(i,1465) * b(i,134)
966          b(i,99) = b(i,99) - lu(i,1464) * b(i,134)
967          b(i,88) = b(i,88) - lu(i,1463) * b(i,134)
968          b(i,36) = b(i,36) - lu(i,1462) * b(i,134)
969          b(i,34) = b(i,34) - lu(i,1461) * b(i,134)
970          b(i,26) = b(i,26) - lu(i,1460) * b(i,134)
971          b(i,133) = b(i,133) * lu(i,1457)
972          b(i,132) = b(i,132) - lu(i,1456) * b(i,133)
973          b(i,131) = b(i,131) - lu(i,1455) * b(i,133)
974          b(i,130) = b(i,130) - lu(i,1454) * b(i,133)
975          b(i,129) = b(i,129) - lu(i,1453) * b(i,133)
976          b(i,128) = b(i,128) - lu(i,1452) * b(i,133)
977          b(i,127) = b(i,127) - lu(i,1451) * b(i,133)
978          b(i,126) = b(i,126) - lu(i,1450) * b(i,133)
979          b(i,125) = b(i,125) - lu(i,1449) * b(i,133)
980          b(i,124) = b(i,124) - lu(i,1448) * b(i,133)
981          b(i,123) = b(i,123) - lu(i,1447) * b(i,133)
982          b(i,122) = b(i,122) - lu(i,1446) * b(i,133)
983          b(i,121) = b(i,121) - lu(i,1445) * b(i,133)
984          b(i,120) = b(i,120) - lu(i,1444) * b(i,133)
985          b(i,119) = b(i,119) - lu(i,1443) * b(i,133)
986          b(i,118) = b(i,118) - lu(i,1442) * b(i,133)
987          b(i,117) = b(i,117) - lu(i,1441) * b(i,133)
988          b(i,108) = b(i,108) - lu(i,1440) * b(i,133)
989          b(i,88) = b(i,88) - lu(i,1439) * b(i,133)
990          b(i,34) = b(i,34) - lu(i,1438) * b(i,133)
991          b(i,132) = b(i,132) * lu(i,1434)
992          b(i,131) = b(i,131) - lu(i,1433) * b(i,132)
993          b(i,130) = b(i,130) - lu(i,1432) * b(i,132)
994          b(i,129) = b(i,129) - lu(i,1431) * b(i,132)
995          b(i,128) = b(i,128) - lu(i,1430) * b(i,132)
996          b(i,127) = b(i,127) - lu(i,1429) * b(i,132)
997          b(i,126) = b(i,126) - lu(i,1428) * b(i,132)
998          b(i,125) = b(i,125) - lu(i,1427) * b(i,132)
999          b(i,124) = b(i,124) - lu(i,1426) * b(i,132)
1000          b(i,123) = b(i,123) - lu(i,1425) * b(i,132)
1001          b(i,122) = b(i,122) - lu(i,1424) * b(i,132)
1002          b(i,121) = b(i,121) - lu(i,1423) * b(i,132)
1003          b(i,120) = b(i,120) - lu(i,1422) * b(i,132)
1004          b(i,119) = b(i,119) - lu(i,1421) * b(i,132)
1005          b(i,118) = b(i,118) - lu(i,1420) * b(i,132)
1006          b(i,116) = b(i,116) - lu(i,1419) * b(i,132)
1007          b(i,115) = b(i,115) - lu(i,1418) * b(i,132)
1008          b(i,114) = b(i,114) - lu(i,1417) * b(i,132)
1009          b(i,113) = b(i,113) - lu(i,1416) * b(i,132)
1010          b(i,112) = b(i,112) - lu(i,1415) * b(i,132)
1011          b(i,111) = b(i,111) - lu(i,1414) * b(i,132)
1012          b(i,110) = b(i,110) - lu(i,1413) * b(i,132)
1013          b(i,109) = b(i,109) - lu(i,1412) * b(i,132)
1014          b(i,107) = b(i,107) - lu(i,1411) * b(i,132)
1015          b(i,106) = b(i,106) - lu(i,1410) * b(i,132)
1016          b(i,105) = b(i,105) - lu(i,1409) * b(i,132)
1017          b(i,104) = b(i,104) - lu(i,1408) * b(i,132)
1018          b(i,103) = b(i,103) - lu(i,1407) * b(i,132)
1019          b(i,102) = b(i,102) - lu(i,1406) * b(i,132)
1020          b(i,101) = b(i,101) - lu(i,1405) * b(i,132)
1021          b(i,99) = b(i,99) - lu(i,1404) * b(i,132)
1022          b(i,98) = b(i,98) - lu(i,1403) * b(i,132)
1023          b(i,97) = b(i,97) - lu(i,1402) * b(i,132)
1024          b(i,95) = b(i,95) - lu(i,1401) * b(i,132)
1025          b(i,94) = b(i,94) - lu(i,1400) * b(i,132)
1026          b(i,81) = b(i,81) - lu(i,1399) * b(i,132)
1027          b(i,73) = b(i,73) - lu(i,1398) * b(i,132)
1028          b(i,49) = b(i,49) - lu(i,1397) * b(i,132)
1029          b(i,47) = b(i,47) - lu(i,1396) * b(i,132)
1030          b(i,40) = b(i,40) - lu(i,1395) * b(i,132)
1031          b(i,39) = b(i,39) - lu(i,1394) * b(i,132)
1032          b(i,131) = b(i,131) * lu(i,1389)
1033          b(i,130) = b(i,130) - lu(i,1388) * b(i,131)
1034          b(i,129) = b(i,129) - lu(i,1387) * b(i,131)
1035          b(i,128) = b(i,128) - lu(i,1386) * b(i,131)
1036          b(i,127) = b(i,127) - lu(i,1385) * b(i,131)
1037          b(i,126) = b(i,126) - lu(i,1384) * b(i,131)
1038          b(i,125) = b(i,125) - lu(i,1383) * b(i,131)
1039          b(i,124) = b(i,124) - lu(i,1382) * b(i,131)
1040          b(i,123) = b(i,123) - lu(i,1381) * b(i,131)
1041          b(i,122) = b(i,122) - lu(i,1380) * b(i,131)
1042          b(i,121) = b(i,121) - lu(i,1379) * b(i,131)
1043          b(i,120) = b(i,120) - lu(i,1378) * b(i,131)
1044          b(i,119) = b(i,119) - lu(i,1377) * b(i,131)
1045          b(i,118) = b(i,118) - lu(i,1376) * b(i,131)
1046          b(i,117) = b(i,117) - lu(i,1375) * b(i,131)
1047          b(i,116) = b(i,116) - lu(i,1374) * b(i,131)
1048          b(i,115) = b(i,115) - lu(i,1373) * b(i,131)
1049          b(i,114) = b(i,114) - lu(i,1372) * b(i,131)
1050          b(i,113) = b(i,113) - lu(i,1371) * b(i,131)
1051          b(i,112) = b(i,112) - lu(i,1370) * b(i,131)
1052          b(i,111) = b(i,111) - lu(i,1369) * b(i,131)
1053          b(i,110) = b(i,110) - lu(i,1368) * b(i,131)
1054          b(i,109) = b(i,109) - lu(i,1367) * b(i,131)
1055          b(i,108) = b(i,108) - lu(i,1366) * b(i,131)
1056          b(i,107) = b(i,107) - lu(i,1365) * b(i,131)
1057          b(i,106) = b(i,106) - lu(i,1364) * b(i,131)
1058          b(i,105) = b(i,105) - lu(i,1363) * b(i,131)
1059          b(i,104) = b(i,104) - lu(i,1362) * b(i,131)
1060          b(i,103) = b(i,103) - lu(i,1361) * b(i,131)
1061          b(i,102) = b(i,102) - lu(i,1360) * b(i,131)
1062          b(i,101) = b(i,101) - lu(i,1359) * b(i,131)
1063          b(i,100) = b(i,100) - lu(i,1358) * b(i,131)
1064          b(i,99) = b(i,99) - lu(i,1357) * b(i,131)
1065          b(i,98) = b(i,98) - lu(i,1356) * b(i,131)
1066          b(i,97) = b(i,97) - lu(i,1355) * b(i,131)
1067          b(i,96) = b(i,96) - lu(i,1354) * b(i,131)
1068          b(i,95) = b(i,95) - lu(i,1353) * b(i,131)
1069          b(i,94) = b(i,94) - lu(i,1352) * b(i,131)
1070          b(i,93) = b(i,93) - lu(i,1351) * b(i,131)
1071          b(i,92) = b(i,92) - lu(i,1350) * b(i,131)
1072          b(i,91) = b(i,91) - lu(i,1349) * b(i,131)
1073          b(i,90) = b(i,90) - lu(i,1348) * b(i,131)
1074          b(i,89) = b(i,89) - lu(i,1347) * b(i,131)
1075          b(i,88) = b(i,88) - lu(i,1346) * b(i,131)
1076          b(i,83) = b(i,83) - lu(i,1345) * b(i,131)
1077          b(i,82) = b(i,82) - lu(i,1344) * b(i,131)
1078          b(i,81) = b(i,81) - lu(i,1343) * b(i,131)
1079          b(i,80) = b(i,80) - lu(i,1342) * b(i,131)
1080          b(i,79) = b(i,79) - lu(i,1341) * b(i,131)
1081          b(i,77) = b(i,77) - lu(i,1340) * b(i,131)
1082          b(i,76) = b(i,76) - lu(i,1339) * b(i,131)
1083          b(i,75) = b(i,75) - lu(i,1338) * b(i,131)
1084          b(i,74) = b(i,74) - lu(i,1337) * b(i,131)
1085          b(i,73) = b(i,73) - lu(i,1336) * b(i,131)
1086          b(i,71) = b(i,71) - lu(i,1335) * b(i,131)
1087          b(i,69) = b(i,69) - lu(i,1334) * b(i,131)
1088          b(i,68) = b(i,68) - lu(i,1333) * b(i,131)
1089          b(i,67) = b(i,67) - lu(i,1332) * b(i,131)
1090          b(i,66) = b(i,66) - lu(i,1331) * b(i,131)
1091          b(i,65) = b(i,65) - lu(i,1330) * b(i,131)
1092          b(i,64) = b(i,64) - lu(i,1329) * b(i,131)
1093          b(i,63) = b(i,63) - lu(i,1328) * b(i,131)
1094          b(i,62) = b(i,62) - lu(i,1327) * b(i,131)
1095          b(i,60) = b(i,60) - lu(i,1326) * b(i,131)
1096          b(i,59) = b(i,59) - lu(i,1325) * b(i,131)
1097          b(i,57) = b(i,57) - lu(i,1324) * b(i,131)
1098          b(i,55) = b(i,55) - lu(i,1323) * b(i,131)
1099          b(i,53) = b(i,53) - lu(i,1322) * b(i,131)
1100          b(i,52) = b(i,52) - lu(i,1321) * b(i,131)
1101          b(i,51) = b(i,51) - lu(i,1320) * b(i,131)
1102          b(i,50) = b(i,50) - lu(i,1319) * b(i,131)
1103          b(i,49) = b(i,49) - lu(i,1318) * b(i,131)
1104          b(i,48) = b(i,48) - lu(i,1317) * b(i,131)
1105          b(i,47) = b(i,47) - lu(i,1316) * b(i,131)
1106          b(i,45) = b(i,45) - lu(i,1315) * b(i,131)
1107          b(i,44) = b(i,44) - lu(i,1314) * b(i,131)
1108          b(i,43) = b(i,43) - lu(i,1313) * b(i,131)
1109          b(i,42) = b(i,42) - lu(i,1312) * b(i,131)
1110          b(i,41) = b(i,41) - lu(i,1311) * b(i,131)
1111          b(i,39) = b(i,39) - lu(i,1310) * b(i,131)
1112          b(i,38) = b(i,38) - lu(i,1309) * b(i,131)
1113          b(i,37) = b(i,37) - lu(i,1308) * b(i,131)
1114          b(i,36) = b(i,36) - lu(i,1307) * b(i,131)
1115          b(i,35) = b(i,35) - lu(i,1306) * b(i,131)
1116          b(i,32) = b(i,32) - lu(i,1305) * b(i,131)
1117          b(i,31) = b(i,31) - lu(i,1304) * b(i,131)
1118          b(i,30) = b(i,30) - lu(i,1303) * b(i,131)
1119          b(i,25) = b(i,25) - lu(i,1302) * b(i,131)
1120          b(i,23) = b(i,23) - lu(i,1301) * b(i,131)
1121          b(i,22) = b(i,22) - lu(i,1300) * b(i,131)
1122          b(i,21) = b(i,21) - lu(i,1299) * b(i,131)
1123          b(i,20) = b(i,20) - lu(i,1298) * b(i,131)
1124          b(i,19) = b(i,19) - lu(i,1297) * b(i,131)
1125          b(i,17) = b(i,17) - lu(i,1296) * b(i,131)
1134         SUBROUTINE lu_slv06_vec(ncol,nb,nz,lu, b)  argument
1141             REAL(KIND=r8), intent(in) :: lu(ncol,nz)  local
1153          b(i,130) = b(i,130) * lu(i,1290)
1154          b(i,129) = b(i,129) - lu(i,1289) * b(i,130)
1155          b(i,128) = b(i,128) - lu(i,1288) * b(i,130)
1156          b(i,127) = b(i,127) - lu(i,1287) * b(i,130)
1157          b(i,126) = b(i,126) - lu(i,1286) * b(i,130)
1158          b(i,125) = b(i,125) - lu(i,1285) * b(i,130)
1159          b(i,124) = b(i,124) - lu(i,1284) * b(i,130)
1160          b(i,123) = b(i,123) - lu(i,1283) * b(i,130)
1161          b(i,122) = b(i,122) - lu(i,1282) * b(i,130)
1162          b(i,121) = b(i,121) - lu(i,1281) * b(i,130)
1163          b(i,120) = b(i,120) - lu(i,1280) * b(i,130)
1164          b(i,119) = b(i,119) - lu(i,1279) * b(i,130)
1165          b(i,118) = b(i,118) - lu(i,1278) * b(i,130)
1166          b(i,117) = b(i,117) - lu(i,1277) * b(i,130)
1167          b(i,116) = b(i,116) - lu(i,1276) * b(i,130)
1168          b(i,115) = b(i,115) - lu(i,1275) * b(i,130)
1169          b(i,114) = b(i,114) - lu(i,1274) * b(i,130)
1170          b(i,109) = b(i,109) - lu(i,1273) * b(i,130)
1171          b(i,105) = b(i,105) - lu(i,1272) * b(i,130)
1172          b(i,103) = b(i,103) - lu(i,1271) * b(i,130)
1173          b(i,100) = b(i,100) - lu(i,1270) * b(i,130)
1174          b(i,99) = b(i,99) - lu(i,1269) * b(i,130)
1175          b(i,92) = b(i,92) - lu(i,1268) * b(i,130)
1176          b(i,84) = b(i,84) - lu(i,1267) * b(i,130)
1177          b(i,81) = b(i,81) - lu(i,1266) * b(i,130)
1178          b(i,71) = b(i,71) - lu(i,1265) * b(i,130)
1179          b(i,70) = b(i,70) - lu(i,1264) * b(i,130)
1180          b(i,66) = b(i,66) - lu(i,1263) * b(i,130)
1181          b(i,60) = b(i,60) - lu(i,1262) * b(i,130)
1182          b(i,57) = b(i,57) - lu(i,1261) * b(i,130)
1183          b(i,40) = b(i,40) - lu(i,1260) * b(i,130)
1184          b(i,31) = b(i,31) - lu(i,1259) * b(i,130)
1185          b(i,129) = b(i,129) * lu(i,1252)
1186          b(i,128) = b(i,128) - lu(i,1251) * b(i,129)
1187          b(i,127) = b(i,127) - lu(i,1250) * b(i,129)
1188          b(i,126) = b(i,126) - lu(i,1249) * b(i,129)
1189          b(i,125) = b(i,125) - lu(i,1248) * b(i,129)
1190          b(i,124) = b(i,124) - lu(i,1247) * b(i,129)
1191          b(i,123) = b(i,123) - lu(i,1246) * b(i,129)
1192          b(i,122) = b(i,122) - lu(i,1245) * b(i,129)
1193          b(i,121) = b(i,121) - lu(i,1244) * b(i,129)
1194          b(i,120) = b(i,120) - lu(i,1243) * b(i,129)
1195          b(i,119) = b(i,119) - lu(i,1242) * b(i,129)
1196          b(i,118) = b(i,118) - lu(i,1241) * b(i,129)
1197          b(i,115) = b(i,115) - lu(i,1240) * b(i,129)
1198          b(i,114) = b(i,114) - lu(i,1239) * b(i,129)
1199          b(i,113) = b(i,113) - lu(i,1238) * b(i,129)
1200          b(i,112) = b(i,112) - lu(i,1237) * b(i,129)
1201          b(i,111) = b(i,111) - lu(i,1236) * b(i,129)
1202          b(i,110) = b(i,110) - lu(i,1235) * b(i,129)
1203          b(i,109) = b(i,109) - lu(i,1234) * b(i,129)
1204          b(i,107) = b(i,107) - lu(i,1233) * b(i,129)
1205          b(i,106) = b(i,106) - lu(i,1232) * b(i,129)
1206          b(i,105) = b(i,105) - lu(i,1231) * b(i,129)
1207          b(i,104) = b(i,104) - lu(i,1230) * b(i,129)
1208          b(i,103) = b(i,103) - lu(i,1229) * b(i,129)
1209          b(i,101) = b(i,101) - lu(i,1228) * b(i,129)
1210          b(i,98) = b(i,98) - lu(i,1227) * b(i,129)
1211          b(i,97) = b(i,97) - lu(i,1226) * b(i,129)
1212          b(i,96) = b(i,96) - lu(i,1225) * b(i,129)
1213          b(i,95) = b(i,95) - lu(i,1224) * b(i,129)
1214          b(i,92) = b(i,92) - lu(i,1223) * b(i,129)
1215          b(i,91) = b(i,91) - lu(i,1222) * b(i,129)
1216          b(i,89) = b(i,89) - lu(i,1221) * b(i,129)
1217          b(i,87) = b(i,87) - lu(i,1220) * b(i,129)
1218          b(i,86) = b(i,86) - lu(i,1219) * b(i,129)
1219          b(i,85) = b(i,85) - lu(i,1218) * b(i,129)
1220          b(i,83) = b(i,83) - lu(i,1217) * b(i,129)
1221          b(i,81) = b(i,81) - lu(i,1216) * b(i,129)
1222          b(i,80) = b(i,80) - lu(i,1215) * b(i,129)
1223          b(i,79) = b(i,79) - lu(i,1214) * b(i,129)
1224          b(i,77) = b(i,77) - lu(i,1213) * b(i,129)
1225          b(i,66) = b(i,66) - lu(i,1212) * b(i,129)
1226          b(i,65) = b(i,65) - lu(i,1211) * b(i,129)
1227          b(i,64) = b(i,64) - lu(i,1210) * b(i,129)
1228          b(i,56) = b(i,56) - lu(i,1209) * b(i,129)
1229          b(i,55) = b(i,55) - lu(i,1208) * b(i,129)
1230          b(i,54) = b(i,54) - lu(i,1207) * b(i,129)
1231          b(i,49) = b(i,49) - lu(i,1206) * b(i,129)
1232          b(i,47) = b(i,47) - lu(i,1205) * b(i,129)
1233          b(i,41) = b(i,41) - lu(i,1204) * b(i,129)
1234          b(i,128) = b(i,128) * lu(i,1196)
1235          b(i,127) = b(i,127) - lu(i,1195) * b(i,128)
1236          b(i,126) = b(i,126) - lu(i,1194) * b(i,128)
1237          b(i,125) = b(i,125) - lu(i,1193) * b(i,128)
1238          b(i,124) = b(i,124) - lu(i,1192) * b(i,128)
1239          b(i,123) = b(i,123) - lu(i,1191) * b(i,128)
1240          b(i,122) = b(i,122) - lu(i,1190) * b(i,128)
1241          b(i,121) = b(i,121) - lu(i,1189) * b(i,128)
1242          b(i,120) = b(i,120) - lu(i,1188) * b(i,128)
1243          b(i,118) = b(i,118) - lu(i,1187) * b(i,128)
1244          b(i,117) = b(i,117) - lu(i,1186) * b(i,128)
1245          b(i,116) = b(i,116) - lu(i,1185) * b(i,128)
1246          b(i,99) = b(i,99) - lu(i,1184) * b(i,128)
1247          b(i,84) = b(i,84) - lu(i,1183) * b(i,128)
1248          b(i,70) = b(i,70) - lu(i,1182) * b(i,128)
1249          b(i,46) = b(i,46) - lu(i,1181) * b(i,128)
1250          b(i,33) = b(i,33) - lu(i,1180) * b(i,128)
1251          b(i,127) = b(i,127) * lu(i,1171)
1252          b(i,126) = b(i,126) - lu(i,1170) * b(i,127)
1253          b(i,125) = b(i,125) - lu(i,1169) * b(i,127)
1254          b(i,124) = b(i,124) - lu(i,1168) * b(i,127)
1255          b(i,123) = b(i,123) - lu(i,1167) * b(i,127)
1256          b(i,122) = b(i,122) - lu(i,1166) * b(i,127)
1257          b(i,121) = b(i,121) - lu(i,1165) * b(i,127)
1258          b(i,120) = b(i,120) - lu(i,1164) * b(i,127)
1259          b(i,119) = b(i,119) - lu(i,1163) * b(i,127)
1260          b(i,118) = b(i,118) - lu(i,1162) * b(i,127)
1261          b(i,117) = b(i,117) - lu(i,1161) * b(i,127)
1262          b(i,108) = b(i,108) - lu(i,1160) * b(i,127)
1263          b(i,126) = b(i,126) * lu(i,1150)
1264          b(i,125) = b(i,125) - lu(i,1149) * b(i,126)
1265          b(i,124) = b(i,124) - lu(i,1148) * b(i,126)
1266          b(i,123) = b(i,123) - lu(i,1147) * b(i,126)
1267          b(i,122) = b(i,122) - lu(i,1146) * b(i,126)
1268          b(i,121) = b(i,121) - lu(i,1145) * b(i,126)
1269          b(i,120) = b(i,120) - lu(i,1144) * b(i,126)
1270          b(i,119) = b(i,119) - lu(i,1143) * b(i,126)
1271          b(i,118) = b(i,118) - lu(i,1142) * b(i,126)
1272          b(i,117) = b(i,117) - lu(i,1141) * b(i,126)
1273          b(i,115) = b(i,115) - lu(i,1140) * b(i,126)
1274          b(i,108) = b(i,108) - lu(i,1139) * b(i,126)
1275          b(i,104) = b(i,104) - lu(i,1138) * b(i,126)
1276          b(i,103) = b(i,103) - lu(i,1137) * b(i,126)
1277          b(i,100) = b(i,100) - lu(i,1136) * b(i,126)
1278          b(i,95) = b(i,95) - lu(i,1135) * b(i,126)
1279          b(i,93) = b(i,93) - lu(i,1134) * b(i,126)
1280          b(i,91) = b(i,91) - lu(i,1133) * b(i,126)
1281          b(i,83) = b(i,83) - lu(i,1132) * b(i,126)
1282          b(i,81) = b(i,81) - lu(i,1131) * b(i,126)
1283          b(i,74) = b(i,74) - lu(i,1130) * b(i,126)
1284          b(i,64) = b(i,64) - lu(i,1129) * b(i,126)
1285          b(i,63) = b(i,63) - lu(i,1128) * b(i,126)
1286          b(i,38) = b(i,38) - lu(i,1127) * b(i,126)
1287          b(i,37) = b(i,37) - lu(i,1126) * b(i,126)
1288          b(i,29) = b(i,29) - lu(i,1125) * b(i,126)
1289          b(i,125) = b(i,125) * lu(i,1114)
1290          b(i,124) = b(i,124) - lu(i,1113) * b(i,125)
1291          b(i,123) = b(i,123) - lu(i,1112) * b(i,125)
1292          b(i,122) = b(i,122) - lu(i,1111) * b(i,125)
1293          b(i,121) = b(i,121) - lu(i,1110) * b(i,125)
1294          b(i,120) = b(i,120) - lu(i,1109) * b(i,125)
1295          b(i,119) = b(i,119) - lu(i,1108) * b(i,125)
1296          b(i,118) = b(i,118) - lu(i,1107) * b(i,125)
1297          b(i,117) = b(i,117) - lu(i,1106) * b(i,125)
1298          b(i,115) = b(i,115) - lu(i,1105) * b(i,125)
1299          b(i,114) = b(i,114) - lu(i,1104) * b(i,125)
1300          b(i,113) = b(i,113) - lu(i,1103) * b(i,125)
1301          b(i,112) = b(i,112) - lu(i,1102) * b(i,125)
1302          b(i,111) = b(i,111) - lu(i,1101) * b(i,125)
1303          b(i,110) = b(i,110) - lu(i,1100) * b(i,125)
1304          b(i,109) = b(i,109) - lu(i,1099) * b(i,125)
1305          b(i,108) = b(i,108) - lu(i,1098) * b(i,125)
1306          b(i,107) = b(i,107) - lu(i,1097) * b(i,125)
1307          b(i,106) = b(i,106) - lu(i,1096) * b(i,125)
1308          b(i,105) = b(i,105) - lu(i,1095) * b(i,125)
1309          b(i,104) = b(i,104) - lu(i,1094) * b(i,125)
1310          b(i,103) = b(i,103) - lu(i,1093) * b(i,125)
1311          b(i,101) = b(i,101) - lu(i,1092) * b(i,125)
1312          b(i,98) = b(i,98) - lu(i,1091) * b(i,125)
1313          b(i,97) = b(i,97) - lu(i,1090) * b(i,125)
1314          b(i,96) = b(i,96) - lu(i,1089) * b(i,125)
1315          b(i,95) = b(i,95) - lu(i,1088) * b(i,125)
1316          b(i,93) = b(i,93) - lu(i,1087) * b(i,125)
1317          b(i,91) = b(i,91) - lu(i,1086) * b(i,125)
1318          b(i,90) = b(i,90) - lu(i,1085) * b(i,125)
1319          b(i,89) = b(i,89) - lu(i,1084) * b(i,125)
1320          b(i,84) = b(i,84) - lu(i,1083) * b(i,125)
1321          b(i,83) = b(i,83) - lu(i,1082) * b(i,125)
1322          b(i,81) = b(i,81) - lu(i,1081) * b(i,125)
1323          b(i,80) = b(i,80) - lu(i,1080) * b(i,125)
1324          b(i,79) = b(i,79) - lu(i,1079) * b(i,125)
1325          b(i,77) = b(i,77) - lu(i,1078) * b(i,125)
1326          b(i,76) = b(i,76) - lu(i,1077) * b(i,125)
1327          b(i,75) = b(i,75) - lu(i,1076) * b(i,125)
1328          b(i,74) = b(i,74) - lu(i,1075) * b(i,125)
1329          b(i,69) = b(i,69) - lu(i,1074) * b(i,125)
1330          b(i,67) = b(i,67) - lu(i,1073) * b(i,125)
1331          b(i,66) = b(i,66) - lu(i,1072) * b(i,125)
1332          b(i,65) = b(i,65) - lu(i,1071) * b(i,125)
1333          b(i,64) = b(i,64) - lu(i,1070) * b(i,125)
1334          b(i,62) = b(i,62) - lu(i,1069) * b(i,125)
1335          b(i,60) = b(i,60) - lu(i,1068) * b(i,125)
1336          b(i,59) = b(i,59) - lu(i,1067) * b(i,125)
1337          b(i,56) = b(i,56) - lu(i,1066) * b(i,125)
1338          b(i,54) = b(i,54) - lu(i,1065) * b(i,125)
1339          b(i,53) = b(i,53) - lu(i,1064) * b(i,125)
1340          b(i,52) = b(i,52) - lu(i,1063) * b(i,125)
1341          b(i,51) = b(i,51) - lu(i,1062) * b(i,125)
1342          b(i,50) = b(i,50) - lu(i,1061) * b(i,125)
1343          b(i,45) = b(i,45) - lu(i,1060) * b(i,125)
1344          b(i,44) = b(i,44) - lu(i,1059) * b(i,125)
1345          b(i,43) = b(i,43) - lu(i,1058) * b(i,125)
1346          b(i,42) = b(i,42) - lu(i,1057) * b(i,125)
1347          b(i,24) = b(i,24) - lu(i,1056) * b(i,125)
1348          b(i,124) = b(i,124) * lu(i,1044)
1349          b(i,123) = b(i,123) - lu(i,1043) * b(i,124)
1350          b(i,122) = b(i,122) - lu(i,1042) * b(i,124)
1351          b(i,121) = b(i,121) - lu(i,1041) * b(i,124)
1352          b(i,120) = b(i,120) - lu(i,1040) * b(i,124)
1353          b(i,119) = b(i,119) - lu(i,1039) * b(i,124)
1354          b(i,118) = b(i,118) - lu(i,1038) * b(i,124)
1355          b(i,117) = b(i,117) - lu(i,1037) * b(i,124)
1356          b(i,116) = b(i,116) - lu(i,1036) * b(i,124)
1357          b(i,100) = b(i,100) - lu(i,1035) * b(i,124)
1358          b(i,99) = b(i,99) - lu(i,1034) * b(i,124)
1359          b(i,93) = b(i,93) - lu(i,1033) * b(i,124)
1360          b(i,46) = b(i,46) - lu(i,1032) * b(i,124)
1361          b(i,33) = b(i,33) - lu(i,1031) * b(i,124)
1362          b(i,29) = b(i,29) - lu(i,1030) * b(i,124)
1363          b(i,18) = b(i,18) - lu(i,1029) * b(i,124)
1372         SUBROUTINE lu_slv07_vec(ncol,nb,nz,lu, b)  argument
1379             REAL(KIND=r8), intent(in) :: lu(ncol,nz)  local
1391          b(i,123) = b(i,123) * lu(i,1016)
1392          b(i,122) = b(i,122) - lu(i,1015) * b(i,123)
1393          b(i,121) = b(i,121) - lu(i,1014) * b(i,123)
1394          b(i,120) = b(i,120) - lu(i,1013) * b(i,123)
1395          b(i,119) = b(i,119) - lu(i,1012) * b(i,123)
1396          b(i,118) = b(i,118) - lu(i,1011) * b(i,123)
1397          b(i,116) = b(i,116) - lu(i,1010) * b(i,123)
1398          b(i,115) = b(i,115) - lu(i,1009) * b(i,123)
1399          b(i,114) = b(i,114) - lu(i,1008) * b(i,123)
1400          b(i,113) = b(i,113) - lu(i,1007) * b(i,123)
1401          b(i,112) = b(i,112) - lu(i,1006) * b(i,123)
1402          b(i,111) = b(i,111) - lu(i,1005) * b(i,123)
1403          b(i,110) = b(i,110) - lu(i,1004) * b(i,123)
1404          b(i,109) = b(i,109) - lu(i,1003) * b(i,123)
1405          b(i,107) = b(i,107) - lu(i,1002) * b(i,123)
1406          b(i,106) = b(i,106) - lu(i,1001) * b(i,123)
1407          b(i,105) = b(i,105) - lu(i,1000) * b(i,123)
1408          b(i,104) = b(i,104) - lu(i,999) * b(i,123)
1409          b(i,103) = b(i,103) - lu(i,998) * b(i,123)
1410          b(i,102) = b(i,102) - lu(i,997) * b(i,123)
1411          b(i,101) = b(i,101) - lu(i,996) * b(i,123)
1412          b(i,99) = b(i,99) - lu(i,995) * b(i,123)
1413          b(i,98) = b(i,98) - lu(i,994) * b(i,123)
1414          b(i,95) = b(i,95) - lu(i,993) * b(i,123)
1415          b(i,94) = b(i,94) - lu(i,992) * b(i,123)
1416          b(i,83) = b(i,83) - lu(i,991) * b(i,123)
1417          b(i,82) = b(i,82) - lu(i,990) * b(i,123)
1418          b(i,75) = b(i,75) - lu(i,989) * b(i,123)
1419          b(i,73) = b(i,73) - lu(i,988) * b(i,123)
1420          b(i,64) = b(i,64) - lu(i,987) * b(i,123)
1421          b(i,63) = b(i,63) - lu(i,986) * b(i,123)
1422          b(i,28) = b(i,28) - lu(i,985) * b(i,123)
1423          b(i,27) = b(i,27) - lu(i,984) * b(i,123)
1424          b(i,122) = b(i,122) * lu(i,970)
1425          b(i,121) = b(i,121) - lu(i,969) * b(i,122)
1426          b(i,120) = b(i,120) - lu(i,968) * b(i,122)
1427          b(i,119) = b(i,119) - lu(i,967) * b(i,122)
1428          b(i,118) = b(i,118) - lu(i,966) * b(i,122)
1429          b(i,117) = b(i,117) - lu(i,965) * b(i,122)
1430          b(i,108) = b(i,108) - lu(i,964) * b(i,122)
1431          b(i,90) = b(i,90) - lu(i,963) * b(i,122)
1432          b(i,88) = b(i,88) - lu(i,962) * b(i,122)
1433          b(i,32) = b(i,32) - lu(i,961) * b(i,122)
1434          b(i,30) = b(i,30) - lu(i,960) * b(i,122)
1435          b(i,28) = b(i,28) - lu(i,959) * b(i,122)
1436          b(i,25) = b(i,25) - lu(i,958) * b(i,122)
1437          b(i,121) = b(i,121) * lu(i,943)
1438          b(i,120) = b(i,120) - lu(i,942) * b(i,121)
1439          b(i,119) = b(i,119) - lu(i,941) * b(i,121)
1440          b(i,118) = b(i,118) - lu(i,940) * b(i,121)
1441          b(i,117) = b(i,117) - lu(i,939) * b(i,121)
1442          b(i,116) = b(i,116) - lu(i,938) * b(i,121)
1443          b(i,108) = b(i,108) - lu(i,937) * b(i,121)
1444          b(i,103) = b(i,103) - lu(i,936) * b(i,121)
1445          b(i,100) = b(i,100) - lu(i,935) * b(i,121)
1446          b(i,99) = b(i,99) - lu(i,934) * b(i,121)
1447          b(i,93) = b(i,93) - lu(i,933) * b(i,121)
1448          b(i,92) = b(i,92) - lu(i,932) * b(i,121)
1449          b(i,90) = b(i,90) - lu(i,931) * b(i,121)
1450          b(i,87) = b(i,87) - lu(i,930) * b(i,121)
1451          b(i,86) = b(i,86) - lu(i,929) * b(i,121)
1452          b(i,85) = b(i,85) - lu(i,928) * b(i,121)
1453          b(i,84) = b(i,84) - lu(i,927) * b(i,121)
1454          b(i,82) = b(i,82) - lu(i,926) * b(i,121)
1455          b(i,78) = b(i,78) - lu(i,925) * b(i,121)
1456          b(i,74) = b(i,74) - lu(i,924) * b(i,121)
1457          b(i,72) = b(i,72) - lu(i,923) * b(i,121)
1458          b(i,70) = b(i,70) - lu(i,922) * b(i,121)
1459          b(i,61) = b(i,61) - lu(i,921) * b(i,121)
1460          b(i,58) = b(i,58) - lu(i,920) * b(i,121)
1461          b(i,48) = b(i,48) - lu(i,919) * b(i,121)
1462          b(i,28) = b(i,28) - lu(i,918) * b(i,121)
1463          b(i,27) = b(i,27) - lu(i,917) * b(i,121)
1464          b(i,120) = b(i,120) * lu(i,903)
1465          b(i,118) = b(i,118) - lu(i,902) * b(i,120)
1466          b(i,116) = b(i,116) - lu(i,901) * b(i,120)
1467          b(i,103) = b(i,103) - lu(i,900) * b(i,120)
1468          b(i,99) = b(i,99) - lu(i,899) * b(i,120)
1469          b(i,95) = b(i,95) - lu(i,898) * b(i,120)
1470          b(i,92) = b(i,92) - lu(i,897) * b(i,120)
1471          b(i,87) = b(i,87) - lu(i,896) * b(i,120)
1472          b(i,86) = b(i,86) - lu(i,895) * b(i,120)
1473          b(i,85) = b(i,85) - lu(i,894) * b(i,120)
1474          b(i,82) = b(i,82) - lu(i,893) * b(i,120)
1475          b(i,78) = b(i,78) - lu(i,892) * b(i,120)
1476          b(i,72) = b(i,72) - lu(i,891) * b(i,120)
1477          b(i,61) = b(i,61) - lu(i,890) * b(i,120)
1478          b(i,58) = b(i,58) - lu(i,889) * b(i,120)
1479          b(i,56) = b(i,56) - lu(i,888) * b(i,120)
1480          b(i,28) = b(i,28) - lu(i,887) * b(i,120)
1481          b(i,27) = b(i,27) - lu(i,886) * b(i,120)
1482          b(i,119) = b(i,119) * lu(i,872)
1483          b(i,115) = b(i,115) - lu(i,871) * b(i,119)
1484          b(i,114) = b(i,114) - lu(i,870) * b(i,119)
1485          b(i,113) = b(i,113) - lu(i,869) * b(i,119)
1486          b(i,112) = b(i,112) - lu(i,868) * b(i,119)
1487          b(i,111) = b(i,111) - lu(i,867) * b(i,119)
1488          b(i,110) = b(i,110) - lu(i,866) * b(i,119)
1489          b(i,109) = b(i,109) - lu(i,865) * b(i,119)
1490          b(i,107) = b(i,107) - lu(i,864) * b(i,119)
1491          b(i,106) = b(i,106) - lu(i,863) * b(i,119)
1492          b(i,105) = b(i,105) - lu(i,862) * b(i,119)
1493          b(i,104) = b(i,104) - lu(i,861) * b(i,119)
1494          b(i,103) = b(i,103) - lu(i,860) * b(i,119)
1495          b(i,96) = b(i,96) - lu(i,859) * b(i,119)
1496          b(i,95) = b(i,95) - lu(i,858) * b(i,119)
1497          b(i,91) = b(i,91) - lu(i,857) * b(i,119)
1498          b(i,81) = b(i,81) - lu(i,856) * b(i,119)
1499          b(i,80) = b(i,80) - lu(i,855) * b(i,119)
1500          b(i,75) = b(i,75) - lu(i,854) * b(i,119)
1501          b(i,68) = b(i,68) - lu(i,853) * b(i,119)
1502          b(i,50) = b(i,50) - lu(i,852) * b(i,119)
1503          b(i,47) = b(i,47) - lu(i,851) * b(i,119)
1504          b(i,35) = b(i,35) - lu(i,850) * b(i,119)
1505          b(i,118) = b(i,118) * lu(i,839)
1506          b(i,103) = b(i,103) - lu(i,838) * b(i,118)
1507          b(i,90) = b(i,90) - lu(i,837) * b(i,118)
1508          b(i,117) = b(i,117) * lu(i,824)
1509          b(i,100) = b(i,100) - lu(i,823) * b(i,117)
1510          b(i,93) = b(i,93) - lu(i,822) * b(i,117)
1511          b(i,84) = b(i,84) - lu(i,821) * b(i,117)
1512          b(i,33) = b(i,33) - lu(i,820) * b(i,117)
1513          b(i,29) = b(i,29) - lu(i,819) * b(i,117)
1514          b(i,116) = b(i,116) * lu(i,805)
1515          b(i,99) = b(i,99) - lu(i,804) * b(i,116)
1516          b(i,82) = b(i,82) - lu(i,803) * b(i,116)
1517          b(i,46) = b(i,46) - lu(i,802) * b(i,116)
1518          b(i,115) = b(i,115) * lu(i,789)
1519          b(i,114) = b(i,114) - lu(i,788) * b(i,115)
1520          b(i,113) = b(i,113) - lu(i,787) * b(i,115)
1521          b(i,112) = b(i,112) - lu(i,786) * b(i,115)
1522          b(i,111) = b(i,111) - lu(i,785) * b(i,115)
1523          b(i,110) = b(i,110) - lu(i,784) * b(i,115)
1524          b(i,109) = b(i,109) - lu(i,783) * b(i,115)
1525          b(i,107) = b(i,107) - lu(i,782) * b(i,115)
1526          b(i,105) = b(i,105) - lu(i,781) * b(i,115)
1527          b(i,103) = b(i,103) - lu(i,780) * b(i,115)
1528          b(i,95) = b(i,95) - lu(i,779) * b(i,115)
1529          b(i,81) = b(i,81) - lu(i,778) * b(i,115)
1530          b(i,75) = b(i,75) - lu(i,777) * b(i,115)
1531          b(i,62) = b(i,62) - lu(i,776) * b(i,115)
1532          b(i,57) = b(i,57) - lu(i,775) * b(i,115)
1533          b(i,47) = b(i,47) - lu(i,774) * b(i,115)
1534          b(i,114) = b(i,114) * lu(i,760)
1535          b(i,109) = b(i,109) - lu(i,759) * b(i,114)
1536          b(i,105) = b(i,105) - lu(i,758) * b(i,114)
1537          b(i,75) = b(i,75) - lu(i,757) * b(i,114)
1538          b(i,71) = b(i,71) - lu(i,756) * b(i,114)
1539          b(i,62) = b(i,62) - lu(i,755) * b(i,114)
1540          b(i,113) = b(i,113) * lu(i,740)
1541          b(i,112) = b(i,112) - lu(i,739) * b(i,113)
1542          b(i,109) = b(i,109) - lu(i,738) * b(i,113)
1543          b(i,105) = b(i,105) - lu(i,737) * b(i,113)
1544          b(i,104) = b(i,104) - lu(i,736) * b(i,113)
1545          b(i,103) = b(i,103) - lu(i,735) * b(i,113)
1546          b(i,102) = b(i,102) - lu(i,734) * b(i,113)
1547          b(i,112) = b(i,112) * lu(i,721)
1548          b(i,110) = b(i,110) - lu(i,720) * b(i,112)
1549          b(i,109) = b(i,109) - lu(i,719) * b(i,112)
1550          b(i,105) = b(i,105) - lu(i,718) * b(i,112)
1551          b(i,103) = b(i,103) - lu(i,717) * b(i,112)
1552          b(i,97) = b(i,97) - lu(i,716) * b(i,112)
1553          b(i,95) = b(i,95) - lu(i,715) * b(i,112)
1554          b(i,68) = b(i,68) - lu(i,714) * b(i,112)
1555          b(i,43) = b(i,43) - lu(i,713) * b(i,112)
1556          b(i,111) = b(i,111) * lu(i,697)
1557          b(i,110) = b(i,110) - lu(i,696) * b(i,111)
1558          b(i,109) = b(i,109) - lu(i,695) * b(i,111)
1559          b(i,107) = b(i,107) - lu(i,694) * b(i,111)
1560          b(i,103) = b(i,103) - lu(i,693) * b(i,111)
1561          b(i,97) = b(i,97) - lu(i,692) * b(i,111)
1562          b(i,69) = b(i,69) - lu(i,691) * b(i,111)
1563          b(i,68) = b(i,68) - lu(i,690) * b(i,111)
1564          b(i,47) = b(i,47) - lu(i,689) * b(i,111)
1565          b(i,110) = b(i,110) * lu(i,677)
1566          b(i,109) = b(i,109) - lu(i,676) * b(i,110)
1567          b(i,105) = b(i,105) - lu(i,675) * b(i,110)
1568          b(i,103) = b(i,103) - lu(i,674) * b(i,110)
1569          b(i,95) = b(i,95) - lu(i,673) * b(i,110)
1570          b(i,81) = b(i,81) - lu(i,672) * b(i,110)
1571          b(i,68) = b(i,68) - lu(i,671) * b(i,110)
1572          b(i,45) = b(i,45) - lu(i,670) * b(i,110)
1573          b(i,109) = b(i,109) * lu(i,662)
1574          b(i,103) = b(i,103) - lu(i,661) * b(i,109)
1575          b(i,108) = b(i,108) * lu(i,650)
1576          b(i,88) = b(i,88) - lu(i,649) * b(i,108)
1577          b(i,34) = b(i,34) - lu(i,648) * b(i,108)
1578          b(i,107) = b(i,107) * lu(i,637)
1579          b(i,103) = b(i,103) - lu(i,636) * b(i,107)
1580          b(i,106) = b(i,106) * lu(i,625)
1581          b(i,105) = b(i,105) - lu(i,624) * b(i,106)
1582          b(i,68) = b(i,68) - lu(i,623) * b(i,106)
1583          b(i,53) = b(i,53) - lu(i,622) * b(i,106)
1584          b(i,105) = b(i,105) * lu(i,616)
1585          b(i,104) = b(i,104) * lu(i,607)
1586          b(i,103) = b(i,103) - lu(i,606) * b(i,104)
1587          b(i,103) = b(i,103) * lu(i,602)
1588          b(i,102) = b(i,102) * lu(i,587)
1589          b(i,89) = b(i,89) - lu(i,586) * b(i,102)
1590          b(i,75) = b(i,75) - lu(i,585) * b(i,102)
1591          b(i,49) = b(i,49) - lu(i,584) * b(i,102)
1600         SUBROUTINE lu_slv08_vec(ncol,nb,nz,lu, b)  argument
1607             REAL(KIND=r8), intent(in) :: lu(ncol,nz)  local
1619          b(i,101) = b(i,101) * lu(i,572)
1620          b(i,97) = b(i,97) - lu(i,571) * b(i,101)
1621          b(i,45) = b(i,45) - lu(i,570) * b(i,101)
1622          b(i,100) = b(i,100) * lu(i,560)
1623          b(i,93) = b(i,93) - lu(i,559) * b(i,100)
1624          b(i,29) = b(i,29) - lu(i,558) * b(i,100)
1625          b(i,99) = b(i,99) * lu(i,552)
1626          b(i,36) = b(i,36) - lu(i,551) * b(i,99)
1627          b(i,98) = b(i,98) * lu(i,540)
1628          b(i,80) = b(i,80) - lu(i,539) * b(i,98)
1629          b(i,59) = b(i,59) - lu(i,538) * b(i,98)
1630          b(i,97) = b(i,97) * lu(i,530)
1631          b(i,47) = b(i,47) - lu(i,529) * b(i,97)
1632          b(i,96) = b(i,96) * lu(i,517)
1633          b(i,80) = b(i,80) - lu(i,516) * b(i,96)
1634          b(i,52) = b(i,52) - lu(i,515) * b(i,96)
1635          b(i,95) = b(i,95) * lu(i,510)
1636          b(i,81) = b(i,81) - lu(i,509) * b(i,95)
1637          b(i,94) = b(i,94) * lu(i,494)
1638          b(i,75) = b(i,75) - lu(i,493) * b(i,94)
1639          b(i,93) = b(i,93) * lu(i,486)
1640          b(i,29) = b(i,29) - lu(i,485) * b(i,93)
1641          b(i,92) = b(i,92) * lu(i,476)
1642          b(i,87) = b(i,87) - lu(i,475) * b(i,92)
1643          b(i,86) = b(i,86) - lu(i,474) * b(i,92)
1644          b(i,85) = b(i,85) - lu(i,473) * b(i,92)
1645          b(i,72) = b(i,72) - lu(i,472) * b(i,92)
1646          b(i,58) = b(i,58) - lu(i,471) * b(i,92)
1647          b(i,91) = b(i,91) * lu(i,462)
1648          b(i,68) = b(i,68) - lu(i,461) * b(i,91)
1649          b(i,44) = b(i,44) - lu(i,460) * b(i,91)
1650          b(i,35) = b(i,35) - lu(i,459) * b(i,91)
1651          b(i,90) = b(i,90) * lu(i,452)
1652          b(i,89) = b(i,89) * lu(i,442)
1653          b(i,67) = b(i,67) - lu(i,441) * b(i,89)
1654          b(i,88) = b(i,88) * lu(i,433)
1655          b(i,34) = b(i,34) - lu(i,432) * b(i,88)
1656          b(i,87) = b(i,87) * lu(i,425)
1657          b(i,86) = b(i,86) - lu(i,424) * b(i,87)
1658          b(i,85) = b(i,85) - lu(i,423) * b(i,87)
1659          b(i,78) = b(i,78) - lu(i,422) * b(i,87)
1660          b(i,61) = b(i,61) - lu(i,421) * b(i,87)
1661          b(i,86) = b(i,86) * lu(i,414)
1662          b(i,61) = b(i,61) - lu(i,413) * b(i,86)
1663          b(i,85) = b(i,85) * lu(i,405)
1664          b(i,84) = b(i,84) * lu(i,397)
1665          b(i,33) = b(i,33) - lu(i,396) * b(i,84)
1666          b(i,83) = b(i,83) * lu(i,388)
1667          b(i,56) = b(i,56) - lu(i,387) * b(i,83)
1668          b(i,24) = b(i,24) - lu(i,386) * b(i,83)
1669          b(i,82) = b(i,82) * lu(i,379)
1670          b(i,81) = b(i,81) * lu(i,375)
1671          b(i,80) = b(i,80) * lu(i,369)
1672          b(i,79) = b(i,79) * lu(i,358)
1673          b(i,77) = b(i,77) - lu(i,357) * b(i,79)
1674          b(i,76) = b(i,76) - lu(i,356) * b(i,79)
1675          b(i,55) = b(i,55) - lu(i,355) * b(i,79)
1676          b(i,49) = b(i,49) - lu(i,354) * b(i,79)
1677          b(i,78) = b(i,78) * lu(i,344)
1678          b(i,72) = b(i,72) - lu(i,343) * b(i,78)
1679          b(i,61) = b(i,61) - lu(i,342) * b(i,78)
1680          b(i,77) = b(i,77) * lu(i,335)
1681          b(i,42) = b(i,42) - lu(i,334) * b(i,77)
1682          b(i,76) = b(i,76) * lu(i,324)
1683          b(i,55) = b(i,55) - lu(i,323) * b(i,76)
1684          b(i,75) = b(i,75) * lu(i,319)
1685          b(i,74) = b(i,74) * lu(i,312)
1686          b(i,73) = b(i,73) * lu(i,303)
1687          b(i,72) = b(i,72) * lu(i,296)
1688          b(i,71) = b(i,71) * lu(i,288)
1689          b(i,70) = b(i,70) * lu(i,280)
1690          b(i,69) = b(i,69) * lu(i,272)
1691          b(i,68) = b(i,68) * lu(i,268)
1692          b(i,67) = b(i,67) * lu(i,260)
1693          b(i,66) = b(i,66) * lu(i,254)
1694          b(i,65) = b(i,65) * lu(i,246)
1695          b(i,51) = b(i,51) - lu(i,245) * b(i,65)
1696          b(i,64) = b(i,64) * lu(i,241)
1697          b(i,63) = b(i,63) * lu(i,233)
1698          b(i,62) = b(i,62) * lu(i,227)
1699          b(i,61) = b(i,61) * lu(i,222)
1700          b(i,60) = b(i,60) * lu(i,215)
1701          b(i,59) = b(i,59) * lu(i,208)
1702          b(i,58) = b(i,58) * lu(i,201)
1703          b(i,57) = b(i,57) * lu(i,194)
1704          b(i,56) = b(i,56) * lu(i,189)
1705          b(i,55) = b(i,55) * lu(i,184)
1706          b(i,54) = b(i,54) * lu(i,178)
1707          b(i,53) = b(i,53) * lu(i,172)
1708          b(i,52) = b(i,52) * lu(i,166)
1709          b(i,51) = b(i,51) * lu(i,160)
1710          b(i,50) = b(i,50) * lu(i,154)
1711          b(i,49) = b(i,49) * lu(i,150)
1712          b(i,48) = b(i,48) * lu(i,142)
1713          b(i,47) = b(i,47) * lu(i,139)
1714          b(i,46) = b(i,46) * lu(i,134)
1715          b(i,45) = b(i,45) * lu(i,130)
1716          b(i,44) = b(i,44) * lu(i,125)
1717          b(i,43) = b(i,43) * lu(i,120)
1718          b(i,42) = b(i,42) * lu(i,115)
1719          b(i,41) = b(i,41) * lu(i,108)
1720          b(i,40) = b(i,40) * lu(i,102)
1721          b(i,39) = b(i,39) * lu(i,96)
1722          b(i,38) = b(i,38) * lu(i,90)
1723          b(i,37) = b(i,37) * lu(i,84)
1724          b(i,36) = b(i,36) * lu(i,80)
1725          b(i,26) = b(i,26) - lu(i,79) * b(i,36)
1726          b(i,35) = b(i,35) * lu(i,75)
1727          b(i,34) = b(i,34) * lu(i,72)
1728          b(i,33) = b(i,33) * lu(i,69)
1729          b(i,32) = b(i,32) * lu(i,65)
1730          b(i,31) = b(i,31) * lu(i,61)
1731          b(i,30) = b(i,30) * lu(i,57)
1732          b(i,29) = b(i,29) * lu(i,55)
1733          b(i,28) = b(i,28) * lu(i,53)
1734          b(i,27) = b(i,27) - lu(i,52) * b(i,28)
1735          b(i,27) = b(i,27) * lu(i,50)
1736          b(i,26) = b(i,26) * lu(i,47)
1737          b(i,25) = b(i,25) * lu(i,44)
1738          b(i,24) = b(i,24) * lu(i,41)
1739          b(i,23) = b(i,23) * lu(i,38)
1740          b(i,22) = b(i,22) * lu(i,33)
1741          b(i,21) = b(i,21) * lu(i,29)
1742          b(i,20) = b(i,20) * lu(i,26)
1743          b(i,19) = b(i,19) * lu(i,23)
1744          b(i,18) = b(i,18) * lu(i,20)
1745          b(i,17) = b(i,17) * lu(i,17)
1746          b(i,16) = b(i,16) * lu(i,16)
1747          b(i,15) = b(i,15) * lu(i,15)
1748          b(i,14) = b(i,14) * lu(i,14)
1749          b(i,13) = b(i,13) * lu(i,13)
1750          b(i,12) = b(i,12) * lu(i,12)
1751          b(i,11) = b(i,11) * lu(i,11)
1752          b(i,10) = b(i,10) * lu(i,10)
1753          b(i,9) = b(i,9) * lu(i,9)
1754          b(i,8) = b(i,8) * lu(i,8)
1755          b(i,7) = b(i,7) * lu(i,7)
1756          b(i,6) = b(i,6) * lu(i,6)
1757          b(i,5) = b(i,5) * lu(i,5)
1758          b(i,4) = b(i,4) * lu(i,4)
1759          b(i,3) = b(i,3) * lu(i,3)
1760          b(i,2) = b(i,2) * lu(i,2)
1761          b(i,1) = b(i,1) * lu(i,1)
1765         SUBROUTINE lu_slv_vec(ncol,nb,nz,lu, b)  argument
1772             REAL(KIND=r8), intent(in) :: lu(ncol,nz)  local
1774       call lu_slv01_vec( ncol,nb,nz,lu, b )
1775       call lu_slv02_vec( ncol,nb,nz,lu, b )
1776       call lu_slv03_vec( ncol,nb,nz,lu, b )
1777       call lu_slv04_vec( ncol,nb,nz,lu, b )
1778       call lu_slv05_vec( ncol,nb,nz,lu, b )
1779       call lu_slv06_vec( ncol,nb,nz,lu, b )
1780       call lu_slv07_vec( ncol,nb,nz,lu, b )
1781       call lu_slv08_vec( ncol,nb,nz,lu, b )