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

Generated by: LCOV version 1.14