LCOV - code coverage report
Current view: top level - chemistry/pp_waccm_ma_noaero - mo_prod_loss.F90 (source / functions) Hit Total Coverage
Test: coverage.info Lines: 170 170 100.0 %
Date: 2025-03-14 01:33:33 Functions: 2 2 100.0 %

          Line data    Source code
       1             :       module mo_prod_loss
       2             :       use shr_kind_mod, only : r8 => shr_kind_r8
       3             :       use chem_mods, only : veclen
       4             :       private
       5             :       public :: exp_prod_loss
       6             :       public :: imp_prod_loss
       7             :       contains
       8       72960 :       subroutine exp_prod_loss( ofl, ofu, prod, loss, y, &
       9       72960 :                                 rxt, het_rates, chnkpnts )
      10             :       use chem_mods, only : gas_pcnst,rxntot,clscnt1
      11             :       implicit none
      12             : !--------------------------------------------------------------------
      13             : ! ... dummy args
      14             : !--------------------------------------------------------------------
      15             :       integer, intent(in) :: ofl, ofu, chnkpnts
      16             :       real(r8), dimension(chnkpnts,max(1,clscnt1)), intent(out) :: &
      17             :             prod, &
      18             :             loss
      19             :       real(r8), intent(in) :: y(chnkpnts,gas_pcnst)
      20             :       real(r8), intent(in) :: rxt(chnkpnts,rxntot)
      21             :       real(r8), intent(in) :: het_rates(chnkpnts,gas_pcnst)
      22             : !--------------------------------------------------------------------
      23             : ! ... local variables
      24             : !--------------------------------------------------------------------
      25             :       integer :: k
      26             : !--------------------------------------------------------------------
      27             : ! ... loss and production for Explicit method
      28             : !--------------------------------------------------------------------
      29    73616640 :       do k = ofl,ofu
      30    73543680 :          loss(k,1) = ( + het_rates(k,5))* y(k,5)
      31    73543680 :          prod(k,1) = 0._r8
      32    73543680 :          loss(k,2) = (rxt(k,180)* y(k,73) + rxt(k,30) + het_rates(k,6))* y(k,6)
      33    73543680 :          prod(k,2) = 0._r8
      34    73543680 :          loss(k,3) = (rxt(k,181)* y(k,73) + rxt(k,31) + het_rates(k,7))* y(k,7)
      35    73543680 :          prod(k,3) = 0._r8
      36    73543680 :          loss(k,4) = (rxt(k,207)* y(k,73) + rxt(k,32) + het_rates(k,8))* y(k,8)
      37    73543680 :          prod(k,4) = 0._r8
      38    73543680 :          loss(k,5) = (rxt(k,182)* y(k,73) + rxt(k,33) + het_rates(k,9))* y(k,9)
      39    73543680 :          prod(k,5) = 0._r8
      40    73543680 :          loss(k,6) = (rxt(k,183)* y(k,73) + rxt(k,34) + het_rates(k,10))* y(k,10)
      41    73543680 :          prod(k,6) = 0._r8
      42    73543680 :          loss(k,7) = (rxt(k,184)* y(k,73) + rxt(k,35) + het_rates(k,11))* y(k,11)
      43    73543680 :          prod(k,7) = 0._r8
      44    73543680 :          loss(k,8) = (rxt(k,185)* y(k,73) + rxt(k,36) + het_rates(k,12))* y(k,12)
      45    73543680 :          prod(k,8) = 0._r8
      46    73543680 :          loss(k,9) = (rxt(k,186)* y(k,73) + rxt(k,37) + het_rates(k,13))* y(k,13)
      47    73543680 :          prod(k,9) = 0._r8
      48             :          loss(k,10) = (rxt(k,218)* y(k,23) +rxt(k,230)* y(k,73) +rxt(k,219)* y(k,77) &
      49    73543680 :                   + rxt(k,38) + het_rates(k,14))* y(k,14)
      50    73543680 :          prod(k,10) = 0._r8
      51             :          loss(k,11) = (rxt(k,220)* y(k,23) +rxt(k,231)* y(k,73) +rxt(k,221)* y(k,77) &
      52    73543680 :                   + rxt(k,39) + het_rates(k,16))* y(k,16)
      53    73543680 :          prod(k,11) = 0._r8
      54    73543680 :          loss(k,12) = (rxt(k,222)* y(k,77) + rxt(k,40) + het_rates(k,17))* y(k,17)
      55    73543680 :          prod(k,12) = 0._r8
      56             :          loss(k,13) = (rxt(k,223)* y(k,23) +rxt(k,224)* y(k,77) + rxt(k,41) &
      57    73543680 :                   + het_rates(k,18))* y(k,18)
      58    73543680 :          prod(k,13) = 0._r8
      59             :          loss(k,14) = (rxt(k,156)* y(k,23) +rxt(k,212)* y(k,34) + (rxt(k,243) + &
      60             :                  rxt(k,244) +rxt(k,245))* y(k,73) +rxt(k,241)* y(k,77) + rxt(k,23) &
      61    73543680 :                   + rxt(k,24) + het_rates(k,21))* y(k,21)
      62    73543680 :          prod(k,14) = 0._r8
      63             :          loss(k,15) = (rxt(k,225)* y(k,23) +rxt(k,208)* y(k,73) +rxt(k,226)* y(k,77) &
      64    73543680 :                   + rxt(k,42) + het_rates(k,22))* y(k,22)
      65    73543680 :          prod(k,15) = 0._r8
      66    73543680 :          loss(k,16) = ( + het_rates(k,28))* y(k,28)
      67    73543680 :          prod(k,16) = 0._r8
      68             :          loss(k,17) = (rxt(k,300)* y(k,78) + rxt(k,25) + rxt(k,61) + het_rates(k,30)) &
      69    73543680 :                  * y(k,30)
      70    73543680 :          prod(k,17) =.440_r8*rxt(k,24)*y(k,21)
      71    73543680 :          loss(k,18) = (rxt(k,209)* y(k,73) + rxt(k,50) + het_rates(k,37))* y(k,37)
      72    73543680 :          prod(k,18) = 0._r8
      73             :          loss(k,19) = (rxt(k,232)* y(k,73) +rxt(k,227)* y(k,77) + rxt(k,52) &
      74    73543680 :                   + het_rates(k,41))* y(k,41)
      75    73543680 :          prod(k,19) = 0._r8
      76             :          loss(k,20) = (rxt(k,233)* y(k,73) +rxt(k,228)* y(k,77) + rxt(k,53) &
      77    73543680 :                   + het_rates(k,42))* y(k,42)
      78    73543680 :          prod(k,20) = 0._r8
      79             :          loss(k,21) = (rxt(k,234)* y(k,73) +rxt(k,229)* y(k,77) + rxt(k,54) &
      80    73543680 :                   + het_rates(k,43))* y(k,43)
      81    73543680 :          prod(k,21) = 0._r8
      82             :          loss(k,22) = ((rxt(k,147) +rxt(k,148))* y(k,73) + rxt(k,12) &
      83    73543680 :                   + het_rates(k,51))* y(k,51)
      84    73543680 :          prod(k,22) = 0._r8
      85    73543680 :          loss(k,23) = ( + rxt(k,60) + het_rates(k,62))* y(k,62)
      86    73616640 :          prod(k,23) = 0._r8
      87             :       end do
      88       72960 :       end subroutine exp_prod_loss
      89     9751000 :       subroutine imp_prod_loss( avec_len, prod, loss, y, &
      90             :                                 rxt, het_rates )
      91             :       use chem_mods, only : gas_pcnst,rxntot,clscnt4
      92             :       implicit none
      93             : !--------------------------------------------------------------------
      94             : ! ... dummy args
      95             : !--------------------------------------------------------------------
      96             :       integer, intent(in) :: avec_len
      97             :       real(r8), dimension(veclen,clscnt4), intent(out) :: &
      98             :             prod, &
      99             :             loss
     100             :       real(r8), intent(in) :: y(veclen,gas_pcnst)
     101             :       real(r8), intent(in) :: rxt(veclen,rxntot)
     102             :       real(r8), intent(in) :: het_rates(veclen,gas_pcnst)
     103             : !--------------------------------------------------------------------
     104             : ! ... local variables
     105             : !--------------------------------------------------------------------
     106             :       integer :: k
     107             : !--------------------------------------------------------------------
     108             : ! ... loss and production for Implicit method
     109             : !--------------------------------------------------------------------
     110   318792954 :       do k = 1,avec_len
     111   309041954 :          loss(k,41) = (rxt(k,191)* y(k,15) +rxt(k,193)* y(k,58) +rxt(k,192)* y(k,68) &
     112   309041954 :                   + het_rates(k,1))* y(k,1)
     113             :          prod(k,41) = (rxt(k,27) +2.000_r8*rxt(k,194)*y(k,3) +rxt(k,195)*y(k,26) + &
     114             :                  rxt(k,196)*y(k,26) +rxt(k,199)*y(k,53) +rxt(k,202)*y(k,56) + &
     115             :                  rxt(k,203)*y(k,77) +rxt(k,253)*y(k,63))*y(k,3) + (rxt(k,181)*y(k,7) + &
     116             :                  rxt(k,207)*y(k,8) +3.000_r8*rxt(k,208)*y(k,22) + &
     117             :                  2.000_r8*rxt(k,209)*y(k,37) +2.000_r8*rxt(k,230)*y(k,14) + &
     118             :                  rxt(k,231)*y(k,16) +rxt(k,210)*y(k,40))*y(k,73) &
     119             :                   + (2.000_r8*rxt(k,219)*y(k,14) +rxt(k,221)*y(k,16) + &
     120             :                  3.000_r8*rxt(k,226)*y(k,22) +rxt(k,205)*y(k,40))*y(k,77) &
     121             :                   + (2.000_r8*rxt(k,218)*y(k,14) +rxt(k,220)*y(k,16) + &
     122             :                  3.000_r8*rxt(k,225)*y(k,22))*y(k,23) + (rxt(k,51) + &
     123             :                  rxt(k,204)*y(k,56))*y(k,40) +rxt(k,26)*y(k,2) +rxt(k,29)*y(k,4) &
     124   309041954 :                   +rxt(k,57)*y(k,48)
     125   309041954 :          loss(k,9) = ( + rxt(k,26) + het_rates(k,2))* y(k,2)
     126             :          prod(k,9) = (rxt(k,276)*y(k,48) +rxt(k,281)*y(k,48))*y(k,44) &
     127   309041954 :                   +rxt(k,197)*y(k,26)*y(k,3)
     128             :          loss(k,57) = (2._r8*rxt(k,194)* y(k,3) + (rxt(k,195) +rxt(k,196) +rxt(k,197)) &
     129             :                  * y(k,26) +rxt(k,199)* y(k,53) +rxt(k,200)* y(k,54) +rxt(k,202) &
     130             :                  * y(k,56) +rxt(k,253)* y(k,63) +rxt(k,198)* y(k,68) +rxt(k,203) &
     131   309041954 :                  * y(k,77) + rxt(k,27) + het_rates(k,3))* y(k,3)
     132             :          prod(k,57) = (rxt(k,28) +rxt(k,201)*y(k,56))*y(k,4) +rxt(k,193)*y(k,58) &
     133   309041954 :                  *y(k,1) +rxt(k,211)*y(k,73)*y(k,40) +rxt(k,206)*y(k,56)*y(k,48)
     134             :          loss(k,20) = (rxt(k,201)* y(k,56) + rxt(k,28) + rxt(k,29) + rxt(k,270) &
     135   309041954 :                   + rxt(k,273) + rxt(k,278) + het_rates(k,4))* y(k,4)
     136   309041954 :          prod(k,20) =rxt(k,200)*y(k,54)*y(k,3)
     137             :          loss(k,51) = (rxt(k,191)* y(k,1) +rxt(k,155)* y(k,23) +rxt(k,235)* y(k,55) &
     138             :                   +rxt(k,236)* y(k,56) +rxt(k,237)* y(k,77) + rxt(k,20) + rxt(k,21) &
     139   309041954 :                   + het_rates(k,15))* y(k,15)
     140             :          prod(k,51) = (rxt(k,162)*y(k,26) +rxt(k,239)*y(k,53))*y(k,19) + (rxt(k,22) + &
     141             :                  .300_r8*rxt(k,240)*y(k,77))*y(k,20) + (rxt(k,244)*y(k,73) + &
     142   309041954 :                  rxt(k,245)*y(k,73))*y(k,21)
     143             :          loss(k,43) = (rxt(k,162)* y(k,26) +rxt(k,239)* y(k,53) +rxt(k,238)* y(k,68) &
     144   309041954 :                   + het_rates(k,19))* y(k,19)
     145             :          prod(k,43) = (rxt(k,156)*y(k,23) +rxt(k,212)*y(k,34) +rxt(k,241)*y(k,77) + &
     146   309041954 :                  rxt(k,243)*y(k,73))*y(k,21) +.700_r8*rxt(k,240)*y(k,77)*y(k,20)
     147   309041954 :          loss(k,16) = (rxt(k,240)* y(k,77) + rxt(k,22) + het_rates(k,20))* y(k,20)
     148   309041954 :          prod(k,16) =rxt(k,238)*y(k,68)*y(k,19)
     149             :          loss(k,48) = (rxt(k,218)* y(k,14) +rxt(k,155)* y(k,15) +rxt(k,220)* y(k,16) &
     150             :                   +rxt(k,223)* y(k,18) +rxt(k,156)* y(k,21) +rxt(k,225)* y(k,22) &
     151             :                   +rxt(k,168)* y(k,27) +rxt(k,157)* y(k,36) +rxt(k,158)* y(k,38) &
     152             :                   +rxt(k,177)* y(k,49) +rxt(k,161)* y(k,58) + (rxt(k,159) +rxt(k,160)) &
     153   309041954 :                  * y(k,68) + het_rates(k,23))* y(k,23)
     154             :          prod(k,48) = (4.000_r8*rxt(k,180)*y(k,6) +rxt(k,181)*y(k,7) + &
     155             :                  2.000_r8*rxt(k,182)*y(k,9) +2.000_r8*rxt(k,183)*y(k,10) + &
     156             :                  2.000_r8*rxt(k,184)*y(k,11) +rxt(k,185)*y(k,12) + &
     157             :                  2.000_r8*rxt(k,186)*y(k,13) +rxt(k,232)*y(k,41) +rxt(k,233)*y(k,42) + &
     158             :                  rxt(k,234)*y(k,43) +rxt(k,187)*y(k,44) +rxt(k,217)*y(k,32))*y(k,73) &
     159             :                   + (rxt(k,45) +rxt(k,162)*y(k,19) +2.000_r8*rxt(k,163)*y(k,26) + &
     160             :                  rxt(k,165)*y(k,26) +rxt(k,167)*y(k,53) +rxt(k,172)*y(k,56) + &
     161             :                  rxt(k,173)*y(k,77) +rxt(k,196)*y(k,3) +rxt(k,254)*y(k,63))*y(k,26) &
     162             :                   + (3.000_r8*rxt(k,222)*y(k,17) +rxt(k,224)*y(k,18) + &
     163             :                  rxt(k,227)*y(k,41) +rxt(k,228)*y(k,42) +rxt(k,229)*y(k,43) + &
     164             :                  rxt(k,176)*y(k,44))*y(k,77) + (rxt(k,55) +rxt(k,175)*y(k,56))*y(k,44) &
     165             :                   +rxt(k,26)*y(k,2) +2.000_r8*rxt(k,43)*y(k,24) +2.000_r8*rxt(k,44) &
     166   309041954 :                  *y(k,25) +rxt(k,46)*y(k,27) +rxt(k,49)*y(k,32) +rxt(k,58)*y(k,49)
     167   309041954 :          loss(k,7) = ( + rxt(k,43) + het_rates(k,24))* y(k,24)
     168             :          prod(k,7) = (rxt(k,269)*y(k,49) +rxt(k,274)*y(k,27) +rxt(k,275)*y(k,49) + &
     169             :                  rxt(k,279)*y(k,27) +rxt(k,280)*y(k,49) +rxt(k,284)*y(k,27))*y(k,44) &
     170   309041954 :                   +rxt(k,168)*y(k,27)*y(k,23) +rxt(k,164)*y(k,26)*y(k,26)
     171   309041954 :          loss(k,2) = ( + rxt(k,44) + rxt(k,190) + het_rates(k,25))* y(k,25)
     172   309041954 :          prod(k,2) =rxt(k,189)*y(k,26)*y(k,26)
     173             :          loss(k,47) = ((rxt(k,195) +rxt(k,196) +rxt(k,197))* y(k,3) +rxt(k,162) &
     174             :                  * y(k,19) + 2._r8*(rxt(k,163) +rxt(k,164) +rxt(k,165) +rxt(k,189)) &
     175             :                  * y(k,26) +rxt(k,167)* y(k,53) +rxt(k,169)* y(k,54) +rxt(k,172) &
     176             :                  * y(k,56) +rxt(k,254)* y(k,63) +rxt(k,166)* y(k,68) + (rxt(k,173) + &
     177   309041954 :                  rxt(k,174))* y(k,77) + rxt(k,45) + het_rates(k,26))* y(k,26)
     178             :          prod(k,47) = (rxt(k,160)*y(k,68) +rxt(k,161)*y(k,58) +rxt(k,177)*y(k,49)) &
     179             :                  *y(k,23) + (rxt(k,47) +rxt(k,170)*y(k,56))*y(k,27) &
     180             :                   + (rxt(k,178)*y(k,56) +rxt(k,179)*y(k,77))*y(k,49) + (rxt(k,59) + &
     181             :                  rxt(k,259)*y(k,63))*y(k,59) +2.000_r8*rxt(k,190)*y(k,25) &
     182   309041954 :                   +rxt(k,188)*y(k,73)*y(k,44)
     183             :          loss(k,37) = (rxt(k,168)* y(k,23) + (rxt(k,274) +rxt(k,279) +rxt(k,284)) &
     184             :                  * y(k,44) +rxt(k,170)* y(k,56) +rxt(k,171)* y(k,77) + rxt(k,46) &
     185             :                   + rxt(k,47) + rxt(k,272) + rxt(k,277) + rxt(k,283) &
     186   309041954 :                   + het_rates(k,27))* y(k,27)
     187   309041954 :          prod(k,37) =rxt(k,169)*y(k,54)*y(k,26)
     188   309041954 :          loss(k,21) = ((rxt(k,242) +rxt(k,246))* y(k,77) + het_rates(k,29))* y(k,29)
     189             :          prod(k,21) = (rxt(k,20) +rxt(k,21) +rxt(k,155)*y(k,23) +rxt(k,191)*y(k,1) + &
     190             :                  rxt(k,235)*y(k,55) +rxt(k,236)*y(k,56) +rxt(k,237)*y(k,77))*y(k,15) &
     191             :                   + (rxt(k,88) +rxt(k,249)*y(k,56) +rxt(k,250)*y(k,77))*y(k,60) &
     192   309041954 :                   +rxt(k,223)*y(k,23)*y(k,18) +rxt(k,300)*y(k,78)*y(k,30)
     193   309041954 :          loss(k,3) = (rxt(k,216)* y(k,73) + rxt(k,48) + het_rates(k,31))* y(k,31)
     194             :          prod(k,3) = (rxt(k,181)*y(k,7) +rxt(k,183)*y(k,10) + &
     195             :                  2.000_r8*rxt(k,184)*y(k,11) +2.000_r8*rxt(k,185)*y(k,12) + &
     196             :                  rxt(k,186)*y(k,13) +rxt(k,207)*y(k,8) +2.000_r8*rxt(k,209)*y(k,37) + &
     197             :                  rxt(k,233)*y(k,42) +rxt(k,234)*y(k,43))*y(k,73) &
     198   309041954 :                   + (rxt(k,228)*y(k,42) +rxt(k,229)*y(k,43))*y(k,77)
     199   309041954 :          loss(k,8) = (rxt(k,217)* y(k,73) + rxt(k,49) + het_rates(k,32))* y(k,32)
     200             :          prod(k,8) = (rxt(k,182)*y(k,9) +rxt(k,183)*y(k,10) +rxt(k,232)*y(k,41)) &
     201   309041954 :                  *y(k,73) +rxt(k,227)*y(k,77)*y(k,41)
     202             :          loss(k,12) = (rxt(k,247)* y(k,55) + (rxt(k,248) +rxt(k,261))* y(k,77) &
     203   309041954 :                   + het_rates(k,33))* y(k,33)
     204   309041954 :          prod(k,12) = 0._r8
     205             :          loss(k,26) = (rxt(k,212)* y(k,21) +rxt(k,213)* y(k,36) +rxt(k,215)* y(k,46) &
     206   309041954 :                   +rxt(k,214)* y(k,81) + het_rates(k,34))* y(k,34)
     207             :          prod(k,26) = (rxt(k,185)*y(k,12) +rxt(k,207)*y(k,8) + &
     208             :                  2.000_r8*rxt(k,216)*y(k,31) +rxt(k,217)*y(k,32))*y(k,73) &
     209   309041954 :                   +2.000_r8*rxt(k,48)*y(k,31) +rxt(k,49)*y(k,32) +rxt(k,56)*y(k,45)
     210             :          loss(k,40) = (rxt(k,116)* y(k,57) +rxt(k,119)* y(k,58) + (rxt(k,113) + &
     211   309041954 :                  rxt(k,114) +rxt(k,115))* y(k,68) + het_rates(k,35))* y(k,35)
     212             :          prod(k,40) = (rxt(k,120)*y(k,36) +rxt(k,123)*y(k,56) +rxt(k,143)*y(k,50) + &
     213             :                  rxt(k,237)*y(k,15) +rxt(k,246)*y(k,29) +rxt(k,250)*y(k,60) + &
     214             :                  rxt(k,255)*y(k,61) +rxt(k,260)*y(k,63))*y(k,77) &
     215             :                   + (rxt(k,94)*y(k,73) +rxt(k,111)*y(k,56) +rxt(k,157)*y(k,23) + &
     216             :                  rxt(k,213)*y(k,34))*y(k,36) + (rxt(k,244)*y(k,21) + &
     217             :                  rxt(k,188)*y(k,44) +rxt(k,211)*y(k,40))*y(k,73) &
     218             :                   + (2.000_r8*rxt(k,2) +rxt(k,3))*y(k,81) +2.000_r8*rxt(k,20)*y(k,15) &
     219             :                   +rxt(k,22)*y(k,20) +rxt(k,51)*y(k,40) +rxt(k,55)*y(k,44) +rxt(k,56) &
     220   309041954 :                  *y(k,45)
     221             :          loss(k,56) = (rxt(k,157)* y(k,23) +rxt(k,213)* y(k,34) +rxt(k,111)* y(k,56) &
     222   309041954 :                   +rxt(k,94)* y(k,73) +rxt(k,120)* y(k,77) + het_rates(k,36))* y(k,36)
     223             :          prod(k,56) =rxt(k,21)*y(k,15) +rxt(k,245)*y(k,73)*y(k,21) +rxt(k,113)*y(k,68) &
     224   309041954 :                  *y(k,35) +rxt(k,1)*y(k,81)
     225             :          loss(k,22) = (rxt(k,158)* y(k,23) +rxt(k,112)* y(k,56) +rxt(k,121)* y(k,77) &
     226   309041954 :                   + rxt(k,4) + het_rates(k,38))* y(k,38)
     227   309041954 :          prod(k,22) =rxt(k,127)*y(k,68)*y(k,68) +rxt(k,126)*y(k,77)*y(k,77)
     228   309041954 :          loss(k,4) = ( + rxt(k,87) + het_rates(k,39))* y(k,39)
     229   309041954 :          prod(k,4) =rxt(k,263)*y(k,81)*y(k,65)
     230             :          loss(k,29) = (rxt(k,204)* y(k,56) + (rxt(k,210) +rxt(k,211))* y(k,73) &
     231   309041954 :                   +rxt(k,205)* y(k,77) + rxt(k,51) + het_rates(k,40))* y(k,40)
     232   309041954 :          prod(k,29) = (rxt(k,191)*y(k,15) +rxt(k,192)*y(k,68))*y(k,1)
     233             :          loss(k,42) = ((rxt(k,274) +rxt(k,279) +rxt(k,284))* y(k,27) + (rxt(k,276) + &
     234             :                  rxt(k,281))* y(k,48) + (rxt(k,269) +rxt(k,275) +rxt(k,280))* y(k,49) &
     235             :                   +rxt(k,175)* y(k,56) + (rxt(k,187) +rxt(k,188))* y(k,73) +rxt(k,176) &
     236   309041954 :                  * y(k,77) + rxt(k,55) + het_rates(k,44))* y(k,44)
     237             :          prod(k,42) = (rxt(k,156)*y(k,21) +rxt(k,218)*y(k,14) +rxt(k,220)*y(k,16) + &
     238             :                  2.000_r8*rxt(k,223)*y(k,18) +rxt(k,225)*y(k,22) +rxt(k,155)*y(k,15) + &
     239             :                  rxt(k,157)*y(k,36) +rxt(k,158)*y(k,38) +rxt(k,159)*y(k,68) + &
     240   309041954 :                  rxt(k,177)*y(k,49))*y(k,23) +rxt(k,174)*y(k,77)*y(k,26)
     241   309041954 :          loss(k,10) = ( + rxt(k,56) + het_rates(k,45))* y(k,45)
     242             :          prod(k,10) = (rxt(k,212)*y(k,21) +rxt(k,213)*y(k,36) +rxt(k,214)*y(k,81) + &
     243   309041954 :                  rxt(k,215)*y(k,46))*y(k,34)
     244             :          loss(k,38) = (rxt(k,215)* y(k,34) +rxt(k,152)* y(k,77) + rxt(k,9) &
     245   309041954 :                   + het_rates(k,46))* y(k,46)
     246             :          prod(k,38) = (rxt(k,272) +rxt(k,277) +rxt(k,283) +rxt(k,274)*y(k,44) + &
     247             :                  rxt(k,279)*y(k,44) +rxt(k,284)*y(k,44))*y(k,27) &
     248             :                   + (2.000_r8*rxt(k,265) +2.000_r8*rxt(k,268) +2.000_r8*rxt(k,271) + &
     249             :                  2.000_r8*rxt(k,282))*y(k,52) + (rxt(k,270) +rxt(k,273) +rxt(k,278)) &
     250             :                  *y(k,4) + (rxt(k,267) +rxt(k,235)*y(k,15) +rxt(k,247)*y(k,33)) &
     251   309041954 :                  *y(k,55) + (.500_r8*rxt(k,266) +rxt(k,151)*y(k,77))*y(k,54)
     252             :          loss(k,17) = (rxt(k,128)* y(k,77) + rxt(k,10) + rxt(k,11) + rxt(k,153) &
     253   309041954 :                   + het_rates(k,47))* y(k,47)
     254   309041954 :          prod(k,17) =rxt(k,149)*y(k,68)*y(k,54)
     255             :          loss(k,27) = ((rxt(k,276) +rxt(k,281))* y(k,44) +rxt(k,206)* y(k,56) &
     256   309041954 :                   + rxt(k,57) + het_rates(k,48))* y(k,48)
     257             :          prod(k,27) = (rxt(k,270) +rxt(k,273) +rxt(k,278))*y(k,4) +rxt(k,198)*y(k,68) &
     258   309041954 :                  *y(k,3)
     259             :          loss(k,28) = (rxt(k,177)* y(k,23) + (rxt(k,269) +rxt(k,275) +rxt(k,280)) &
     260             :                  * y(k,44) +rxt(k,178)* y(k,56) +rxt(k,179)* y(k,77) + rxt(k,58) &
     261   309041954 :                   + het_rates(k,49))* y(k,49)
     262             :          prod(k,28) = (rxt(k,272) +rxt(k,277) +rxt(k,283) +rxt(k,171)*y(k,77))*y(k,27) &
     263   309041954 :                   +rxt(k,166)*y(k,68)*y(k,26)
     264             :          loss(k,35) = (rxt(k,131)* y(k,53) + (rxt(k,132) +rxt(k,133) +rxt(k,134)) &
     265             :                  * y(k,54) +rxt(k,135)* y(k,57) +rxt(k,297)* y(k,76) +rxt(k,143) &
     266   309041954 :                  * y(k,77) + rxt(k,62) + het_rates(k,50))* y(k,50)
     267             :          prod(k,35) = (rxt(k,129)*y(k,69) +rxt(k,294)*y(k,72))*y(k,56) &
     268             :                   + (.200_r8*rxt(k,288)*y(k,71) +1.100_r8*rxt(k,290)*y(k,70))*y(k,67) &
     269   309041954 :                   +rxt(k,15)*y(k,53) +rxt(k,295)*y(k,72)*y(k,57) +rxt(k,301)*y(k,78)
     270             :          loss(k,13) = ( + rxt(k,13) + rxt(k,14) + rxt(k,154) + rxt(k,265) + rxt(k,268) &
     271   309041954 :                   + rxt(k,271) + rxt(k,282) + het_rates(k,52))* y(k,52)
     272   309041954 :          prod(k,13) =rxt(k,150)*y(k,55)*y(k,54)
     273             :          loss(k,55) = (rxt(k,199)* y(k,3) +rxt(k,239)* y(k,19) +rxt(k,167)* y(k,26) &
     274             :                   +rxt(k,131)* y(k,50) +rxt(k,140)* y(k,55) +rxt(k,146)* y(k,56) &
     275             :                   +rxt(k,145)* y(k,58) +rxt(k,144)* y(k,68) +rxt(k,299)* y(k,76) &
     276   309041954 :                   + rxt(k,15) + rxt(k,16) + het_rates(k,53))* y(k,53)
     277             :          prod(k,55) = (rxt(k,17) +.500_r8*rxt(k,266) +2.000_r8*rxt(k,133)*y(k,50) + &
     278             :                  rxt(k,136)*y(k,56) +rxt(k,256)*y(k,63))*y(k,54) &
     279             :                   + (rxt(k,135)*y(k,57) +rxt(k,143)*y(k,77))*y(k,50) &
     280             :                   +2.000_r8*rxt(k,147)*y(k,73)*y(k,51) +rxt(k,14)*y(k,52) +rxt(k,19) &
     281             :                  *y(k,55) +rxt(k,130)*y(k,69)*y(k,57) +rxt(k,298)*y(k,76) +rxt(k,311) &
     282   309041954 :                  *y(k,80)
     283             :          loss(k,53) = (rxt(k,200)* y(k,3) +rxt(k,169)* y(k,26) + (rxt(k,132) + &
     284             :                  rxt(k,133) +rxt(k,134))* y(k,50) +rxt(k,150)* y(k,55) + (rxt(k,136) + &
     285             :                  rxt(k,138))* y(k,56) +rxt(k,137)* y(k,58) +rxt(k,256)* y(k,63) &
     286             :                   +rxt(k,149)* y(k,68) +rxt(k,151)* y(k,77) + rxt(k,17) + rxt(k,266) &
     287   309041954 :                   + het_rates(k,54))* y(k,54)
     288             :          prod(k,53) = (2.000_r8*rxt(k,140)*y(k,55) +rxt(k,144)*y(k,68) + &
     289             :                  rxt(k,145)*y(k,58) +rxt(k,146)*y(k,56) +rxt(k,167)*y(k,26) + &
     290             :                  rxt(k,199)*y(k,3) +rxt(k,239)*y(k,19))*y(k,53) + (rxt(k,18) + &
     291             :                  rxt(k,139)*y(k,68) +rxt(k,141)*y(k,56) +rxt(k,142)*y(k,77))*y(k,55) &
     292             :                   + (rxt(k,11) +rxt(k,153) +rxt(k,128)*y(k,77))*y(k,47) + (rxt(k,13) + &
     293             :                  rxt(k,154))*y(k,52) +rxt(k,28)*y(k,4) +rxt(k,47)*y(k,27) +rxt(k,9) &
     294   309041954 :                  *y(k,46)
     295             :          loss(k,44) = (rxt(k,235)* y(k,15) +rxt(k,247)* y(k,33) +rxt(k,140)* y(k,53) &
     296             :                   +rxt(k,150)* y(k,54) +rxt(k,141)* y(k,56) +rxt(k,139)* y(k,68) &
     297             :                   +rxt(k,142)* y(k,77) + rxt(k,18) + rxt(k,19) + rxt(k,267) &
     298   309041954 :                   + het_rates(k,55))* y(k,55)
     299             :          prod(k,44) = (rxt(k,46) +rxt(k,168)*y(k,23) +rxt(k,170)*y(k,56) + &
     300             :                  rxt(k,171)*y(k,77))*y(k,27) + (rxt(k,13) +rxt(k,14) +rxt(k,154)) &
     301             :                  *y(k,52) + (rxt(k,29) +rxt(k,201)*y(k,56))*y(k,4) &
     302             :                   + (rxt(k,152)*y(k,77) +rxt(k,215)*y(k,34))*y(k,46) &
     303             :                   + (rxt(k,137)*y(k,58) +rxt(k,138)*y(k,56))*y(k,54) +rxt(k,10) &
     304   309041954 :                  *y(k,47)
     305             :          loss(k,46) = (rxt(k,202)* y(k,3) +rxt(k,201)* y(k,4) +rxt(k,236)* y(k,15) &
     306             :                   +rxt(k,172)* y(k,26) +rxt(k,170)* y(k,27) +rxt(k,111)* y(k,36) &
     307             :                   +rxt(k,112)* y(k,38) +rxt(k,204)* y(k,40) +rxt(k,175)* y(k,44) &
     308             :                   +rxt(k,206)* y(k,48) +rxt(k,178)* y(k,49) +rxt(k,146)* y(k,53) &
     309             :                   + (rxt(k,136) +rxt(k,138))* y(k,54) +rxt(k,141)* y(k,55) &
     310             :                   + 2._r8*rxt(k,109)* y(k,56) +rxt(k,110)* y(k,57) +rxt(k,108) &
     311             :                  * y(k,58) +rxt(k,249)* y(k,60) +rxt(k,117)* y(k,68) + (rxt(k,292) + &
     312             :                  rxt(k,293))* y(k,70) +rxt(k,294)* y(k,72) +rxt(k,123)* y(k,77) &
     313             :                   + rxt(k,71) + rxt(k,72) + rxt(k,73) + rxt(k,74) + rxt(k,75) &
     314   309041954 :                   + rxt(k,76) + het_rates(k,56))* y(k,56)
     315             :          prod(k,46) = (2.000_r8*rxt(k,5) +rxt(k,6) +rxt(k,77) +rxt(k,79) +rxt(k,81) + &
     316             :                  2.000_r8*rxt(k,82) +2.000_r8*rxt(k,83) +rxt(k,84) +rxt(k,85) + &
     317             :                  rxt(k,86) +rxt(k,97)*y(k,73) +rxt(k,98)*y(k,73) +rxt(k,135)*y(k,50) + &
     318             :                  rxt(k,251)*y(k,61) +rxt(k,257)*y(k,63) +rxt(k,296)*y(k,72) + &
     319             :                  rxt(k,303)*y(k,78) +rxt(k,307)*y(k,79))*y(k,57) &
     320             :                   + (rxt(k,131)*y(k,53) +rxt(k,132)*y(k,54) +rxt(k,297)*y(k,76)) &
     321             :                  *y(k,50) + (rxt(k,288)*y(k,71) +1.150_r8*rxt(k,289)*y(k,76))*y(k,67) &
     322             :                   +rxt(k,27)*y(k,3) +rxt(k,45)*y(k,26) +rxt(k,115)*y(k,68)*y(k,35) &
     323             :                   +rxt(k,14)*y(k,52) +rxt(k,15)*y(k,53) +rxt(k,17)*y(k,54) +rxt(k,18) &
     324             :                  *y(k,55) +rxt(k,8)*y(k,58) +rxt(k,59)*y(k,59) +rxt(k,89)*y(k,63) &
     325             :                   +rxt(k,90)*y(k,64) +rxt(k,91)*y(k,65) +rxt(k,302)*y(k,78)*y(k,69) &
     326             :                   +rxt(k,96)*y(k,73) +rxt(k,125)*y(k,77)*y(k,77) +rxt(k,305)*y(k,79) &
     327   309041954 :                   +rxt(k,310)*y(k,80) +rxt(k,2)*y(k,81)
     328             :          loss(k,49) = (rxt(k,116)* y(k,35) +rxt(k,135)* y(k,50) +rxt(k,110)* y(k,56) &
     329             :                   +rxt(k,251)* y(k,61) +rxt(k,257)* y(k,63) +rxt(k,130)* y(k,69) &
     330             :                   +rxt(k,291)* y(k,70) + (rxt(k,295) +rxt(k,296))* y(k,72) +rxt(k,97) &
     331             :                  * y(k,73) +rxt(k,102)* y(k,74) +rxt(k,303)* y(k,78) +rxt(k,307) &
     332             :                  * y(k,79) + rxt(k,5) + rxt(k,6) + rxt(k,77) + rxt(k,78) + rxt(k,79) &
     333             :                   + rxt(k,80) + rxt(k,81) + rxt(k,82) + rxt(k,83) + rxt(k,84) &
     334   309041954 :                   + rxt(k,85) + rxt(k,86) + het_rates(k,57))* y(k,57)
     335             :          prod(k,49) = (rxt(k,8) +2.000_r8*rxt(k,99)*y(k,73) + &
     336             :                  2.000_r8*rxt(k,108)*y(k,56) +2.000_r8*rxt(k,118)*y(k,68) + &
     337             :                  rxt(k,119)*y(k,35) +rxt(k,124)*y(k,77) +rxt(k,137)*y(k,54) + &
     338             :                  rxt(k,145)*y(k,53) +rxt(k,161)*y(k,23) +rxt(k,193)*y(k,1) + &
     339             :                  rxt(k,252)*y(k,61) +rxt(k,258)*y(k,63))*y(k,58) &
     340             :                   + (rxt(k,113)*y(k,35) +rxt(k,117)*y(k,56) +rxt(k,122)*y(k,77) + &
     341             :                  rxt(k,127)*y(k,68) +rxt(k,139)*y(k,55) +rxt(k,159)*y(k,23) + &
     342             :                  rxt(k,166)*y(k,26) +rxt(k,192)*y(k,1) +rxt(k,198)*y(k,3) + &
     343             :                  rxt(k,238)*y(k,19))*y(k,68) + (rxt(k,101)*y(k,74) + &
     344             :                  rxt(k,109)*y(k,56) +rxt(k,123)*y(k,77) +rxt(k,136)*y(k,54) + &
     345             :                  rxt(k,141)*y(k,55) +rxt(k,172)*y(k,26) +rxt(k,202)*y(k,3))*y(k,56) &
     346             :                   + (rxt(k,163)*y(k,26) +rxt(k,164)*y(k,26) +rxt(k,174)*y(k,77) + &
     347             :                  rxt(k,196)*y(k,3) +rxt(k,197)*y(k,3))*y(k,26) + (rxt(k,92) + &
     348             :                  rxt(k,100) +2.000_r8*rxt(k,102)*y(k,57))*y(k,74) +rxt(k,194)*y(k,3) &
     349             :                  *y(k,3) +rxt(k,128)*y(k,77)*y(k,47) +rxt(k,134)*y(k,54)*y(k,50) &
     350             :                   +rxt(k,148)*y(k,73)*y(k,51) +rxt(k,299)*y(k,76)*y(k,53) +rxt(k,19) &
     351   309041954 :                  *y(k,55) +rxt(k,93)*y(k,75)
     352             :          loss(k,52) = (rxt(k,193)* y(k,1) +rxt(k,161)* y(k,23) +rxt(k,119)* y(k,35) &
     353             :                   +rxt(k,145)* y(k,53) +rxt(k,137)* y(k,54) +rxt(k,108)* y(k,56) &
     354             :                   +rxt(k,252)* y(k,61) +rxt(k,258)* y(k,63) +rxt(k,118)* y(k,68) &
     355             :                   +rxt(k,99)* y(k,73) +rxt(k,124)* y(k,77) + rxt(k,7) + rxt(k,8) &
     356   309041954 :                   + het_rates(k,58))* y(k,58)
     357   309041954 :          prod(k,52) =rxt(k,110)*y(k,57)*y(k,56)
     358   309041954 :          loss(k,14) = (rxt(k,259)* y(k,63) + rxt(k,59) + het_rates(k,59))* y(k,59)
     359   309041954 :          prod(k,14) = (rxt(k,165)*y(k,26) +rxt(k,195)*y(k,3))*y(k,26)
     360             :          loss(k,15) = (rxt(k,249)* y(k,56) +rxt(k,250)* y(k,77) + rxt(k,88) &
     361   309041954 :                   + het_rates(k,60))* y(k,60)
     362   309041954 :          prod(k,15) = 0._r8
     363             :          loss(k,24) = (rxt(k,251)* y(k,57) +rxt(k,252)* y(k,58) +rxt(k,255)* y(k,77) &
     364   309041954 :                   + het_rates(k,61))* y(k,61)
     365   309041954 :          prod(k,24) =rxt(k,88)*y(k,60) +rxt(k,89)*y(k,63)
     366             :          loss(k,39) = (rxt(k,253)* y(k,3) +rxt(k,254)* y(k,26) +rxt(k,256)* y(k,54) &
     367             :                   +rxt(k,257)* y(k,57) +rxt(k,258)* y(k,58) +rxt(k,259)* y(k,59) &
     368   309041954 :                   +rxt(k,260)* y(k,77) + rxt(k,89) + het_rates(k,63))* y(k,63)
     369             :          prod(k,39) = (rxt(k,251)*y(k,57) +rxt(k,252)*y(k,58) +rxt(k,255)*y(k,77)) &
     370   309041954 :                  *y(k,61) +rxt(k,249)*y(k,60)*y(k,56) +rxt(k,90)*y(k,64)
     371   309041954 :          loss(k,36) = (rxt(k,262)* y(k,77) + rxt(k,90) + het_rates(k,64))* y(k,64)
     372             :          prod(k,36) = (rxt(k,253)*y(k,3) +rxt(k,254)*y(k,26) +rxt(k,256)*y(k,54) + &
     373             :                  rxt(k,257)*y(k,57) +rxt(k,258)*y(k,58) +rxt(k,259)*y(k,59) + &
     374             :                  rxt(k,260)*y(k,77))*y(k,63) + (rxt(k,247)*y(k,55) + &
     375             :                  rxt(k,248)*y(k,77) +.500_r8*rxt(k,261)*y(k,77))*y(k,33) &
     376   309041954 :                   +rxt(k,250)*y(k,77)*y(k,60) +rxt(k,91)*y(k,65)
     377   309041954 :          loss(k,11) = (rxt(k,263)* y(k,81) + rxt(k,91) + het_rates(k,65))* y(k,65)
     378   309041954 :          prod(k,11) =rxt(k,87)*y(k,39) +rxt(k,262)*y(k,77)*y(k,64)
     379   309041954 :          loss(k,1) = ( + het_rates(k,66))* y(k,66)
     380   309041954 :          prod(k,1) = 0._r8
     381             :          loss(k,33) = (rxt(k,290)* y(k,70) +rxt(k,288)* y(k,71) +rxt(k,289)* y(k,76) &
     382   309041954 :                   + het_rates(k,67))* y(k,67)
     383             :          prod(k,33) = (rxt(k,77) +rxt(k,78) +rxt(k,79) +rxt(k,80) +rxt(k,81) + &
     384             :                  rxt(k,84) +rxt(k,85) +rxt(k,86))*y(k,57) + (rxt(k,71) +rxt(k,72) + &
     385             :                  rxt(k,73) +rxt(k,74) +rxt(k,75) +rxt(k,76))*y(k,56) +rxt(k,62) &
     386   309041954 :                  *y(k,50) +rxt(k,16)*y(k,53)
     387             :          loss(k,45) = (rxt(k,192)* y(k,1) +rxt(k,198)* y(k,3) +rxt(k,238)* y(k,19) &
     388             :                   + (rxt(k,159) +rxt(k,160))* y(k,23) +rxt(k,166)* y(k,26) &
     389             :                   + (rxt(k,113) +rxt(k,114) +rxt(k,115))* y(k,35) +rxt(k,144)* y(k,53) &
     390             :                   +rxt(k,149)* y(k,54) +rxt(k,139)* y(k,55) +rxt(k,117)* y(k,56) &
     391             :                   +rxt(k,118)* y(k,58) + 2._r8*rxt(k,127)* y(k,68) +rxt(k,122) &
     392   309041954 :                  * y(k,77) + rxt(k,264) + het_rates(k,68))* y(k,68)
     393             :          prod(k,45) = (rxt(k,221)*y(k,16) +rxt(k,224)*y(k,18) +rxt(k,121)*y(k,38) + &
     394             :                  rxt(k,124)*y(k,58) +rxt(k,142)*y(k,55) +rxt(k,173)*y(k,26) + &
     395             :                  rxt(k,203)*y(k,3) +rxt(k,242)*y(k,29) +.500_r8*rxt(k,261)*y(k,33) + &
     396             :                  rxt(k,262)*y(k,64))*y(k,77) + (rxt(k,155)*y(k,23) + &
     397             :                  rxt(k,191)*y(k,1) +rxt(k,235)*y(k,55) +rxt(k,236)*y(k,56))*y(k,15) &
     398             :                   + (rxt(k,220)*y(k,16) +rxt(k,223)*y(k,18) +rxt(k,158)*y(k,38)) &
     399             :                  *y(k,23) + (rxt(k,162)*y(k,26) +rxt(k,239)*y(k,53))*y(k,19) &
     400             :                   + (rxt(k,11) +rxt(k,153))*y(k,47) +rxt(k,244)*y(k,73)*y(k,21) &
     401   309041954 :                   +rxt(k,116)*y(k,57)*y(k,35) +rxt(k,112)*y(k,56)*y(k,38)
     402             :          loss(k,34) = (rxt(k,129)* y(k,56) +rxt(k,130)* y(k,57) +rxt(k,302)* y(k,78) &
     403   309041954 :                   + het_rates(k,69))* y(k,69)
     404             :          prod(k,34) = (.800_r8*rxt(k,288)*y(k,71) +.900_r8*rxt(k,290)*y(k,70))*y(k,67) &
     405   309041954 :                   +rxt(k,292)*y(k,70)*y(k,56)
     406             :          loss(k,23) = ((rxt(k,292) +rxt(k,293))* y(k,56) +rxt(k,291)* y(k,57) &
     407   309041954 :                   +rxt(k,290)* y(k,67) + het_rates(k,70))* y(k,70)
     408   309041954 :          prod(k,23) =rxt(k,305)*y(k,79) +rxt(k,310)*y(k,80)
     409   309041954 :          loss(k,25) = (rxt(k,288)* y(k,67) + het_rates(k,71))* y(k,71)
     410             :          prod(k,25) = (rxt(k,298) +rxt(k,297)*y(k,50) +rxt(k,299)*y(k,53))*y(k,76) &
     411             :                   +rxt(k,16)*y(k,53) +rxt(k,292)*y(k,70)*y(k,56) +rxt(k,296)*y(k,72) &
     412   309041954 :                  *y(k,57) +rxt(k,301)*y(k,78)
     413             :          loss(k,30) = (rxt(k,294)* y(k,56) + (rxt(k,295) +rxt(k,296))* y(k,57) &
     414   309041954 :                   + het_rates(k,72))* y(k,72)
     415   309041954 :          prod(k,30) =rxt(k,62)*y(k,50) +rxt(k,302)*y(k,78)*y(k,69) +rxt(k,311)*y(k,80)
     416             :          loss(k,50) = (rxt(k,180)* y(k,6) +rxt(k,181)* y(k,7) +rxt(k,207)* y(k,8) &
     417             :                   +rxt(k,182)* y(k,9) +rxt(k,183)* y(k,10) +rxt(k,184)* y(k,11) &
     418             :                   +rxt(k,185)* y(k,12) +rxt(k,186)* y(k,13) +rxt(k,230)* y(k,14) &
     419             :                   +rxt(k,231)* y(k,16) + (rxt(k,243) +rxt(k,244) +rxt(k,245))* y(k,21) &
     420             :                   +rxt(k,208)* y(k,22) +rxt(k,216)* y(k,31) +rxt(k,217)* y(k,32) &
     421             :                   +rxt(k,94)* y(k,36) +rxt(k,209)* y(k,37) + (rxt(k,210) +rxt(k,211)) &
     422             :                  * y(k,40) +rxt(k,232)* y(k,41) +rxt(k,233)* y(k,42) +rxt(k,234) &
     423             :                  * y(k,43) + (rxt(k,187) +rxt(k,188))* y(k,44) + (rxt(k,147) + &
     424             :                  rxt(k,148))* y(k,51) + (rxt(k,97) +rxt(k,98))* y(k,57) +rxt(k,99) &
     425   309041954 :                  * y(k,58) +rxt(k,95)* y(k,81) + rxt(k,96) + het_rates(k,73))* y(k,73)
     426             :          prod(k,50) = (rxt(k,6) +rxt(k,130)*y(k,69))*y(k,57) +rxt(k,7)*y(k,58) &
     427   309041954 :                   +.850_r8*rxt(k,289)*y(k,76)*y(k,67) +rxt(k,1)*y(k,81)
     428             :          loss(k,5) = (rxt(k,101)* y(k,56) +rxt(k,102)* y(k,57) + rxt(k,92) &
     429   309041954 :                   + rxt(k,100) + het_rates(k,74))* y(k,74)
     430             :          prod(k,5) = (rxt(k,104) +rxt(k,103)*y(k,30) +rxt(k,105)*y(k,56) + &
     431   309041954 :                  rxt(k,106)*y(k,57) +rxt(k,107)*y(k,58))*y(k,75) +rxt(k,7)*y(k,58)
     432             :          loss(k,6) = (rxt(k,103)* y(k,30) +rxt(k,105)* y(k,56) +rxt(k,106)* y(k,57) &
     433             :                   +rxt(k,107)* y(k,58) + rxt(k,93) + rxt(k,104) + het_rates(k,75)) &
     434   309041954 :                  * y(k,75)
     435   309041954 :          prod(k,6) =rxt(k,97)*y(k,73)*y(k,57)
     436             :          loss(k,32) = (rxt(k,297)* y(k,50) +rxt(k,299)* y(k,53) +rxt(k,289)* y(k,67) &
     437   309041954 :                   + rxt(k,298) + het_rates(k,76))* y(k,76)
     438             :          prod(k,32) = (rxt(k,78) +rxt(k,80) +rxt(k,291)*y(k,70) +rxt(k,295)*y(k,72) + &
     439             :                  rxt(k,303)*y(k,78) +rxt(k,307)*y(k,79))*y(k,57) +rxt(k,300)*y(k,78) &
     440   309041954 :                  *y(k,30)
     441             :          loss(k,54) = (rxt(k,203)* y(k,3) +rxt(k,219)* y(k,14) +rxt(k,237)* y(k,15) &
     442             :                   +rxt(k,221)* y(k,16) +rxt(k,222)* y(k,17) +rxt(k,224)* y(k,18) &
     443             :                   +rxt(k,240)* y(k,20) +rxt(k,241)* y(k,21) +rxt(k,226)* y(k,22) &
     444             :                   + (rxt(k,173) +rxt(k,174))* y(k,26) +rxt(k,171)* y(k,27) &
     445             :                   + (rxt(k,242) +rxt(k,246))* y(k,29) + (rxt(k,248) +rxt(k,261)) &
     446             :                  * y(k,33) +rxt(k,120)* y(k,36) +rxt(k,121)* y(k,38) +rxt(k,205) &
     447             :                  * y(k,40) +rxt(k,227)* y(k,41) +rxt(k,228)* y(k,42) +rxt(k,229) &
     448             :                  * y(k,43) +rxt(k,176)* y(k,44) +rxt(k,152)* y(k,46) +rxt(k,128) &
     449             :                  * y(k,47) +rxt(k,179)* y(k,49) +rxt(k,143)* y(k,50) +rxt(k,151) &
     450             :                  * y(k,54) +rxt(k,142)* y(k,55) +rxt(k,123)* y(k,56) +rxt(k,124) &
     451             :                  * y(k,58) +rxt(k,250)* y(k,60) +rxt(k,255)* y(k,61) +rxt(k,260) &
     452             :                  * y(k,63) +rxt(k,262)* y(k,64) +rxt(k,122)* y(k,68) &
     453             :                   + 2._r8*(rxt(k,125) +rxt(k,126))* y(k,77) + het_rates(k,77)) &
     454   309041954 :                  * y(k,77)
     455             :          prod(k,54) = (rxt(k,111)*y(k,36) +rxt(k,112)*y(k,38) +rxt(k,117)*y(k,68) + &
     456             :                  rxt(k,175)*y(k,44) +rxt(k,178)*y(k,49) +rxt(k,204)*y(k,40) + &
     457             :                  rxt(k,206)*y(k,48) +rxt(k,236)*y(k,15))*y(k,56) &
     458             :                   + (2.000_r8*rxt(k,114)*y(k,35) +rxt(k,118)*y(k,58) + &
     459             :                  rxt(k,139)*y(k,55) +rxt(k,144)*y(k,53) +rxt(k,160)*y(k,23))*y(k,68) &
     460             :                   + (rxt(k,243)*y(k,21) +rxt(k,94)*y(k,36) + &
     461             :                  2.000_r8*rxt(k,95)*y(k,81) +rxt(k,187)*y(k,44) +rxt(k,210)*y(k,40)) &
     462             :                  *y(k,73) + (rxt(k,22) +.300_r8*rxt(k,240)*y(k,77))*y(k,20) &
     463             :                   + (rxt(k,3) +rxt(k,214)*y(k,34))*y(k,81) +rxt(k,119)*y(k,58)*y(k,35) &
     464             :                   +2.000_r8*rxt(k,4)*y(k,38) +rxt(k,9)*y(k,46) +rxt(k,10)*y(k,47) &
     465   309041954 :                   +rxt(k,57)*y(k,48) +rxt(k,58)*y(k,49) +.500_r8*rxt(k,266)*y(k,54)
     466             :          loss(k,31) = (rxt(k,300)* y(k,30) +rxt(k,303)* y(k,57) +rxt(k,302)* y(k,69) &
     467   309041954 :                   + rxt(k,301) + het_rates(k,78))* y(k,78)
     468             :          prod(k,31) = (rxt(k,73) +rxt(k,74) +rxt(k,293)*y(k,70) +rxt(k,294)*y(k,72) + &
     469             :                  rxt(k,306)*y(k,79) +rxt(k,312)*y(k,80))*y(k,56) + (rxt(k,79) + &
     470             :                  rxt(k,81))*y(k,57) + (rxt(k,304)*y(k,79) +rxt(k,309)*y(k,80))*y(k,67) &
     471   309041954 :                   +rxt(k,286)*y(k,79) +rxt(k,285)*y(k,80)
     472             :          loss(k,19) = (rxt(k,306)* y(k,56) +rxt(k,307)* y(k,57) +rxt(k,304)* y(k,67) &
     473   309041954 :                   + rxt(k,286) + rxt(k,305) + het_rates(k,79))* y(k,79)
     474             :          prod(k,19) = (rxt(k,75) +rxt(k,76))*y(k,56) + (rxt(k,85) +rxt(k,86))*y(k,57) &
     475   309041954 :                   + (rxt(k,287) +rxt(k,308)*y(k,67))*y(k,80)
     476             :          loss(k,18) = (rxt(k,312)* y(k,56) + (rxt(k,308) +rxt(k,309))* y(k,67) &
     477             :                   + rxt(k,285) + rxt(k,287) + rxt(k,310) + rxt(k,311) &
     478   309041954 :                   + het_rates(k,80))* y(k,80)
     479   309041954 :          prod(k,18) = (rxt(k,71) +rxt(k,72))*y(k,56) + (rxt(k,77) +rxt(k,84))*y(k,57)
     480             :          loss(k,58) = (rxt(k,214)* y(k,34) +rxt(k,263)* y(k,65) +rxt(k,95)* y(k,73) &
     481   309041954 :                   + rxt(k,1) + rxt(k,2) + rxt(k,3) + het_rates(k,81))* y(k,81)
     482             :          prod(k,58) = (rxt(k,219)*y(k,14) +rxt(k,221)*y(k,16) +rxt(k,222)*y(k,17) + &
     483             :                  rxt(k,224)*y(k,18) +rxt(k,229)*y(k,43) +rxt(k,241)*y(k,21) + &
     484             :                  rxt(k,120)*y(k,36) +rxt(k,121)*y(k,38) +rxt(k,122)*y(k,68) + &
     485             :                  rxt(k,125)*y(k,77) +rxt(k,128)*y(k,47) +rxt(k,152)*y(k,46) + &
     486             :                  rxt(k,176)*y(k,44) +rxt(k,179)*y(k,49) +rxt(k,205)*y(k,40) + &
     487             :                  rxt(k,237)*y(k,15) +rxt(k,240)*y(k,20))*y(k,77) &
     488             :                   + (rxt(k,269)*y(k,49) +rxt(k,275)*y(k,49) +rxt(k,276)*y(k,48) + &
     489             :                  rxt(k,280)*y(k,49) +rxt(k,281)*y(k,48))*y(k,44) + (rxt(k,264) + &
     490   318792954 :                  rxt(k,115)*y(k,35))*y(k,68) +rxt(k,87)*y(k,39)
     491             :       end do
     492     9751000 :       end subroutine imp_prod_loss
     493             :       end module mo_prod_loss

Generated by: LCOV version 1.14