LCOV - code coverage report
Current view: top level - chemistry/pp_waccm_tsmlt_mam5 - mo_lu_solve.F90 (source / functions) Hit Total Coverage
Test: coverage.info Lines: 2499 2499 100.0 %
Date: 2025-03-14 01:23:43 Functions: 13 13 100.0 %

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

Generated by: LCOV version 1.14