LCOV - code coverage report
Current view: top level - chemistry/pp_waccm_mad_mam5 - mo_prod_loss.F90 (source / functions) Hit Total Coverage
Test: coverage.info Lines: 283 283 100.0 %
Date: 2025-03-14 01:26:08 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       21888 :       subroutine exp_prod_loss( ofl, ofu, prod, loss, y, &
       9             :                                 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       21888 :       end subroutine exp_prod_loss
      23     3585597 :       subroutine imp_prod_loss( avec_len, prod, loss, y, &
      24             :                                 rxt, het_rates )
      25             :       use chem_mods, only : gas_pcnst,rxntot,clscnt4
      26             :       implicit none
      27             : !--------------------------------------------------------------------
      28             : ! ... dummy args
      29             : !--------------------------------------------------------------------
      30             :       integer, intent(in) :: avec_len
      31             :       real(r8), dimension(veclen,clscnt4), intent(out) :: &
      32             :             prod, &
      33             :             loss
      34             :       real(r8), intent(in) :: y(veclen,gas_pcnst)
      35             :       real(r8), intent(in) :: rxt(veclen,rxntot)
      36             :       real(r8), intent(in) :: het_rates(veclen,gas_pcnst)
      37             : !--------------------------------------------------------------------
      38             : ! ... local variables
      39             : !--------------------------------------------------------------------
      40             :       integer :: k
      41             : !--------------------------------------------------------------------
      42             : ! ... loss and production for Implicit method
      43             : !--------------------------------------------------------------------
      44   117656469 :       do k = 1,avec_len
      45   114070872 :          loss(k,1) = ( + het_rates(k,1))* y(k,1)
      46   114070872 :          prod(k,1) = 0._r8
      47   114070872 :          loss(k,2) = ( + het_rates(k,2))* y(k,2)
      48   114070872 :          prod(k,2) = 0._r8
      49             :          loss(k,99) = (rxt(k,486)* y(k,17) +rxt(k,488)* y(k,73) +rxt(k,487)* y(k,101) &
      50   114070872 :                   + het_rates(k,3))* y(k,3)
      51             :          prod(k,99) = (rxt(k,28) +2.000_r8*rxt(k,489)*y(k,5) +rxt(k,490)*y(k,28) + &
      52             :                  rxt(k,491)*y(k,28) +rxt(k,494)*y(k,62) +rxt(k,497)*y(k,70) + &
      53             :                  rxt(k,498)*y(k,133) +rxt(k,548)*y(k,80))*y(k,5) &
      54             :                   + (rxt(k,476)*y(k,9) +rxt(k,502)*y(k,10) + &
      55             :                  3.000_r8*rxt(k,503)*y(k,24) +2.000_r8*rxt(k,504)*y(k,42) + &
      56             :                  rxt(k,505)*y(k,45) +2.000_r8*rxt(k,525)*y(k,16) +rxt(k,526)*y(k,18)) &
      57             :                  *y(k,71) + (rxt(k,500)*y(k,45) +2.000_r8*rxt(k,514)*y(k,16) + &
      58             :                  rxt(k,516)*y(k,18) +3.000_r8*rxt(k,521)*y(k,24))*y(k,133) &
      59             :                   + (2.000_r8*rxt(k,513)*y(k,16) +rxt(k,515)*y(k,18) + &
      60             :                  3.000_r8*rxt(k,520)*y(k,24))*y(k,25) + (rxt(k,52) + &
      61             :                  rxt(k,499)*y(k,70))*y(k,45) +rxt(k,27)*y(k,4) +rxt(k,30)*y(k,6) &
      62             :                   +rxt(k,32)*y(k,9) +rxt(k,33)*y(k,10) +2.000_r8*rxt(k,39)*y(k,16) &
      63             :                   +rxt(k,40)*y(k,18) +3.000_r8*rxt(k,43)*y(k,24) +2.000_r8*rxt(k,51) &
      64   114070872 :                  *y(k,42) +rxt(k,58)*y(k,53)
      65   114070872 :          loss(k,43) = ( + rxt(k,27) + het_rates(k,4))* y(k,4)
      66             :          prod(k,43) = (rxt(k,570)*y(k,53) +rxt(k,575)*y(k,53))*y(k,49) &
      67   114070872 :                   +rxt(k,492)*y(k,28)*y(k,5)
      68             :          loss(k,106) = (2._r8*rxt(k,489)* y(k,5) + (rxt(k,490) +rxt(k,491) + &
      69             :                  rxt(k,492))* y(k,28) +rxt(k,494)* y(k,62) +rxt(k,495)* y(k,63) &
      70             :                   +rxt(k,497)* y(k,70) +rxt(k,548)* y(k,80) +rxt(k,493)* y(k,101) &
      71   114070872 :                   +rxt(k,498)* y(k,133) + rxt(k,28) + het_rates(k,5))* y(k,5)
      72             :          prod(k,106) = (rxt(k,29) +rxt(k,496)*y(k,70))*y(k,6) +rxt(k,488)*y(k,73) &
      73   114070872 :                  *y(k,3) +rxt(k,506)*y(k,71)*y(k,45) +rxt(k,501)*y(k,70)*y(k,53)
      74             :          loss(k,63) = (rxt(k,496)* y(k,70) + rxt(k,29) + rxt(k,30) + rxt(k,564) &
      75   114070872 :                   + rxt(k,567) + rxt(k,572) + het_rates(k,6))* y(k,6)
      76   114070872 :          prod(k,63) =rxt(k,495)*y(k,63)*y(k,5)
      77   114070872 :          loss(k,3) = ( + het_rates(k,7))* y(k,7)
      78   114070872 :          prod(k,3) = 0._r8
      79   114070872 :          loss(k,27) = (rxt(k,475)* y(k,71) + rxt(k,31) + het_rates(k,8))* y(k,8)
      80   114070872 :          prod(k,27) = 0._r8
      81   114070872 :          loss(k,34) = (rxt(k,476)* y(k,71) + rxt(k,32) + het_rates(k,9))* y(k,9)
      82   114070872 :          prod(k,34) = 0._r8
      83   114070872 :          loss(k,35) = (rxt(k,502)* y(k,71) + rxt(k,33) + het_rates(k,10))* y(k,10)
      84   114070872 :          prod(k,35) = 0._r8
      85   114070872 :          loss(k,29) = (rxt(k,477)* y(k,71) + rxt(k,34) + het_rates(k,11))* y(k,11)
      86   114070872 :          prod(k,29) = 0._r8
      87   114070872 :          loss(k,36) = (rxt(k,478)* y(k,71) + rxt(k,35) + het_rates(k,12))* y(k,12)
      88   114070872 :          prod(k,36) = 0._r8
      89   114070872 :          loss(k,30) = (rxt(k,479)* y(k,71) + rxt(k,36) + het_rates(k,13))* y(k,13)
      90   114070872 :          prod(k,30) = 0._r8
      91   114070872 :          loss(k,37) = (rxt(k,480)* y(k,71) + rxt(k,37) + het_rates(k,14))* y(k,14)
      92   114070872 :          prod(k,37) = 0._r8
      93   114070872 :          loss(k,31) = (rxt(k,481)* y(k,71) + rxt(k,38) + het_rates(k,15))* y(k,15)
      94   114070872 :          prod(k,31) = 0._r8
      95             :          loss(k,60) = (rxt(k,513)* y(k,25) +rxt(k,525)* y(k,71) +rxt(k,514)* y(k,133) &
      96   114070872 :                   + rxt(k,39) + het_rates(k,16))* y(k,16)
      97   114070872 :          prod(k,60) = 0._r8
      98             :          loss(k,107) = (rxt(k,486)* y(k,3) +rxt(k,450)* y(k,25) +rxt(k,530)* y(k,64) &
      99             :                   +rxt(k,531)* y(k,70) +rxt(k,532)* y(k,133) + rxt(k,21) + rxt(k,22) &
     100   114070872 :                   + het_rates(k,17))* y(k,17)
     101             :          prod(k,107) = (.180_r8*rxt(k,25) +rxt(k,538)*y(k,71) +rxt(k,539)*y(k,71)) &
     102             :                  *y(k,23) + (rxt(k,457)*y(k,28) +rxt(k,534)*y(k,62))*y(k,21) &
     103   114070872 :                   + (rxt(k,23) +.300_r8*rxt(k,535)*y(k,133))*y(k,22)
     104             :          loss(k,65) = (rxt(k,515)* y(k,25) +rxt(k,526)* y(k,71) +rxt(k,516)* y(k,133) &
     105   114070872 :                   + rxt(k,40) + het_rates(k,18))* y(k,18)
     106   114070872 :          prod(k,65) = 0._r8
     107   114070872 :          loss(k,32) = (rxt(k,517)* y(k,133) + rxt(k,41) + het_rates(k,19))* y(k,19)
     108   114070872 :          prod(k,32) = 0._r8
     109             :          loss(k,54) = (rxt(k,518)* y(k,25) +rxt(k,519)* y(k,133) + rxt(k,42) &
     110   114070872 :                   + het_rates(k,20))* y(k,20)
     111   114070872 :          prod(k,54) = 0._r8
     112             :          loss(k,89) = (rxt(k,457)* y(k,28) +rxt(k,534)* y(k,62) +rxt(k,533)* y(k,101) &
     113   114070872 :                   + het_rates(k,21))* y(k,21)
     114             :          prod(k,89) = (rxt(k,24) +rxt(k,451)*y(k,25) +rxt(k,507)*y(k,39) + &
     115             :                  rxt(k,536)*y(k,133) +rxt(k,537)*y(k,71))*y(k,23) +rxt(k,40)*y(k,18) &
     116   114070872 :                   +rxt(k,42)*y(k,20) +.700_r8*rxt(k,535)*y(k,133)*y(k,22)
     117   114070872 :          loss(k,55) = (rxt(k,535)* y(k,133) + rxt(k,23) + het_rates(k,22))* y(k,22)
     118   114070872 :          prod(k,55) =rxt(k,533)*y(k,101)*y(k,21)
     119             :          loss(k,88) = (rxt(k,451)* y(k,25) +rxt(k,507)* y(k,39) + (rxt(k,537) + &
     120             :                  rxt(k,538) +rxt(k,539))* y(k,71) +rxt(k,536)* y(k,133) + rxt(k,24) &
     121   114070872 :                   + rxt(k,25) + het_rates(k,23))* y(k,23)
     122   114070872 :          prod(k,88) = 0._r8
     123             :          loss(k,56) = (rxt(k,520)* y(k,25) +rxt(k,503)* y(k,71) +rxt(k,521)* y(k,133) &
     124   114070872 :                   + rxt(k,43) + het_rates(k,24))* y(k,24)
     125   114070872 :          prod(k,56) = 0._r8
     126             :          loss(k,125) = (rxt(k,513)* y(k,16) +rxt(k,450)* y(k,17) +rxt(k,515)* y(k,18) &
     127             :                   +rxt(k,518)* y(k,20) +rxt(k,451)* y(k,23) +rxt(k,520)* y(k,24) &
     128             :                   +rxt(k,463)* y(k,29) +rxt(k,452)* y(k,41) +rxt(k,453)* y(k,43) &
     129             :                   +rxt(k,472)* y(k,54) +rxt(k,456)* y(k,73) + (rxt(k,114) +rxt(k,115)) &
     130             :                  * y(k,94) +rxt(k,127)* y(k,97) + (rxt(k,454) +rxt(k,455))* y(k,101) &
     131             :                   +rxt(k,285)* y(k,111) +rxt(k,314)* y(k,127) +rxt(k,341)* y(k,134) &
     132   114070872 :                   +rxt(k,350)* y(k,135) + het_rates(k,25))* y(k,25)
     133             :          prod(k,125) = (4.000_r8*rxt(k,475)*y(k,8) +rxt(k,476)*y(k,9) + &
     134             :                  2.000_r8*rxt(k,477)*y(k,11) +2.000_r8*rxt(k,478)*y(k,12) + &
     135             :                  2.000_r8*rxt(k,479)*y(k,13) +rxt(k,480)*y(k,14) + &
     136             :                  2.000_r8*rxt(k,481)*y(k,15) +rxt(k,482)*y(k,49) +rxt(k,512)*y(k,34) + &
     137             :                  rxt(k,527)*y(k,46) +rxt(k,528)*y(k,47) +rxt(k,529)*y(k,48))*y(k,71) &
     138             :                   + (rxt(k,46) +rxt(k,457)*y(k,21) +2.000_r8*rxt(k,458)*y(k,28) + &
     139             :                  rxt(k,460)*y(k,28) +rxt(k,462)*y(k,62) +rxt(k,467)*y(k,70) + &
     140             :                  rxt(k,468)*y(k,133) +rxt(k,491)*y(k,5) +rxt(k,549)*y(k,80))*y(k,28) &
     141             :                   + (rxt(k,110)*y(k,63) +rxt(k,147)*y(k,105) +rxt(k,154)*y(k,128) + &
     142             :                  rxt(k,177)*y(k,107) +rxt(k,195)*y(k,103) +rxt(k,212)*y(k,122) + &
     143             :                  rxt(k,230)*y(k,119) +rxt(k,248)*y(k,118))*y(k,90) &
     144             :                   + (rxt(k,159)*y(k,128) +rxt(k,164)*y(k,105) +rxt(k,182)*y(k,107) + &
     145             :                  rxt(k,199)*y(k,103) +rxt(k,217)*y(k,122) +rxt(k,234)*y(k,119) + &
     146             :                  rxt(k,252)*y(k,118))*y(k,91) + (rxt(k,169)*y(k,107) + &
     147             :                  rxt(k,181)*y(k,105) +rxt(k,187)*y(k,103) +rxt(k,205)*y(k,122) + &
     148             :                  rxt(k,222)*y(k,119) +rxt(k,240)*y(k,118) +rxt(k,257)*y(k,128)) &
     149             :                  *y(k,92) + (rxt(k,471)*y(k,49) +3.000_r8*rxt(k,517)*y(k,19) + &
     150             :                  rxt(k,519)*y(k,20) +rxt(k,522)*y(k,46) +rxt(k,523)*y(k,47) + &
     151             :                  rxt(k,524)*y(k,48))*y(k,133) + (rxt(k,56) +rxt(k,470)*y(k,70)) &
     152             :                  *y(k,49) +rxt(k,27)*y(k,4) +4.000_r8*rxt(k,31)*y(k,8) +rxt(k,32) &
     153             :                  *y(k,9) +2.000_r8*rxt(k,34)*y(k,11) +2.000_r8*rxt(k,35)*y(k,12) &
     154             :                   +2.000_r8*rxt(k,36)*y(k,13) +rxt(k,37)*y(k,14) +2.000_r8*rxt(k,38) &
     155             :                  *y(k,15) +3.000_r8*rxt(k,41)*y(k,19) +rxt(k,42)*y(k,20) &
     156             :                   +2.000_r8*rxt(k,44)*y(k,26) +2.000_r8*rxt(k,45)*y(k,27) +rxt(k,47) &
     157             :                  *y(k,29) +rxt(k,50)*y(k,34) +rxt(k,53)*y(k,46) +rxt(k,54)*y(k,47) &
     158   114070872 :                   +rxt(k,55)*y(k,48) +rxt(k,59)*y(k,54) +rxt(k,111)*y(k,93)*y(k,62)
     159   114070872 :          loss(k,38) = ( + rxt(k,44) + het_rates(k,26))* y(k,26)
     160             :          prod(k,38) = (rxt(k,563)*y(k,54) +rxt(k,568)*y(k,29) +rxt(k,569)*y(k,54) + &
     161             :                  rxt(k,573)*y(k,29) +rxt(k,574)*y(k,54) +rxt(k,578)*y(k,29))*y(k,49) &
     162   114070872 :                   +rxt(k,463)*y(k,29)*y(k,25) +rxt(k,459)*y(k,28)*y(k,28)
     163   114070872 :          loss(k,26) = ( + rxt(k,45) + rxt(k,485) + het_rates(k,27))* y(k,27)
     164   114070872 :          prod(k,26) =rxt(k,484)*y(k,28)*y(k,28)
     165             :          loss(k,117) = ((rxt(k,490) +rxt(k,491) +rxt(k,492))* y(k,5) +rxt(k,457) &
     166             :                  * y(k,21) + 2._r8*(rxt(k,458) +rxt(k,459) +rxt(k,460) +rxt(k,484)) &
     167             :                  * y(k,28) +rxt(k,462)* y(k,62) +rxt(k,464)* y(k,63) +rxt(k,467) &
     168             :                  * y(k,70) +rxt(k,549)* y(k,80) +rxt(k,116)* y(k,94) +rxt(k,128) &
     169             :                  * y(k,97) +rxt(k,461)* y(k,101) +rxt(k,286)* y(k,111) +rxt(k,315) &
     170             :                  * y(k,127) + (rxt(k,468) +rxt(k,469))* y(k,133) +rxt(k,342)* y(k,134) &
     171   114070872 :                   +rxt(k,351)* y(k,135) + rxt(k,46) + het_rates(k,28))* y(k,28)
     172             :          prod(k,117) = (rxt(k,455)*y(k,101) +rxt(k,456)*y(k,73) +rxt(k,472)*y(k,54)) &
     173             :                  *y(k,25) + (rxt(k,48) +rxt(k,465)*y(k,70))*y(k,29) &
     174             :                   + (rxt(k,473)*y(k,70) +rxt(k,474)*y(k,133))*y(k,54) + (rxt(k,60) + &
     175             :                  rxt(k,554)*y(k,80))*y(k,74) +2.000_r8*rxt(k,485)*y(k,27) &
     176   114070872 :                   +rxt(k,483)*y(k,71)*y(k,49)
     177             :          loss(k,83) = (rxt(k,463)* y(k,25) + (rxt(k,568) +rxt(k,573) +rxt(k,578)) &
     178             :                  * y(k,49) +rxt(k,465)* y(k,70) +rxt(k,466)* y(k,133) + rxt(k,47) &
     179             :                   + rxt(k,48) + rxt(k,566) + rxt(k,571) + rxt(k,577) &
     180   114070872 :                   + het_rates(k,29))* y(k,29)
     181   114070872 :          prod(k,83) =rxt(k,464)*y(k,63)*y(k,28)
     182   114070872 :          loss(k,4) = ( + het_rates(k,30))* y(k,30)
     183   114070872 :          prod(k,4) = 0._r8
     184   114070872 :          loss(k,69) = (rxt(k,540)* y(k,133) + het_rates(k,31))* y(k,31)
     185             :          prod(k,69) = (rxt(k,21) +rxt(k,22) +rxt(k,450)*y(k,25) +rxt(k,486)*y(k,3) + &
     186             :                  rxt(k,530)*y(k,64) +rxt(k,531)*y(k,70) +rxt(k,532)*y(k,133))*y(k,17) &
     187             :                   + (rxt(k,26) +rxt(k,62) +rxt(k,594)*y(k,136))*y(k,32) &
     188             :                   + (rxt(k,101) +rxt(k,543)*y(k,70) +rxt(k,544)*y(k,133))*y(k,75) &
     189   114070872 :                   +rxt(k,518)*y(k,25)*y(k,20) +.380_r8*rxt(k,25)*y(k,23)
     190             :          loss(k,121) = (rxt(k,368)* y(k,118) +rxt(k,312)* y(k,123) +rxt(k,316) &
     191             :                  * y(k,127) +rxt(k,330)* y(k,130) +rxt(k,335)* y(k,131) +rxt(k,343) &
     192             :                  * y(k,134) +rxt(k,352)* y(k,135) +rxt(k,594)* y(k,136) + rxt(k,26) &
     193   114070872 :                   + rxt(k,62) + het_rates(k,32))* y(k,32)
     194             :          prod(k,121) = (rxt(k,63) +rxt(k,114)*y(k,25) +rxt(k,115)*y(k,25) + &
     195             :                  rxt(k,116)*y(k,28) +rxt(k,117)*y(k,40) +rxt(k,124)*y(k,51) + &
     196             :                  rxt(k,125)*y(k,70) +rxt(k,126)*y(k,72) +rxt(k,168)*y(k,107) + &
     197             :                  rxt(k,170)*y(k,105) +rxt(k,186)*y(k,103) +rxt(k,204)*y(k,122) + &
     198             :                  rxt(k,221)*y(k,119) +rxt(k,239)*y(k,118) +rxt(k,256)*y(k,128) + &
     199             :                  rxt(k,258)*y(k,105) +rxt(k,265)*y(k,107) +rxt(k,280)*y(k,62) + &
     200             :                  rxt(k,281)*y(k,63))*y(k,94) + (rxt(k,120)*y(k,63) + &
     201             :                  rxt(k,121)*y(k,63) +rxt(k,122)*y(k,62) +rxt(k,123)*y(k,62) + &
     202             :                  rxt(k,155)*y(k,128) +rxt(k,158)*y(k,105) +rxt(k,178)*y(k,107) + &
     203             :                  rxt(k,196)*y(k,103) +rxt(k,213)*y(k,122) +rxt(k,231)*y(k,119) + &
     204             :                  rxt(k,249)*y(k,118) +rxt(k,260)*y(k,105) +rxt(k,261)*y(k,107)) &
     205             :                  *y(k,96) + (rxt(k,65) +rxt(k,127)*y(k,25) +rxt(k,128)*y(k,28) + &
     206             :                  rxt(k,130)*y(k,49) +rxt(k,132)*y(k,73) +rxt(k,151)*y(k,128) + &
     207             :                  rxt(k,174)*y(k,107) +rxt(k,191)*y(k,103) +rxt(k,209)*y(k,122) + &
     208             :                  rxt(k,225)*y(k,105) +rxt(k,227)*y(k,119) +rxt(k,244)*y(k,118)) &
     209             :                  *y(k,97) + (rxt(k,153)*y(k,128) +rxt(k,176)*y(k,107) + &
     210             :                  rxt(k,194)*y(k,103) +rxt(k,211)*y(k,122) +rxt(k,229)*y(k,119) + &
     211             :                  rxt(k,246)*y(k,118) +rxt(k,247)*y(k,105) +rxt(k,259)*y(k,107) + &
     212             :                  rxt(k,271)*y(k,105))*y(k,95) + (rxt(k,149)*y(k,128) + &
     213             :                  rxt(k,172)*y(k,107) +rxt(k,189)*y(k,103) +rxt(k,203)*y(k,105) + &
     214             :                  rxt(k,207)*y(k,122) +rxt(k,224)*y(k,119) +rxt(k,242)*y(k,118)) &
     215             :                  *y(k,100) + (rxt(k,369) +rxt(k,306)*y(k,98) +rxt(k,307)*y(k,139)) &
     216   114070872 :                  *y(k,121) +.440_r8*rxt(k,25)*y(k,23) +rxt(k,540)*y(k,133)*y(k,31)
     217   114070872 :          loss(k,39) = (rxt(k,511)* y(k,71) + rxt(k,49) + het_rates(k,33))* y(k,33)
     218             :          prod(k,39) = (rxt(k,476)*y(k,9) +rxt(k,478)*y(k,12) + &
     219             :                  2.000_r8*rxt(k,479)*y(k,13) +2.000_r8*rxt(k,480)*y(k,14) + &
     220             :                  rxt(k,481)*y(k,15) +rxt(k,502)*y(k,10) +2.000_r8*rxt(k,504)*y(k,42) + &
     221             :                  rxt(k,528)*y(k,47) +rxt(k,529)*y(k,48))*y(k,71) + (rxt(k,54) + &
     222             :                  rxt(k,523)*y(k,133))*y(k,47) + (rxt(k,55) +rxt(k,524)*y(k,133)) &
     223             :                  *y(k,48) +rxt(k,32)*y(k,9) +rxt(k,33)*y(k,10) +rxt(k,35)*y(k,12) &
     224             :                   +2.000_r8*rxt(k,36)*y(k,13) +2.000_r8*rxt(k,37)*y(k,14) +rxt(k,38) &
     225   114070872 :                  *y(k,15) +2.000_r8*rxt(k,51)*y(k,42)
     226   114070872 :          loss(k,41) = (rxt(k,512)* y(k,71) + rxt(k,50) + het_rates(k,34))* y(k,34)
     227             :          prod(k,41) = (rxt(k,53) +rxt(k,522)*y(k,133) +rxt(k,527)*y(k,71))*y(k,46) &
     228             :                   + (rxt(k,34) +rxt(k,477)*y(k,71))*y(k,11) + (rxt(k,35) + &
     229   114070872 :                  rxt(k,478)*y(k,71))*y(k,12)
     230             :          loss(k,50) = (rxt(k,541)* y(k,64) + (rxt(k,542) +rxt(k,556))* y(k,133) &
     231   114070872 :                   + het_rates(k,35))* y(k,35)
     232   114070872 :          prod(k,50) = 0._r8
     233   114070872 :          loss(k,5) = ( + het_rates(k,36))* y(k,36)
     234   114070872 :          prod(k,5) = 0._r8
     235   114070872 :          loss(k,6) = ( + het_rates(k,37))* y(k,37)
     236   114070872 :          prod(k,6) = 0._r8
     237   114070872 :          loss(k,7) = ( + het_rates(k,38))* y(k,38)
     238   114070872 :          prod(k,7) = 0._r8
     239             :          loss(k,81) = (rxt(k,507)* y(k,23) +rxt(k,508)* y(k,41) +rxt(k,510)* y(k,51) &
     240   114070872 :                   +rxt(k,509)* y(k,139) + het_rates(k,39))* y(k,39)
     241             :          prod(k,81) = (rxt(k,480)*y(k,14) +rxt(k,502)*y(k,10) + &
     242             :                  2.000_r8*rxt(k,511)*y(k,33) +rxt(k,512)*y(k,34))*y(k,71) +rxt(k,33) &
     243             :                  *y(k,10) +rxt(k,37)*y(k,14) +2.000_r8*rxt(k,49)*y(k,33) +rxt(k,50) &
     244   114070872 :                  *y(k,34) +rxt(k,57)*y(k,50)
     245             :          loss(k,132) = (rxt(k,411)* y(k,72) +rxt(k,414)* y(k,73) +rxt(k,105)* y(k,90) &
     246             :                   +rxt(k,117)* y(k,94) +rxt(k,129)* y(k,97) + (rxt(k,408) + &
     247             :                  rxt(k,409) +rxt(k,410))* y(k,101) +rxt(k,287)* y(k,111) +rxt(k,309) &
     248             :                  * y(k,122) +rxt(k,317)* y(k,127) +rxt(k,331)* y(k,130) +rxt(k,344) &
     249   114070872 :                  * y(k,134) + het_rates(k,40))* y(k,40)
     250             :          prod(k,132) = (rxt(k,141)*y(k,98) +rxt(k,147)*y(k,90) +rxt(k,158)*y(k,96) + &
     251             :                  rxt(k,162)*y(k,112) +rxt(k,163)*y(k,116) +rxt(k,164)*y(k,91) + &
     252             :                  rxt(k,165)*y(k,114) +rxt(k,166)*y(k,111) +rxt(k,170)*y(k,94) + &
     253             :                  rxt(k,181)*y(k,92) +rxt(k,203)*y(k,100) +rxt(k,214)*y(k,127) + &
     254             :                  rxt(k,225)*y(k,97) +rxt(k,236)*y(k,115) +rxt(k,247)*y(k,95) + &
     255             :                  rxt(k,258)*y(k,94) +rxt(k,260)*y(k,96) +rxt(k,262)*y(k,115) + &
     256             :                  rxt(k,271)*y(k,95))*y(k,105) + (rxt(k,144)*y(k,98) + &
     257             :                  rxt(k,168)*y(k,94) +rxt(k,169)*y(k,92) +rxt(k,172)*y(k,100) + &
     258             :                  rxt(k,173)*y(k,127) +rxt(k,174)*y(k,97) +rxt(k,175)*y(k,115) + &
     259             :                  rxt(k,176)*y(k,95) +rxt(k,177)*y(k,90) +rxt(k,178)*y(k,96) + &
     260             :                  rxt(k,179)*y(k,112) +rxt(k,180)*y(k,116) +rxt(k,182)*y(k,91) + &
     261             :                  rxt(k,183)*y(k,114) +rxt(k,184)*y(k,111) +rxt(k,259)*y(k,95) + &
     262             :                  rxt(k,261)*y(k,96) +rxt(k,263)*y(k,115) +rxt(k,265)*y(k,94))*y(k,107) &
     263             :                   + (rxt(k,186)*y(k,94) +rxt(k,187)*y(k,92) +rxt(k,189)*y(k,100) + &
     264             :                  rxt(k,190)*y(k,127) +rxt(k,191)*y(k,97) +rxt(k,193)*y(k,115) + &
     265             :                  rxt(k,194)*y(k,95) +rxt(k,195)*y(k,90) +rxt(k,196)*y(k,96) + &
     266             :                  rxt(k,197)*y(k,112) +rxt(k,198)*y(k,116) +rxt(k,199)*y(k,91) + &
     267             :                  rxt(k,200)*y(k,114) +rxt(k,201)*y(k,111) +rxt(k,383)*y(k,98)) &
     268             :                  *y(k,103) + (rxt(k,415)*y(k,41) +rxt(k,418)*y(k,70) + &
     269             :                  rxt(k,438)*y(k,56) +rxt(k,532)*y(k,17) +rxt(k,544)*y(k,75) + &
     270             :                  rxt(k,550)*y(k,78) +rxt(k,555)*y(k,80))*y(k,133) &
     271             :                   + (rxt(k,354)*y(k,135) +rxt(k,389)*y(k,71) +rxt(k,406)*y(k,70) + &
     272             :                  rxt(k,452)*y(k,25) +rxt(k,508)*y(k,39))*y(k,41) + (rxt(k,24) + &
     273             :                  .330_r8*rxt(k,25) +rxt(k,538)*y(k,71))*y(k,23) &
     274             :                   + (rxt(k,138)*y(k,99) +rxt(k,381)*y(k,108) +rxt(k,382)*y(k,102)) &
     275             :                  *y(k,98) + (rxt(k,52) +rxt(k,506)*y(k,71))*y(k,45) + (rxt(k,56) + &
     276             :                  rxt(k,483)*y(k,71))*y(k,49) + (2.000_r8*rxt(k,2) +rxt(k,3))*y(k,139) &
     277   114070872 :                   +2.000_r8*rxt(k,21)*y(k,17) +rxt(k,23)*y(k,22) +rxt(k,57)*y(k,50)
     278             :          loss(k,103) = (rxt(k,452)* y(k,25) +rxt(k,508)* y(k,39) +rxt(k,406)* y(k,70) &
     279             :                   +rxt(k,389)* y(k,71) +rxt(k,415)* y(k,133) + (rxt(k,353) + &
     280   114070872 :                  rxt(k,354))* y(k,135) + het_rates(k,41))* y(k,41)
     281             :          prod(k,103) = (1.440_r8*rxt(k,25) +rxt(k,539)*y(k,71))*y(k,23) +rxt(k,22) &
     282   114070872 :                  *y(k,17) +rxt(k,408)*y(k,101)*y(k,40) +rxt(k,1)*y(k,139)
     283   114070872 :          loss(k,28) = (rxt(k,504)* y(k,71) + rxt(k,51) + het_rates(k,42))* y(k,42)
     284   114070872 :          prod(k,28) = 0._r8
     285             :          loss(k,64) = (rxt(k,453)* y(k,25) +rxt(k,407)* y(k,70) +rxt(k,416)* y(k,133) &
     286   114070872 :                   + rxt(k,4) + het_rates(k,43))* y(k,43)
     287             :          prod(k,64) = (.500_r8*rxt(k,558) +rxt(k,422)*y(k,101))*y(k,101) &
     288   114070872 :                   +rxt(k,421)*y(k,133)*y(k,133)
     289   114070872 :          loss(k,33) = ( + rxt(k,100) + het_rates(k,44))* y(k,44)
     290   114070872 :          prod(k,33) =rxt(k,557)*y(k,139)*y(k,82)
     291             :          loss(k,79) = (rxt(k,499)* y(k,70) + (rxt(k,505) +rxt(k,506))* y(k,71) &
     292   114070872 :                   +rxt(k,500)* y(k,133) + rxt(k,52) + het_rates(k,45))* y(k,45)
     293   114070872 :          prod(k,79) = (rxt(k,486)*y(k,17) +rxt(k,487)*y(k,101))*y(k,3)
     294             :          loss(k,40) = (rxt(k,527)* y(k,71) +rxt(k,522)* y(k,133) + rxt(k,53) &
     295   114070872 :                   + het_rates(k,46))* y(k,46)
     296   114070872 :          prod(k,40) = 0._r8
     297             :          loss(k,42) = (rxt(k,528)* y(k,71) +rxt(k,523)* y(k,133) + rxt(k,54) &
     298   114070872 :                   + het_rates(k,47))* y(k,47)
     299   114070872 :          prod(k,42) = 0._r8
     300             :          loss(k,48) = (rxt(k,529)* y(k,71) +rxt(k,524)* y(k,133) + rxt(k,55) &
     301   114070872 :                   + het_rates(k,48))* y(k,48)
     302   114070872 :          prod(k,48) = 0._r8
     303             :          loss(k,138) = ((rxt(k,568) +rxt(k,573) +rxt(k,578))* y(k,29) + (rxt(k,570) + &
     304             :                  rxt(k,575))* y(k,53) + (rxt(k,563) +rxt(k,569) +rxt(k,574))* y(k,54) &
     305             :                   +rxt(k,470)* y(k,70) + (rxt(k,482) +rxt(k,483))* y(k,71) +rxt(k,108) &
     306             :                  * y(k,90) +rxt(k,106)* y(k,91) +rxt(k,130)* y(k,97) +rxt(k,289) &
     307             :                  * y(k,111) + (rxt(k,276) +rxt(k,298))* y(k,113) +rxt(k,318)* y(k,127) &
     308             :                   +rxt(k,471)* y(k,133) +rxt(k,345)* y(k,134) +rxt(k,356)* y(k,135) &
     309   114070872 :                   + rxt(k,56) + het_rates(k,49))* y(k,49)
     310             :          prod(k,138) = (rxt(k,450)*y(k,17) +rxt(k,451)*y(k,23) +rxt(k,452)*y(k,41) + &
     311             :                  rxt(k,453)*y(k,43) +rxt(k,454)*y(k,101) +rxt(k,472)*y(k,54) + &
     312             :                  rxt(k,513)*y(k,16) +rxt(k,515)*y(k,18) +2.000_r8*rxt(k,518)*y(k,20) + &
     313             :                  rxt(k,520)*y(k,24))*y(k,25) + (rxt(k,386) +rxt(k,169)*y(k,107) + &
     314             :                  rxt(k,181)*y(k,105) +rxt(k,187)*y(k,103) +rxt(k,205)*y(k,122) + &
     315             :                  rxt(k,222)*y(k,119) +rxt(k,240)*y(k,118) +rxt(k,257)*y(k,128) + &
     316             :                  2.000_r8*rxt(k,267)*y(k,105) +2.000_r8*rxt(k,268)*y(k,107))*y(k,92) &
     317             :                   + (rxt(k,157)*y(k,128) +rxt(k,163)*y(k,105) +rxt(k,180)*y(k,107) + &
     318             :                  rxt(k,198)*y(k,103) +rxt(k,216)*y(k,122) +rxt(k,233)*y(k,119) + &
     319             :                  rxt(k,251)*y(k,118) +rxt(k,299)*y(k,51))*y(k,116) &
     320             :                   + (rxt(k,105)*y(k,40) +rxt(k,109)*y(k,51))*y(k,90) &
     321   114070872 :                   +rxt(k,469)*y(k,133)*y(k,28)
     322   114070872 :          loss(k,49) = ( + rxt(k,57) + het_rates(k,50))* y(k,50)
     323             :          prod(k,49) = (rxt(k,507)*y(k,23) +rxt(k,508)*y(k,41) +rxt(k,509)*y(k,139) + &
     324   114070872 :                  rxt(k,510)*y(k,51))*y(k,39)
     325             :          loss(k,126) = (rxt(k,510)* y(k,39) +rxt(k,109)* y(k,90) +rxt(k,124)* y(k,94) &
     326             :                   +rxt(k,290)* y(k,111) +rxt(k,300)* y(k,113) +rxt(k,295)* y(k,115) &
     327             :                   +rxt(k,299)* y(k,116) +rxt(k,319)* y(k,127) +rxt(k,447)* y(k,133) &
     328   114070872 :                   +rxt(k,357)* y(k,135) + rxt(k,9) + het_rates(k,51))* y(k,51)
     329             :          prod(k,126) = (rxt(k,275) +2.000_r8*rxt(k,146)*y(k,105) + &
     330             :                  2.000_r8*rxt(k,167)*y(k,107) +2.000_r8*rxt(k,185)*y(k,103) + &
     331             :                  rxt(k,202)*y(k,122) +rxt(k,220)*y(k,119) +rxt(k,238)*y(k,118) + &
     332             :                  rxt(k,255)*y(k,128) +2.000_r8*rxt(k,269)*y(k,105) + &
     333             :                  2.000_r8*rxt(k,270)*y(k,107))*y(k,117) + (2.000_r8*rxt(k,559) + &
     334             :                  2.000_r8*rxt(k,562) +2.000_r8*rxt(k,565) +2.000_r8*rxt(k,576) + &
     335             :                  rxt(k,142)*y(k,105) +rxt(k,145)*y(k,107) +rxt(k,293)*y(k,114) + &
     336             :                  rxt(k,297)*y(k,115))*y(k,58) + (rxt(k,566) +rxt(k,571) +rxt(k,577) + &
     337             :                  rxt(k,568)*y(k,49) +rxt(k,573)*y(k,49) +rxt(k,578)*y(k,49))*y(k,29) &
     338             :                   + (rxt(k,171)*y(k,107) +rxt(k,188)*y(k,103) +rxt(k,192)*y(k,105) + &
     339             :                  rxt(k,264)*y(k,105) +rxt(k,266)*y(k,107) +rxt(k,298)*y(k,49)) &
     340             :                  *y(k,113) + (rxt(k,564) +rxt(k,567) +rxt(k,572))*y(k,6) &
     341             :                   + (rxt(k,561) +rxt(k,530)*y(k,17) +rxt(k,541)*y(k,35))*y(k,64) &
     342             :                   + (.500_r8*rxt(k,560) +rxt(k,446)*y(k,133))*y(k,63) &
     343   114070872 :                   + (rxt(k,140)*y(k,104) +rxt(k,143)*y(k,106))*y(k,139)
     344             :          loss(k,59) = (rxt(k,423)* y(k,133) + rxt(k,10) + rxt(k,11) + rxt(k,448) &
     345   114070872 :                   + het_rates(k,52))* y(k,52)
     346   114070872 :          prod(k,59) =rxt(k,444)*y(k,101)*y(k,63)
     347             :          loss(k,77) = ((rxt(k,570) +rxt(k,575))* y(k,49) +rxt(k,501)* y(k,70) &
     348   114070872 :                   + rxt(k,58) + het_rates(k,53))* y(k,53)
     349             :          prod(k,77) = (rxt(k,564) +rxt(k,567) +rxt(k,572))*y(k,6) +rxt(k,493)*y(k,101) &
     350   114070872 :                  *y(k,5)
     351             :          loss(k,80) = (rxt(k,472)* y(k,25) + (rxt(k,563) +rxt(k,569) +rxt(k,574)) &
     352             :                  * y(k,49) +rxt(k,473)* y(k,70) +rxt(k,474)* y(k,133) + rxt(k,59) &
     353   114070872 :                   + het_rates(k,54))* y(k,54)
     354             :          prod(k,80) = (rxt(k,566) +rxt(k,571) +rxt(k,577) +rxt(k,466)*y(k,133)) &
     355   114070872 :                  *y(k,29) +rxt(k,461)*y(k,101)*y(k,28)
     356   114070872 :          loss(k,73) = (rxt(k,340)* y(k,133) + rxt(k,12) + het_rates(k,55))* y(k,55)
     357             :          prod(k,73) = (rxt(k,289)*y(k,49) +rxt(k,290)*y(k,51))*y(k,111) &
     358   114070872 :                   +rxt(k,349)*y(k,133)*y(k,62) +rxt(k,305)*y(k,139)*y(k,120)
     359             :          loss(k,90) = (rxt(k,426)* y(k,62) + (rxt(k,427) +rxt(k,428) +rxt(k,429)) &
     360             :                  * y(k,63) +rxt(k,430)* y(k,72) +rxt(k,591)* y(k,128) +rxt(k,438) &
     361   114070872 :                  * y(k,133) + rxt(k,66) + het_rates(k,56))* y(k,56)
     362             :          prod(k,90) = (rxt(k,424)*y(k,109) +rxt(k,588)*y(k,124))*y(k,70) &
     363             :                   + (.200_r8*rxt(k,582)*y(k,118) +1.100_r8*rxt(k,584)*y(k,110)) &
     364             :                  *y(k,98) +rxt(k,17)*y(k,62) +rxt(k,589)*y(k,124)*y(k,72) +rxt(k,595) &
     365   114070872 :                  *y(k,136)
     366             :          loss(k,44) = ((rxt(k,442) +rxt(k,443))* y(k,71) + rxt(k,13) &
     367   114070872 :                   + het_rates(k,57))* y(k,57)
     368   114070872 :          prod(k,44) =rxt(k,427)*y(k,63)*y(k,56)
     369             :          loss(k,87) = (rxt(k,142)* y(k,105) +rxt(k,145)* y(k,107) +rxt(k,293) &
     370             :                  * y(k,114) +rxt(k,297)* y(k,115) + rxt(k,14) + rxt(k,15) + rxt(k,449) &
     371             :                   + rxt(k,559) + rxt(k,562) + rxt(k,565) + rxt(k,576) &
     372   114070872 :                   + het_rates(k,58))* y(k,58)
     373   114070872 :          prod(k,87) =rxt(k,445)*y(k,64)*y(k,63)
     374   114070872 :          loss(k,8) = ( + het_rates(k,59))* y(k,59)
     375   114070872 :          prod(k,8) = 0._r8
     376   114070872 :          loss(k,9) = ( + het_rates(k,60))* y(k,60)
     377   114070872 :          prod(k,9) = 0._r8
     378   114070872 :          loss(k,10) = ( + het_rates(k,61))* y(k,61)
     379   114070872 :          prod(k,10) = 0._r8
     380             :          loss(k,137) = (rxt(k,494)* y(k,5) +rxt(k,534)* y(k,21) +rxt(k,462)* y(k,28) &
     381             :                   +rxt(k,426)* y(k,56) +rxt(k,435)* y(k,64) +rxt(k,441)* y(k,70) &
     382             :                   +rxt(k,440)* y(k,73) + (rxt(k,111) +rxt(k,112))* y(k,93) +rxt(k,280) &
     383             :                  * y(k,94) + (rxt(k,122) +rxt(k,123))* y(k,96) +rxt(k,439)* y(k,101) &
     384             :                   +rxt(k,593)* y(k,128) + (rxt(k,272) +rxt(k,279))* y(k,130) &
     385             :                   +rxt(k,349)* y(k,133) +rxt(k,136)* y(k,135) + rxt(k,16) + rxt(k,17) &
     386   114070872 :                   + het_rates(k,62))* y(k,62)
     387             :          prod(k,137) = (rxt(k,202)*y(k,117) +rxt(k,204)*y(k,94) +rxt(k,205)*y(k,92) + &
     388             :                  rxt(k,206)*y(k,113) +rxt(k,207)*y(k,100) +rxt(k,208)*y(k,127) + &
     389             :                  rxt(k,209)*y(k,97) +rxt(k,210)*y(k,115) +rxt(k,211)*y(k,95) + &
     390             :                  rxt(k,212)*y(k,90) +rxt(k,213)*y(k,96) +rxt(k,215)*y(k,112) + &
     391             :                  rxt(k,216)*y(k,116) +rxt(k,217)*y(k,91) +rxt(k,218)*y(k,114) + &
     392             :                  rxt(k,219)*y(k,111) +rxt(k,308)*y(k,98) +rxt(k,309)*y(k,40))*y(k,122) &
     393             :                   + (rxt(k,220)*y(k,117) +rxt(k,221)*y(k,94) +rxt(k,222)*y(k,92) + &
     394             :                  rxt(k,223)*y(k,113) +rxt(k,224)*y(k,100) +rxt(k,226)*y(k,127) + &
     395             :                  rxt(k,227)*y(k,97) +rxt(k,228)*y(k,115) +rxt(k,229)*y(k,95) + &
     396             :                  rxt(k,230)*y(k,90) +rxt(k,231)*y(k,96) +rxt(k,232)*y(k,112) + &
     397             :                  rxt(k,233)*y(k,116) +rxt(k,234)*y(k,91) +rxt(k,235)*y(k,114) + &
     398             :                  rxt(k,237)*y(k,111) +rxt(k,303)*y(k,98))*y(k,119) &
     399             :                   + (rxt(k,238)*y(k,117) +rxt(k,239)*y(k,94) +rxt(k,240)*y(k,92) + &
     400             :                  rxt(k,241)*y(k,113) +rxt(k,242)*y(k,100) +rxt(k,243)*y(k,127) + &
     401             :                  rxt(k,244)*y(k,97) +rxt(k,245)*y(k,115) +rxt(k,246)*y(k,95) + &
     402             :                  rxt(k,248)*y(k,90) +rxt(k,249)*y(k,96) +rxt(k,250)*y(k,112) + &
     403             :                  rxt(k,251)*y(k,116) +rxt(k,252)*y(k,91) +rxt(k,253)*y(k,114) + &
     404             :                  rxt(k,254)*y(k,111))*y(k,118) + (rxt(k,18) +.500_r8*rxt(k,560) + &
     405             :                  rxt(k,291)*y(k,111) +2.000_r8*rxt(k,428)*y(k,56) + &
     406             :                  rxt(k,431)*y(k,70) +rxt(k,551)*y(k,80))*y(k,63) &
     407             :                   + (rxt(k,304)*y(k,120) +rxt(k,306)*y(k,121) +rxt(k,384)*y(k,123)) &
     408             :                  *y(k,98) + (rxt(k,430)*y(k,72) +rxt(k,438)*y(k,133))*y(k,56) &
     409             :                   +rxt(k,287)*y(k,111)*y(k,40) +rxt(k,12)*y(k,55) &
     410             :                   +2.000_r8*rxt(k,442)*y(k,71)*y(k,57) +rxt(k,15)*y(k,58) +rxt(k,20) &
     411             :                  *y(k,64) +rxt(k,425)*y(k,109)*y(k,72) +rxt(k,592)*y(k,128) &
     412   114070872 :                   +rxt(k,605)*y(k,138)
     413             :          loss(k,133) = (rxt(k,495)* y(k,5) +rxt(k,464)* y(k,28) + (rxt(k,427) + &
     414             :                  rxt(k,428) +rxt(k,429))* y(k,56) +rxt(k,445)* y(k,64) + (rxt(k,431) + &
     415             :                  rxt(k,433))* y(k,70) +rxt(k,432)* y(k,73) +rxt(k,551)* y(k,80) &
     416             :                   +rxt(k,110)* y(k,90) +rxt(k,281)* y(k,94) + (rxt(k,120) +rxt(k,121)) &
     417             :                  * y(k,96) +rxt(k,444)* y(k,101) +rxt(k,291)* y(k,111) +rxt(k,320) &
     418             :                  * y(k,127) + (rxt(k,277) +rxt(k,278))* y(k,130) +rxt(k,446)* y(k,133) &
     419             :                   +rxt(k,346)* y(k,134) +rxt(k,359)* y(k,135) + rxt(k,18) + rxt(k,560) &
     420   114070872 :                   + het_rates(k,63))* y(k,63)
     421             :          prod(k,133) = (rxt(k,112)*y(k,93) +rxt(k,136)*y(k,135) + &
     422             :                  2.000_r8*rxt(k,435)*y(k,64) +rxt(k,439)*y(k,101) + &
     423             :                  rxt(k,440)*y(k,73) +rxt(k,441)*y(k,70) +rxt(k,462)*y(k,28) + &
     424             :                  rxt(k,494)*y(k,5) +rxt(k,534)*y(k,21))*y(k,62) + (rxt(k,75) + &
     425             :                  rxt(k,161)*y(k,128) +rxt(k,166)*y(k,105) +rxt(k,184)*y(k,107) + &
     426             :                  rxt(k,201)*y(k,103) +rxt(k,219)*y(k,122) +rxt(k,237)*y(k,119) + &
     427             :                  rxt(k,254)*y(k,118) +rxt(k,285)*y(k,25))*y(k,111) &
     428             :                   + (rxt(k,156)*y(k,128) +rxt(k,162)*y(k,105) +rxt(k,179)*y(k,107) + &
     429             :                  rxt(k,197)*y(k,103) +rxt(k,215)*y(k,122) +rxt(k,232)*y(k,119) + &
     430             :                  rxt(k,250)*y(k,118))*y(k,112) + (rxt(k,19) +rxt(k,434)*y(k,101) + &
     431             :                  rxt(k,436)*y(k,70) +rxt(k,437)*y(k,133))*y(k,64) + (rxt(k,11) + &
     432             :                  rxt(k,448) +rxt(k,423)*y(k,133))*y(k,52) + (rxt(k,14) +rxt(k,449)) &
     433             :                  *y(k,58) + (rxt(k,311)*y(k,122) +rxt(k,340)*y(k,55))*y(k,133) &
     434   114070872 :                   +rxt(k,29)*y(k,6) +rxt(k,48)*y(k,29) +rxt(k,9)*y(k,51)
     435             :          loss(k,131) = (rxt(k,530)* y(k,17) +rxt(k,541)* y(k,35) +rxt(k,435)* y(k,62) &
     436             :                   +rxt(k,445)* y(k,63) +rxt(k,436)* y(k,70) +rxt(k,434)* y(k,101) &
     437             :                   +rxt(k,437)* y(k,133) + rxt(k,19) + rxt(k,20) + rxt(k,561) &
     438   114070872 :                   + het_rates(k,64))* y(k,64)
     439             :          prod(k,131) = (rxt(k,152)*y(k,128) +rxt(k,175)*y(k,107) + &
     440             :                  rxt(k,193)*y(k,103) +rxt(k,210)*y(k,122) +rxt(k,228)*y(k,119) + &
     441             :                  rxt(k,236)*y(k,105) +rxt(k,245)*y(k,118) +rxt(k,262)*y(k,105) + &
     442             :                  rxt(k,263)*y(k,107))*y(k,115) + (rxt(k,160)*y(k,128) + &
     443             :                  rxt(k,165)*y(k,105) +rxt(k,183)*y(k,107) +rxt(k,200)*y(k,103) + &
     444             :                  rxt(k,218)*y(k,122) +rxt(k,235)*y(k,119) +rxt(k,253)*y(k,118)) &
     445             :                  *y(k,114) + (rxt(k,157)*y(k,128) +rxt(k,163)*y(k,105) + &
     446             :                  rxt(k,180)*y(k,107) +rxt(k,198)*y(k,103) +rxt(k,216)*y(k,122) + &
     447             :                  rxt(k,233)*y(k,119) +rxt(k,251)*y(k,118))*y(k,116) + (rxt(k,76) + &
     448             :                  rxt(k,148)*y(k,128) +rxt(k,206)*y(k,122) +rxt(k,223)*y(k,119) + &
     449             :                  rxt(k,241)*y(k,118))*y(k,113) + (rxt(k,47) +rxt(k,463)*y(k,25) + &
     450             :                  rxt(k,465)*y(k,70) +rxt(k,466)*y(k,133))*y(k,29) &
     451             :                   + (rxt(k,202)*y(k,122) +rxt(k,220)*y(k,119) +rxt(k,238)*y(k,118) + &
     452             :                  rxt(k,255)*y(k,128))*y(k,117) + (rxt(k,14) +rxt(k,15) +rxt(k,449)) &
     453             :                  *y(k,58) + (rxt(k,30) +rxt(k,496)*y(k,70))*y(k,6) &
     454             :                   + (rxt(k,447)*y(k,133) +rxt(k,510)*y(k,39))*y(k,51) &
     455             :                   + (rxt(k,432)*y(k,73) +rxt(k,433)*y(k,70))*y(k,63) &
     456             :                   +rxt(k,286)*y(k,111)*y(k,28) +rxt(k,10)*y(k,52) +rxt(k,310)*y(k,122) &
     457   114070872 :                  *y(k,101)
     458   114070872 :          loss(k,11) = ( + het_rates(k,65))* y(k,65)
     459   114070872 :          prod(k,11) = 0._r8
     460   114070872 :          loss(k,12) = ( + het_rates(k,66))* y(k,66)
     461   114070872 :          prod(k,12) = 0._r8
     462   114070872 :          loss(k,13) = ( + het_rates(k,67))* y(k,67)
     463   114070872 :          prod(k,13) = 0._r8
     464   114070872 :          loss(k,14) = ( + het_rates(k,68))* y(k,68)
     465   114070872 :          prod(k,14) = 0._r8
     466   114070872 :          loss(k,15) = ( + het_rates(k,69))* y(k,69)
     467   114070872 :          prod(k,15) = 0._r8
     468             :          loss(k,123) = (rxt(k,497)* y(k,5) +rxt(k,496)* y(k,6) +rxt(k,531)* y(k,17) &
     469             :                   +rxt(k,467)* y(k,28) +rxt(k,465)* y(k,29) +rxt(k,406)* y(k,41) &
     470             :                   +rxt(k,407)* y(k,43) +rxt(k,499)* y(k,45) +rxt(k,470)* y(k,49) &
     471             :                   +rxt(k,501)* y(k,53) +rxt(k,473)* y(k,54) +rxt(k,441)* y(k,62) &
     472             :                   + (rxt(k,431) +rxt(k,433))* y(k,63) +rxt(k,436)* y(k,64) &
     473             :                   + 2._r8*rxt(k,404)* y(k,70) +rxt(k,405)* y(k,72) +rxt(k,403) &
     474             :                  * y(k,73) +rxt(k,543)* y(k,75) +rxt(k,113)* y(k,93) +rxt(k,125) &
     475             :                  * y(k,94) +rxt(k,131)* y(k,97) +rxt(k,412)* y(k,101) + (rxt(k,586) + &
     476             :                  rxt(k,587))* y(k,110) +rxt(k,301)* y(k,113) +rxt(k,588)* y(k,124) &
     477             :                   + (rxt(k,324) +rxt(k,325))* y(k,127) + (rxt(k,333) +rxt(k,334)) &
     478             :                  * y(k,130) +rxt(k,336)* y(k,131) +rxt(k,338)* y(k,132) +rxt(k,418) &
     479             :                  * y(k,133) +rxt(k,347)* y(k,134) +rxt(k,360)* y(k,135) + rxt(k,77) &
     480             :                   + rxt(k,78) + rxt(k,79) + rxt(k,80) + rxt(k,81) + rxt(k,82) &
     481   114070872 :                   + het_rates(k,70))* y(k,70)
     482             :          prod(k,123) = (2.000_r8*rxt(k,5) +rxt(k,6) +rxt(k,83) +rxt(k,85) +rxt(k,87) + &
     483             :                  2.000_r8*rxt(k,88) +2.000_r8*rxt(k,89) +rxt(k,90) +rxt(k,91) + &
     484             :                  rxt(k,92) +rxt(k,392)*y(k,71) +rxt(k,393)*y(k,71) + &
     485             :                  rxt(k,430)*y(k,56) +rxt(k,545)*y(k,78) +rxt(k,552)*y(k,80) + &
     486             :                  rxt(k,590)*y(k,124) +rxt(k,597)*y(k,136) +rxt(k,601)*y(k,137)) &
     487             :                  *y(k,72) + (rxt(k,114)*y(k,25) +rxt(k,168)*y(k,107) + &
     488             :                  rxt(k,170)*y(k,105) +rxt(k,186)*y(k,103) +rxt(k,204)*y(k,122) + &
     489             :                  rxt(k,221)*y(k,119) +rxt(k,239)*y(k,118) +rxt(k,256)*y(k,128) + &
     490             :                  rxt(k,258)*y(k,105) +rxt(k,265)*y(k,107))*y(k,94) &
     491             :                   + (rxt(k,153)*y(k,128) +rxt(k,176)*y(k,107) +rxt(k,194)*y(k,103) + &
     492             :                  rxt(k,211)*y(k,122) +rxt(k,229)*y(k,119) +rxt(k,246)*y(k,118) + &
     493             :                  rxt(k,247)*y(k,105) +rxt(k,259)*y(k,107) +rxt(k,271)*y(k,105)) &
     494             :                  *y(k,95) + (rxt(k,155)*y(k,128) +rxt(k,158)*y(k,105) + &
     495             :                  rxt(k,178)*y(k,107) +rxt(k,196)*y(k,103) +rxt(k,213)*y(k,122) + &
     496             :                  rxt(k,231)*y(k,119) +rxt(k,249)*y(k,118) +rxt(k,260)*y(k,105) + &
     497             :                  rxt(k,261)*y(k,107))*y(k,96) + (rxt(k,99) +rxt(k,358) + &
     498             :                  rxt(k,350)*y(k,25) +rxt(k,359)*y(k,63) +rxt(k,363)*y(k,73))*y(k,135) &
     499             :                   + (rxt(k,426)*y(k,62) +rxt(k,427)*y(k,63) +rxt(k,591)*y(k,128)) &
     500             :                  *y(k,56) + (rxt(k,26) +rxt(k,62))*y(k,32) + (rxt(k,17) + &
     501             :                  rxt(k,272)*y(k,130))*y(k,62) + (rxt(k,582)*y(k,118) + &
     502             :                  1.150_r8*rxt(k,583)*y(k,128))*y(k,98) +rxt(k,28)*y(k,5) &
     503             :                   +.180_r8*rxt(k,25)*y(k,23) +rxt(k,46)*y(k,28) +rxt(k,410)*y(k,101) &
     504             :                  *y(k,40) +rxt(k,15)*y(k,58) +rxt(k,18)*y(k,63) +rxt(k,19)*y(k,64) &
     505             :                   +rxt(k,391)*y(k,71) +rxt(k,8)*y(k,73) +rxt(k,60)*y(k,74) +rxt(k,102) &
     506             :                  *y(k,80) +rxt(k,103)*y(k,81) +rxt(k,104)*y(k,82) +rxt(k,596)*y(k,136) &
     507             :                  *y(k,109) +rxt(k,420)*y(k,133)*y(k,133) +rxt(k,599)*y(k,137) &
     508   114070872 :                   +rxt(k,604)*y(k,138) +rxt(k,2)*y(k,139)
     509             :          loss(k,108) = (rxt(k,475)* y(k,8) +rxt(k,476)* y(k,9) +rxt(k,502)* y(k,10) &
     510             :                   +rxt(k,477)* y(k,11) +rxt(k,478)* y(k,12) +rxt(k,479)* y(k,13) &
     511             :                   +rxt(k,480)* y(k,14) +rxt(k,481)* y(k,15) +rxt(k,525)* y(k,16) &
     512             :                   +rxt(k,526)* y(k,18) + (rxt(k,537) +rxt(k,538) +rxt(k,539))* y(k,23) &
     513             :                   +rxt(k,503)* y(k,24) +rxt(k,511)* y(k,33) +rxt(k,512)* y(k,34) &
     514             :                   +rxt(k,389)* y(k,41) +rxt(k,504)* y(k,42) + (rxt(k,505) +rxt(k,506)) &
     515             :                  * y(k,45) +rxt(k,527)* y(k,46) +rxt(k,528)* y(k,47) +rxt(k,529) &
     516             :                  * y(k,48) + (rxt(k,482) +rxt(k,483))* y(k,49) + (rxt(k,442) + &
     517             :                  rxt(k,443))* y(k,57) + (rxt(k,392) +rxt(k,393))* y(k,72) +rxt(k,394) &
     518             :                  * y(k,73) +rxt(k,390)* y(k,139) + rxt(k,391) + het_rates(k,71)) &
     519   114070872 :                  * y(k,71)
     520             :          prod(k,108) = (rxt(k,6) +rxt(k,425)*y(k,109))*y(k,72) +rxt(k,13)*y(k,57) &
     521             :                   +rxt(k,7)*y(k,73) +.850_r8*rxt(k,583)*y(k,128)*y(k,98) +rxt(k,1) &
     522   114070872 :                  *y(k,139)
     523             :          loss(k,124) = (rxt(k,411)* y(k,40) +rxt(k,430)* y(k,56) +rxt(k,405)* y(k,70) &
     524             :                   +rxt(k,392)* y(k,71) +rxt(k,545)* y(k,78) +rxt(k,552)* y(k,80) &
     525             :                   +rxt(k,126)* y(k,94) + (rxt(k,133) +rxt(k,135))* y(k,98) +rxt(k,425) &
     526             :                  * y(k,109) +rxt(k,585)* y(k,110) + (rxt(k,589) +rxt(k,590))* y(k,124) &
     527             :                   +rxt(k,397)* y(k,125) +rxt(k,322)* y(k,127) +rxt(k,364)* y(k,128) &
     528             :                   +rxt(k,362)* y(k,135) +rxt(k,597)* y(k,136) +rxt(k,601)* y(k,137) &
     529             :                   + rxt(k,5) + rxt(k,6) + rxt(k,83) + rxt(k,84) + rxt(k,85) &
     530             :                   + rxt(k,86) + rxt(k,87) + rxt(k,88) + rxt(k,89) + rxt(k,90) &
     531   114070872 :                   + rxt(k,91) + rxt(k,92) + het_rates(k,72))* y(k,72)
     532             :          prod(k,124) = (rxt(k,8) +rxt(k,132)*y(k,97) +rxt(k,134)*y(k,98) + &
     533             :                  rxt(k,292)*y(k,111) +2.000_r8*rxt(k,302)*y(k,113) + &
     534             :                  rxt(k,323)*y(k,127) +3.000_r8*rxt(k,332)*y(k,130) + &
     535             :                  2.000_r8*rxt(k,394)*y(k,71) +2.000_r8*rxt(k,403)*y(k,70) + &
     536             :                  2.000_r8*rxt(k,413)*y(k,101) +rxt(k,414)*y(k,40) + &
     537             :                  rxt(k,419)*y(k,133) +rxt(k,432)*y(k,63) +rxt(k,440)*y(k,62) + &
     538             :                  rxt(k,456)*y(k,25) +rxt(k,488)*y(k,3) +rxt(k,547)*y(k,78) + &
     539             :                  rxt(k,553)*y(k,80))*y(k,73) + (rxt(k,113)*y(k,93) + &
     540             :                  rxt(k,131)*y(k,97) +rxt(k,301)*y(k,113) +rxt(k,325)*y(k,127) + &
     541             :                  2.000_r8*rxt(k,333)*y(k,130) +rxt(k,334)*y(k,130) + &
     542             :                  rxt(k,336)*y(k,131) +rxt(k,360)*y(k,135) +rxt(k,396)*y(k,125) + &
     543             :                  rxt(k,404)*y(k,70) +rxt(k,412)*y(k,101) +rxt(k,418)*y(k,133) + &
     544             :                  rxt(k,431)*y(k,63) +rxt(k,436)*y(k,64) +rxt(k,467)*y(k,28) + &
     545             :                  rxt(k,497)*y(k,5))*y(k,70) + (rxt(k,148)*y(k,113) + &
     546             :                  rxt(k,149)*y(k,100) +2.000_r8*rxt(k,151)*y(k,97) + &
     547             :                  rxt(k,152)*y(k,115) +rxt(k,153)*y(k,95) +rxt(k,154)*y(k,90) + &
     548             :                  rxt(k,155)*y(k,96) +rxt(k,156)*y(k,112) +rxt(k,157)*y(k,116) + &
     549             :                  rxt(k,159)*y(k,91) +rxt(k,160)*y(k,114) +rxt(k,161)*y(k,111) + &
     550             :                  rxt(k,255)*y(k,117) +rxt(k,256)*y(k,94) +rxt(k,257)*y(k,92) + &
     551             :                  rxt(k,593)*y(k,62))*y(k,128) + (rxt(k,93) +rxt(k,137) + &
     552             :                  rxt(k,173)*y(k,107) +rxt(k,190)*y(k,103) +rxt(k,208)*y(k,122) + &
     553             :                  rxt(k,214)*y(k,105) +rxt(k,226)*y(k,119) +rxt(k,243)*y(k,118) + &
     554             :                  rxt(k,314)*y(k,25) +rxt(k,315)*y(k,28) +rxt(k,320)*y(k,63) + &
     555             :                  2.000_r8*rxt(k,321)*y(k,125))*y(k,127) + (rxt(k,116)*y(k,94) + &
     556             :                  rxt(k,128)*y(k,97) +rxt(k,351)*y(k,135) +rxt(k,458)*y(k,28) + &
     557             :                  rxt(k,459)*y(k,28) +rxt(k,461)*y(k,101) +rxt(k,469)*y(k,133) + &
     558             :                  rxt(k,491)*y(k,5) +rxt(k,492)*y(k,5))*y(k,28) + (rxt(k,408)*y(k,40) + &
     559             :                  rxt(k,417)*y(k,133) +rxt(k,422)*y(k,101) +rxt(k,434)*y(k,64) + &
     560             :                  rxt(k,454)*y(k,25) +rxt(k,487)*y(k,3) +rxt(k,493)*y(k,5) + &
     561             :                  rxt(k,533)*y(k,21))*y(k,101) + (rxt(k,127)*y(k,25) + &
     562             :                  rxt(k,174)*y(k,107) +rxt(k,191)*y(k,103) +rxt(k,209)*y(k,122) + &
     563             :                  rxt(k,225)*y(k,105) +rxt(k,227)*y(k,119) +rxt(k,244)*y(k,118)) &
     564             :                  *y(k,97) + (rxt(k,95) +rxt(k,277)*y(k,63) +rxt(k,279)*y(k,62) + &
     565             :                  rxt(k,330)*y(k,32) +rxt(k,331)*y(k,40))*y(k,130) + (rxt(k,387) + &
     566             :                  rxt(k,395) +2.000_r8*rxt(k,339)*y(k,132) + &
     567             :                  2.000_r8*rxt(k,397)*y(k,72))*y(k,125) + (rxt(k,326)*y(k,98) + &
     568             :                  rxt(k,327)*y(k,139) +rxt(k,328)*y(k,139))*y(k,129) + (rxt(k,97) + &
     569             :                  rxt(k,335)*y(k,32))*y(k,131) + (rxt(k,337)*y(k,139) + &
     570             :                  2.000_r8*rxt(k,380)*y(k,98))*y(k,132) +rxt(k,489)*y(k,5)*y(k,5) &
     571             :                   +rxt(k,423)*y(k,133)*y(k,52) +rxt(k,429)*y(k,63)*y(k,56) &
     572   114070872 :                   +rxt(k,443)*y(k,71)*y(k,57) +rxt(k,20)*y(k,64) +rxt(k,388)*y(k,126)
     573             :          loss(k,136) = (rxt(k,488)* y(k,3) +rxt(k,456)* y(k,25) +rxt(k,414)* y(k,40) &
     574             :                   +rxt(k,440)* y(k,62) +rxt(k,432)* y(k,63) +rxt(k,403)* y(k,70) &
     575             :                   +rxt(k,394)* y(k,71) +rxt(k,547)* y(k,78) +rxt(k,553)* y(k,80) &
     576             :                   +rxt(k,132)* y(k,97) +rxt(k,134)* y(k,98) +rxt(k,413)* y(k,101) &
     577             :                   +rxt(k,292)* y(k,111) +rxt(k,302)* y(k,113) +rxt(k,323)* y(k,127) &
     578             :                   +rxt(k,332)* y(k,130) +rxt(k,419)* y(k,133) +rxt(k,348)* y(k,134) &
     579             :                   +rxt(k,363)* y(k,135) + rxt(k,7) + rxt(k,8) + het_rates(k,73)) &
     580   114070872 :                  * y(k,73)
     581             :          prod(k,136) = (rxt(k,324)*y(k,127) +rxt(k,338)*y(k,132) +rxt(k,405)*y(k,72)) &
     582             :                  *y(k,70) + (rxt(k,96) +rxt(k,278)*y(k,63))*y(k,130) &
     583   114070872 :                   +rxt(k,361)*y(k,135)*y(k,125)
     584   114070872 :          loss(k,52) = (rxt(k,554)* y(k,80) + rxt(k,60) + het_rates(k,74))* y(k,74)
     585   114070872 :          prod(k,52) = (rxt(k,460)*y(k,28) +rxt(k,490)*y(k,5))*y(k,28)
     586             :          loss(k,53) = (rxt(k,543)* y(k,70) +rxt(k,544)* y(k,133) + rxt(k,101) &
     587   114070872 :                   + het_rates(k,75))* y(k,75)
     588   114070872 :          prod(k,53) = 0._r8
     589   114070872 :          loss(k,16) = ( + het_rates(k,76))* y(k,76)
     590   114070872 :          prod(k,16) = 0._r8
     591   114070872 :          loss(k,17) = ( + het_rates(k,77))* y(k,77)
     592   114070872 :          prod(k,17) = 0._r8
     593             :          loss(k,72) = (rxt(k,545)* y(k,72) +rxt(k,547)* y(k,73) +rxt(k,550)* y(k,133) &
     594   114070872 :                   + het_rates(k,78))* y(k,78)
     595   114070872 :          prod(k,72) =rxt(k,101)*y(k,75) +rxt(k,102)*y(k,80)
     596   114070872 :          loss(k,18) = ( + rxt(k,61) + het_rates(k,79))* y(k,79)
     597   114070872 :          prod(k,18) = 0._r8
     598             :          loss(k,92) = (rxt(k,548)* y(k,5) +rxt(k,549)* y(k,28) +rxt(k,551)* y(k,63) &
     599             :                   +rxt(k,552)* y(k,72) +rxt(k,553)* y(k,73) +rxt(k,554)* y(k,74) &
     600   114070872 :                   +rxt(k,555)* y(k,133) + rxt(k,102) + het_rates(k,80))* y(k,80)
     601             :          prod(k,92) = (rxt(k,545)*y(k,72) +rxt(k,547)*y(k,73) +rxt(k,550)*y(k,133)) &
     602   114070872 :                  *y(k,78) +rxt(k,543)*y(k,75)*y(k,70) +rxt(k,103)*y(k,81)
     603   114070872 :          loss(k,82) = (rxt(k,546)* y(k,133) + rxt(k,103) + het_rates(k,81))* y(k,81)
     604             :          prod(k,82) = (rxt(k,548)*y(k,5) +rxt(k,549)*y(k,28) +rxt(k,551)*y(k,63) + &
     605             :                  rxt(k,552)*y(k,72) +rxt(k,553)*y(k,73) +rxt(k,554)*y(k,74) + &
     606             :                  rxt(k,555)*y(k,133))*y(k,80) + (rxt(k,541)*y(k,64) + &
     607             :                  rxt(k,542)*y(k,133) +.500_r8*rxt(k,556)*y(k,133))*y(k,35) &
     608   114070872 :                   +rxt(k,544)*y(k,133)*y(k,75) +rxt(k,104)*y(k,82)
     609   114070872 :          loss(k,45) = (rxt(k,557)* y(k,139) + rxt(k,104) + het_rates(k,82))* y(k,82)
     610   114070872 :          prod(k,45) =rxt(k,100)*y(k,44) +rxt(k,546)*y(k,133)*y(k,81)
     611   114070872 :          loss(k,19) = ( + het_rates(k,83))* y(k,83)
     612   114070872 :          prod(k,19) = 0._r8
     613   114070872 :          loss(k,20) = ( + het_rates(k,84))* y(k,84)
     614   114070872 :          prod(k,20) = 0._r8
     615   114070872 :          loss(k,21) = ( + het_rates(k,85))* y(k,85)
     616   114070872 :          prod(k,21) = 0._r8
     617   114070872 :          loss(k,22) = ( + het_rates(k,86))* y(k,86)
     618   114070872 :          prod(k,22) = 0._r8
     619   114070872 :          loss(k,23) = ( + het_rates(k,87))* y(k,87)
     620   114070872 :          prod(k,23) = 0._r8
     621   114070872 :          loss(k,24) = ( + het_rates(k,88))* y(k,88)
     622   114070872 :          prod(k,24) = 0._r8
     623   114070872 :          loss(k,25) = ( + het_rates(k,89))* y(k,89)
     624   114070872 :          prod(k,25) = 0._r8
     625             :          loss(k,109) = (rxt(k,105)* y(k,40) +rxt(k,108)* y(k,49) +rxt(k,109)* y(k,51) &
     626             :                   +rxt(k,110)* y(k,63) +rxt(k,195)* y(k,103) +rxt(k,147)* y(k,105) &
     627             :                   +rxt(k,177)* y(k,107) +rxt(k,248)* y(k,118) +rxt(k,230)* y(k,119) &
     628             :                   +rxt(k,212)* y(k,122) +rxt(k,154)* y(k,128) +rxt(k,107)* y(k,139) &
     629   114070872 :                   + het_rates(k,90))* y(k,90)
     630             :          prod(k,109) = (rxt(k,114)*y(k,94) +rxt(k,127)*y(k,97) +rxt(k,285)*y(k,111) + &
     631             :                  rxt(k,314)*y(k,127) +rxt(k,341)*y(k,134) +rxt(k,350)*y(k,135)) &
     632             :                  *y(k,25) + (rxt(k,130)*y(k,97) +rxt(k,289)*y(k,111) + &
     633             :                  rxt(k,298)*y(k,113) +rxt(k,318)*y(k,127) +rxt(k,345)*y(k,134) + &
     634             :                  rxt(k,356)*y(k,135))*y(k,49) + (rxt(k,116)*y(k,94) + &
     635             :                  rxt(k,286)*y(k,111) +rxt(k,351)*y(k,135))*y(k,28) &
     636             :                   + (rxt(k,112)*y(k,62) +rxt(k,113)*y(k,70))*y(k,93) +rxt(k,385) &
     637   114070872 :                  *y(k,91) +rxt(k,386)*y(k,92)
     638             :          loss(k,96) = (rxt(k,106)* y(k,49) +rxt(k,199)* y(k,103) +rxt(k,164)* y(k,105) &
     639             :                   +rxt(k,182)* y(k,107) +rxt(k,252)* y(k,118) +rxt(k,234)* y(k,119) &
     640             :                   +rxt(k,217)* y(k,122) +rxt(k,159)* y(k,128) + rxt(k,385) &
     641   114070872 :                   + het_rates(k,91))* y(k,91)
     642   114070872 :          prod(k,96) =rxt(k,107)*y(k,139)*y(k,90)
     643             :          loss(k,95) = (rxt(k,187)* y(k,103) + (rxt(k,181) +rxt(k,267))* y(k,105) &
     644             :                   + (rxt(k,169) +rxt(k,268))* y(k,107) +rxt(k,240)* y(k,118) &
     645             :                   +rxt(k,222)* y(k,119) +rxt(k,205)* y(k,122) +rxt(k,257)* y(k,128) &
     646   114070872 :                   + rxt(k,386) + het_rates(k,92))* y(k,92)
     647   114070872 :          prod(k,95) = (rxt(k,106)*y(k,91) +rxt(k,108)*y(k,90))*y(k,49)
     648             :          loss(k,84) = ((rxt(k,111) +rxt(k,112))* y(k,62) +rxt(k,113)* y(k,70) &
     649   114070872 :                   + het_rates(k,93))* y(k,93)
     650             :          prod(k,84) = (rxt(k,128)*y(k,97) +rxt(k,315)*y(k,127) +rxt(k,342)*y(k,134)) &
     651   114070872 :                  *y(k,28) +rxt(k,115)*y(k,94)*y(k,25)
     652             :          loss(k,115) = ((rxt(k,114) +rxt(k,115))* y(k,25) +rxt(k,116)* y(k,28) &
     653             :                   +rxt(k,117)* y(k,40) +rxt(k,124)* y(k,51) +rxt(k,280)* y(k,62) &
     654             :                   +rxt(k,281)* y(k,63) +rxt(k,125)* y(k,70) +rxt(k,126)* y(k,72) &
     655             :                   +rxt(k,186)* y(k,103) + (rxt(k,170) +rxt(k,258))* y(k,105) &
     656             :                   + (rxt(k,168) +rxt(k,265))* y(k,107) +rxt(k,239)* y(k,118) &
     657             :                   +rxt(k,221)* y(k,119) +rxt(k,204)* y(k,122) +rxt(k,256)* y(k,128) &
     658   114070872 :                   +rxt(k,119)* y(k,139) + rxt(k,63) + het_rates(k,94))* y(k,94)
     659             :          prod(k,115) = (rxt(k,330)*y(k,130) +rxt(k,352)*y(k,135))*y(k,32) &
     660             :                   + (rxt(k,64) +rxt(k,283))*y(k,96) + (rxt(k,129)*y(k,40) + &
     661   114070872 :                  rxt(k,131)*y(k,70))*y(k,97)
     662             :          loss(k,94) = (rxt(k,194)* y(k,103) + (rxt(k,247) +rxt(k,271))* y(k,105) &
     663             :                   + (rxt(k,176) +rxt(k,259))* y(k,107) +rxt(k,246)* y(k,118) &
     664             :                   +rxt(k,229)* y(k,119) +rxt(k,211)* y(k,122) +rxt(k,153)* y(k,128) &
     665   114070872 :                   + rxt(k,284) + het_rates(k,95))* y(k,95)
     666   114070872 :          prod(k,94) =rxt(k,118)*y(k,139)*y(k,96)
     667             :          loss(k,104) = ((rxt(k,122) +rxt(k,123))* y(k,62) + (rxt(k,120) +rxt(k,121)) &
     668             :                  * y(k,63) +rxt(k,196)* y(k,103) + (rxt(k,158) +rxt(k,260))* y(k,105) &
     669             :                   + (rxt(k,178) +rxt(k,261))* y(k,107) +rxt(k,249)* y(k,118) &
     670             :                   +rxt(k,231)* y(k,119) +rxt(k,213)* y(k,122) +rxt(k,155)* y(k,128) &
     671             :                   +rxt(k,118)* y(k,139) + rxt(k,64) + rxt(k,283) + het_rates(k,96)) &
     672   114070872 :                  * y(k,96)
     673   114070872 :          prod(k,104) =rxt(k,119)*y(k,139)*y(k,94) +rxt(k,284)*y(k,95)
     674             :          loss(k,110) = (rxt(k,127)* y(k,25) +rxt(k,128)* y(k,28) +rxt(k,129)* y(k,40) &
     675             :                   +rxt(k,130)* y(k,49) +rxt(k,131)* y(k,70) +rxt(k,132)* y(k,73) &
     676             :                   +rxt(k,191)* y(k,103) +rxt(k,225)* y(k,105) +rxt(k,174)* y(k,107) &
     677             :                   +rxt(k,244)* y(k,118) +rxt(k,227)* y(k,119) +rxt(k,209)* y(k,122) &
     678   114070872 :                   +rxt(k,151)* y(k,128) + rxt(k,65) + het_rates(k,97))* y(k,97)
     679   114070872 :          prod(k,110) = (rxt(k,316)*y(k,127) +rxt(k,335)*y(k,131))*y(k,32)
     680             :          loss(k,116) = ((rxt(k,133) +rxt(k,135))* y(k,72) +rxt(k,134)* y(k,73) &
     681             :                   +rxt(k,138)* y(k,99) +rxt(k,382)* y(k,102) +rxt(k,383)* y(k,103) &
     682             :                   +rxt(k,141)* y(k,105) +rxt(k,144)* y(k,107) +rxt(k,381)* y(k,108) &
     683             :                   +rxt(k,584)* y(k,110) +rxt(k,582)* y(k,118) +rxt(k,303)* y(k,119) &
     684             :                   +rxt(k,304)* y(k,120) +rxt(k,306)* y(k,121) +rxt(k,308)* y(k,122) &
     685             :                   +rxt(k,384)* y(k,123) +rxt(k,583)* y(k,128) +rxt(k,326)* y(k,129) &
     686   114070872 :                   +rxt(k,380)* y(k,132) + het_rates(k,98))* y(k,98)
     687             :          prod(k,116) = (rxt(k,77) +rxt(k,78) +rxt(k,79) +rxt(k,80) +rxt(k,81) + &
     688             :                  rxt(k,82) +rxt(k,324)*y(k,127) +rxt(k,333)*y(k,130) + &
     689             :                  rxt(k,347)*y(k,134) +rxt(k,360)*y(k,135))*y(k,70) + (rxt(k,83) + &
     690             :                  rxt(k,84) +rxt(k,85) +rxt(k,86) +rxt(k,87) +rxt(k,90) +rxt(k,91) + &
     691             :                  rxt(k,92))*y(k,72) + (rxt(k,99) +rxt(k,358) +rxt(k,136)*y(k,62) + &
     692             :                  rxt(k,353)*y(k,41) +rxt(k,361)*y(k,125))*y(k,135) + (rxt(k,93) + &
     693             :                  rxt(k,137) +rxt(k,317)*y(k,40) +rxt(k,321)*y(k,125))*y(k,127) &
     694             :                   + (rxt(k,105)*y(k,90) +rxt(k,344)*y(k,134))*y(k,40) + (rxt(k,96) + &
     695             :                  rxt(k,332)*y(k,73))*y(k,130) +rxt(k,66)*y(k,56) +rxt(k,16)*y(k,62) &
     696   114070872 :                   +rxt(k,75)*y(k,111) +rxt(k,76)*y(k,113) +rxt(k,98)*y(k,134)
     697             :          loss(k,57) = (rxt(k,138)* y(k,98) +rxt(k,139)* y(k,139) + het_rates(k,99)) &
     698   114070872 :                  * y(k,99)
     699   114070872 :          prod(k,57) =rxt(k,327)*y(k,139)*y(k,129)
     700             :          loss(k,93) = (rxt(k,189)* y(k,103) +rxt(k,203)* y(k,105) +rxt(k,172) &
     701             :                  * y(k,107) +rxt(k,242)* y(k,118) +rxt(k,224)* y(k,119) +rxt(k,207) &
     702   114070872 :                  * y(k,122) +rxt(k,149)* y(k,128) + het_rates(k,100))* y(k,100)
     703   114070872 :          prod(k,93) =rxt(k,343)*y(k,134)*y(k,32)
     704             :          loss(k,112) = (rxt(k,487)* y(k,3) +rxt(k,493)* y(k,5) +rxt(k,533)* y(k,21) &
     705             :                   + (rxt(k,454) +rxt(k,455))* y(k,25) +rxt(k,461)* y(k,28) &
     706             :                   + (rxt(k,408) +rxt(k,409) +rxt(k,410))* y(k,40) +rxt(k,439)* y(k,62) &
     707             :                   +rxt(k,444)* y(k,63) +rxt(k,434)* y(k,64) +rxt(k,412)* y(k,70) &
     708             :                   +rxt(k,413)* y(k,73) + 2._r8*rxt(k,422)* y(k,101) +rxt(k,310) &
     709             :                  * y(k,122) +rxt(k,417)* y(k,133) + rxt(k,558) + het_rates(k,101)) &
     710   114070872 :                  * y(k,101)
     711             :          prod(k,112) = (rxt(k,416)*y(k,43) +rxt(k,419)*y(k,73) +rxt(k,437)*y(k,64) + &
     712             :                  rxt(k,468)*y(k,28) +rxt(k,498)*y(k,5) +rxt(k,516)*y(k,18) + &
     713             :                  rxt(k,519)*y(k,20) +rxt(k,540)*y(k,31) +rxt(k,546)*y(k,81) + &
     714             :                  .500_r8*rxt(k,556)*y(k,35))*y(k,133) + (rxt(k,450)*y(k,25) + &
     715             :                  rxt(k,486)*y(k,3) +rxt(k,530)*y(k,64) +rxt(k,531)*y(k,70))*y(k,17) &
     716             :                   + (rxt(k,453)*y(k,43) +rxt(k,515)*y(k,18) +rxt(k,518)*y(k,20)) &
     717             :                  *y(k,25) + (rxt(k,317)*y(k,40) +rxt(k,318)*y(k,49) + &
     718             :                  rxt(k,319)*y(k,51))*y(k,127) + (rxt(k,457)*y(k,28) + &
     719             :                  rxt(k,534)*y(k,62))*y(k,21) + (rxt(k,11) +rxt(k,448))*y(k,52) &
     720             :                   + (rxt(k,347)*y(k,134) +rxt(k,407)*y(k,43))*y(k,70) &
     721             :                   +rxt(k,538)*y(k,71)*y(k,23) +rxt(k,411)*y(k,72)*y(k,40) &
     722   114070872 :                   +rxt(k,130)*y(k,97)*y(k,49)
     723             :          loss(k,68) = (rxt(k,382)* y(k,98) +rxt(k,374)* y(k,139) + rxt(k,373) &
     724   114070872 :                   + het_rates(k,102))* y(k,102)
     725             :          prod(k,68) = (rxt(k,139)*y(k,99) +rxt(k,372)*y(k,108))*y(k,139) +rxt(k,375) &
     726   114070872 :                  *y(k,103)
     727             :          loss(k,120) = (rxt(k,195)* y(k,90) +rxt(k,199)* y(k,91) +rxt(k,187)* y(k,92) &
     728             :                   +rxt(k,186)* y(k,94) +rxt(k,194)* y(k,95) +rxt(k,196)* y(k,96) &
     729             :                   +rxt(k,191)* y(k,97) +rxt(k,383)* y(k,98) +rxt(k,189)* y(k,100) &
     730             :                   +rxt(k,201)* y(k,111) +rxt(k,197)* y(k,112) +rxt(k,188)* y(k,113) &
     731             :                   +rxt(k,200)* y(k,114) +rxt(k,193)* y(k,115) +rxt(k,198)* y(k,116) &
     732             :                   +rxt(k,185)* y(k,117) +rxt(k,190)* y(k,127) +rxt(k,376)* y(k,139) &
     733   114070872 :                   + rxt(k,375) + het_rates(k,103))* y(k,103)
     734             :          prod(k,120) = (rxt(k,305)*y(k,120) +rxt(k,374)*y(k,102))*y(k,139) +rxt(k,377) &
     735   114070872 :                  *y(k,105)
     736   114070872 :          loss(k,46) = (rxt(k,140)* y(k,139) + het_rates(k,104))* y(k,104)
     737   114070872 :          prod(k,46) =rxt(k,142)*y(k,105)*y(k,58)
     738             :          loss(k,130) = (rxt(k,142)* y(k,58) +rxt(k,147)* y(k,90) +rxt(k,164)* y(k,91) &
     739             :                   + (rxt(k,181) +rxt(k,267))* y(k,92) + (rxt(k,170) +rxt(k,258)) &
     740             :                  * y(k,94) + (rxt(k,247) +rxt(k,271))* y(k,95) + (rxt(k,158) + &
     741             :                  rxt(k,260))* y(k,96) +rxt(k,225)* y(k,97) +rxt(k,141)* y(k,98) &
     742             :                   +rxt(k,203)* y(k,100) +rxt(k,166)* y(k,111) +rxt(k,162)* y(k,112) &
     743             :                   + (rxt(k,192) +rxt(k,264))* y(k,113) +rxt(k,165)* y(k,114) &
     744             :                   + (rxt(k,236) +rxt(k,262))* y(k,115) +rxt(k,163)* y(k,116) &
     745             :                   + (rxt(k,146) +rxt(k,269))* y(k,117) +rxt(k,214)* y(k,127) &
     746   114070872 :                   +rxt(k,378)* y(k,139) + rxt(k,377) + het_rates(k,105))* y(k,105)
     747             :          prod(k,130) = (rxt(k,140)*y(k,104) +rxt(k,376)*y(k,103))*y(k,139) +rxt(k,379) &
     748   114070872 :                  *y(k,107)
     749   114070872 :          loss(k,47) = (rxt(k,143)* y(k,139) + het_rates(k,106))* y(k,106)
     750   114070872 :          prod(k,47) =rxt(k,145)*y(k,107)*y(k,58)
     751             :          loss(k,129) = (rxt(k,145)* y(k,58) +rxt(k,177)* y(k,90) +rxt(k,182)* y(k,91) &
     752             :                   + (rxt(k,169) +rxt(k,268))* y(k,92) + (rxt(k,168) +rxt(k,265)) &
     753             :                  * y(k,94) + (rxt(k,176) +rxt(k,259))* y(k,95) + (rxt(k,178) + &
     754             :                  rxt(k,261))* y(k,96) +rxt(k,174)* y(k,97) +rxt(k,144)* y(k,98) &
     755             :                   +rxt(k,172)* y(k,100) +rxt(k,184)* y(k,111) +rxt(k,179)* y(k,112) &
     756             :                   + (rxt(k,171) +rxt(k,266))* y(k,113) +rxt(k,183)* y(k,114) &
     757             :                   + (rxt(k,175) +rxt(k,263))* y(k,115) +rxt(k,180)* y(k,116) &
     758             :                   + (rxt(k,167) +rxt(k,270))* y(k,117) +rxt(k,173)* y(k,127) &
     759   114070872 :                   + rxt(k,379) + het_rates(k,107))* y(k,107)
     760   114070872 :          prod(k,129) = (rxt(k,143)*y(k,106) +rxt(k,378)*y(k,105))*y(k,139)
     761             :          loss(k,74) = (rxt(k,381)* y(k,98) +rxt(k,372)* y(k,139) + het_rates(k,108)) &
     762   114070872 :                  * y(k,108)
     763             :          prod(k,74) = (rxt(k,309)*y(k,40) +rxt(k,310)*y(k,101) +rxt(k,311)*y(k,133)) &
     764   114070872 :                  *y(k,122) +rxt(k,373)*y(k,102) +rxt(k,328)*y(k,139)*y(k,129)
     765             :          loss(k,76) = (rxt(k,424)* y(k,70) +rxt(k,425)* y(k,72) +rxt(k,596)* y(k,136) &
     766   114070872 :                   + het_rates(k,109))* y(k,109)
     767             :          prod(k,76) = (.800_r8*rxt(k,582)*y(k,118) +.900_r8*rxt(k,584)*y(k,110)) &
     768   114070872 :                  *y(k,98) +rxt(k,586)*y(k,110)*y(k,70)
     769             :          loss(k,66) = ((rxt(k,586) +rxt(k,587))* y(k,70) +rxt(k,585)* y(k,72) &
     770   114070872 :                   +rxt(k,584)* y(k,98) + het_rates(k,110))* y(k,110)
     771   114070872 :          prod(k,66) =rxt(k,599)*y(k,137) +rxt(k,604)*y(k,138)
     772             :          loss(k,113) = (rxt(k,285)* y(k,25) +rxt(k,286)* y(k,28) +rxt(k,287)* y(k,40) &
     773             :                   +rxt(k,289)* y(k,49) +rxt(k,290)* y(k,51) +rxt(k,291)* y(k,63) &
     774             :                   +rxt(k,292)* y(k,73) +rxt(k,201)* y(k,103) +rxt(k,166)* y(k,105) &
     775             :                   +rxt(k,184)* y(k,107) +rxt(k,254)* y(k,118) +rxt(k,237)* y(k,119) &
     776             :                   +rxt(k,219)* y(k,122) +rxt(k,161)* y(k,128) +rxt(k,288)* y(k,139) &
     777   114070872 :                   + rxt(k,75) + het_rates(k,111))* y(k,111)
     778             :          prod(k,113) = (rxt(k,110)*y(k,90) +rxt(k,278)*y(k,130) +rxt(k,320)*y(k,127) + &
     779             :                  rxt(k,346)*y(k,134) +rxt(k,359)*y(k,135))*y(k,63) &
     780             :                   + (rxt(k,111)*y(k,93) +rxt(k,123)*y(k,96) +rxt(k,279)*y(k,130) + &
     781             :                  rxt(k,280)*y(k,94))*y(k,62) + (rxt(k,301)*y(k,70) + &
     782   114070872 :                  rxt(k,302)*y(k,73))*y(k,113) +rxt(k,273)*y(k,112)
     783             :          loss(k,97) = (rxt(k,197)* y(k,103) +rxt(k,162)* y(k,105) +rxt(k,179) &
     784             :                  * y(k,107) +rxt(k,250)* y(k,118) +rxt(k,232)* y(k,119) +rxt(k,215) &
     785             :                  * y(k,122) +rxt(k,156)* y(k,128) + rxt(k,273) + het_rates(k,112)) &
     786   114070872 :                  * y(k,112)
     787   114070872 :          prod(k,97) =rxt(k,122)*y(k,96)*y(k,62) +rxt(k,288)*y(k,139)*y(k,111)
     788             :          loss(k,111) = ((rxt(k,276) +rxt(k,298))* y(k,49) +rxt(k,300)* y(k,51) &
     789             :                   +rxt(k,301)* y(k,70) +rxt(k,302)* y(k,73) +rxt(k,188)* y(k,103) &
     790             :                   + (rxt(k,192) +rxt(k,264))* y(k,105) + (rxt(k,171) +rxt(k,266)) &
     791             :                  * y(k,107) +rxt(k,241)* y(k,118) +rxt(k,223)* y(k,119) +rxt(k,206) &
     792             :                  * y(k,122) +rxt(k,148)* y(k,128) +rxt(k,296)* y(k,139) + rxt(k,76) &
     793   114070872 :                   + het_rates(k,113))* y(k,113)
     794             :          prod(k,111) = (rxt(k,109)*y(k,90) +rxt(k,124)*y(k,94) +rxt(k,290)*y(k,111) + &
     795             :                  rxt(k,319)*y(k,127) +rxt(k,357)*y(k,135))*y(k,51) &
     796             :                   + (rxt(k,120)*y(k,96) +rxt(k,277)*y(k,130) +rxt(k,281)*y(k,94) + &
     797             :                  rxt(k,291)*y(k,111))*y(k,63) +rxt(k,272)*y(k,130)*y(k,62) &
     798             :                   +rxt(k,292)*y(k,111)*y(k,73) +rxt(k,282)*y(k,115) +rxt(k,275) &
     799   114070872 :                  *y(k,117)
     800             :          loss(k,102) = (rxt(k,293)* y(k,58) +rxt(k,200)* y(k,103) +rxt(k,165) &
     801             :                  * y(k,105) +rxt(k,183)* y(k,107) +rxt(k,253)* y(k,118) +rxt(k,235) &
     802             :                  * y(k,119) +rxt(k,218)* y(k,122) +rxt(k,160)* y(k,128) + rxt(k,274) &
     803   114070872 :                   + het_rates(k,114))* y(k,114)
     804   114070872 :          prod(k,102) =rxt(k,294)*y(k,139)*y(k,115)
     805             :          loss(k,105) = (rxt(k,295)* y(k,51) +rxt(k,297)* y(k,58) +rxt(k,193)* y(k,103) &
     806             :                   + (rxt(k,236) +rxt(k,262))* y(k,105) + (rxt(k,175) +rxt(k,263)) &
     807             :                  * y(k,107) +rxt(k,245)* y(k,118) +rxt(k,228)* y(k,119) +rxt(k,210) &
     808             :                  * y(k,122) +rxt(k,152)* y(k,128) +rxt(k,294)* y(k,139) + rxt(k,282) &
     809   114070872 :                   + het_rates(k,115))* y(k,115)
     810             :          prod(k,105) =rxt(k,121)*y(k,96)*y(k,63) +rxt(k,296)*y(k,139)*y(k,113) &
     811   114070872 :                   +rxt(k,274)*y(k,114)
     812             :          loss(k,98) = (rxt(k,299)* y(k,51) +rxt(k,198)* y(k,103) +rxt(k,163)* y(k,105) &
     813             :                   +rxt(k,180)* y(k,107) +rxt(k,251)* y(k,118) +rxt(k,233)* y(k,119) &
     814             :                   +rxt(k,216)* y(k,122) +rxt(k,157)* y(k,128) + het_rates(k,116)) &
     815   114070872 :                  * y(k,116)
     816   114070872 :          prod(k,98) =rxt(k,276)*y(k,113)*y(k,49)
     817             :          loss(k,101) = (rxt(k,185)* y(k,103) + (rxt(k,146) +rxt(k,269))* y(k,105) &
     818             :                   + (rxt(k,167) +rxt(k,270))* y(k,107) +rxt(k,238)* y(k,118) &
     819             :                   +rxt(k,220)* y(k,119) +rxt(k,202)* y(k,122) +rxt(k,255)* y(k,128) &
     820   114070872 :                   + rxt(k,275) + het_rates(k,117))* y(k,117)
     821             :          prod(k,101) = (rxt(k,295)*y(k,115) +rxt(k,299)*y(k,116) +rxt(k,300)*y(k,113)) &
     822   114070872 :                  *y(k,51) + (rxt(k,293)*y(k,114) +rxt(k,297)*y(k,115))*y(k,58)
     823             :          loss(k,118) = (rxt(k,368)* y(k,32) +rxt(k,248)* y(k,90) +rxt(k,252)* y(k,91) &
     824             :                   +rxt(k,240)* y(k,92) +rxt(k,239)* y(k,94) +rxt(k,246)* y(k,95) &
     825             :                   +rxt(k,249)* y(k,96) +rxt(k,244)* y(k,97) +rxt(k,582)* y(k,98) &
     826             :                   +rxt(k,242)* y(k,100) +rxt(k,254)* y(k,111) +rxt(k,250)* y(k,112) &
     827             :                   +rxt(k,241)* y(k,113) +rxt(k,253)* y(k,114) +rxt(k,245)* y(k,115) &
     828             :                   +rxt(k,251)* y(k,116) +rxt(k,238)* y(k,117) +rxt(k,243)* y(k,127) &
     829   114070872 :                   +rxt(k,365)* y(k,139) + rxt(k,370) + het_rates(k,118))* y(k,118)
     830             :          prod(k,118) = (rxt(k,592) +rxt(k,591)*y(k,56) +rxt(k,593)*y(k,62))*y(k,128) &
     831             :                   +rxt(k,16)*y(k,62) +rxt(k,586)*y(k,110)*y(k,70) +rxt(k,590)*y(k,124) &
     832             :                  *y(k,72) +rxt(k,369)*y(k,121) +rxt(k,371)*y(k,123) +rxt(k,595) &
     833   114070872 :                  *y(k,136)
     834             :          loss(k,119) = (rxt(k,230)* y(k,90) +rxt(k,234)* y(k,91) +rxt(k,222)* y(k,92) &
     835             :                   +rxt(k,221)* y(k,94) +rxt(k,229)* y(k,95) +rxt(k,231)* y(k,96) &
     836             :                   +rxt(k,227)* y(k,97) +rxt(k,303)* y(k,98) +rxt(k,224)* y(k,100) &
     837             :                   +rxt(k,237)* y(k,111) +rxt(k,232)* y(k,112) +rxt(k,223)* y(k,113) &
     838             :                   +rxt(k,235)* y(k,114) +rxt(k,228)* y(k,115) +rxt(k,233)* y(k,116) &
     839             :                   +rxt(k,220)* y(k,117) +rxt(k,226)* y(k,127) +rxt(k,367)* y(k,139) &
     840   114070872 :                   + het_rates(k,119))* y(k,119)
     841   114070872 :          prod(k,119) =rxt(k,366)*y(k,139)*y(k,122)
     842             :          loss(k,58) = (rxt(k,304)* y(k,98) +rxt(k,305)* y(k,139) + het_rates(k,120)) &
     843   114070872 :                  * y(k,120)
     844   114070872 :          prod(k,58) =rxt(k,367)*y(k,139)*y(k,119)
     845             :          loss(k,71) = (rxt(k,306)* y(k,98) +rxt(k,307)* y(k,139) + rxt(k,369) &
     846   114070872 :                   + het_rates(k,121))* y(k,121)
     847   114070872 :          prod(k,71) = (rxt(k,312)*y(k,123) +rxt(k,368)*y(k,118))*y(k,32)
     848             :          loss(k,122) = (rxt(k,309)* y(k,40) +rxt(k,212)* y(k,90) +rxt(k,217)* y(k,91) &
     849             :                   +rxt(k,205)* y(k,92) +rxt(k,204)* y(k,94) +rxt(k,211)* y(k,95) &
     850             :                   +rxt(k,213)* y(k,96) +rxt(k,209)* y(k,97) +rxt(k,308)* y(k,98) &
     851             :                   +rxt(k,207)* y(k,100) +rxt(k,310)* y(k,101) +rxt(k,219)* y(k,111) &
     852             :                   +rxt(k,215)* y(k,112) +rxt(k,206)* y(k,113) +rxt(k,218)* y(k,114) &
     853             :                   +rxt(k,210)* y(k,115) +rxt(k,216)* y(k,116) +rxt(k,202)* y(k,117) &
     854             :                   +rxt(k,208)* y(k,127) +rxt(k,311)* y(k,133) +rxt(k,366)* y(k,139) &
     855   114070872 :                   + het_rates(k,122))* y(k,122)
     856             :          prod(k,122) = (rxt(k,307)*y(k,121) +rxt(k,313)*y(k,123) +rxt(k,365)*y(k,118)) &
     857   114070872 :                  *y(k,139)
     858             :          loss(k,70) = (rxt(k,312)* y(k,32) +rxt(k,384)* y(k,98) +rxt(k,313)* y(k,139) &
     859   114070872 :                   + rxt(k,371) + het_rates(k,123))* y(k,123)
     860   114070872 :          prod(k,70) =rxt(k,370)*y(k,118)
     861             :          loss(k,78) = (rxt(k,588)* y(k,70) + (rxt(k,589) +rxt(k,590))* y(k,72) &
     862   114070872 :                   + het_rates(k,124))* y(k,124)
     863             :          prod(k,78) =rxt(k,66)*y(k,56) +rxt(k,596)*y(k,136)*y(k,109) +rxt(k,605) &
     864   114070872 :                  *y(k,138)
     865             :          loss(k,85) = (rxt(k,396)* y(k,70) +rxt(k,397)* y(k,72) +rxt(k,321)* y(k,127) &
     866             :                   +rxt(k,339)* y(k,132) +rxt(k,361)* y(k,135) + rxt(k,387) &
     867   114070872 :                   + rxt(k,395) + het_rates(k,125))* y(k,125)
     868             :          prod(k,85) = (rxt(k,399) +rxt(k,398)*y(k,32) +rxt(k,400)*y(k,70) + &
     869   114070872 :                  rxt(k,401)*y(k,72) +rxt(k,402)*y(k,73))*y(k,126) +rxt(k,7)*y(k,73)
     870             :          loss(k,51) = (rxt(k,398)* y(k,32) +rxt(k,400)* y(k,70) +rxt(k,401)* y(k,72) &
     871             :                   +rxt(k,402)* y(k,73) + rxt(k,388) + rxt(k,399) + het_rates(k,126)) &
     872   114070872 :                  * y(k,126)
     873   114070872 :          prod(k,51) =rxt(k,392)*y(k,72)*y(k,71)
     874             :          loss(k,127) = (rxt(k,314)* y(k,25) +rxt(k,315)* y(k,28) +rxt(k,316)* y(k,32) &
     875             :                   +rxt(k,317)* y(k,40) +rxt(k,318)* y(k,49) +rxt(k,319)* y(k,51) &
     876             :                   +rxt(k,320)* y(k,63) + (rxt(k,324) +rxt(k,325))* y(k,70) +rxt(k,322) &
     877             :                  * y(k,72) +rxt(k,323)* y(k,73) +rxt(k,190)* y(k,103) +rxt(k,214) &
     878             :                  * y(k,105) +rxt(k,173)* y(k,107) +rxt(k,243)* y(k,118) +rxt(k,226) &
     879             :                  * y(k,119) +rxt(k,208)* y(k,122) +rxt(k,321)* y(k,125) +rxt(k,150) &
     880   114070872 :                  * y(k,128) + rxt(k,93) + rxt(k,137) + het_rates(k,127))* y(k,127)
     881             :          prod(k,127) = (rxt(k,125)*y(k,94) +rxt(k,334)*y(k,130))*y(k,70) &
     882             :                   + (rxt(k,133)*y(k,98) +rxt(k,135)*y(k,98))*y(k,72) +rxt(k,65) &
     883   114070872 :                  *y(k,97) +rxt(k,97)*y(k,131)
     884             :          loss(k,128) = (rxt(k,591)* y(k,56) +rxt(k,593)* y(k,62) +rxt(k,364)* y(k,72) &
     885             :                   +rxt(k,154)* y(k,90) +rxt(k,159)* y(k,91) +rxt(k,257)* y(k,92) &
     886             :                   +rxt(k,256)* y(k,94) +rxt(k,153)* y(k,95) +rxt(k,155)* y(k,96) &
     887             :                   +rxt(k,151)* y(k,97) +rxt(k,583)* y(k,98) +rxt(k,149)* y(k,100) &
     888             :                   +rxt(k,161)* y(k,111) +rxt(k,156)* y(k,112) +rxt(k,148)* y(k,113) &
     889             :                   +rxt(k,160)* y(k,114) +rxt(k,152)* y(k,115) +rxt(k,157)* y(k,116) &
     890             :                   +rxt(k,255)* y(k,117) +rxt(k,150)* y(k,127) +rxt(k,329)* y(k,139) &
     891   114070872 :                   + rxt(k,592) + het_rates(k,128))* y(k,128)
     892             :          prod(k,128) = (rxt(k,84) +rxt(k,86) +rxt(k,585)*y(k,110) + &
     893             :                  rxt(k,589)*y(k,124) +rxt(k,597)*y(k,136) +rxt(k,601)*y(k,137)) &
     894             :                  *y(k,72) + (rxt(k,338)*y(k,70) +rxt(k,339)*y(k,125))*y(k,132) &
     895             :                   +rxt(k,594)*y(k,136)*y(k,32) +2.000_r8*rxt(k,150)*y(k,128)*y(k,127) &
     896   114070872 :                   +rxt(k,94)*y(k,129)
     897             :          loss(k,75) = (rxt(k,326)* y(k,98) + (rxt(k,327) +rxt(k,328))* y(k,139) &
     898   114070872 :                   + rxt(k,94) + het_rates(k,129))* y(k,129)
     899   114070872 :          prod(k,75) = (rxt(k,329)*y(k,128) +rxt(k,337)*y(k,132))*y(k,139)
     900             :          loss(k,100) = (rxt(k,330)* y(k,32) +rxt(k,331)* y(k,40) + (rxt(k,272) + &
     901             :                  rxt(k,279))* y(k,62) + (rxt(k,277) +rxt(k,278))* y(k,63) &
     902             :                   + (rxt(k,333) +rxt(k,334))* y(k,70) +rxt(k,332)* y(k,73) + rxt(k,95) &
     903   114070872 :                   + rxt(k,96) + het_rates(k,130))* y(k,130)
     904             :          prod(k,100) = (rxt(k,132)*y(k,97) +rxt(k,323)*y(k,127) +rxt(k,348)*y(k,134) + &
     905             :                  rxt(k,363)*y(k,135))*y(k,73) + (rxt(k,126)*y(k,94) + &
     906   114070872 :                  rxt(k,362)*y(k,135))*y(k,72) +rxt(k,336)*y(k,131)*y(k,70)
     907             :          loss(k,67) = (rxt(k,335)* y(k,32) +rxt(k,336)* y(k,70) + rxt(k,97) &
     908   114070872 :                   + het_rates(k,131))* y(k,131)
     909   114070872 :          prod(k,67) =rxt(k,322)*y(k,127)*y(k,72)
     910             :          loss(k,91) = (rxt(k,338)* y(k,70) +rxt(k,380)* y(k,98) +rxt(k,339)* y(k,125) &
     911   114070872 :                   +rxt(k,337)* y(k,139) + het_rates(k,132))* y(k,132)
     912   114070872 :          prod(k,91) =rxt(k,364)*y(k,128)*y(k,72)
     913             :          loss(k,114) = (rxt(k,498)* y(k,5) +rxt(k,514)* y(k,16) +rxt(k,532)* y(k,17) &
     914             :                   +rxt(k,516)* y(k,18) +rxt(k,517)* y(k,19) +rxt(k,519)* y(k,20) &
     915             :                   +rxt(k,535)* y(k,22) +rxt(k,536)* y(k,23) +rxt(k,521)* y(k,24) &
     916             :                   + (rxt(k,468) +rxt(k,469))* y(k,28) +rxt(k,466)* y(k,29) +rxt(k,540) &
     917             :                  * y(k,31) + (rxt(k,542) +rxt(k,556))* y(k,35) +rxt(k,415)* y(k,41) &
     918             :                   +rxt(k,416)* y(k,43) +rxt(k,500)* y(k,45) +rxt(k,522)* y(k,46) &
     919             :                   +rxt(k,523)* y(k,47) +rxt(k,524)* y(k,48) +rxt(k,471)* y(k,49) &
     920             :                   +rxt(k,447)* y(k,51) +rxt(k,423)* y(k,52) +rxt(k,474)* y(k,54) &
     921             :                   +rxt(k,340)* y(k,55) +rxt(k,438)* y(k,56) +rxt(k,349)* y(k,62) &
     922             :                   +rxt(k,446)* y(k,63) +rxt(k,437)* y(k,64) +rxt(k,418)* y(k,70) &
     923             :                   +rxt(k,419)* y(k,73) +rxt(k,544)* y(k,75) +rxt(k,550)* y(k,78) &
     924             :                   +rxt(k,555)* y(k,80) +rxt(k,546)* y(k,81) +rxt(k,417)* y(k,101) &
     925             :                   +rxt(k,311)* y(k,122) + 2._r8*(rxt(k,420) +rxt(k,421))* y(k,133) &
     926   114070872 :                   + het_rates(k,133))* y(k,133)
     927             :          prod(k,114) = (rxt(k,406)*y(k,41) +rxt(k,407)*y(k,43) +rxt(k,412)*y(k,101) + &
     928             :                  rxt(k,470)*y(k,49) +rxt(k,473)*y(k,54) +rxt(k,499)*y(k,45) + &
     929             :                  rxt(k,501)*y(k,53) +rxt(k,531)*y(k,17))*y(k,70) &
     930             :                   + (rxt(k,149)*y(k,128) +rxt(k,172)*y(k,107) +rxt(k,189)*y(k,103) + &
     931             :                  rxt(k,203)*y(k,105) +rxt(k,207)*y(k,122) +rxt(k,224)*y(k,119) + &
     932             :                  rxt(k,242)*y(k,118))*y(k,100) + (rxt(k,3) +rxt(k,139)*y(k,99) + &
     933             :                  rxt(k,328)*y(k,129) +rxt(k,355)*y(k,135) + &
     934             :                  2.000_r8*rxt(k,390)*y(k,71) +rxt(k,509)*y(k,39))*y(k,139) &
     935             :                   + (2.000_r8*rxt(k,409)*y(k,40) +rxt(k,413)*y(k,73) + &
     936             :                  rxt(k,434)*y(k,64) +rxt(k,439)*y(k,62) +rxt(k,455)*y(k,25))*y(k,101) &
     937             :                   + (rxt(k,98) +rxt(k,341)*y(k,25) +rxt(k,342)*y(k,28) + &
     938             :                  rxt(k,346)*y(k,63) +rxt(k,348)*y(k,73))*y(k,134) &
     939             :                   + (rxt(k,389)*y(k,41) +rxt(k,482)*y(k,49) +rxt(k,505)*y(k,45) + &
     940             :                  rxt(k,537)*y(k,23))*y(k,71) + (rxt(k,9) +rxt(k,124)*y(k,94) + &
     941             :                  rxt(k,357)*y(k,135))*y(k,51) + (rxt(k,23) + &
     942             :                  .300_r8*rxt(k,535)*y(k,133))*y(k,22) + (rxt(k,129)*y(k,97) + &
     943             :                  rxt(k,414)*y(k,73))*y(k,40) +.330_r8*rxt(k,25)*y(k,23) &
     944             :                   +2.000_r8*rxt(k,4)*y(k,43) +rxt(k,356)*y(k,135)*y(k,49) +rxt(k,10) &
     945             :                  *y(k,52) +rxt(k,58)*y(k,53) +rxt(k,59)*y(k,54) +rxt(k,12)*y(k,55) &
     946   114070872 :                   +.500_r8*rxt(k,560)*y(k,63) +rxt(k,138)*y(k,99)*y(k,98)
     947             :          loss(k,134) = (rxt(k,341)* y(k,25) +rxt(k,342)* y(k,28) +rxt(k,343)* y(k,32) &
     948             :                   +rxt(k,344)* y(k,40) +rxt(k,345)* y(k,49) +rxt(k,346)* y(k,63) &
     949             :                   +rxt(k,347)* y(k,70) +rxt(k,348)* y(k,73) + rxt(k,98) &
     950   114070872 :                   + het_rates(k,134))* y(k,134)
     951             :          prod(k,134) = (rxt(k,117)*y(k,94) +rxt(k,287)*y(k,111) +rxt(k,331)*y(k,130)) &
     952   114070872 :                  *y(k,40) + (rxt(k,354)*y(k,41) +rxt(k,355)*y(k,139))*y(k,135)
     953             :          loss(k,135) = (rxt(k,350)* y(k,25) +rxt(k,351)* y(k,28) +rxt(k,352)* y(k,32) &
     954             :                   + (rxt(k,353) +rxt(k,354))* y(k,41) +rxt(k,356)* y(k,49) +rxt(k,357) &
     955             :                  * y(k,51) +rxt(k,136)* y(k,62) +rxt(k,359)* y(k,63) +rxt(k,360) &
     956             :                  * y(k,70) +rxt(k,362)* y(k,72) +rxt(k,363)* y(k,73) +rxt(k,361) &
     957             :                  * y(k,125) +rxt(k,355)* y(k,139) + rxt(k,99) + rxt(k,358) &
     958   114070872 :                   + het_rates(k,135))* y(k,135)
     959             :          prod(k,135) =rxt(k,325)*y(k,127)*y(k,70) +rxt(k,134)*y(k,98)*y(k,73) &
     960   114070872 :                   +rxt(k,63)*y(k,94) +rxt(k,95)*y(k,130)
     961             :          loss(k,86) = (rxt(k,594)* y(k,32) +rxt(k,597)* y(k,72) +rxt(k,596)* y(k,109) &
     962   114070872 :                   + rxt(k,595) + het_rates(k,136))* y(k,136)
     963             :          prod(k,86) = (rxt(k,77) +rxt(k,81) +rxt(k,587)*y(k,110) + &
     964             :                  rxt(k,588)*y(k,124) +rxt(k,600)*y(k,137) +rxt(k,606)*y(k,138)) &
     965             :                  *y(k,70) + (rxt(k,85) +rxt(k,87))*y(k,72) + (rxt(k,598)*y(k,137) + &
     966             :                  rxt(k,603)*y(k,138))*y(k,98) +rxt(k,580)*y(k,137) +rxt(k,579) &
     967   114070872 :                  *y(k,138)
     968             :          loss(k,62) = (rxt(k,600)* y(k,70) +rxt(k,601)* y(k,72) +rxt(k,598)* y(k,98) &
     969   114070872 :                   + rxt(k,580) + rxt(k,599) + het_rates(k,137))* y(k,137)
     970             :          prod(k,62) = (rxt(k,78) +rxt(k,82))*y(k,70) + (rxt(k,83) +rxt(k,92))*y(k,72) &
     971   114070872 :                   + (rxt(k,581) +rxt(k,602)*y(k,98))*y(k,138)
     972             :          loss(k,61) = (rxt(k,606)* y(k,70) + (rxt(k,602) +rxt(k,603))* y(k,98) &
     973             :                   + rxt(k,579) + rxt(k,581) + rxt(k,604) + rxt(k,605) &
     974   114070872 :                   + het_rates(k,138))* y(k,138)
     975   114070872 :          prod(k,61) = (rxt(k,79) +rxt(k,80))*y(k,70) + (rxt(k,90) +rxt(k,91))*y(k,72)
     976             :          loss(k,139) = (rxt(k,509)* y(k,39) +rxt(k,390)* y(k,71) +rxt(k,557)* y(k,82) &
     977             :                   +rxt(k,107)* y(k,90) +rxt(k,119)* y(k,94) +rxt(k,118)* y(k,96) &
     978             :                   +rxt(k,139)* y(k,99) +rxt(k,374)* y(k,102) +rxt(k,376)* y(k,103) &
     979             :                   +rxt(k,140)* y(k,104) +rxt(k,378)* y(k,105) +rxt(k,143)* y(k,106) &
     980             :                   +rxt(k,372)* y(k,108) +rxt(k,288)* y(k,111) +rxt(k,296)* y(k,113) &
     981             :                   +rxt(k,294)* y(k,115) +rxt(k,365)* y(k,118) +rxt(k,367)* y(k,119) &
     982             :                   +rxt(k,305)* y(k,120) +rxt(k,307)* y(k,121) +rxt(k,366)* y(k,122) &
     983             :                   +rxt(k,313)* y(k,123) +rxt(k,329)* y(k,128) + (rxt(k,327) + &
     984             :                  rxt(k,328))* y(k,129) +rxt(k,337)* y(k,132) +rxt(k,355)* y(k,135) &
     985   114070872 :                   + rxt(k,1) + rxt(k,2) + rxt(k,3) + het_rates(k,139))* y(k,139)
     986             :          prod(k,139) = (rxt(k,377) +4.000_r8*rxt(k,141)*y(k,98) + &
     987             :                  4.000_r8*rxt(k,146)*y(k,117) +4.000_r8*rxt(k,147)*y(k,90) + &
     988             :                  5.000_r8*rxt(k,158)*y(k,96) +5.000_r8*rxt(k,162)*y(k,112) + &
     989             :                  4.000_r8*rxt(k,163)*y(k,116) +5.000_r8*rxt(k,164)*y(k,91) + &
     990             :                  6.000_r8*rxt(k,165)*y(k,114) +4.000_r8*rxt(k,166)*y(k,111) + &
     991             :                  4.000_r8*rxt(k,170)*y(k,94) +4.000_r8*rxt(k,181)*y(k,92) + &
     992             :                  4.000_r8*rxt(k,192)*y(k,113) +4.000_r8*rxt(k,203)*y(k,100) + &
     993             :                  4.000_r8*rxt(k,214)*y(k,127) +4.000_r8*rxt(k,225)*y(k,97) + &
     994             :                  5.000_r8*rxt(k,236)*y(k,115) +6.000_r8*rxt(k,247)*y(k,95) + &
     995             :                  4.000_r8*rxt(k,258)*y(k,94) +5.000_r8*rxt(k,260)*y(k,96) + &
     996             :                  5.000_r8*rxt(k,262)*y(k,115) +4.000_r8*rxt(k,264)*y(k,113) + &
     997             :                  4.000_r8*rxt(k,267)*y(k,92) +4.000_r8*rxt(k,269)*y(k,117) + &
     998             :                  6.000_r8*rxt(k,271)*y(k,95))*y(k,105) + (rxt(k,379) + &
     999             :                  5.000_r8*rxt(k,144)*y(k,98) +5.000_r8*rxt(k,167)*y(k,117) + &
    1000             :                  5.000_r8*rxt(k,168)*y(k,94) +5.000_r8*rxt(k,169)*y(k,92) + &
    1001             :                  5.000_r8*rxt(k,171)*y(k,113) +5.000_r8*rxt(k,172)*y(k,100) + &
    1002             :                  5.000_r8*rxt(k,173)*y(k,127) +5.000_r8*rxt(k,174)*y(k,97) + &
    1003             :                  6.000_r8*rxt(k,175)*y(k,115) +7.000_r8*rxt(k,176)*y(k,95) + &
    1004             :                  5.000_r8*rxt(k,177)*y(k,90) +6.000_r8*rxt(k,178)*y(k,96) + &
    1005             :                  6.000_r8*rxt(k,179)*y(k,112) +5.000_r8*rxt(k,180)*y(k,116) + &
    1006             :                  6.000_r8*rxt(k,182)*y(k,91) +7.000_r8*rxt(k,183)*y(k,114) + &
    1007             :                  5.000_r8*rxt(k,184)*y(k,111) +7.000_r8*rxt(k,259)*y(k,95) + &
    1008             :                  6.000_r8*rxt(k,261)*y(k,96) +6.000_r8*rxt(k,263)*y(k,115) + &
    1009             :                  5.000_r8*rxt(k,265)*y(k,94) +5.000_r8*rxt(k,266)*y(k,113) + &
    1010             :                  5.000_r8*rxt(k,268)*y(k,92) +5.000_r8*rxt(k,270)*y(k,117))*y(k,107) &
    1011             :                   + (rxt(k,375) +3.000_r8*rxt(k,185)*y(k,117) + &
    1012             :                  3.000_r8*rxt(k,186)*y(k,94) +3.000_r8*rxt(k,187)*y(k,92) + &
    1013             :                  3.000_r8*rxt(k,188)*y(k,113) +3.000_r8*rxt(k,189)*y(k,100) + &
    1014             :                  3.000_r8*rxt(k,190)*y(k,127) +3.000_r8*rxt(k,191)*y(k,97) + &
    1015             :                  4.000_r8*rxt(k,193)*y(k,115) +5.000_r8*rxt(k,194)*y(k,95) + &
    1016             :                  3.000_r8*rxt(k,195)*y(k,90) +4.000_r8*rxt(k,196)*y(k,96) + &
    1017             :                  4.000_r8*rxt(k,197)*y(k,112) +3.000_r8*rxt(k,198)*y(k,116) + &
    1018             :                  4.000_r8*rxt(k,199)*y(k,91) +5.000_r8*rxt(k,200)*y(k,114) + &
    1019             :                  3.000_r8*rxt(k,201)*y(k,111) +3.000_r8*rxt(k,383)*y(k,98))*y(k,103) &
    1020             :                   + (rxt(k,340)*y(k,55) +rxt(k,415)*y(k,41) +rxt(k,416)*y(k,43) + &
    1021             :                  rxt(k,417)*y(k,101) +rxt(k,420)*y(k,133) +rxt(k,423)*y(k,52) + &
    1022             :                  rxt(k,447)*y(k,51) +rxt(k,471)*y(k,49) +rxt(k,474)*y(k,54) + &
    1023             :                  rxt(k,500)*y(k,45) +rxt(k,514)*y(k,16) +rxt(k,516)*y(k,18) + &
    1024             :                  rxt(k,517)*y(k,19) +rxt(k,519)*y(k,20) +rxt(k,524)*y(k,48) + &
    1025             :                  rxt(k,532)*y(k,17) +rxt(k,535)*y(k,22) +rxt(k,536)*y(k,23))*y(k,133) &
    1026             :                   + (2.000_r8*rxt(k,220)*y(k,117) +2.000_r8*rxt(k,221)*y(k,94) + &
    1027             :                  2.000_r8*rxt(k,222)*y(k,92) +2.000_r8*rxt(k,223)*y(k,113) + &
    1028             :                  2.000_r8*rxt(k,224)*y(k,100) +2.000_r8*rxt(k,226)*y(k,127) + &
    1029             :                  2.000_r8*rxt(k,227)*y(k,97) +3.000_r8*rxt(k,228)*y(k,115) + &
    1030             :                  4.000_r8*rxt(k,229)*y(k,95) +2.000_r8*rxt(k,230)*y(k,90) + &
    1031             :                  3.000_r8*rxt(k,231)*y(k,96) +3.000_r8*rxt(k,232)*y(k,112) + &
    1032             :                  2.000_r8*rxt(k,233)*y(k,116) +3.000_r8*rxt(k,234)*y(k,91) + &
    1033             :                  4.000_r8*rxt(k,235)*y(k,114) +2.000_r8*rxt(k,237)*y(k,111) + &
    1034             :                  2.000_r8*rxt(k,303)*y(k,98))*y(k,119) + (rxt(k,202)*y(k,117) + &
    1035             :                  rxt(k,204)*y(k,94) +rxt(k,205)*y(k,92) +rxt(k,206)*y(k,113) + &
    1036             :                  rxt(k,207)*y(k,100) +rxt(k,208)*y(k,127) +rxt(k,209)*y(k,97) + &
    1037             :                  2.000_r8*rxt(k,210)*y(k,115) +3.000_r8*rxt(k,211)*y(k,95) + &
    1038             :                  rxt(k,212)*y(k,90) +2.000_r8*rxt(k,213)*y(k,96) + &
    1039             :                  2.000_r8*rxt(k,215)*y(k,112) +rxt(k,216)*y(k,116) + &
    1040             :                  2.000_r8*rxt(k,217)*y(k,91) +3.000_r8*rxt(k,218)*y(k,114) + &
    1041             :                  rxt(k,219)*y(k,111) +rxt(k,308)*y(k,98))*y(k,122) &
    1042             :                   + (rxt(k,106)*y(k,91) +rxt(k,345)*y(k,134) +rxt(k,563)*y(k,54) + &
    1043             :                  rxt(k,569)*y(k,54) +rxt(k,570)*y(k,53) +rxt(k,574)*y(k,54) + &
    1044             :                  rxt(k,575)*y(k,53))*y(k,49) + (rxt(k,64) +rxt(k,283) + &
    1045             :                  rxt(k,120)*y(k,63) +rxt(k,123)*y(k,62) +rxt(k,155)*y(k,128) + &
    1046             :                  rxt(k,249)*y(k,118))*y(k,96) + (rxt(k,138)*y(k,99) + &
    1047             :                  3.000_r8*rxt(k,304)*y(k,120) +rxt(k,326)*y(k,129) + &
    1048             :                  rxt(k,381)*y(k,108) +2.000_r8*rxt(k,382)*y(k,102))*y(k,98) &
    1049             :                   + (rxt(k,245)*y(k,115) +2.000_r8*rxt(k,246)*y(k,95) + &
    1050             :                  rxt(k,250)*y(k,112) +rxt(k,252)*y(k,91) + &
    1051             :                  2.000_r8*rxt(k,253)*y(k,114))*y(k,118) + (rxt(k,152)*y(k,115) + &
    1052             :                  2.000_r8*rxt(k,153)*y(k,95) +rxt(k,156)*y(k,112) + &
    1053             :                  rxt(k,159)*y(k,91) +2.000_r8*rxt(k,160)*y(k,114))*y(k,128) &
    1054             :                   + (rxt(k,344)*y(k,134) +rxt(k,410)*y(k,101))*y(k,40) + (rxt(k,274) + &
    1055             :                  rxt(k,293)*y(k,58))*y(k,114) + (rxt(k,282) +rxt(k,295)*y(k,51)) &
    1056             :                  *y(k,115) +.050_r8*rxt(k,25)*y(k,23) +rxt(k,353)*y(k,135)*y(k,41) &
    1057             :                   +rxt(k,100)*y(k,44) +rxt(k,385)*y(k,91) +rxt(k,284)*y(k,95) &
    1058   117656469 :                   +rxt(k,373)*y(k,102) +rxt(k,273)*y(k,112) +rxt(k,94)*y(k,129)
    1059             :       end do
    1060     3585597 :       end subroutine imp_prod_loss
    1061             :       end module mo_prod_loss

Generated by: LCOV version 1.14