Line data Source code
1 : module mo_prod_loss
2 : use shr_kind_mod, only : r8 => shr_kind_r8
3 : use chem_mods, only : veclen
4 : private
5 : public :: exp_prod_loss
6 : public :: imp_prod_loss
7 : contains
8 72960 : subroutine exp_prod_loss( ofl, ofu, prod, loss, y, &
9 72960 : rxt, het_rates, chnkpnts )
10 : use chem_mods, only : gas_pcnst,rxntot,clscnt1
11 : implicit none
12 : !--------------------------------------------------------------------
13 : ! ... dummy args
14 : !--------------------------------------------------------------------
15 : integer, intent(in) :: ofl, ofu, chnkpnts
16 : real(r8), dimension(chnkpnts,max(1,clscnt1)), intent(out) :: &
17 : prod, &
18 : loss
19 : real(r8), intent(in) :: y(chnkpnts,gas_pcnst)
20 : real(r8), intent(in) :: rxt(chnkpnts,rxntot)
21 : real(r8), intent(in) :: het_rates(chnkpnts,gas_pcnst)
22 : !--------------------------------------------------------------------
23 : ! ... local variables
24 : !--------------------------------------------------------------------
25 : integer :: k
26 : !--------------------------------------------------------------------
27 : ! ... loss and production for Explicit method
28 : !--------------------------------------------------------------------
29 73616640 : do k = ofl,ofu
30 73543680 : loss(k,1) = ( + het_rates(k,189))* y(k,189)
31 73543680 : prod(k,1) = 0._r8
32 73543680 : loss(k,2) = ( + het_rates(k,190))* y(k,190)
33 73616640 : prod(k,2) = 0._r8
34 : end do
35 72960 : end subroutine exp_prod_loss
36 9949925 : subroutine imp_prod_loss( avec_len, prod, loss, y, &
37 : rxt, het_rates )
38 : use chem_mods, only : gas_pcnst,rxntot,clscnt4
39 : implicit none
40 : !--------------------------------------------------------------------
41 : ! ... dummy args
42 : !--------------------------------------------------------------------
43 : integer, intent(in) :: avec_len
44 : real(r8), dimension(veclen,clscnt4), intent(out) :: &
45 : prod, &
46 : loss
47 : real(r8), intent(in) :: y(veclen,gas_pcnst)
48 : real(r8), intent(in) :: rxt(veclen,rxntot)
49 : real(r8), intent(in) :: het_rates(veclen,gas_pcnst)
50 : !--------------------------------------------------------------------
51 : ! ... local variables
52 : !--------------------------------------------------------------------
53 : integer :: k
54 : !--------------------------------------------------------------------
55 : ! ... loss and production for Implicit method
56 : !--------------------------------------------------------------------
57 325174377 : do k = 1,avec_len
58 315224452 : loss(k,155) = (rxt(k,396)* y(k,228) + rxt(k,20) + het_rates(k,1))* y(k,1)
59 315224452 : prod(k,155) =rxt(k,399)*y(k,192)*y(k,124)
60 315224452 : loss(k,159) = (rxt(k,400)* y(k,228) + rxt(k,21) + het_rates(k,2))* y(k,2)
61 315224452 : prod(k,159) =rxt(k,397)*y(k,192)*y(k,90)
62 315224452 : loss(k,1) = ( + het_rates(k,3))* y(k,3)
63 315224452 : prod(k,1) = 0._r8
64 315224452 : loss(k,2) = ( + het_rates(k,4))* y(k,4)
65 315224452 : prod(k,2) = 0._r8
66 315224452 : loss(k,3) = ( + het_rates(k,5))* y(k,5)
67 315224452 : prod(k,3) = 0._r8
68 : loss(k,188) = (rxt(k,479)* y(k,126) +rxt(k,480)* y(k,136) +rxt(k,481) &
69 315224452 : * y(k,228) + het_rates(k,6))* y(k,6)
70 315224452 : prod(k,188) = 0._r8
71 315224452 : loss(k,71) = (rxt(k,438)* y(k,228) + het_rates(k,7))* y(k,7)
72 315224452 : prod(k,71) = 0._r8
73 315224452 : loss(k,122) = (rxt(k,441)* y(k,228) + rxt(k,22) + het_rates(k,8))* y(k,8)
74 315224452 : prod(k,122) =rxt(k,439)*y(k,194)*y(k,90)
75 315224452 : loss(k,72) = ( + rxt(k,23) + het_rates(k,9))* y(k,9)
76 315224452 : prod(k,72) =.120_r8*rxt(k,438)*y(k,228)*y(k,7)
77 315224452 : loss(k,116) = ( + rxt(k,24) + het_rates(k,10))* y(k,10)
78 : prod(k,116) = (.100_r8*rxt(k,480)*y(k,6) +.100_r8*rxt(k,483)*y(k,110)) &
79 315224452 : *y(k,136)
80 315224452 : loss(k,130) = ( + rxt(k,25) + het_rates(k,11))* y(k,11)
81 : prod(k,130) = (.500_r8*rxt(k,440)*y(k,194) +.200_r8*rxt(k,467)*y(k,235) + &
82 : .060_r8*rxt(k,473)*y(k,238))*y(k,124) +.500_r8*rxt(k,22)*y(k,8) &
83 : +rxt(k,23)*y(k,9) +.200_r8*rxt(k,71)*y(k,182) +.060_r8*rxt(k,73) &
84 315224452 : *y(k,186)
85 315224452 : loss(k,98) = ( + rxt(k,26) + het_rates(k,12))* y(k,12)
86 : prod(k,98) = (.200_r8*rxt(k,467)*y(k,235) +.200_r8*rxt(k,473)*y(k,238)) &
87 315224452 : *y(k,124) +.200_r8*rxt(k,71)*y(k,182) +.200_r8*rxt(k,73)*y(k,186)
88 315224452 : loss(k,149) = ( + rxt(k,27) + het_rates(k,13))* y(k,13)
89 : prod(k,149) = (.200_r8*rxt(k,467)*y(k,235) +.150_r8*rxt(k,473)*y(k,238)) &
90 : *y(k,124) +rxt(k,47)*y(k,95) +rxt(k,57)*y(k,116) +.200_r8*rxt(k,71) &
91 315224452 : *y(k,182) +.150_r8*rxt(k,73)*y(k,186)
92 315224452 : loss(k,108) = ( + rxt(k,28) + het_rates(k,14))* y(k,14)
93 315224452 : prod(k,108) =.210_r8*rxt(k,473)*y(k,238)*y(k,124) +.210_r8*rxt(k,73)*y(k,186)
94 315224452 : loss(k,85) = (rxt(k,401)* y(k,228) + het_rates(k,15))* y(k,15)
95 : prod(k,85) = (.050_r8*rxt(k,480)*y(k,6) +.050_r8*rxt(k,483)*y(k,110)) &
96 315224452 : *y(k,136)
97 : loss(k,113) = (rxt(k,367)* y(k,126) +rxt(k,368)* y(k,228) + het_rates(k,16)) &
98 315224452 : * y(k,16)
99 315224452 : prod(k,113) = 0._r8
100 : loss(k,221) = (rxt(k,250)* y(k,42) +rxt(k,251)* y(k,90) +rxt(k,252)* y(k,136) &
101 315224452 : + het_rates(k,17))* y(k,17)
102 : prod(k,221) = (rxt(k,76) +2.000_r8*rxt(k,253)*y(k,19) +rxt(k,254)*y(k,59) + &
103 : rxt(k,255)*y(k,59) +rxt(k,258)*y(k,124) +rxt(k,261)*y(k,134) + &
104 : rxt(k,262)*y(k,228) +rxt(k,509)*y(k,152))*y(k,19) &
105 : + (rxt(k,240)*y(k,34) +rxt(k,266)*y(k,35) + &
106 : 3.000_r8*rxt(k,267)*y(k,55) +2.000_r8*rxt(k,268)*y(k,78) + &
107 : rxt(k,269)*y(k,81) +2.000_r8*rxt(k,289)*y(k,41) +rxt(k,290)*y(k,43)) &
108 : *y(k,224) + (rxt(k,264)*y(k,81) +2.000_r8*rxt(k,278)*y(k,41) + &
109 : rxt(k,280)*y(k,43) +3.000_r8*rxt(k,285)*y(k,55))*y(k,228) &
110 : + (2.000_r8*rxt(k,277)*y(k,41) +rxt(k,279)*y(k,43) + &
111 : 3.000_r8*rxt(k,284)*y(k,55))*y(k,56) + (rxt(k,100) + &
112 : rxt(k,263)*y(k,134))*y(k,81) +rxt(k,75)*y(k,18) +rxt(k,78)*y(k,20) &
113 : +rxt(k,80)*y(k,34) +rxt(k,81)*y(k,35) +2.000_r8*rxt(k,87)*y(k,41) &
114 : +rxt(k,88)*y(k,43) +3.000_r8*rxt(k,91)*y(k,55) +2.000_r8*rxt(k,99) &
115 315224452 : *y(k,78) +rxt(k,106)*y(k,92)
116 315224452 : loss(k,86) = ( + rxt(k,75) + het_rates(k,18))* y(k,18)
117 : prod(k,86) = (rxt(k,570)*y(k,92) +rxt(k,575)*y(k,92))*y(k,85) &
118 315224452 : +rxt(k,256)*y(k,59)*y(k,19)
119 : loss(k,223) = (2._r8*rxt(k,253)* y(k,19) + (rxt(k,254) +rxt(k,255) + &
120 : rxt(k,256))* y(k,59) +rxt(k,257)* y(k,90) +rxt(k,258)* y(k,124) &
121 : +rxt(k,259)* y(k,125) +rxt(k,261)* y(k,134) +rxt(k,509)* y(k,152) &
122 315224452 : +rxt(k,262)* y(k,228) + rxt(k,76) + het_rates(k,19))* y(k,19)
123 : prod(k,223) = (rxt(k,77) +rxt(k,260)*y(k,134))*y(k,20) +rxt(k,252)*y(k,136) &
124 315224452 : *y(k,17) +rxt(k,270)*y(k,224)*y(k,81) +rxt(k,265)*y(k,134)*y(k,92)
125 : loss(k,145) = (rxt(k,260)* y(k,134) + rxt(k,77) + rxt(k,78) + rxt(k,564) &
126 315224452 : + rxt(k,567) + rxt(k,572) + het_rates(k,20))* y(k,20)
127 315224452 : prod(k,145) =rxt(k,259)*y(k,125)*y(k,19)
128 315224452 : loss(k,4) = ( + het_rates(k,21))* y(k,21)
129 315224452 : prod(k,4) = 0._r8
130 315224452 : loss(k,88) = (rxt(k,442)* y(k,228) + het_rates(k,22))* y(k,22)
131 315224452 : prod(k,88) =rxt(k,29)*y(k,23) +rxt(k,445)*y(k,196)*y(k,124)
132 315224452 : loss(k,111) = (rxt(k,444)* y(k,228) + rxt(k,29) + het_rates(k,23))* y(k,23)
133 315224452 : prod(k,111) =rxt(k,443)*y(k,196)*y(k,90)
134 : loss(k,100) = (rxt(k,315)* y(k,56) +rxt(k,316)* y(k,228) + het_rates(k,24)) &
135 315224452 : * y(k,24)
136 315224452 : prod(k,100) = 0._r8
137 : loss(k,143) = (rxt(k,317)* y(k,56) +rxt(k,318)* y(k,136) +rxt(k,343) &
138 315224452 : * y(k,228) + het_rates(k,25))* y(k,25)
139 315224452 : prod(k,143) = 0._r8
140 315224452 : loss(k,95) = (rxt(k,323)* y(k,228) + het_rates(k,26))* y(k,26)
141 : prod(k,95) = (.400_r8*rxt(k,319)*y(k,197) +.200_r8*rxt(k,320)*y(k,201)) &
142 315224452 : *y(k,197)
143 315224452 : loss(k,112) = (rxt(k,324)* y(k,228) + rxt(k,30) + het_rates(k,27))* y(k,27)
144 315224452 : prod(k,112) =rxt(k,321)*y(k,197)*y(k,90)
145 : loss(k,102) = (rxt(k,325)* y(k,56) +rxt(k,326)* y(k,228) + het_rates(k,28)) &
146 315224452 : * y(k,28)
147 315224452 : prod(k,102) = 0._r8
148 : loss(k,199) = (rxt(k,346)* y(k,126) +rxt(k,347)* y(k,136) +rxt(k,365) &
149 315224452 : * y(k,228) + het_rates(k,29))* y(k,29)
150 315224452 : prod(k,199) =.130_r8*rxt(k,425)*y(k,136)*y(k,99) +.700_r8*rxt(k,56)*y(k,111)
151 315224452 : loss(k,121) = (rxt(k,351)* y(k,228) + rxt(k,31) + het_rates(k,30))* y(k,30)
152 315224452 : prod(k,121) =rxt(k,349)*y(k,198)*y(k,90)
153 : loss(k,103) = (rxt(k,355)* y(k,56) +rxt(k,352)* y(k,228) + het_rates(k,31)) &
154 315224452 : * y(k,31)
155 315224452 : prod(k,103) = 0._r8
156 315224452 : loss(k,96) = (rxt(k,448)* y(k,228) + rxt(k,32) + het_rates(k,32))* y(k,32)
157 315224452 : prod(k,96) =rxt(k,446)*y(k,199)*y(k,90)
158 315224452 : loss(k,56) = (rxt(k,239)* y(k,224) + rxt(k,79) + het_rates(k,33))* y(k,33)
159 315224452 : prod(k,56) = 0._r8
160 315224452 : loss(k,67) = (rxt(k,240)* y(k,224) + rxt(k,80) + het_rates(k,34))* y(k,34)
161 315224452 : prod(k,67) = 0._r8
162 315224452 : loss(k,68) = (rxt(k,266)* y(k,224) + rxt(k,81) + het_rates(k,35))* y(k,35)
163 315224452 : prod(k,68) = 0._r8
164 315224452 : loss(k,59) = (rxt(k,241)* y(k,224) + rxt(k,82) + het_rates(k,36))* y(k,36)
165 315224452 : prod(k,59) = 0._r8
166 315224452 : loss(k,69) = (rxt(k,242)* y(k,224) + rxt(k,83) + het_rates(k,37))* y(k,37)
167 315224452 : prod(k,69) = 0._r8
168 315224452 : loss(k,60) = (rxt(k,243)* y(k,224) + rxt(k,84) + het_rates(k,38))* y(k,38)
169 315224452 : prod(k,60) = 0._r8
170 315224452 : loss(k,70) = (rxt(k,244)* y(k,224) + rxt(k,85) + het_rates(k,39))* y(k,39)
171 315224452 : prod(k,70) = 0._r8
172 315224452 : loss(k,61) = (rxt(k,245)* y(k,224) + rxt(k,86) + het_rates(k,40))* y(k,40)
173 315224452 : prod(k,61) = 0._r8
174 : loss(k,133) = (rxt(k,277)* y(k,56) +rxt(k,289)* y(k,224) +rxt(k,278) &
175 315224452 : * y(k,228) + rxt(k,87) + het_rates(k,41))* y(k,41)
176 315224452 : prod(k,133) = 0._r8
177 : loss(k,229) = (rxt(k,250)* y(k,17) +rxt(k,214)* y(k,56) +rxt(k,294)* y(k,90) &
178 : +rxt(k,295)* y(k,126) +rxt(k,296)* y(k,134) +rxt(k,297)* y(k,228) &
179 315224452 : + rxt(k,33) + rxt(k,34) + het_rates(k,42))* y(k,42)
180 : prod(k,229) = (rxt(k,221)*y(k,59) +2.000_r8*rxt(k,298)*y(k,201) + &
181 : rxt(k,299)*y(k,201) +rxt(k,301)*y(k,124) + &
182 : .700_r8*rxt(k,320)*y(k,197) +rxt(k,331)*y(k,200) + &
183 : rxt(k,348)*y(k,198) +.800_r8*rxt(k,361)*y(k,232) + &
184 : .880_r8*rxt(k,373)*y(k,213) +2.000_r8*rxt(k,382)*y(k,215) + &
185 : 1.500_r8*rxt(k,406)*y(k,208) +.750_r8*rxt(k,411)*y(k,209) + &
186 : .800_r8*rxt(k,420)*y(k,210) +.800_r8*rxt(k,431)*y(k,237) + &
187 : .750_r8*rxt(k,485)*y(k,223) +.930_r8*rxt(k,490)*y(k,233) + &
188 : .950_r8*rxt(k,495)*y(k,234))*y(k,201) &
189 : + (.500_r8*rxt(k,337)*y(k,206) +rxt(k,359)*y(k,231) + &
190 : rxt(k,363)*y(k,232) +.500_r8*rxt(k,369)*y(k,204) + &
191 : .250_r8*rxt(k,376)*y(k,213) +rxt(k,385)*y(k,215) + &
192 : .100_r8*rxt(k,398)*y(k,192) +.920_r8*rxt(k,408)*y(k,208) + &
193 : .250_r8*rxt(k,433)*y(k,237) +.340_r8*rxt(k,492)*y(k,233) + &
194 : .320_r8*rxt(k,497)*y(k,234))*y(k,124) + (rxt(k,302)*y(k,52) + &
195 : .300_r8*rxt(k,303)*y(k,53) +.500_r8*rxt(k,335)*y(k,51) + &
196 : .800_r8*rxt(k,340)*y(k,74) +rxt(k,342)*y(k,141) + &
197 : .500_r8*rxt(k,391)*y(k,109) +.400_r8*rxt(k,396)*y(k,1) + &
198 : .300_r8*rxt(k,416)*y(k,100) +.680_r8*rxt(k,501)*y(k,181))*y(k,228) &
199 : + (rxt(k,318)*y(k,25) +.500_r8*rxt(k,347)*y(k,29) + &
200 : .120_r8*rxt(k,378)*y(k,105) +.600_r8*rxt(k,392)*y(k,111) + &
201 : .910_r8*rxt(k,425)*y(k,99) +.340_r8*rxt(k,480)*y(k,6) + &
202 : .340_r8*rxt(k,483)*y(k,110))*y(k,136) + (.500_r8*rxt(k,367)*y(k,16) + &
203 : .250_r8*rxt(k,375)*y(k,213) +rxt(k,386)*y(k,215) + &
204 : rxt(k,409)*y(k,208))*y(k,126) + (.250_r8*rxt(k,372)*y(k,213) + &
205 : rxt(k,381)*y(k,215) +rxt(k,405)*y(k,208) + &
206 : .250_r8*rxt(k,430)*y(k,237))*y(k,200) + (.180_r8*rxt(k,40) + &
207 : rxt(k,311)*y(k,224) +rxt(k,312)*y(k,224))*y(k,54) &
208 : + (.150_r8*rxt(k,362)*y(k,232) +.450_r8*rxt(k,383)*y(k,215))*y(k,90) &
209 : +.100_r8*rxt(k,20)*y(k,1) +.100_r8*rxt(k,21)*y(k,2) +rxt(k,39) &
210 : *y(k,53) +rxt(k,44)*y(k,74) +.330_r8*rxt(k,46)*y(k,94) +rxt(k,48) &
211 : *y(k,96) +rxt(k,50)*y(k,103) +1.340_r8*rxt(k,51)*y(k,105) +rxt(k,58) &
212 : *y(k,127) +rxt(k,63)*y(k,148) +rxt(k,64)*y(k,149) +.375_r8*rxt(k,66) &
213 : *y(k,177) +.400_r8*rxt(k,68)*y(k,179) +.680_r8*rxt(k,70)*y(k,181) &
214 : +2.000_r8*rxt(k,338)*y(k,205) +rxt(k,308)*y(k,207) &
215 315224452 : +2.000_r8*rxt(k,384)*y(k,215)*y(k,215)
216 : loss(k,150) = (rxt(k,279)* y(k,56) +rxt(k,290)* y(k,224) +rxt(k,280) &
217 315224452 : * y(k,228) + rxt(k,88) + het_rates(k,43))* y(k,43)
218 315224452 : prod(k,150) = 0._r8
219 315224452 : loss(k,62) = (rxt(k,281)* y(k,228) + rxt(k,89) + het_rates(k,44))* y(k,44)
220 315224452 : prod(k,62) = 0._r8
221 : loss(k,200) = (rxt(k,327)* y(k,126) +rxt(k,328)* y(k,228) + rxt(k,35) &
222 315224452 : + het_rates(k,45))* y(k,45)
223 : prod(k,200) = (rxt(k,322)*y(k,197) +.270_r8*rxt(k,350)*y(k,198) + &
224 : rxt(k,359)*y(k,231) +rxt(k,369)*y(k,204) +rxt(k,388)*y(k,217) + &
225 : .400_r8*rxt(k,398)*y(k,192))*y(k,124) + (rxt(k,323)*y(k,26) + &
226 : .500_r8*rxt(k,324)*y(k,27) +.800_r8*rxt(k,396)*y(k,1))*y(k,228) &
227 : + (.500_r8*rxt(k,347)*y(k,29) +.100_r8*rxt(k,392)*y(k,111))*y(k,136) &
228 : + (1.600_r8*rxt(k,319)*y(k,197) +.800_r8*rxt(k,320)*y(k,201)) &
229 : *y(k,197) +.400_r8*rxt(k,20)*y(k,1) +.400_r8*rxt(k,21)*y(k,2) &
230 : +rxt(k,367)*y(k,126)*y(k,16) +rxt(k,30)*y(k,27) &
231 : +.200_r8*rxt(k,387)*y(k,217)*y(k,90) +.330_r8*rxt(k,46)*y(k,94) &
232 315224452 : +rxt(k,54)*y(k,108) +rxt(k,63)*y(k,148)
233 : loss(k,115) = (rxt(k,282)* y(k,56) +rxt(k,283)* y(k,228) + rxt(k,90) &
234 315224452 : + het_rates(k,46))* y(k,46)
235 315224452 : prod(k,115) = 0._r8
236 315224452 : loss(k,57) = (rxt(k,329)* y(k,228) + het_rates(k,47))* y(k,47)
237 315224452 : prod(k,57) = 0._r8
238 315224452 : loss(k,195) = (rxt(k,366)* y(k,228) + rxt(k,36) + het_rates(k,48))* y(k,48)
239 : prod(k,195) = (.820_r8*rxt(k,350)*y(k,198) +.500_r8*rxt(k,369)*y(k,204) + &
240 : .250_r8*rxt(k,398)*y(k,192) +.270_r8*rxt(k,492)*y(k,233) + &
241 : .040_r8*rxt(k,497)*y(k,234))*y(k,124) &
242 : + (.820_r8*rxt(k,348)*y(k,198) +.150_r8*rxt(k,490)*y(k,233) + &
243 : .025_r8*rxt(k,495)*y(k,234))*y(k,201) + (.250_r8*rxt(k,20) + &
244 : .800_r8*rxt(k,396)*y(k,228))*y(k,1) + (.520_r8*rxt(k,480)*y(k,6) + &
245 : .520_r8*rxt(k,483)*y(k,110))*y(k,136) + (.500_r8*rxt(k,70) + &
246 : .500_r8*rxt(k,501)*y(k,228))*y(k,181) +.250_r8*rxt(k,21)*y(k,2) &
247 : +.500_r8*rxt(k,367)*y(k,126)*y(k,16) +.820_r8*rxt(k,31)*y(k,30) &
248 : +.170_r8*rxt(k,46)*y(k,94) +.300_r8*rxt(k,66)*y(k,177) &
249 315224452 : +.050_r8*rxt(k,68)*y(k,179)
250 : loss(k,210) = (rxt(k,353)* y(k,126) +rxt(k,354)* y(k,228) + rxt(k,37) &
251 315224452 : + het_rates(k,49))* y(k,49)
252 : prod(k,210) = (.250_r8*rxt(k,376)*y(k,213) +.050_r8*rxt(k,414)*y(k,209) + &
253 : .250_r8*rxt(k,433)*y(k,237) +.170_r8*rxt(k,451)*y(k,202) + &
254 : .170_r8*rxt(k,457)*y(k,216) +.400_r8*rxt(k,467)*y(k,235) + &
255 : .540_r8*rxt(k,473)*y(k,238) +.510_r8*rxt(k,476)*y(k,240))*y(k,124) &
256 : + (.250_r8*rxt(k,375)*y(k,213) +.050_r8*rxt(k,415)*y(k,209) + &
257 : .250_r8*rxt(k,434)*y(k,237))*y(k,126) &
258 : + (.500_r8*rxt(k,361)*y(k,232) +.240_r8*rxt(k,373)*y(k,213) + &
259 : .100_r8*rxt(k,431)*y(k,237))*y(k,201) &
260 : + (.070_r8*rxt(k,450)*y(k,202) +.070_r8*rxt(k,456)*y(k,216))*y(k,90) &
261 : + (.880_r8*rxt(k,378)*y(k,105) +.500_r8*rxt(k,392)*y(k,111)) &
262 : *y(k,136) + (.250_r8*rxt(k,372)*y(k,213) + &
263 : .250_r8*rxt(k,430)*y(k,237))*y(k,200) + (rxt(k,356)*y(k,96) + &
264 : rxt(k,357)*y(k,127))*y(k,228) +.180_r8*rxt(k,24)*y(k,10) +rxt(k,28) &
265 : *y(k,14) +.400_r8*rxt(k,71)*y(k,182) +.540_r8*rxt(k,73)*y(k,186) &
266 315224452 : +.510_r8*rxt(k,74)*y(k,188)
267 315224452 : loss(k,161) = (rxt(k,334)* y(k,228) + het_rates(k,50))* y(k,50)
268 : prod(k,161) = (.150_r8*rxt(k,332)*y(k,200) +.150_r8*rxt(k,383)*y(k,215)) &
269 : *y(k,90) +.120_r8*rxt(k,347)*y(k,136)*y(k,29) &
270 315224452 : +.100_r8*rxt(k,331)*y(k,201)*y(k,200)
271 315224452 : loss(k,152) = (rxt(k,335)* y(k,228) + rxt(k,38) + het_rates(k,51))* y(k,51)
272 : prod(k,152) = (.400_r8*rxt(k,332)*y(k,200) +.400_r8*rxt(k,383)*y(k,215)) &
273 315224452 : *y(k,90)
274 315224452 : loss(k,173) = (rxt(k,302)* y(k,228) + het_rates(k,52))* y(k,52)
275 : prod(k,173) = (rxt(k,299)*y(k,201) +.300_r8*rxt(k,320)*y(k,197) + &
276 : .500_r8*rxt(k,361)*y(k,232) +.250_r8*rxt(k,373)*y(k,213) + &
277 : .250_r8*rxt(k,406)*y(k,208) +.250_r8*rxt(k,411)*y(k,209) + &
278 : .200_r8*rxt(k,420)*y(k,210) +.300_r8*rxt(k,431)*y(k,237) + &
279 : .250_r8*rxt(k,485)*y(k,223) +.250_r8*rxt(k,490)*y(k,233) + &
280 315224452 : .250_r8*rxt(k,495)*y(k,234))*y(k,201)
281 315224452 : loss(k,128) = (rxt(k,303)* y(k,228) + rxt(k,39) + het_rates(k,53))* y(k,53)
282 315224452 : prod(k,128) =rxt(k,300)*y(k,201)*y(k,90)
283 : loss(k,236) = (rxt(k,215)* y(k,56) +rxt(k,271)* y(k,73) + (rxt(k,310) + &
284 : rxt(k,311) +rxt(k,312))* y(k,224) +rxt(k,304)* y(k,228) + rxt(k,40) &
285 315224452 : + rxt(k,41) + het_rates(k,54))* y(k,54)
286 315224452 : prod(k,236) =.100_r8*rxt(k,347)*y(k,136)*y(k,29)
287 : loss(k,117) = (rxt(k,284)* y(k,56) +rxt(k,267)* y(k,224) +rxt(k,285) &
288 315224452 : * y(k,228) + rxt(k,91) + het_rates(k,55))* y(k,55)
289 315224452 : prod(k,117) = 0._r8
290 : loss(k,234) = (rxt(k,325)* y(k,28) +rxt(k,355)* y(k,31) +rxt(k,277)* y(k,41) &
291 : +rxt(k,214)* y(k,42) +rxt(k,279)* y(k,43) +rxt(k,282)* y(k,46) &
292 : +rxt(k,215)* y(k,54) +rxt(k,284)* y(k,55) +rxt(k,227)* y(k,60) &
293 : +rxt(k,216)* y(k,77) +rxt(k,217)* y(k,79) + (rxt(k,218) +rxt(k,219)) &
294 : * y(k,90) +rxt(k,236)* y(k,93) +rxt(k,220)* y(k,136) &
295 315224452 : + het_rates(k,56))* y(k,56)
296 : prod(k,234) = (4.000_r8*rxt(k,239)*y(k,33) +rxt(k,240)*y(k,34) + &
297 : 2.000_r8*rxt(k,241)*y(k,36) +2.000_r8*rxt(k,242)*y(k,37) + &
298 : 2.000_r8*rxt(k,243)*y(k,38) +rxt(k,244)*y(k,39) + &
299 : 2.000_r8*rxt(k,245)*y(k,40) +rxt(k,246)*y(k,85) +rxt(k,276)*y(k,65) + &
300 : rxt(k,291)*y(k,82) +rxt(k,292)*y(k,83) +rxt(k,293)*y(k,84))*y(k,224) &
301 : + (rxt(k,94) +rxt(k,221)*y(k,201) +2.000_r8*rxt(k,222)*y(k,59) + &
302 : rxt(k,224)*y(k,59) +rxt(k,226)*y(k,124) +rxt(k,231)*y(k,134) + &
303 : rxt(k,232)*y(k,228) +rxt(k,255)*y(k,19) +rxt(k,510)*y(k,152))*y(k,59) &
304 : + (rxt(k,235)*y(k,85) +3.000_r8*rxt(k,281)*y(k,44) + &
305 : rxt(k,283)*y(k,46) +rxt(k,286)*y(k,82) +rxt(k,287)*y(k,83) + &
306 : rxt(k,288)*y(k,84))*y(k,228) + (rxt(k,104) +rxt(k,234)*y(k,134)) &
307 : *y(k,85) +rxt(k,75)*y(k,18) +4.000_r8*rxt(k,79)*y(k,33) +rxt(k,80) &
308 : *y(k,34) +2.000_r8*rxt(k,82)*y(k,36) +2.000_r8*rxt(k,83)*y(k,37) &
309 : +2.000_r8*rxt(k,84)*y(k,38) +rxt(k,85)*y(k,39) +2.000_r8*rxt(k,86) &
310 : *y(k,40) +3.000_r8*rxt(k,89)*y(k,44) +rxt(k,90)*y(k,46) &
311 : +2.000_r8*rxt(k,92)*y(k,57) +2.000_r8*rxt(k,93)*y(k,58) +rxt(k,95) &
312 : *y(k,60) +rxt(k,98)*y(k,65) +rxt(k,101)*y(k,82) +rxt(k,102)*y(k,83) &
313 315224452 : +rxt(k,103)*y(k,84) +rxt(k,107)*y(k,93)
314 315224452 : loss(k,74) = ( + rxt(k,92) + het_rates(k,57))* y(k,57)
315 : prod(k,74) = (rxt(k,563)*y(k,93) +rxt(k,568)*y(k,60) +rxt(k,569)*y(k,93) + &
316 : rxt(k,573)*y(k,60) +rxt(k,574)*y(k,93) +rxt(k,578)*y(k,60))*y(k,85) &
317 315224452 : +rxt(k,227)*y(k,60)*y(k,56) +rxt(k,223)*y(k,59)*y(k,59)
318 315224452 : loss(k,54) = ( + rxt(k,93) + rxt(k,249) + het_rates(k,58))* y(k,58)
319 315224452 : prod(k,54) =rxt(k,248)*y(k,59)*y(k,59)
320 : loss(k,224) = ((rxt(k,254) +rxt(k,255) +rxt(k,256))* y(k,19) &
321 : + 2._r8*(rxt(k,222) +rxt(k,223) +rxt(k,224) +rxt(k,248))* y(k,59) &
322 : +rxt(k,225)* y(k,90) +rxt(k,226)* y(k,124) +rxt(k,228)* y(k,125) &
323 : +rxt(k,231)* y(k,134) +rxt(k,510)* y(k,152) +rxt(k,221)* y(k,201) &
324 : + (rxt(k,232) +rxt(k,233))* y(k,228) + rxt(k,94) + het_rates(k,59)) &
325 315224452 : * y(k,59)
326 : prod(k,224) = (rxt(k,219)*y(k,90) +rxt(k,220)*y(k,136) +rxt(k,236)*y(k,93)) &
327 : *y(k,56) + (rxt(k,96) +rxt(k,229)*y(k,134))*y(k,60) &
328 : + (rxt(k,237)*y(k,134) +rxt(k,238)*y(k,228))*y(k,93) + (rxt(k,108) + &
329 : rxt(k,515)*y(k,152))*y(k,138) +2.000_r8*rxt(k,249)*y(k,58) &
330 315224452 : +rxt(k,247)*y(k,224)*y(k,85)
331 : loss(k,187) = (rxt(k,227)* y(k,56) + (rxt(k,568) +rxt(k,573) +rxt(k,578)) &
332 : * y(k,85) +rxt(k,229)* y(k,134) +rxt(k,230)* y(k,228) + rxt(k,95) &
333 : + rxt(k,96) + rxt(k,566) + rxt(k,571) + rxt(k,577) &
334 315224452 : + het_rates(k,60))* y(k,60)
335 315224452 : prod(k,187) =rxt(k,228)*y(k,125)*y(k,59)
336 315224452 : loss(k,5) = ( + het_rates(k,61))* y(k,61)
337 315224452 : prod(k,5) = 0._r8
338 315224452 : loss(k,201) = (rxt(k,314)* y(k,228) + het_rates(k,62))* y(k,62)
339 : prod(k,201) = (rxt(k,33) +rxt(k,34) +rxt(k,214)*y(k,56) +rxt(k,250)*y(k,17) + &
340 : rxt(k,295)*y(k,126) +rxt(k,296)*y(k,134) +rxt(k,297)*y(k,228)) &
341 : *y(k,42) + (.630_r8*rxt(k,318)*y(k,25) +.560_r8*rxt(k,347)*y(k,29) + &
342 : .650_r8*rxt(k,378)*y(k,105) +.560_r8*rxt(k,392)*y(k,111) + &
343 : .620_r8*rxt(k,425)*y(k,99) +.230_r8*rxt(k,480)*y(k,6) + &
344 : .230_r8*rxt(k,483)*y(k,110))*y(k,136) &
345 : + (.220_r8*rxt(k,376)*y(k,213) +.250_r8*rxt(k,433)*y(k,237) + &
346 : .170_r8*rxt(k,451)*y(k,202) +.400_r8*rxt(k,454)*y(k,214) + &
347 : .350_r8*rxt(k,457)*y(k,216) +.225_r8*rxt(k,492)*y(k,233))*y(k,124) &
348 : + (.350_r8*rxt(k,316)*y(k,24) +rxt(k,341)*y(k,75) + &
349 : rxt(k,354)*y(k,49) +.700_r8*rxt(k,501)*y(k,181) +rxt(k,505)*y(k,139)) &
350 : *y(k,228) + (rxt(k,42) +rxt(k,110) +rxt(k,591)*y(k,229))*y(k,63) &
351 : + (.070_r8*rxt(k,450)*y(k,202) +.160_r8*rxt(k,453)*y(k,214) + &
352 : .140_r8*rxt(k,456)*y(k,216))*y(k,90) + (rxt(k,353)*y(k,49) + &
353 : .220_r8*rxt(k,375)*y(k,213) +.500_r8*rxt(k,434)*y(k,237))*y(k,126) &
354 : + (.110_r8*rxt(k,373)*y(k,213) +.200_r8*rxt(k,431)*y(k,237) + &
355 : .125_r8*rxt(k,490)*y(k,233))*y(k,201) + (rxt(k,137) + &
356 : rxt(k,504)*y(k,134))*y(k,139) + (.220_r8*rxt(k,372)*y(k,213) + &
357 : .250_r8*rxt(k,430)*y(k,237))*y(k,200) +1.500_r8*rxt(k,23)*y(k,9) &
358 : +.450_r8*rxt(k,24)*y(k,10) +.600_r8*rxt(k,27)*y(k,13) +rxt(k,28) &
359 : *y(k,14) +rxt(k,35)*y(k,45) +rxt(k,282)*y(k,56)*y(k,46) +rxt(k,37) &
360 : *y(k,49) +.380_r8*rxt(k,40)*y(k,54) +rxt(k,44)*y(k,74) &
361 : +2.000_r8*rxt(k,45)*y(k,75) +.330_r8*rxt(k,46)*y(k,94) &
362 : +1.340_r8*rxt(k,52)*y(k,105) +.700_r8*rxt(k,56)*y(k,111) &
363 : +1.500_r8*rxt(k,65)*y(k,176) +.250_r8*rxt(k,66)*y(k,177) +rxt(k,69) &
364 315224452 : *y(k,180) +1.700_r8*rxt(k,70)*y(k,181)
365 : loss(k,217) = (rxt(k,591)* y(k,229) + rxt(k,42) + rxt(k,110) &
366 315224452 : + het_rates(k,63))* y(k,63)
367 : prod(k,217) = (rxt(k,306)*y(k,87) +rxt(k,314)*y(k,62) +rxt(k,334)*y(k,50) + &
368 : .500_r8*rxt(k,335)*y(k,51) +.800_r8*rxt(k,340)*y(k,74) + &
369 : rxt(k,341)*y(k,75) +.500_r8*rxt(k,391)*y(k,109) + &
370 : 1.800_r8*rxt(k,501)*y(k,181))*y(k,228) &
371 : + (2.000_r8*rxt(k,330)*y(k,200) +.900_r8*rxt(k,331)*y(k,201) + &
372 : rxt(k,333)*y(k,124) +2.000_r8*rxt(k,381)*y(k,215) + &
373 : rxt(k,405)*y(k,208) +rxt(k,430)*y(k,237))*y(k,200) &
374 : + (.200_r8*rxt(k,347)*y(k,29) +.100_r8*rxt(k,392)*y(k,111) + &
375 : .270_r8*rxt(k,480)*y(k,6) +.270_r8*rxt(k,483)*y(k,110))*y(k,136) &
376 : + (rxt(k,382)*y(k,201) +.450_r8*rxt(k,383)*y(k,90) + &
377 : 2.000_r8*rxt(k,384)*y(k,215))*y(k,215) &
378 : + (.500_r8*rxt(k,490)*y(k,201) +.900_r8*rxt(k,492)*y(k,124)) &
379 : *y(k,233) +rxt(k,38)*y(k,51) +.440_r8*rxt(k,40)*y(k,54) &
380 : +.400_r8*rxt(k,61)*y(k,141) +rxt(k,66)*y(k,177) +.800_r8*rxt(k,70) &
381 315224452 : *y(k,181)
382 315224452 : loss(k,89) = (rxt(k,275)* y(k,224) + rxt(k,97) + het_rates(k,64))* y(k,64)
383 : prod(k,89) = (rxt(k,240)*y(k,34) +rxt(k,242)*y(k,37) + &
384 : 2.000_r8*rxt(k,243)*y(k,38) +2.000_r8*rxt(k,244)*y(k,39) + &
385 : rxt(k,245)*y(k,40) +rxt(k,266)*y(k,35) +2.000_r8*rxt(k,268)*y(k,78) + &
386 : rxt(k,292)*y(k,83) +rxt(k,293)*y(k,84))*y(k,224) + (rxt(k,102) + &
387 : rxt(k,287)*y(k,228))*y(k,83) + (rxt(k,103) +rxt(k,288)*y(k,228)) &
388 : *y(k,84) +rxt(k,80)*y(k,34) +rxt(k,81)*y(k,35) +rxt(k,83)*y(k,37) &
389 : +2.000_r8*rxt(k,84)*y(k,38) +2.000_r8*rxt(k,85)*y(k,39) +rxt(k,86) &
390 315224452 : *y(k,40) +2.000_r8*rxt(k,99)*y(k,78)
391 315224452 : loss(k,91) = (rxt(k,276)* y(k,224) + rxt(k,98) + het_rates(k,65))* y(k,65)
392 : prod(k,91) = (rxt(k,101) +rxt(k,286)*y(k,228) +rxt(k,291)*y(k,224))*y(k,82) &
393 : + (rxt(k,82) +rxt(k,241)*y(k,224))*y(k,36) + (rxt(k,83) + &
394 315224452 : rxt(k,242)*y(k,224))*y(k,37)
395 315224452 : loss(k,79) = (rxt(k,449)* y(k,228) + het_rates(k,66))* y(k,66)
396 315224452 : prod(k,79) =.180_r8*rxt(k,469)*y(k,228)*y(k,183)
397 : loss(k,104) = (rxt(k,502)* y(k,126) + (rxt(k,503) +rxt(k,517))* y(k,228) &
398 315224452 : + het_rates(k,67))* y(k,67)
399 315224452 : prod(k,104) = 0._r8
400 315224452 : loss(k,6) = ( + het_rates(k,68))* y(k,68)
401 315224452 : prod(k,6) = 0._r8
402 315224452 : loss(k,7) = ( + het_rates(k,69))* y(k,69)
403 315224452 : prod(k,7) = 0._r8
404 315224452 : loss(k,8) = ( + het_rates(k,70))* y(k,70)
405 315224452 : prod(k,8) = 0._r8
406 315224452 : loss(k,9) = ( + rxt(k,594) + het_rates(k,71))* y(k,71)
407 315224452 : prod(k,9) = 0._r8
408 315224452 : loss(k,63) = ( + rxt(k,43) + het_rates(k,72))* y(k,72)
409 315224452 : prod(k,63) =rxt(k,336)*y(k,206)*y(k,90)
410 : loss(k,184) = (rxt(k,271)* y(k,54) +rxt(k,272)* y(k,77) +rxt(k,274)* y(k,89) &
411 315224452 : +rxt(k,273)* y(k,241) + het_rates(k,73))* y(k,73)
412 : prod(k,184) = (rxt(k,244)*y(k,39) +rxt(k,266)*y(k,35) + &
413 : 2.000_r8*rxt(k,275)*y(k,64) +rxt(k,276)*y(k,65))*y(k,224) +rxt(k,81) &
414 : *y(k,35) +rxt(k,85)*y(k,39) +2.000_r8*rxt(k,97)*y(k,64) +rxt(k,98) &
415 315224452 : *y(k,65) +rxt(k,105)*y(k,88)
416 315224452 : loss(k,202) = (rxt(k,340)* y(k,228) + rxt(k,44) + het_rates(k,74))* y(k,74)
417 : prod(k,202) = (.530_r8*rxt(k,376)*y(k,213) +.050_r8*rxt(k,414)*y(k,209) + &
418 : .250_r8*rxt(k,433)*y(k,237) +.225_r8*rxt(k,492)*y(k,233))*y(k,124) &
419 : + (.530_r8*rxt(k,375)*y(k,213) +.050_r8*rxt(k,415)*y(k,209) + &
420 : .250_r8*rxt(k,434)*y(k,237))*y(k,126) &
421 : + (.260_r8*rxt(k,373)*y(k,213) +.100_r8*rxt(k,431)*y(k,237) + &
422 : .125_r8*rxt(k,490)*y(k,233))*y(k,201) &
423 : + (.700_r8*rxt(k,416)*y(k,100) +.500_r8*rxt(k,417)*y(k,101) + &
424 : rxt(k,428)*y(k,115))*y(k,228) + (.530_r8*rxt(k,372)*y(k,213) + &
425 : .250_r8*rxt(k,430)*y(k,237))*y(k,200) +.330_r8*rxt(k,46)*y(k,94) &
426 315224452 : +rxt(k,339)*y(k,205)*y(k,135) +.250_r8*rxt(k,66)*y(k,177)
427 : loss(k,192) = (rxt(k,341)* y(k,228) + rxt(k,45) + rxt(k,520) &
428 315224452 : + het_rates(k,75))* y(k,75)
429 : prod(k,192) = (.050_r8*rxt(k,414)*y(k,209) +.250_r8*rxt(k,433)*y(k,237) + &
430 : rxt(k,440)*y(k,194) +.400_r8*rxt(k,454)*y(k,214) + &
431 : .170_r8*rxt(k,457)*y(k,216) +.700_r8*rxt(k,460)*y(k,230) + &
432 : .600_r8*rxt(k,467)*y(k,235) +.340_r8*rxt(k,473)*y(k,238) + &
433 : .170_r8*rxt(k,476)*y(k,240))*y(k,124) + (.650_r8*rxt(k,316)*y(k,24) + &
434 : .200_r8*rxt(k,340)*y(k,74) +rxt(k,429)*y(k,116))*y(k,228) &
435 : + (.250_r8*rxt(k,430)*y(k,200) +.100_r8*rxt(k,431)*y(k,201) + &
436 : .250_r8*rxt(k,434)*y(k,126))*y(k,237) &
437 : + (.160_r8*rxt(k,453)*y(k,214) +.070_r8*rxt(k,456)*y(k,216))*y(k,90) &
438 : +rxt(k,22)*y(k,8) +.130_r8*rxt(k,24)*y(k,10) &
439 : +.050_r8*rxt(k,415)*y(k,209)*y(k,126) +.700_r8*rxt(k,62)*y(k,145) &
440 : +.600_r8*rxt(k,71)*y(k,182) +.340_r8*rxt(k,73)*y(k,186) &
441 315224452 : +.170_r8*rxt(k,74)*y(k,188)
442 : loss(k,233) = ((rxt(k,172) +rxt(k,173) +rxt(k,174))* y(k,90) +rxt(k,175) &
443 315224452 : * y(k,135) +rxt(k,178)* y(k,136) + het_rates(k,76))* y(k,76)
444 : prod(k,233) = (rxt(k,179)*y(k,77) +rxt(k,182)*y(k,134) +rxt(k,202)*y(k,112) + &
445 : rxt(k,297)*y(k,42) +rxt(k,505)*y(k,139) +rxt(k,511)*y(k,150) + &
446 : rxt(k,516)*y(k,152))*y(k,228) + (rxt(k,153)*y(k,224) + &
447 : rxt(k,170)*y(k,134) +rxt(k,216)*y(k,56) +rxt(k,272)*y(k,73))*y(k,77) &
448 : + (.330_r8*rxt(k,40) +rxt(k,41) +rxt(k,311)*y(k,224))*y(k,54) &
449 : + (rxt(k,100) +rxt(k,270)*y(k,224))*y(k,81) + (rxt(k,104) + &
450 : rxt(k,247)*y(k,224))*y(k,85) + (2.000_r8*rxt(k,2) +rxt(k,3))*y(k,241) &
451 315224452 : +2.000_r8*rxt(k,34)*y(k,42) +rxt(k,39)*y(k,53) +rxt(k,105)*y(k,88)
452 : loss(k,218) = (rxt(k,216)* y(k,56) +rxt(k,272)* y(k,73) +rxt(k,170)* y(k,134) &
453 : +rxt(k,153)* y(k,224) +rxt(k,179)* y(k,228) + het_rates(k,77)) &
454 315224452 : * y(k,77)
455 : prod(k,218) = (1.440_r8*rxt(k,40) +rxt(k,312)*y(k,224))*y(k,54) +rxt(k,33) &
456 315224452 : *y(k,42) +rxt(k,172)*y(k,90)*y(k,76) +rxt(k,1)*y(k,241)
457 315224452 : loss(k,58) = (rxt(k,268)* y(k,224) + rxt(k,99) + het_rates(k,78))* y(k,78)
458 315224452 : prod(k,58) = 0._r8
459 : loss(k,147) = (rxt(k,217)* y(k,56) +rxt(k,171)* y(k,134) +rxt(k,180) &
460 315224452 : * y(k,228) + rxt(k,4) + het_rates(k,79))* y(k,79)
461 315224452 : prod(k,147) =rxt(k,186)*y(k,90)*y(k,90) +rxt(k,185)*y(k,228)*y(k,228)
462 315224452 : loss(k,64) = ( + rxt(k,136) + het_rates(k,80))* y(k,80)
463 315224452 : prod(k,64) =rxt(k,518)*y(k,241)*y(k,154)
464 : loss(k,174) = (rxt(k,263)* y(k,134) + (rxt(k,269) +rxt(k,270))* y(k,224) &
465 315224452 : +rxt(k,264)* y(k,228) + rxt(k,100) + het_rates(k,81))* y(k,81)
466 315224452 : prod(k,174) = (rxt(k,250)*y(k,42) +rxt(k,251)*y(k,90))*y(k,17)
467 : loss(k,90) = (rxt(k,291)* y(k,224) +rxt(k,286)* y(k,228) + rxt(k,101) &
468 315224452 : + het_rates(k,82))* y(k,82)
469 315224452 : prod(k,90) = 0._r8
470 : loss(k,92) = (rxt(k,292)* y(k,224) +rxt(k,287)* y(k,228) + rxt(k,102) &
471 315224452 : + het_rates(k,83))* y(k,83)
472 315224452 : prod(k,92) = 0._r8
473 : loss(k,101) = (rxt(k,293)* y(k,224) +rxt(k,288)* y(k,228) + rxt(k,103) &
474 315224452 : + het_rates(k,84))* y(k,84)
475 315224452 : prod(k,101) = 0._r8
476 : loss(k,220) = ((rxt(k,568) +rxt(k,573) +rxt(k,578))* y(k,60) + (rxt(k,570) + &
477 : rxt(k,575))* y(k,92) + (rxt(k,563) +rxt(k,569) +rxt(k,574))* y(k,93) &
478 : +rxt(k,234)* y(k,134) + (rxt(k,246) +rxt(k,247))* y(k,224) &
479 315224452 : +rxt(k,235)* y(k,228) + rxt(k,104) + het_rates(k,85))* y(k,85)
480 : prod(k,220) = (rxt(k,214)*y(k,42) +rxt(k,215)*y(k,54) +rxt(k,216)*y(k,77) + &
481 : rxt(k,217)*y(k,79) +rxt(k,218)*y(k,90) +rxt(k,236)*y(k,93) + &
482 : rxt(k,277)*y(k,41) +rxt(k,279)*y(k,43) +2.000_r8*rxt(k,282)*y(k,46) + &
483 : rxt(k,284)*y(k,55) +rxt(k,325)*y(k,28) +rxt(k,355)*y(k,31))*y(k,56) &
484 315224452 : +rxt(k,233)*y(k,228)*y(k,59)
485 : loss(k,77) = (rxt(k,313)* y(k,224) +rxt(k,305)* y(k,228) + het_rates(k,86)) &
486 315224452 : * y(k,86)
487 315224452 : prod(k,77) = 0._r8
488 315224452 : loss(k,183) = (rxt(k,306)* y(k,228) + het_rates(k,87))* y(k,87)
489 : prod(k,183) = (.370_r8*rxt(k,318)*y(k,25) +.120_r8*rxt(k,347)*y(k,29) + &
490 : .330_r8*rxt(k,378)*y(k,105) +.120_r8*rxt(k,392)*y(k,111) + &
491 : .110_r8*rxt(k,425)*y(k,99) +.050_r8*rxt(k,480)*y(k,6) + &
492 : .050_r8*rxt(k,483)*y(k,110))*y(k,136) + (rxt(k,307)*y(k,90) + &
493 315224452 : rxt(k,309)*y(k,124))*y(k,207) +.350_r8*rxt(k,316)*y(k,228)*y(k,24)
494 315224452 : loss(k,99) = ( + rxt(k,105) + het_rates(k,88))* y(k,88)
495 : prod(k,99) = (rxt(k,271)*y(k,54) +rxt(k,272)*y(k,77) +rxt(k,273)*y(k,241) + &
496 315224452 : rxt(k,274)*y(k,89))*y(k,73)
497 : loss(k,231) = (rxt(k,274)* y(k,73) +rxt(k,211)* y(k,228) + rxt(k,9) &
498 315224452 : + het_rates(k,89))* y(k,89)
499 : prod(k,231) = (rxt(k,566) +rxt(k,571) +rxt(k,577) +rxt(k,568)*y(k,85) + &
500 : rxt(k,573)*y(k,85) +rxt(k,578)*y(k,85))*y(k,60) + (rxt(k,530) + &
501 : rxt(k,295)*y(k,42) +rxt(k,327)*y(k,45) +rxt(k,353)*y(k,49) + &
502 : rxt(k,502)*y(k,67))*y(k,126) + (2.000_r8*rxt(k,525) + &
503 : 2.000_r8*rxt(k,562) +2.000_r8*rxt(k,565) +2.000_r8*rxt(k,576)) &
504 : *y(k,114) + (rxt(k,564) +rxt(k,567) +rxt(k,572))*y(k,20) &
505 : + (.500_r8*rxt(k,529) +rxt(k,210)*y(k,228))*y(k,125) +rxt(k,522) &
506 : *y(k,94) +rxt(k,523)*y(k,100) +rxt(k,524)*y(k,101) +rxt(k,526) &
507 : *y(k,115) +rxt(k,527)*y(k,116) +rxt(k,531)*y(k,128) +rxt(k,532) &
508 315224452 : *y(k,140) +rxt(k,533)*y(k,178)
509 : loss(k,237) = (rxt(k,251)* y(k,17) +rxt(k,257)* y(k,19) +rxt(k,294)* y(k,42) &
510 : + (rxt(k,218) +rxt(k,219))* y(k,56) +rxt(k,225)* y(k,59) &
511 : + (rxt(k,172) +rxt(k,173) +rxt(k,174))* y(k,76) + 2._r8*rxt(k,186) &
512 : * y(k,90) +rxt(k,203)* y(k,124) +rxt(k,208)* y(k,125) +rxt(k,198) &
513 : * y(k,126) +rxt(k,176)* y(k,134) +rxt(k,177)* y(k,136) +rxt(k,436) &
514 : * y(k,191) +rxt(k,397)* y(k,192) +rxt(k,439)* y(k,194) +rxt(k,443) &
515 : * y(k,196) +rxt(k,321)* y(k,197) +rxt(k,349)* y(k,198) +rxt(k,446) &
516 : * y(k,199) +rxt(k,332)* y(k,200) +rxt(k,300)* y(k,201) +rxt(k,450) &
517 : * y(k,202) +rxt(k,336)* y(k,206) +rxt(k,307)* y(k,207) +rxt(k,407) &
518 : * y(k,208) +rxt(k,412)* y(k,209) +rxt(k,421)* y(k,210) +rxt(k,374) &
519 : * y(k,213) +rxt(k,453)* y(k,214) +rxt(k,383)* y(k,215) +rxt(k,456) &
520 : * y(k,216) +rxt(k,387)* y(k,217) +rxt(k,486)* y(k,223) +rxt(k,181) &
521 : * y(k,228) +rxt(k,459)* y(k,230) +rxt(k,358)* y(k,231) +rxt(k,362) &
522 : * y(k,232) +rxt(k,491)* y(k,233) +rxt(k,496)* y(k,234) +rxt(k,466) &
523 : * y(k,235) +rxt(k,432)* y(k,237) +rxt(k,472)* y(k,238) +rxt(k,475) &
524 315224452 : * y(k,240) + rxt(k,521) + het_rates(k,90))* y(k,90)
525 : prod(k,237) = (rxt(k,180)*y(k,79) +rxt(k,183)*y(k,136) +rxt(k,201)*y(k,126) + &
526 : rxt(k,232)*y(k,59) +rxt(k,262)*y(k,19) +rxt(k,280)*y(k,43) + &
527 : rxt(k,283)*y(k,46) +rxt(k,302)*y(k,52) +rxt(k,305)*y(k,86) + &
528 : rxt(k,306)*y(k,87) +rxt(k,314)*y(k,62) +.350_r8*rxt(k,316)*y(k,24) + &
529 : rxt(k,323)*y(k,26) +rxt(k,329)*y(k,47) +rxt(k,340)*y(k,74) + &
530 : rxt(k,341)*y(k,75) +rxt(k,356)*y(k,96) +rxt(k,371)*y(k,94) + &
531 : .200_r8*rxt(k,380)*y(k,106) +.500_r8*rxt(k,391)*y(k,109) + &
532 : .300_r8*rxt(k,416)*y(k,100) +rxt(k,417)*y(k,101) + &
533 : rxt(k,424)*y(k,102) +rxt(k,428)*y(k,115) +rxt(k,429)*y(k,116) + &
534 : .650_r8*rxt(k,438)*y(k,7) +.730_r8*rxt(k,449)*y(k,66) + &
535 : .800_r8*rxt(k,461)*y(k,144) +.280_r8*rxt(k,469)*y(k,183) + &
536 : .380_r8*rxt(k,471)*y(k,185) +.630_r8*rxt(k,477)*y(k,187) + &
537 : .200_r8*rxt(k,501)*y(k,181) +rxt(k,507)*y(k,153) + &
538 : .500_r8*rxt(k,517)*y(k,67))*y(k,228) + (rxt(k,301)*y(k,201) + &
539 : rxt(k,309)*y(k,207) +rxt(k,322)*y(k,197) + &
540 : .250_r8*rxt(k,337)*y(k,206) +rxt(k,350)*y(k,198) + &
541 : rxt(k,359)*y(k,231) +rxt(k,369)*y(k,204) + &
542 : .470_r8*rxt(k,376)*y(k,213) +rxt(k,398)*y(k,192) + &
543 : .920_r8*rxt(k,408)*y(k,208) +.920_r8*rxt(k,414)*y(k,209) + &
544 : rxt(k,422)*y(k,210) +rxt(k,433)*y(k,237) +rxt(k,440)*y(k,194) + &
545 : rxt(k,445)*y(k,196) +.170_r8*rxt(k,451)*y(k,202) + &
546 : .400_r8*rxt(k,454)*y(k,214) +.830_r8*rxt(k,457)*y(k,216) + &
547 : rxt(k,460)*y(k,230) +rxt(k,467)*y(k,235) +rxt(k,473)*y(k,238) + &
548 : rxt(k,476)*y(k,240) +.900_r8*rxt(k,492)*y(k,233) + &
549 : .800_r8*rxt(k,497)*y(k,234))*y(k,124) + (rxt(k,221)*y(k,59) + &
550 : 2.000_r8*rxt(k,298)*y(k,201) +rxt(k,320)*y(k,197) + &
551 : .900_r8*rxt(k,331)*y(k,200) +rxt(k,348)*y(k,198) + &
552 : .300_r8*rxt(k,361)*y(k,232) +.730_r8*rxt(k,373)*y(k,213) + &
553 : rxt(k,382)*y(k,215) +rxt(k,406)*y(k,208) +rxt(k,411)*y(k,209) + &
554 : 1.200_r8*rxt(k,420)*y(k,210) +.800_r8*rxt(k,431)*y(k,237) + &
555 : .500_r8*rxt(k,485)*y(k,223) +rxt(k,490)*y(k,233) + &
556 : rxt(k,495)*y(k,234))*y(k,201) + (.130_r8*rxt(k,318)*y(k,25) + &
557 : .280_r8*rxt(k,347)*y(k,29) +.140_r8*rxt(k,378)*y(k,105) + &
558 : .280_r8*rxt(k,392)*y(k,111) +.370_r8*rxt(k,425)*y(k,99) + &
559 : .570_r8*rxt(k,480)*y(k,6) +.570_r8*rxt(k,483)*y(k,110))*y(k,136) &
560 : + (rxt(k,295)*y(k,42) +.470_r8*rxt(k,375)*y(k,213) + &
561 : rxt(k,409)*y(k,208) +rxt(k,415)*y(k,209) +rxt(k,423)*y(k,210) + &
562 : rxt(k,434)*y(k,237))*y(k,126) + (.470_r8*rxt(k,372)*y(k,213) + &
563 : rxt(k,405)*y(k,208) +rxt(k,410)*y(k,209) +rxt(k,419)*y(k,210) + &
564 : rxt(k,430)*y(k,237))*y(k,200) + (rxt(k,214)*y(k,42) + &
565 : rxt(k,217)*y(k,79) +rxt(k,279)*y(k,43) +rxt(k,282)*y(k,46))*y(k,56) &
566 : + (.070_r8*rxt(k,450)*y(k,202) +.160_r8*rxt(k,453)*y(k,214) + &
567 : .330_r8*rxt(k,456)*y(k,216))*y(k,90) + (rxt(k,250)*y(k,17) + &
568 : rxt(k,296)*y(k,134))*y(k,42) + (rxt(k,11) +rxt(k,212))*y(k,91) &
569 : + (1.340_r8*rxt(k,51) +.660_r8*rxt(k,52))*y(k,105) &
570 : + (rxt(k,175)*y(k,76) +rxt(k,339)*y(k,205))*y(k,135) +rxt(k,20) &
571 : *y(k,1) +.900_r8*rxt(k,21)*y(k,2) +rxt(k,22)*y(k,8) &
572 : +1.500_r8*rxt(k,23)*y(k,9) +.560_r8*rxt(k,24)*y(k,10) +rxt(k,25) &
573 : *y(k,11) +.600_r8*rxt(k,26)*y(k,12) +.600_r8*rxt(k,27)*y(k,13) &
574 : +rxt(k,28)*y(k,14) +rxt(k,29)*y(k,23) +rxt(k,30)*y(k,27) +rxt(k,31) &
575 : *y(k,30) +rxt(k,35)*y(k,45) +rxt(k,37)*y(k,49) +rxt(k,311)*y(k,224) &
576 : *y(k,54) +2.000_r8*rxt(k,44)*y(k,74) +2.000_r8*rxt(k,45)*y(k,75) &
577 : +rxt(k,171)*y(k,134)*y(k,79) +.670_r8*rxt(k,46)*y(k,94) +rxt(k,47) &
578 : *y(k,95) +rxt(k,48)*y(k,96) +rxt(k,49)*y(k,102) +rxt(k,50)*y(k,103) &
579 : +rxt(k,57)*y(k,116) +rxt(k,62)*y(k,145) +rxt(k,63)*y(k,148) &
580 : +rxt(k,65)*y(k,176) +rxt(k,66)*y(k,177) +rxt(k,67)*y(k,178) &
581 : +rxt(k,68)*y(k,179) +rxt(k,69)*y(k,180) +1.200_r8*rxt(k,70)*y(k,181) &
582 : +rxt(k,71)*y(k,182) +rxt(k,73)*y(k,186) +rxt(k,74)*y(k,188) &
583 : +1.200_r8*rxt(k,319)*y(k,197)*y(k,197) +rxt(k,338)*y(k,205) &
584 315224452 : +rxt(k,308)*y(k,207) +rxt(k,413)*y(k,209)
585 : loss(k,132) = (rxt(k,187)* y(k,228) + rxt(k,10) + rxt(k,11) + rxt(k,212) &
586 315224452 : + het_rates(k,91))* y(k,91)
587 315224452 : prod(k,132) =rxt(k,208)*y(k,125)*y(k,90)
588 : loss(k,169) = ((rxt(k,570) +rxt(k,575))* y(k,85) +rxt(k,265)* y(k,134) &
589 315224452 : + rxt(k,106) + het_rates(k,92))* y(k,92)
590 : prod(k,169) = (rxt(k,564) +rxt(k,567) +rxt(k,572))*y(k,20) &
591 315224452 : +rxt(k,257)*y(k,90)*y(k,19)
592 : loss(k,175) = (rxt(k,236)* y(k,56) + (rxt(k,563) +rxt(k,569) +rxt(k,574)) &
593 : * y(k,85) +rxt(k,237)* y(k,134) +rxt(k,238)* y(k,228) + rxt(k,107) &
594 315224452 : + het_rates(k,93))* y(k,93)
595 : prod(k,175) = (rxt(k,566) +rxt(k,571) +rxt(k,577) +rxt(k,230)*y(k,228)) &
596 315224452 : *y(k,60) +rxt(k,225)*y(k,90)*y(k,59)
597 : loss(k,204) = (rxt(k,371)* y(k,228) + rxt(k,46) + rxt(k,522) &
598 315224452 : + het_rates(k,94))* y(k,94)
599 : prod(k,204) = (rxt(k,370)*y(k,204) +rxt(k,377)*y(k,213))*y(k,124) &
600 : + (.300_r8*rxt(k,416)*y(k,100) +.500_r8*rxt(k,417)*y(k,101)) &
601 315224452 : *y(k,228)
602 315224452 : loss(k,84) = (rxt(k,402)* y(k,228) + rxt(k,47) + het_rates(k,95))* y(k,95)
603 315224452 : prod(k,84) =rxt(k,413)*y(k,209)
604 315224452 : loss(k,205) = (rxt(k,356)* y(k,228) + rxt(k,48) + het_rates(k,96))* y(k,96)
605 : prod(k,205) = (.220_r8*rxt(k,372)*y(k,200) +.230_r8*rxt(k,373)*y(k,201) + &
606 : .220_r8*rxt(k,375)*y(k,126) +.220_r8*rxt(k,376)*y(k,124))*y(k,213) &
607 : + (.500_r8*rxt(k,360)*y(k,148) +.500_r8*rxt(k,391)*y(k,109) + &
608 : .700_r8*rxt(k,416)*y(k,100) +.500_r8*rxt(k,417)*y(k,101))*y(k,228) &
609 : + (.250_r8*rxt(k,430)*y(k,200) +.100_r8*rxt(k,431)*y(k,201) + &
610 : .250_r8*rxt(k,433)*y(k,124) +.250_r8*rxt(k,434)*y(k,126))*y(k,237) &
611 : + (.050_r8*rxt(k,414)*y(k,124) +.050_r8*rxt(k,415)*y(k,126)) &
612 : *y(k,209) +.170_r8*rxt(k,46)*y(k,94) +.200_r8*rxt(k,361)*y(k,232) &
613 315224452 : *y(k,201)
614 315224452 : loss(k,107) = (rxt(k,403)* y(k,228) + het_rates(k,97))* y(k,97)
615 : prod(k,107) = (rxt(k,410)*y(k,200) +.750_r8*rxt(k,411)*y(k,201) + &
616 315224452 : .870_r8*rxt(k,414)*y(k,124) +.950_r8*rxt(k,415)*y(k,126))*y(k,209)
617 315224452 : loss(k,65) = (rxt(k,404)* y(k,228) + het_rates(k,98))* y(k,98)
618 315224452 : prod(k,65) =.600_r8*rxt(k,427)*y(k,228)*y(k,103)
619 : loss(k,179) = (rxt(k,418)* y(k,126) +rxt(k,425)* y(k,136) +rxt(k,426) &
620 315224452 : * y(k,228) + het_rates(k,99))* y(k,99)
621 315224452 : prod(k,179) = 0._r8
622 : loss(k,148) = (rxt(k,416)* y(k,228) + rxt(k,523) + het_rates(k,100)) &
623 315224452 : * y(k,100)
624 315224452 : prod(k,148) =.080_r8*rxt(k,408)*y(k,208)*y(k,124)
625 : loss(k,144) = (rxt(k,417)* y(k,228) + rxt(k,524) + het_rates(k,101)) &
626 315224452 : * y(k,101)
627 315224452 : prod(k,144) =.080_r8*rxt(k,414)*y(k,209)*y(k,124)
628 315224452 : loss(k,120) = (rxt(k,424)* y(k,228) + rxt(k,49) + het_rates(k,102))* y(k,102)
629 315224452 : prod(k,120) =rxt(k,421)*y(k,210)*y(k,90)
630 315224452 : loss(k,162) = (rxt(k,427)* y(k,228) + rxt(k,50) + het_rates(k,103))* y(k,103)
631 : prod(k,162) = (rxt(k,407)*y(k,208) +rxt(k,412)*y(k,209))*y(k,90) +rxt(k,49) &
632 315224452 : *y(k,102)
633 315224452 : loss(k,50) = (rxt(k,549)* y(k,228) + het_rates(k,104))* y(k,104)
634 315224452 : prod(k,50) = 0._r8
635 : loss(k,209) = (rxt(k,378)* y(k,136) +rxt(k,379)* y(k,228) + rxt(k,51) &
636 315224452 : + rxt(k,52) + het_rates(k,105))* y(k,105)
637 : prod(k,209) = (.390_r8*rxt(k,405)*y(k,200) +.310_r8*rxt(k,406)*y(k,201) + &
638 : .360_r8*rxt(k,408)*y(k,124) +.400_r8*rxt(k,409)*y(k,126))*y(k,208) &
639 315224452 : +.300_r8*rxt(k,425)*y(k,136)*y(k,99) +.300_r8*rxt(k,50)*y(k,103)
640 315224452 : loss(k,105) = (rxt(k,380)* y(k,228) + het_rates(k,106))* y(k,106)
641 315224452 : prod(k,105) =rxt(k,374)*y(k,213)*y(k,90)
642 315224452 : loss(k,138) = (rxt(k,389)* y(k,228) + rxt(k,53) + het_rates(k,107))* y(k,107)
643 : prod(k,138) =.800_r8*rxt(k,20)*y(k,1) +.800_r8*rxt(k,21)*y(k,2) &
644 315224452 : +.800_r8*rxt(k,398)*y(k,192)*y(k,124)
645 315224452 : loss(k,106) = (rxt(k,390)* y(k,228) + rxt(k,54) + het_rates(k,108))* y(k,108)
646 315224452 : prod(k,106) =.800_r8*rxt(k,387)*y(k,217)*y(k,90)
647 : loss(k,151) = (rxt(k,391)* y(k,228) + rxt(k,55) + rxt(k,395) &
648 315224452 : + het_rates(k,109))* y(k,109)
649 315224452 : prod(k,151) =rxt(k,394)*y(k,215)*y(k,125)
650 : loss(k,190) = (rxt(k,482)* y(k,126) +rxt(k,483)* y(k,136) +rxt(k,484) &
651 315224452 : * y(k,228) + het_rates(k,110))* y(k,110)
652 315224452 : prod(k,190) = 0._r8
653 : loss(k,214) = (rxt(k,392)* y(k,136) +rxt(k,393)* y(k,228) + rxt(k,56) &
654 315224452 : + het_rates(k,111))* y(k,111)
655 : prod(k,214) = (.610_r8*rxt(k,405)*y(k,200) +.440_r8*rxt(k,406)*y(k,201) + &
656 : .560_r8*rxt(k,408)*y(k,124) +.600_r8*rxt(k,409)*y(k,126))*y(k,208) &
657 315224452 : +.200_r8*rxt(k,425)*y(k,136)*y(k,99) +.700_r8*rxt(k,50)*y(k,103)
658 : loss(k,189) = (rxt(k,190)* y(k,124) + (rxt(k,191) +rxt(k,192) +rxt(k,193)) &
659 : * y(k,125) +rxt(k,194)* y(k,135) +rxt(k,588)* y(k,227) +rxt(k,202) &
660 315224452 : * y(k,228) + rxt(k,111) + het_rates(k,112))* y(k,112)
661 : prod(k,189) = (rxt(k,188)*y(k,219) +rxt(k,585)*y(k,222))*y(k,134) &
662 : + (.200_r8*rxt(k,579)*y(k,221) +1.100_r8*rxt(k,581)*y(k,220)) &
663 : *y(k,203) +rxt(k,15)*y(k,124) +rxt(k,586)*y(k,222)*y(k,135) &
664 315224452 : +rxt(k,592)*y(k,229)
665 : loss(k,93) = ((rxt(k,206) +rxt(k,207))* y(k,224) + rxt(k,12) &
666 315224452 : + het_rates(k,113))* y(k,113)
667 315224452 : prod(k,93) =rxt(k,191)*y(k,125)*y(k,112)
668 : loss(k,97) = ( + rxt(k,13) + rxt(k,14) + rxt(k,213) + rxt(k,525) + rxt(k,562) &
669 315224452 : + rxt(k,565) + rxt(k,576) + het_rates(k,114))* y(k,114)
670 315224452 : prod(k,97) =rxt(k,209)*y(k,126)*y(k,125)
671 : loss(k,124) = (rxt(k,428)* y(k,228) + rxt(k,526) + het_rates(k,115)) &
672 315224452 : * y(k,115)
673 315224452 : prod(k,124) =.200_r8*rxt(k,420)*y(k,210)*y(k,201)
674 : loss(k,196) = (rxt(k,429)* y(k,228) + rxt(k,57) + rxt(k,527) &
675 315224452 : + het_rates(k,116))* y(k,116)
676 : prod(k,196) = (rxt(k,419)*y(k,200) +.800_r8*rxt(k,420)*y(k,201) + &
677 315224452 : rxt(k,422)*y(k,124) +rxt(k,423)*y(k,126))*y(k,210)
678 315224452 : loss(k,10) = ( + het_rates(k,117))* y(k,117)
679 315224452 : prod(k,10) = 0._r8
680 315224452 : loss(k,11) = ( + het_rates(k,118))* y(k,118)
681 315224452 : prod(k,11) = 0._r8
682 315224452 : loss(k,12) = ( + het_rates(k,119))* y(k,119)
683 315224452 : prod(k,12) = 0._r8
684 315224452 : loss(k,55) = (rxt(k,519)* y(k,228) + het_rates(k,120))* y(k,120)
685 315224452 : prod(k,55) = 0._r8
686 315224452 : loss(k,13) = ( + rxt(k,528) + het_rates(k,121))* y(k,121)
687 315224452 : prod(k,13) = 0._r8
688 315224452 : loss(k,14) = ( + rxt(k,596) + het_rates(k,122))* y(k,122)
689 315224452 : prod(k,14) = 0._r8
690 315224452 : loss(k,15) = ( + rxt(k,595) + het_rates(k,123))* y(k,123)
691 315224452 : prod(k,15) = 0._r8
692 : loss(k,227) = (rxt(k,258)* y(k,19) +rxt(k,226)* y(k,59) +rxt(k,203)* y(k,90) &
693 : +rxt(k,190)* y(k,112) +rxt(k,199)* y(k,126) +rxt(k,205)* y(k,134) &
694 : +rxt(k,204)* y(k,136) +rxt(k,437)* y(k,191) + (rxt(k,398) + &
695 : rxt(k,399))* y(k,192) +rxt(k,440)* y(k,194) +rxt(k,445)* y(k,196) &
696 : +rxt(k,322)* y(k,197) +rxt(k,350)* y(k,198) +rxt(k,447)* y(k,199) &
697 : +rxt(k,333)* y(k,200) +rxt(k,301)* y(k,201) +rxt(k,451)* y(k,202) &
698 : + (rxt(k,369) +rxt(k,370))* y(k,204) +rxt(k,337)* y(k,206) &
699 : +rxt(k,309)* y(k,207) +rxt(k,408)* y(k,208) +rxt(k,414)* y(k,209) &
700 : +rxt(k,422)* y(k,210) + (rxt(k,376) +rxt(k,377))* y(k,213) &
701 : +rxt(k,454)* y(k,214) +rxt(k,385)* y(k,215) +rxt(k,457)* y(k,216) &
702 : +rxt(k,388)* y(k,217) +rxt(k,487)* y(k,223) +rxt(k,590)* y(k,227) &
703 : +rxt(k,460)* y(k,230) +rxt(k,359)* y(k,231) +rxt(k,363)* y(k,232) &
704 : +rxt(k,492)* y(k,233) +rxt(k,497)* y(k,234) +rxt(k,467)* y(k,235) &
705 : +rxt(k,433)* y(k,237) +rxt(k,473)* y(k,238) +rxt(k,476)* y(k,240) &
706 315224452 : + rxt(k,15) + rxt(k,16) + het_rates(k,124))* y(k,124)
707 : prod(k,227) = (rxt(k,17) +.500_r8*rxt(k,529) +2.000_r8*rxt(k,192)*y(k,112) + &
708 : rxt(k,195)*y(k,134) +rxt(k,512)*y(k,152))*y(k,125) &
709 : + (rxt(k,194)*y(k,135) +rxt(k,202)*y(k,228))*y(k,112) &
710 : +2.000_r8*rxt(k,206)*y(k,224)*y(k,113) +rxt(k,14)*y(k,114) &
711 : +rxt(k,19)*y(k,126) +rxt(k,189)*y(k,219)*y(k,135) +rxt(k,589) &
712 315224452 : *y(k,227)
713 : loss(k,228) = (rxt(k,259)* y(k,19) +rxt(k,228)* y(k,59) +rxt(k,208)* y(k,90) &
714 : + (rxt(k,191) +rxt(k,192) +rxt(k,193))* y(k,112) +rxt(k,209) &
715 : * y(k,126) + (rxt(k,195) +rxt(k,197))* y(k,134) +rxt(k,196)* y(k,136) &
716 : +rxt(k,462)* y(k,143) +rxt(k,512)* y(k,152) +rxt(k,465)* y(k,191) &
717 : +rxt(k,344)* y(k,200) +rxt(k,452)* y(k,202) +rxt(k,455)* y(k,214) &
718 : +rxt(k,394)* y(k,215) +rxt(k,458)* y(k,216) +rxt(k,210)* y(k,228) &
719 315224452 : + rxt(k,17) + rxt(k,529) + het_rates(k,125))* y(k,125)
720 : prod(k,228) = (2.000_r8*rxt(k,199)*y(k,126) +rxt(k,203)*y(k,90) + &
721 : rxt(k,204)*y(k,136) +rxt(k,205)*y(k,134) +rxt(k,226)*y(k,59) + &
722 : rxt(k,258)*y(k,19) +rxt(k,301)*y(k,201) +rxt(k,309)*y(k,207) + &
723 : rxt(k,322)*y(k,197) +rxt(k,333)*y(k,200) +rxt(k,337)*y(k,206) + &
724 : rxt(k,350)*y(k,198) +rxt(k,359)*y(k,231) +rxt(k,363)*y(k,232) + &
725 : rxt(k,369)*y(k,204) +rxt(k,376)*y(k,213) +rxt(k,385)*y(k,215) + &
726 : rxt(k,388)*y(k,217) +rxt(k,398)*y(k,192) + &
727 : .920_r8*rxt(k,408)*y(k,208) +.920_r8*rxt(k,414)*y(k,209) + &
728 : rxt(k,422)*y(k,210) +rxt(k,433)*y(k,237) +rxt(k,437)*y(k,191) + &
729 : rxt(k,440)*y(k,194) +rxt(k,445)*y(k,196) +rxt(k,447)*y(k,199) + &
730 : rxt(k,451)*y(k,202) +rxt(k,454)*y(k,214) +rxt(k,457)*y(k,216) + &
731 : rxt(k,460)*y(k,230) +rxt(k,467)*y(k,235) +rxt(k,473)*y(k,238) + &
732 : rxt(k,476)*y(k,240) +1.600_r8*rxt(k,487)*y(k,223) + &
733 : .900_r8*rxt(k,492)*y(k,233) +.800_r8*rxt(k,497)*y(k,234))*y(k,124) &
734 : + (rxt(k,18) +rxt(k,198)*y(k,90) +rxt(k,200)*y(k,134) + &
735 : rxt(k,201)*y(k,228) +rxt(k,367)*y(k,16) +rxt(k,375)*y(k,213) + &
736 : rxt(k,386)*y(k,215) +rxt(k,409)*y(k,208) +rxt(k,415)*y(k,209) + &
737 : rxt(k,423)*y(k,210) +rxt(k,434)*y(k,237) + &
738 : 2.000_r8*rxt(k,488)*y(k,223))*y(k,126) + (rxt(k,187)*y(k,91) + &
739 : rxt(k,357)*y(k,127) +rxt(k,396)*y(k,1) +.700_r8*rxt(k,416)*y(k,100) + &
740 : rxt(k,494)*y(k,178))*y(k,228) + (rxt(k,11) +rxt(k,212))*y(k,91) &
741 : + (rxt(k,55) +rxt(k,395))*y(k,109) + (rxt(k,13) +rxt(k,213)) &
742 : *y(k,114) + (.600_r8*rxt(k,61) +rxt(k,345))*y(k,141) +rxt(k,20) &
743 : *y(k,1) +rxt(k,77)*y(k,20) +rxt(k,96)*y(k,60) +rxt(k,9)*y(k,89) &
744 : +rxt(k,46)*y(k,94) +rxt(k,49)*y(k,102) +rxt(k,57)*y(k,116) &
745 : +rxt(k,58)*y(k,127) +rxt(k,59)*y(k,128) +rxt(k,60)*y(k,140) &
746 : +rxt(k,470)*y(k,142) +rxt(k,67)*y(k,178) &
747 315224452 : +.500_r8*rxt(k,485)*y(k,223)*y(k,201)
748 : loss(k,225) = (rxt(k,479)* y(k,6) +rxt(k,367)* y(k,16) +rxt(k,346)* y(k,29) &
749 : +rxt(k,295)* y(k,42) +rxt(k,327)* y(k,45) +rxt(k,353)* y(k,49) &
750 : +rxt(k,502)* y(k,67) +rxt(k,198)* y(k,90) +rxt(k,418)* y(k,99) &
751 : +rxt(k,482)* y(k,110) +rxt(k,199)* y(k,124) +rxt(k,209)* y(k,125) &
752 : +rxt(k,200)* y(k,134) +rxt(k,499)* y(k,180) +rxt(k,409)* y(k,208) &
753 : +rxt(k,415)* y(k,209) +rxt(k,423)* y(k,210) +rxt(k,375)* y(k,213) &
754 : +rxt(k,386)* y(k,215) +rxt(k,488)* y(k,223) +rxt(k,201)* y(k,228) &
755 : +rxt(k,434)* y(k,237) + rxt(k,18) + rxt(k,19) + rxt(k,530) &
756 315224452 : + het_rates(k,126))* y(k,126)
757 : prod(k,225) = (rxt(k,95) +rxt(k,227)*y(k,56) +rxt(k,229)*y(k,134) + &
758 : rxt(k,230)*y(k,228))*y(k,60) + (rxt(k,13) +rxt(k,14) +rxt(k,213)) &
759 : *y(k,114) + (rxt(k,211)*y(k,89) +rxt(k,342)*y(k,141) + &
760 : .500_r8*rxt(k,391)*y(k,109))*y(k,228) + (rxt(k,78) + &
761 : rxt(k,260)*y(k,134))*y(k,20) + (rxt(k,196)*y(k,136) + &
762 : rxt(k,197)*y(k,134))*y(k,125) +rxt(k,274)*y(k,89)*y(k,73) +rxt(k,10) &
763 315224452 : *y(k,91) +.400_r8*rxt(k,61)*y(k,141)
764 315224452 : loss(k,180) = (rxt(k,357)* y(k,228) + rxt(k,58) + het_rates(k,127))* y(k,127)
765 : prod(k,180) = (.500_r8*rxt(k,417)*y(k,101) +rxt(k,424)*y(k,102) + &
766 : rxt(k,428)*y(k,115) +rxt(k,429)*y(k,116))*y(k,228) &
767 315224452 : +rxt(k,346)*y(k,126)*y(k,29)
768 : loss(k,118) = (rxt(k,489)* y(k,228) + rxt(k,59) + rxt(k,531) &
769 315224452 : + het_rates(k,128))* y(k,128)
770 315224452 : prod(k,118) =rxt(k,486)*y(k,223)*y(k,90)
771 315224452 : loss(k,16) = ( + het_rates(k,129))* y(k,129)
772 315224452 : prod(k,16) = 0._r8
773 315224452 : loss(k,17) = ( + het_rates(k,130))* y(k,130)
774 315224452 : prod(k,17) = 0._r8
775 315224452 : loss(k,18) = ( + het_rates(k,131))* y(k,131)
776 315224452 : prod(k,18) = 0._r8
777 315224452 : loss(k,19) = ( + het_rates(k,132))* y(k,132)
778 315224452 : prod(k,19) = 0._r8
779 315224452 : loss(k,20) = ( + het_rates(k,133))* y(k,133)
780 315224452 : prod(k,20) = 0._r8
781 : loss(k,232) = (rxt(k,261)* y(k,19) +rxt(k,260)* y(k,20) +rxt(k,296)* y(k,42) &
782 : +rxt(k,231)* y(k,59) +rxt(k,229)* y(k,60) +rxt(k,170)* y(k,77) &
783 : +rxt(k,171)* y(k,79) +rxt(k,263)* y(k,81) +rxt(k,234)* y(k,85) &
784 : +rxt(k,176)* y(k,90) +rxt(k,265)* y(k,92) +rxt(k,237)* y(k,93) &
785 : +rxt(k,205)* y(k,124) + (rxt(k,195) +rxt(k,197))* y(k,125) &
786 : +rxt(k,200)* y(k,126) + 2._r8*rxt(k,168)* y(k,134) +rxt(k,169) &
787 : * y(k,135) +rxt(k,167)* y(k,136) +rxt(k,504)* y(k,139) &
788 : + (rxt(k,583) +rxt(k,584))* y(k,220) +rxt(k,585)* y(k,222) &
789 : +rxt(k,182)* y(k,228) + rxt(k,120) + rxt(k,121) + rxt(k,122) &
790 315224452 : + rxt(k,123) + rxt(k,124) + rxt(k,125) + het_rates(k,134))* y(k,134)
791 : prod(k,232) = (2.000_r8*rxt(k,5) +rxt(k,6) +rxt(k,127) +rxt(k,128) + &
792 : rxt(k,129) +rxt(k,131) +rxt(k,132) +rxt(k,133) +2.000_r8*rxt(k,134) + &
793 : 2.000_r8*rxt(k,135) +rxt(k,156)*y(k,224) +rxt(k,157)*y(k,224) + &
794 : rxt(k,194)*y(k,112) +rxt(k,506)*y(k,150) +rxt(k,513)*y(k,152) + &
795 : rxt(k,587)*y(k,222) +rxt(k,593)*y(k,229))*y(k,135) &
796 : + (rxt(k,190)*y(k,124) +rxt(k,191)*y(k,125) +rxt(k,588)*y(k,227)) &
797 : *y(k,112) + (rxt(k,42) +rxt(k,110))*y(k,63) + (rxt(k,579)*y(k,221) + &
798 : 1.150_r8*rxt(k,580)*y(k,227))*y(k,203) +rxt(k,76)*y(k,19) &
799 : +.180_r8*rxt(k,40)*y(k,54) +rxt(k,94)*y(k,59) +rxt(k,174)*y(k,90) &
800 : *y(k,76) +rxt(k,14)*y(k,114) +rxt(k,15)*y(k,124) +rxt(k,17)*y(k,125) &
801 : +rxt(k,18)*y(k,126) +rxt(k,8)*y(k,136) +rxt(k,108)*y(k,138) &
802 : +rxt(k,138)*y(k,152) +rxt(k,139)*y(k,153) +rxt(k,140)*y(k,154) &
803 : +rxt(k,155)*y(k,224) +rxt(k,184)*y(k,228)*y(k,228) +rxt(k,2) &
804 315224452 : *y(k,241)
805 : loss(k,222) = (rxt(k,175)* y(k,76) +rxt(k,194)* y(k,112) +rxt(k,169) &
806 : * y(k,134) +rxt(k,506)* y(k,150) +rxt(k,513)* y(k,152) +rxt(k,339) &
807 : * y(k,205) +rxt(k,189)* y(k,219) +rxt(k,582)* y(k,220) &
808 : + (rxt(k,586) +rxt(k,587))* y(k,222) +rxt(k,156)* y(k,224) &
809 : +rxt(k,161)* y(k,225) +rxt(k,593)* y(k,229) + rxt(k,5) + rxt(k,6) &
810 : + rxt(k,126) + rxt(k,127) + rxt(k,128) + rxt(k,129) + rxt(k,130) &
811 : + rxt(k,131) + rxt(k,132) + rxt(k,133) + rxt(k,134) + rxt(k,135) &
812 315224452 : + het_rates(k,135))* y(k,135)
813 : prod(k,222) = (rxt(k,172)*y(k,76) +rxt(k,176)*y(k,134) + &
814 : 2.000_r8*rxt(k,177)*y(k,136) +rxt(k,181)*y(k,228) + &
815 : rxt(k,186)*y(k,90) +rxt(k,198)*y(k,126) +rxt(k,218)*y(k,56) + &
816 : rxt(k,225)*y(k,59) +rxt(k,251)*y(k,17) +rxt(k,257)*y(k,19) + &
817 : rxt(k,300)*y(k,201) +rxt(k,321)*y(k,197) +rxt(k,349)*y(k,198) + &
818 : rxt(k,358)*y(k,231))*y(k,90) + (rxt(k,8) + &
819 : 2.000_r8*rxt(k,158)*y(k,224) +2.000_r8*rxt(k,167)*y(k,134) + &
820 : rxt(k,178)*y(k,76) +rxt(k,183)*y(k,228) +rxt(k,196)*y(k,125) + &
821 : rxt(k,204)*y(k,124) +rxt(k,220)*y(k,56) +rxt(k,252)*y(k,17) + &
822 : rxt(k,508)*y(k,150) +rxt(k,514)*y(k,152))*y(k,136) &
823 : + (rxt(k,160)*y(k,225) +rxt(k,168)*y(k,134) +rxt(k,182)*y(k,228) + &
824 : rxt(k,195)*y(k,125) +rxt(k,200)*y(k,126) +rxt(k,231)*y(k,59) + &
825 : rxt(k,261)*y(k,19))*y(k,134) + (rxt(k,222)*y(k,59) + &
826 : rxt(k,223)*y(k,59) +rxt(k,233)*y(k,228) +rxt(k,255)*y(k,19) + &
827 : rxt(k,256)*y(k,19))*y(k,59) + (rxt(k,151) +rxt(k,159) + &
828 : 2.000_r8*rxt(k,161)*y(k,135))*y(k,225) +rxt(k,253)*y(k,19)*y(k,19) &
829 : +rxt(k,187)*y(k,228)*y(k,91) +rxt(k,193)*y(k,125)*y(k,112) &
830 : +rxt(k,207)*y(k,224)*y(k,113) +rxt(k,590)*y(k,227)*y(k,124) &
831 315224452 : +rxt(k,19)*y(k,126) +rxt(k,152)*y(k,226)
832 : loss(k,235) = (rxt(k,480)* y(k,6) +rxt(k,252)* y(k,17) +rxt(k,318)* y(k,25) &
833 : +rxt(k,347)* y(k,29) +rxt(k,220)* y(k,56) +rxt(k,178)* y(k,76) &
834 : +rxt(k,177)* y(k,90) +rxt(k,425)* y(k,99) +rxt(k,378)* y(k,105) &
835 : +rxt(k,483)* y(k,110) +rxt(k,392)* y(k,111) +rxt(k,204)* y(k,124) &
836 : +rxt(k,196)* y(k,125) +rxt(k,167)* y(k,134) +rxt(k,463)* y(k,143) &
837 : +rxt(k,508)* y(k,150) +rxt(k,514)* y(k,152) +rxt(k,158)* y(k,224) &
838 : +rxt(k,183)* y(k,228) + rxt(k,7) + rxt(k,8) + het_rates(k,136)) &
839 315224452 : * y(k,136)
840 : prod(k,235) = (.150_r8*rxt(k,332)*y(k,200) +.150_r8*rxt(k,383)*y(k,215)) &
841 315224452 : *y(k,90) +rxt(k,169)*y(k,135)*y(k,134)
842 315224452 : loss(k,21) = ( + het_rates(k,137))* y(k,137)
843 315224452 : prod(k,21) = 0._r8
844 : loss(k,109) = (rxt(k,515)* y(k,152) + rxt(k,108) + het_rates(k,138)) &
845 315224452 : * y(k,138)
846 315224452 : prod(k,109) = (rxt(k,224)*y(k,59) +rxt(k,254)*y(k,19))*y(k,59)
847 : loss(k,114) = (rxt(k,504)* y(k,134) +rxt(k,505)* y(k,228) + rxt(k,137) &
848 315224452 : + het_rates(k,139))* y(k,139)
849 315224452 : prod(k,114) = 0._r8
850 315224452 : loss(k,87) = ( + rxt(k,60) + rxt(k,532) + het_rates(k,140))* y(k,140)
851 315224452 : prod(k,87) =rxt(k,371)*y(k,228)*y(k,94) +.100_r8*rxt(k,492)*y(k,233)*y(k,124)
852 : loss(k,141) = (rxt(k,342)* y(k,228) + rxt(k,61) + rxt(k,345) &
853 315224452 : + het_rates(k,141))* y(k,141)
854 315224452 : prod(k,141) =rxt(k,344)*y(k,200)*y(k,125)
855 315224452 : loss(k,66) = ( + rxt(k,470) + het_rates(k,142))* y(k,142)
856 315224452 : prod(k,66) =rxt(k,465)*y(k,191)*y(k,125)
857 : loss(k,131) = (rxt(k,462)* y(k,125) +rxt(k,463)* y(k,136) + het_rates(k,143)) &
858 315224452 : * y(k,143)
859 : prod(k,131) = (.070_r8*rxt(k,449)*y(k,66) +.060_r8*rxt(k,461)*y(k,144) + &
860 : .070_r8*rxt(k,477)*y(k,187))*y(k,228) +rxt(k,32)*y(k,32) &
861 315224452 : +rxt(k,447)*y(k,199)*y(k,124)
862 315224452 : loss(k,73) = (rxt(k,461)* y(k,228) + het_rates(k,144))* y(k,144)
863 315224452 : prod(k,73) =.530_r8*rxt(k,438)*y(k,228)*y(k,7)
864 315224452 : loss(k,110) = (rxt(k,464)* y(k,228) + rxt(k,62) + het_rates(k,145))* y(k,145)
865 315224452 : prod(k,110) =rxt(k,459)*y(k,230)*y(k,90)
866 315224452 : loss(k,22) = ( + het_rates(k,146))* y(k,146)
867 315224452 : prod(k,22) = 0._r8
868 315224452 : loss(k,23) = ( + het_rates(k,147))* y(k,147)
869 315224452 : prod(k,23) = 0._r8
870 315224452 : loss(k,142) = (rxt(k,360)* y(k,228) + rxt(k,63) + het_rates(k,148))* y(k,148)
871 315224452 : prod(k,142) =rxt(k,358)*y(k,231)*y(k,90)
872 315224452 : loss(k,119) = (rxt(k,364)* y(k,228) + rxt(k,64) + het_rates(k,149))* y(k,149)
873 315224452 : prod(k,119) =.850_r8*rxt(k,362)*y(k,232)*y(k,90)
874 : loss(k,163) = (rxt(k,506)* y(k,135) +rxt(k,508)* y(k,136) +rxt(k,511) &
875 315224452 : * y(k,228) + het_rates(k,150))* y(k,150)
876 315224452 : prod(k,163) =rxt(k,137)*y(k,139) +rxt(k,138)*y(k,152)
877 315224452 : loss(k,24) = ( + rxt(k,109) + het_rates(k,151))* y(k,151)
878 315224452 : prod(k,24) = 0._r8
879 : loss(k,219) = (rxt(k,509)* y(k,19) +rxt(k,510)* y(k,59) +rxt(k,512)* y(k,125) &
880 : +rxt(k,513)* y(k,135) +rxt(k,514)* y(k,136) +rxt(k,515)* y(k,138) &
881 315224452 : +rxt(k,516)* y(k,228) + rxt(k,138) + het_rates(k,152))* y(k,152)
882 : prod(k,219) = (rxt(k,506)*y(k,135) +rxt(k,508)*y(k,136) +rxt(k,511)*y(k,228)) &
883 315224452 : *y(k,150) +rxt(k,504)*y(k,139)*y(k,134) +rxt(k,139)*y(k,153)
884 : loss(k,186) = (rxt(k,507)* y(k,228) + rxt(k,139) + het_rates(k,153)) &
885 315224452 : * y(k,153)
886 : prod(k,186) = (rxt(k,509)*y(k,19) +rxt(k,510)*y(k,59) +rxt(k,512)*y(k,125) + &
887 : rxt(k,513)*y(k,135) +rxt(k,514)*y(k,136) +rxt(k,515)*y(k,138) + &
888 : rxt(k,516)*y(k,228))*y(k,152) + (rxt(k,502)*y(k,126) + &
889 : rxt(k,503)*y(k,228) +.500_r8*rxt(k,517)*y(k,228))*y(k,67) &
890 315224452 : +rxt(k,505)*y(k,228)*y(k,139) +rxt(k,140)*y(k,154)
891 315224452 : loss(k,94) = (rxt(k,518)* y(k,241) + rxt(k,140) + het_rates(k,154))* y(k,154)
892 315224452 : prod(k,94) =rxt(k,136)*y(k,80) +rxt(k,507)*y(k,228)*y(k,153)
893 315224452 : loss(k,25) = ( + het_rates(k,155))* y(k,155)
894 315224452 : prod(k,25) = 0._r8
895 315224452 : loss(k,26) = ( + het_rates(k,156))* y(k,156)
896 315224452 : prod(k,26) = 0._r8
897 315224452 : loss(k,27) = ( + het_rates(k,157))* y(k,157)
898 315224452 : prod(k,27) = 0._r8
899 315224452 : loss(k,28) = ( + het_rates(k,158))* y(k,158)
900 315224452 : prod(k,28) = 0._r8
901 315224452 : loss(k,29) = ( + rxt(k,141) + het_rates(k,159))* y(k,159)
902 315224452 : prod(k,29) = 0._r8
903 315224452 : loss(k,30) = ( + rxt(k,142) + het_rates(k,160))* y(k,160)
904 315224452 : prod(k,30) = 0._r8
905 315224452 : loss(k,31) = ( + rxt(k,143) + het_rates(k,161))* y(k,161)
906 315224452 : prod(k,31) = 0._r8
907 315224452 : loss(k,32) = ( + rxt(k,144) + het_rates(k,162))* y(k,162)
908 315224452 : prod(k,32) = 0._r8
909 315224452 : loss(k,33) = ( + rxt(k,145) + het_rates(k,163))* y(k,163)
910 315224452 : prod(k,33) = 0._r8
911 315224452 : loss(k,34) = ( + rxt(k,146) + het_rates(k,164))* y(k,164)
912 315224452 : prod(k,34) = 0._r8
913 315224452 : loss(k,35) = ( + rxt(k,147) + het_rates(k,165))* y(k,165)
914 315224452 : prod(k,35) = 0._r8
915 315224452 : loss(k,36) = ( + rxt(k,148) + het_rates(k,166))* y(k,166)
916 315224452 : prod(k,36) = 0._r8
917 315224452 : loss(k,37) = ( + rxt(k,149) + het_rates(k,167))* y(k,167)
918 315224452 : prod(k,37) = 0._r8
919 315224452 : loss(k,38) = ( + rxt(k,150) + het_rates(k,168))* y(k,168)
920 315224452 : prod(k,38) = 0._r8
921 315224452 : loss(k,39) = ( + het_rates(k,169))* y(k,169)
922 : prod(k,39) = (.2202005_r8*rxt(k,535)*y(k,193) + &
923 : .0023005_r8*rxt(k,540)*y(k,195) +.0031005_r8*rxt(k,543)*y(k,211) + &
924 : .2381005_r8*rxt(k,547)*y(k,212) +.0508005_r8*rxt(k,551)*y(k,218) + &
925 : .1364005_r8*rxt(k,557)*y(k,236) +.1677005_r8*rxt(k,560)*y(k,239)) &
926 : *y(k,90) + (.1279005_r8*rxt(k,536)*y(k,193) + &
927 : .0097005_r8*rxt(k,541)*y(k,195) +.0003005_r8*rxt(k,544)*y(k,211) + &
928 : .1056005_r8*rxt(k,548)*y(k,212) +.0245005_r8*rxt(k,552)*y(k,218) + &
929 : .0154005_r8*rxt(k,558)*y(k,236) +.0063005_r8*rxt(k,561)*y(k,239)) &
930 : *y(k,124) + (.2202005_r8*rxt(k,537)*y(k,6) + &
931 : .0508005_r8*rxt(k,553)*y(k,110))*y(k,136) +rxt(k,520)*y(k,75) &
932 315224452 : +.5931005_r8*rxt(k,555)*y(k,228)*y(k,175)
933 315224452 : loss(k,40) = ( + het_rates(k,170))* y(k,170)
934 : prod(k,40) = (.2067005_r8*rxt(k,535)*y(k,193) + &
935 : .0008005_r8*rxt(k,540)*y(k,195) +.0035005_r8*rxt(k,543)*y(k,211) + &
936 : .1308005_r8*rxt(k,547)*y(k,212) +.1149005_r8*rxt(k,551)*y(k,218) + &
937 : .0101005_r8*rxt(k,557)*y(k,236) +.0174005_r8*rxt(k,560)*y(k,239)) &
938 : *y(k,90) + (.1792005_r8*rxt(k,536)*y(k,193) + &
939 : .0034005_r8*rxt(k,541)*y(k,195) +.0003005_r8*rxt(k,544)*y(k,211) + &
940 : .1026005_r8*rxt(k,548)*y(k,212) +.0082005_r8*rxt(k,552)*y(k,218) + &
941 : .0452005_r8*rxt(k,558)*y(k,236) +.0237005_r8*rxt(k,561)*y(k,239)) &
942 : *y(k,124) + (.2067005_r8*rxt(k,537)*y(k,6) + &
943 : .1149005_r8*rxt(k,553)*y(k,110))*y(k,136) &
944 315224452 : +.1534005_r8*rxt(k,555)*y(k,228)*y(k,175)
945 315224452 : loss(k,41) = ( + het_rates(k,171))* y(k,171)
946 : prod(k,41) = (.0653005_r8*rxt(k,535)*y(k,193) + &
947 : .0843005_r8*rxt(k,540)*y(k,195) +.0003005_r8*rxt(k,543)*y(k,211) + &
948 : .0348005_r8*rxt(k,547)*y(k,212) +.0348005_r8*rxt(k,551)*y(k,218) + &
949 : .0763005_r8*rxt(k,557)*y(k,236) +.086_r8*rxt(k,560)*y(k,239))*y(k,90) &
950 : + (.0676005_r8*rxt(k,536)*y(k,193) + &
951 : .1579005_r8*rxt(k,541)*y(k,195) +.0073005_r8*rxt(k,544)*y(k,211) + &
952 : .0521005_r8*rxt(k,548)*y(k,212) +.0772005_r8*rxt(k,552)*y(k,218) + &
953 : .0966005_r8*rxt(k,558)*y(k,236) +.0025005_r8*rxt(k,561)*y(k,239)) &
954 : *y(k,124) + (.0653005_r8*rxt(k,537)*y(k,6) + &
955 : .0348005_r8*rxt(k,553)*y(k,110))*y(k,136) &
956 315224452 : +.0459005_r8*rxt(k,555)*y(k,228)*y(k,175)
957 315224452 : loss(k,42) = ( + het_rates(k,172))* y(k,172)
958 : prod(k,42) = (.1284005_r8*rxt(k,535)*y(k,193) + &
959 : .0443005_r8*rxt(k,540)*y(k,195) +.0271005_r8*rxt(k,543)*y(k,211) + &
960 : .0076005_r8*rxt(k,547)*y(k,212) +.0554005_r8*rxt(k,551)*y(k,218) + &
961 : .2157005_r8*rxt(k,557)*y(k,236) +.0512005_r8*rxt(k,560)*y(k,239)) &
962 : *y(k,90) + (.079_r8*rxt(k,536)*y(k,193) + &
963 : .0059005_r8*rxt(k,541)*y(k,195) +.0057005_r8*rxt(k,544)*y(k,211) + &
964 : .0143005_r8*rxt(k,548)*y(k,212) +.0332005_r8*rxt(k,552)*y(k,218) + &
965 : .0073005_r8*rxt(k,558)*y(k,236) +.011_r8*rxt(k,561)*y(k,239)) &
966 : *y(k,124) + (.1749305_r8*rxt(k,534)*y(k,6) + &
967 : .0590245_r8*rxt(k,542)*y(k,99) +.1749305_r8*rxt(k,550)*y(k,110)) &
968 : *y(k,126) + (.1284005_r8*rxt(k,537)*y(k,6) + &
969 : .0033005_r8*rxt(k,545)*y(k,99) +.0554005_r8*rxt(k,553)*y(k,110)) &
970 315224452 : *y(k,136) +.0085005_r8*rxt(k,555)*y(k,228)*y(k,175)
971 315224452 : loss(k,43) = ( + het_rates(k,173))* y(k,173)
972 : prod(k,43) = (.114_r8*rxt(k,535)*y(k,193) +.1621005_r8*rxt(k,540)*y(k,195) + &
973 : .0474005_r8*rxt(k,543)*y(k,211) +.0113005_r8*rxt(k,547)*y(k,212) + &
974 : .1278005_r8*rxt(k,551)*y(k,218) +.0738005_r8*rxt(k,557)*y(k,236) + &
975 : .1598005_r8*rxt(k,560)*y(k,239))*y(k,90) &
976 : + (.1254005_r8*rxt(k,536)*y(k,193) + &
977 : .0536005_r8*rxt(k,541)*y(k,195) +.0623005_r8*rxt(k,544)*y(k,211) + &
978 : .0166005_r8*rxt(k,548)*y(k,212) +.130_r8*rxt(k,552)*y(k,218) + &
979 : .238_r8*rxt(k,558)*y(k,236) +.1185005_r8*rxt(k,561)*y(k,239)) &
980 : *y(k,124) + (.5901905_r8*rxt(k,534)*y(k,6) + &
981 : .0250245_r8*rxt(k,542)*y(k,99) +.5901905_r8*rxt(k,550)*y(k,110)) &
982 : *y(k,126) + (.114_r8*rxt(k,537)*y(k,6) + &
983 : .1278005_r8*rxt(k,553)*y(k,110))*y(k,136) &
984 315224452 : +.0128005_r8*rxt(k,555)*y(k,228)*y(k,175)
985 315224452 : loss(k,44) = ( + rxt(k,597) + het_rates(k,174))* y(k,174)
986 315224452 : prod(k,44) = 0._r8
987 315224452 : loss(k,45) = (rxt(k,555)* y(k,228) + het_rates(k,175))* y(k,175)
988 315224452 : prod(k,45) = 0._r8
989 315224452 : loss(k,80) = ( + rxt(k,65) + het_rates(k,176))* y(k,176)
990 : prod(k,80) = (.100_r8*rxt(k,469)*y(k,183) +.230_r8*rxt(k,471)*y(k,185)) &
991 315224452 : *y(k,228)
992 315224452 : loss(k,156) = (rxt(k,493)* y(k,228) + rxt(k,66) + het_rates(k,177))* y(k,177)
993 315224452 : prod(k,156) =rxt(k,491)*y(k,233)*y(k,90)
994 : loss(k,153) = (rxt(k,494)* y(k,228) + rxt(k,67) + rxt(k,533) &
995 315224452 : + het_rates(k,178))* y(k,178)
996 : prod(k,153) = (.200_r8*rxt(k,487)*y(k,223) +.200_r8*rxt(k,497)*y(k,234)) &
997 315224452 : *y(k,124) +.500_r8*rxt(k,485)*y(k,223)*y(k,201)
998 315224452 : loss(k,134) = (rxt(k,498)* y(k,228) + rxt(k,68) + het_rates(k,179))* y(k,179)
999 315224452 : prod(k,134) =rxt(k,496)*y(k,234)*y(k,90)
1000 : loss(k,194) = (rxt(k,499)* y(k,126) +rxt(k,500)* y(k,228) + rxt(k,69) &
1001 315224452 : + het_rates(k,180))* y(k,180)
1002 : prod(k,194) = (.500_r8*rxt(k,485)*y(k,201) +.800_r8*rxt(k,487)*y(k,124) + &
1003 : rxt(k,488)*y(k,126))*y(k,223) + (.330_r8*rxt(k,480)*y(k,6) + &
1004 : .330_r8*rxt(k,483)*y(k,110))*y(k,136) + (rxt(k,67) + &
1005 : rxt(k,494)*y(k,228))*y(k,178) + (rxt(k,495)*y(k,201) + &
1006 : .800_r8*rxt(k,497)*y(k,124))*y(k,234) +rxt(k,59)*y(k,128) +rxt(k,68) &
1007 315224452 : *y(k,179)
1008 315224452 : loss(k,198) = (rxt(k,501)* y(k,228) + rxt(k,70) + het_rates(k,181))* y(k,181)
1009 : prod(k,198) = (.300_r8*rxt(k,480)*y(k,6) +.300_r8*rxt(k,483)*y(k,110)) &
1010 : *y(k,136) + (rxt(k,490)*y(k,201) +.900_r8*rxt(k,492)*y(k,124)) &
1011 315224452 : *y(k,233) +rxt(k,66)*y(k,177) +rxt(k,69)*y(k,180)
1012 315224452 : loss(k,157) = (rxt(k,468)* y(k,228) + rxt(k,71) + het_rates(k,182))* y(k,182)
1013 315224452 : prod(k,157) =rxt(k,466)*y(k,235)*y(k,90)
1014 315224452 : loss(k,78) = (rxt(k,469)* y(k,228) + het_rates(k,183))* y(k,183)
1015 315224452 : prod(k,78) = 0._r8
1016 315224452 : loss(k,81) = (rxt(k,435)* y(k,228) + rxt(k,72) + het_rates(k,184))* y(k,184)
1017 315224452 : prod(k,81) =rxt(k,432)*y(k,237)*y(k,90)
1018 315224452 : loss(k,82) = (rxt(k,471)* y(k,228) + het_rates(k,185))* y(k,185)
1019 315224452 : prod(k,82) = 0._r8
1020 315224452 : loss(k,164) = (rxt(k,474)* y(k,228) + rxt(k,73) + het_rates(k,186))* y(k,186)
1021 315224452 : prod(k,164) =rxt(k,472)*y(k,238)*y(k,90)
1022 315224452 : loss(k,83) = (rxt(k,477)* y(k,228) + het_rates(k,187))* y(k,187)
1023 315224452 : prod(k,83) =.150_r8*rxt(k,471)*y(k,228)*y(k,185)
1024 315224452 : loss(k,123) = (rxt(k,478)* y(k,228) + rxt(k,74) + het_rates(k,188))* y(k,188)
1025 315224452 : prod(k,123) =rxt(k,475)*y(k,240)*y(k,90)
1026 : loss(k,139) = (rxt(k,436)* y(k,90) +rxt(k,437)* y(k,124) +rxt(k,465) &
1027 315224452 : * y(k,125) + het_rates(k,191))* y(k,191)
1028 315224452 : prod(k,139) =rxt(k,442)*y(k,228)*y(k,22) +rxt(k,470)*y(k,142)
1029 : loss(k,191) = (rxt(k,397)* y(k,90) + (rxt(k,398) +rxt(k,399))* y(k,124) &
1030 315224452 : + het_rates(k,192))* y(k,192)
1031 315224452 : prod(k,191) = (rxt(k,400)*y(k,2) +rxt(k,401)*y(k,15))*y(k,228)
1032 : loss(k,46) = (rxt(k,535)* y(k,90) +rxt(k,536)* y(k,124) + het_rates(k,193)) &
1033 315224452 : * y(k,193)
1034 315224452 : prod(k,46) =rxt(k,538)*y(k,228)*y(k,6)
1035 : loss(k,135) = (rxt(k,439)* y(k,90) +rxt(k,440)* y(k,124) + het_rates(k,194)) &
1036 315224452 : * y(k,194)
1037 315224452 : prod(k,135) = (.350_r8*rxt(k,438)*y(k,7) +rxt(k,441)*y(k,8))*y(k,228)
1038 : loss(k,47) = (rxt(k,540)* y(k,90) +rxt(k,541)* y(k,124) + het_rates(k,195)) &
1039 315224452 : * y(k,195)
1040 315224452 : prod(k,47) =rxt(k,539)*y(k,228)*y(k,7)
1041 : loss(k,125) = (rxt(k,443)* y(k,90) +rxt(k,445)* y(k,124) + het_rates(k,196)) &
1042 315224452 : * y(k,196)
1043 : prod(k,125) = (rxt(k,444)*y(k,23) +.070_r8*rxt(k,469)*y(k,183) + &
1044 315224452 : .060_r8*rxt(k,471)*y(k,185))*y(k,228)
1045 : loss(k,181) = (rxt(k,321)* y(k,90) +rxt(k,322)* y(k,124) + 2._r8*rxt(k,319) &
1046 315224452 : * y(k,197) +rxt(k,320)* y(k,201) + het_rates(k,197))* y(k,197)
1047 : prod(k,181) = (rxt(k,325)*y(k,56) +rxt(k,326)*y(k,228))*y(k,28) &
1048 315224452 : +.500_r8*rxt(k,324)*y(k,228)*y(k,27) +rxt(k,53)*y(k,107)
1049 : loss(k,185) = (rxt(k,349)* y(k,90) +rxt(k,350)* y(k,124) +rxt(k,348) &
1050 315224452 : * y(k,201) + het_rates(k,198))* y(k,198)
1051 : prod(k,185) = (rxt(k,352)*y(k,228) +rxt(k,355)*y(k,56))*y(k,31) &
1052 315224452 : +rxt(k,351)*y(k,228)*y(k,30)
1053 : loss(k,154) = (rxt(k,446)* y(k,90) +rxt(k,447)* y(k,124) + het_rates(k,199)) &
1054 315224452 : * y(k,199)
1055 : prod(k,154) = (.400_r8*rxt(k,436)*y(k,90) +rxt(k,437)*y(k,124))*y(k,191) &
1056 315224452 : +rxt(k,448)*y(k,228)*y(k,32) +rxt(k,463)*y(k,143)*y(k,136)
1057 : loss(k,216) = (rxt(k,332)* y(k,90) +rxt(k,333)* y(k,124) +rxt(k,344) &
1058 : * y(k,125) + 2._r8*rxt(k,330)* y(k,200) +rxt(k,331)* y(k,201) &
1059 : +rxt(k,405)* y(k,208) +rxt(k,410)* y(k,209) +rxt(k,419)* y(k,210) &
1060 : +rxt(k,372)* y(k,213) +rxt(k,430)* y(k,237) + het_rates(k,200)) &
1061 315224452 : * y(k,200)
1062 : prod(k,216) = (.100_r8*rxt(k,378)*y(k,105) +.280_r8*rxt(k,392)*y(k,111) + &
1063 : .080_r8*rxt(k,425)*y(k,99) +.060_r8*rxt(k,480)*y(k,6) + &
1064 : .060_r8*rxt(k,483)*y(k,110))*y(k,136) + (rxt(k,382)*y(k,201) + &
1065 : .450_r8*rxt(k,383)*y(k,90) +2.000_r8*rxt(k,384)*y(k,215) + &
1066 : rxt(k,385)*y(k,124) +rxt(k,386)*y(k,126))*y(k,215) &
1067 : + (.530_r8*rxt(k,372)*y(k,200) +.260_r8*rxt(k,373)*y(k,201) + &
1068 : .530_r8*rxt(k,375)*y(k,126) +.530_r8*rxt(k,376)*y(k,124))*y(k,213) &
1069 : + (rxt(k,328)*y(k,45) +.500_r8*rxt(k,335)*y(k,51) + &
1070 : rxt(k,354)*y(k,49) +.650_r8*rxt(k,501)*y(k,181))*y(k,228) &
1071 : + (.300_r8*rxt(k,361)*y(k,201) +.150_r8*rxt(k,362)*y(k,90) + &
1072 : rxt(k,363)*y(k,124))*y(k,232) + (rxt(k,37) +rxt(k,353)*y(k,126)) &
1073 : *y(k,49) + (.600_r8*rxt(k,61) +rxt(k,345))*y(k,141) &
1074 : + (.200_r8*rxt(k,387)*y(k,90) +rxt(k,388)*y(k,124))*y(k,217) &
1075 : +.130_r8*rxt(k,24)*y(k,10) +rxt(k,28)*y(k,14) +rxt(k,327)*y(k,126) &
1076 : *y(k,45) +rxt(k,36)*y(k,48) +.330_r8*rxt(k,46)*y(k,94) +rxt(k,48) &
1077 : *y(k,96) +1.340_r8*rxt(k,51)*y(k,105) +rxt(k,53)*y(k,107) +rxt(k,54) &
1078 : *y(k,108) +.300_r8*rxt(k,56)*y(k,111) +rxt(k,58)*y(k,127) +rxt(k,64) &
1079 315224452 : *y(k,149) +.500_r8*rxt(k,65)*y(k,176) +.650_r8*rxt(k,70)*y(k,181)
1080 : loss(k,238) = (rxt(k,221)* y(k,59) +rxt(k,300)* y(k,90) +rxt(k,301)* y(k,124) &
1081 : +rxt(k,320)* y(k,197) +rxt(k,348)* y(k,198) +rxt(k,331)* y(k,200) &
1082 : + 2._r8*(rxt(k,298) +rxt(k,299))* y(k,201) +rxt(k,406)* y(k,208) &
1083 : +rxt(k,411)* y(k,209) +rxt(k,420)* y(k,210) +rxt(k,373)* y(k,213) &
1084 : +rxt(k,382)* y(k,215) +rxt(k,485)* y(k,223) +rxt(k,361)* y(k,232) &
1085 : +rxt(k,490)* y(k,233) +rxt(k,495)* y(k,234) +rxt(k,431)* y(k,237) &
1086 315224452 : + het_rates(k,201))* y(k,201)
1087 : prod(k,238) = (2.000_r8*rxt(k,330)*y(k,200) +.900_r8*rxt(k,331)*y(k,201) + &
1088 : .450_r8*rxt(k,332)*y(k,90) +rxt(k,333)*y(k,124) + &
1089 : rxt(k,372)*y(k,213) +rxt(k,381)*y(k,215) +rxt(k,405)*y(k,208) + &
1090 : rxt(k,410)*y(k,209) +rxt(k,419)*y(k,210) +rxt(k,430)*y(k,237)) &
1091 : *y(k,200) + (rxt(k,41) +rxt(k,215)*y(k,56) +rxt(k,271)*y(k,73) + &
1092 : rxt(k,304)*y(k,228) +rxt(k,310)*y(k,224))*y(k,54) &
1093 : + (.330_r8*rxt(k,450)*y(k,202) +.070_r8*rxt(k,456)*y(k,216))*y(k,90) &
1094 : + (.830_r8*rxt(k,451)*y(k,202) +.170_r8*rxt(k,457)*y(k,216)) &
1095 : *y(k,124) + (.280_r8*rxt(k,347)*y(k,29) +.050_r8*rxt(k,425)*y(k,99)) &
1096 : *y(k,136) + (.700_r8*rxt(k,303)*y(k,53) +rxt(k,334)*y(k,50))*y(k,228) &
1097 : +rxt(k,88)*y(k,43) +rxt(k,35)*y(k,45) +rxt(k,90)*y(k,46) +rxt(k,36) &
1098 : *y(k,48) +rxt(k,38)*y(k,51) +.300_r8*rxt(k,56)*y(k,111) &
1099 315224452 : +.400_r8*rxt(k,61)*y(k,141)
1100 : loss(k,168) = (rxt(k,450)* y(k,90) +rxt(k,451)* y(k,124) +rxt(k,452) &
1101 315224452 : * y(k,125) + het_rates(k,202))* y(k,202)
1102 315224452 : prod(k,168) =.600_r8*rxt(k,26)*y(k,12)
1103 : loss(k,178) = (rxt(k,581)* y(k,220) +rxt(k,579)* y(k,221) +rxt(k,580) &
1104 315224452 : * y(k,227) + het_rates(k,203))* y(k,203)
1105 : prod(k,178) = (rxt(k,126) +rxt(k,127) +rxt(k,128) +rxt(k,129) +rxt(k,130) + &
1106 : rxt(k,131) +rxt(k,132) +rxt(k,133))*y(k,135) + (rxt(k,120) + &
1107 : rxt(k,121) +rxt(k,122) +rxt(k,123) +rxt(k,124) +rxt(k,125))*y(k,134) &
1108 315224452 : +rxt(k,111)*y(k,112) +rxt(k,16)*y(k,124)
1109 : loss(k,146) = ((rxt(k,369) +rxt(k,370))* y(k,124) + het_rates(k,204)) &
1110 315224452 : * y(k,204)
1111 315224452 : prod(k,146) =rxt(k,368)*y(k,228)*y(k,16)
1112 : loss(k,129) = (rxt(k,339)* y(k,135) + rxt(k,338) + het_rates(k,205)) &
1113 315224452 : * y(k,205)
1114 315224452 : prod(k,129) =rxt(k,43)*y(k,72) +.750_r8*rxt(k,337)*y(k,206)*y(k,124)
1115 : loss(k,170) = (rxt(k,336)* y(k,90) +rxt(k,337)* y(k,124) + het_rates(k,206)) &
1116 315224452 : * y(k,206)
1117 315224452 : prod(k,170) =rxt(k,343)*y(k,228)*y(k,25)
1118 : loss(k,126) = (rxt(k,307)* y(k,90) +rxt(k,309)* y(k,124) + rxt(k,308) &
1119 315224452 : + het_rates(k,207))* y(k,207)
1120 315224452 : prod(k,126) =rxt(k,294)*y(k,90)*y(k,42)
1121 : loss(k,211) = (rxt(k,407)* y(k,90) +rxt(k,408)* y(k,124) +rxt(k,409) &
1122 : * y(k,126) +rxt(k,405)* y(k,200) +rxt(k,406)* y(k,201) &
1123 315224452 : + het_rates(k,208))* y(k,208)
1124 315224452 : prod(k,211) =.600_r8*rxt(k,426)*y(k,228)*y(k,99)
1125 : loss(k,212) = (rxt(k,412)* y(k,90) +rxt(k,414)* y(k,124) +rxt(k,415) &
1126 : * y(k,126) +rxt(k,410)* y(k,200) +rxt(k,411)* y(k,201) + rxt(k,413) &
1127 315224452 : + het_rates(k,209))* y(k,209)
1128 315224452 : prod(k,212) =.400_r8*rxt(k,426)*y(k,228)*y(k,99)
1129 : loss(k,208) = (rxt(k,421)* y(k,90) +rxt(k,422)* y(k,124) +rxt(k,423) &
1130 : * y(k,126) +rxt(k,419)* y(k,200) +rxt(k,420)* y(k,201) &
1131 315224452 : + het_rates(k,210))* y(k,210)
1132 315224452 : prod(k,208) =rxt(k,418)*y(k,126)*y(k,99)
1133 : loss(k,48) = (rxt(k,543)* y(k,90) +rxt(k,544)* y(k,124) + het_rates(k,211)) &
1134 315224452 : * y(k,211)
1135 315224452 : prod(k,48) =rxt(k,546)*y(k,228)*y(k,99)
1136 : loss(k,49) = (rxt(k,547)* y(k,90) +rxt(k,548)* y(k,124) + het_rates(k,212)) &
1137 315224452 : * y(k,212)
1138 315224452 : prod(k,49) =rxt(k,549)*y(k,228)*y(k,104)
1139 : loss(k,213) = (rxt(k,374)* y(k,90) + (rxt(k,376) +rxt(k,377))* y(k,124) &
1140 : +rxt(k,375)* y(k,126) +rxt(k,372)* y(k,200) +rxt(k,373)* y(k,201) &
1141 315224452 : + het_rates(k,213))* y(k,213)
1142 : prod(k,213) = (.500_r8*rxt(k,379)*y(k,105) +.200_r8*rxt(k,380)*y(k,106) + &
1143 315224452 : rxt(k,393)*y(k,111))*y(k,228)
1144 : loss(k,165) = (rxt(k,453)* y(k,90) +rxt(k,454)* y(k,124) +rxt(k,455) &
1145 315224452 : * y(k,125) + het_rates(k,214))* y(k,214)
1146 315224452 : prod(k,165) =.600_r8*rxt(k,25)*y(k,11)
1147 : loss(k,215) = (rxt(k,383)* y(k,90) +rxt(k,385)* y(k,124) +rxt(k,394) &
1148 : * y(k,125) +rxt(k,386)* y(k,126) +rxt(k,381)* y(k,200) +rxt(k,382) &
1149 315224452 : * y(k,201) + 2._r8*rxt(k,384)* y(k,215) + het_rates(k,215))* y(k,215)
1150 : prod(k,215) = (.660_r8*rxt(k,51) +.500_r8*rxt(k,379)*y(k,228))*y(k,105) &
1151 : + (rxt(k,55) +rxt(k,395))*y(k,109) +.500_r8*rxt(k,380)*y(k,228) &
1152 315224452 : *y(k,106)
1153 : loss(k,182) = (rxt(k,456)* y(k,90) +rxt(k,457)* y(k,124) +rxt(k,458) &
1154 315224452 : * y(k,125) + het_rates(k,216))* y(k,216)
1155 315224452 : prod(k,182) =.600_r8*rxt(k,27)*y(k,13)
1156 : loss(k,160) = (rxt(k,387)* y(k,90) +rxt(k,388)* y(k,124) + het_rates(k,217)) &
1157 315224452 : * y(k,217)
1158 315224452 : prod(k,160) = (rxt(k,389)*y(k,107) +rxt(k,390)*y(k,108))*y(k,228)
1159 : loss(k,51) = (rxt(k,551)* y(k,90) +rxt(k,552)* y(k,124) + het_rates(k,218)) &
1160 315224452 : * y(k,218)
1161 315224452 : prod(k,51) =rxt(k,554)*y(k,228)*y(k,110)
1162 : loss(k,140) = (rxt(k,188)* y(k,134) +rxt(k,189)* y(k,135) + het_rates(k,219)) &
1163 315224452 : * y(k,219)
1164 : prod(k,140) = (.800_r8*rxt(k,579)*y(k,221) +.900_r8*rxt(k,581)*y(k,220)) &
1165 315224452 : *y(k,203) +rxt(k,583)*y(k,220)*y(k,134)
1166 : loss(k,158) = ((rxt(k,583) +rxt(k,584))* y(k,134) +rxt(k,582)* y(k,135) &
1167 315224452 : +rxt(k,581)* y(k,203) + het_rates(k,220))* y(k,220)
1168 315224452 : prod(k,158) = 0._r8
1169 315224452 : loss(k,176) = (rxt(k,579)* y(k,203) + het_rates(k,221))* y(k,221)
1170 : prod(k,176) = (rxt(k,589) +rxt(k,588)*y(k,112) +rxt(k,590)*y(k,124))*y(k,227) &
1171 : +rxt(k,16)*y(k,124) +rxt(k,583)*y(k,220)*y(k,134) &
1172 315224452 : +rxt(k,587)*y(k,222)*y(k,135) +rxt(k,592)*y(k,229)
1173 : loss(k,136) = (rxt(k,585)* y(k,134) + (rxt(k,586) +rxt(k,587))* y(k,135) &
1174 315224452 : + het_rates(k,222))* y(k,222)
1175 315224452 : prod(k,136) =rxt(k,111)*y(k,112)
1176 : loss(k,197) = (rxt(k,486)* y(k,90) +rxt(k,487)* y(k,124) +rxt(k,488) &
1177 315224452 : * y(k,126) +rxt(k,485)* y(k,201) + het_rates(k,223))* y(k,223)
1178 : prod(k,197) = (rxt(k,479)*y(k,6) +rxt(k,482)*y(k,110) + &
1179 315224452 : .500_r8*rxt(k,499)*y(k,180))*y(k,126) +rxt(k,489)*y(k,228)*y(k,128)
1180 : loss(k,230) = (rxt(k,239)* y(k,33) +rxt(k,240)* y(k,34) +rxt(k,266)* y(k,35) &
1181 : +rxt(k,241)* y(k,36) +rxt(k,242)* y(k,37) +rxt(k,243)* y(k,38) &
1182 : +rxt(k,244)* y(k,39) +rxt(k,245)* y(k,40) +rxt(k,289)* y(k,41) &
1183 : +rxt(k,290)* y(k,43) + (rxt(k,310) +rxt(k,311) +rxt(k,312))* y(k,54) &
1184 : +rxt(k,267)* y(k,55) +rxt(k,275)* y(k,64) +rxt(k,276)* y(k,65) &
1185 : +rxt(k,153)* y(k,77) +rxt(k,268)* y(k,78) + (rxt(k,269) +rxt(k,270)) &
1186 : * y(k,81) +rxt(k,291)* y(k,82) +rxt(k,292)* y(k,83) +rxt(k,293) &
1187 : * y(k,84) + (rxt(k,246) +rxt(k,247))* y(k,85) +rxt(k,313)* y(k,86) &
1188 : + (rxt(k,206) +rxt(k,207))* y(k,113) + (rxt(k,156) +rxt(k,157)) &
1189 : * y(k,135) +rxt(k,158)* y(k,136) +rxt(k,154)* y(k,241) + rxt(k,155) &
1190 315224452 : + het_rates(k,224))* y(k,224)
1191 : prod(k,230) = (rxt(k,6) +rxt(k,189)*y(k,219))*y(k,135) +rxt(k,12)*y(k,113) &
1192 : +rxt(k,7)*y(k,136) +.850_r8*rxt(k,580)*y(k,227)*y(k,203) +rxt(k,1) &
1193 315224452 : *y(k,241)
1194 : loss(k,75) = (rxt(k,160)* y(k,134) +rxt(k,161)* y(k,135) + rxt(k,151) &
1195 315224452 : + rxt(k,159) + het_rates(k,225))* y(k,225)
1196 : prod(k,75) = (rxt(k,163) +rxt(k,162)*y(k,63) +rxt(k,164)*y(k,134) + &
1197 315224452 : rxt(k,165)*y(k,135) +rxt(k,166)*y(k,136))*y(k,226) +rxt(k,7)*y(k,136)
1198 : loss(k,76) = (rxt(k,162)* y(k,63) +rxt(k,164)* y(k,134) +rxt(k,165)* y(k,135) &
1199 : +rxt(k,166)* y(k,136) + rxt(k,152) + rxt(k,163) + het_rates(k,226)) &
1200 315224452 : * y(k,226)
1201 315224452 : prod(k,76) =rxt(k,156)*y(k,224)*y(k,135)
1202 : loss(k,177) = (rxt(k,588)* y(k,112) +rxt(k,590)* y(k,124) +rxt(k,580) &
1203 315224452 : * y(k,203) + rxt(k,589) + het_rates(k,227))* y(k,227)
1204 : prod(k,177) = (rxt(k,126) +rxt(k,130) +rxt(k,582)*y(k,220) + &
1205 : rxt(k,586)*y(k,222) +rxt(k,593)*y(k,229))*y(k,135) &
1206 315224452 : +rxt(k,591)*y(k,229)*y(k,63)
1207 : loss(k,226) = (rxt(k,396)* y(k,1) +rxt(k,400)* y(k,2) +rxt(k,481)* y(k,6) &
1208 : +rxt(k,438)* y(k,7) +rxt(k,441)* y(k,8) +rxt(k,401)* y(k,15) &
1209 : +rxt(k,368)* y(k,16) +rxt(k,262)* y(k,19) +rxt(k,442)* y(k,22) &
1210 : +rxt(k,444)* y(k,23) +rxt(k,316)* y(k,24) +rxt(k,343)* y(k,25) &
1211 : +rxt(k,323)* y(k,26) +rxt(k,324)* y(k,27) +rxt(k,326)* y(k,28) &
1212 : +rxt(k,365)* y(k,29) +rxt(k,351)* y(k,30) +rxt(k,352)* y(k,31) &
1213 : +rxt(k,448)* y(k,32) +rxt(k,278)* y(k,41) +rxt(k,297)* y(k,42) &
1214 : +rxt(k,280)* y(k,43) +rxt(k,281)* y(k,44) +rxt(k,328)* y(k,45) &
1215 : +rxt(k,283)* y(k,46) +rxt(k,329)* y(k,47) +rxt(k,366)* y(k,48) &
1216 : +rxt(k,354)* y(k,49) +rxt(k,334)* y(k,50) +rxt(k,335)* y(k,51) &
1217 : +rxt(k,302)* y(k,52) +rxt(k,303)* y(k,53) +rxt(k,304)* y(k,54) &
1218 : +rxt(k,285)* y(k,55) + (rxt(k,232) +rxt(k,233))* y(k,59) +rxt(k,230) &
1219 : * y(k,60) +rxt(k,314)* y(k,62) +rxt(k,449)* y(k,66) + (rxt(k,503) + &
1220 : rxt(k,517))* y(k,67) +rxt(k,340)* y(k,74) +rxt(k,341)* y(k,75) &
1221 : +rxt(k,179)* y(k,77) +rxt(k,180)* y(k,79) +rxt(k,264)* y(k,81) &
1222 : +rxt(k,286)* y(k,82) +rxt(k,287)* y(k,83) +rxt(k,288)* y(k,84) &
1223 : +rxt(k,235)* y(k,85) +rxt(k,305)* y(k,86) +rxt(k,306)* y(k,87) &
1224 : +rxt(k,211)* y(k,89) +rxt(k,181)* y(k,90) +rxt(k,187)* y(k,91) &
1225 : +rxt(k,238)* y(k,93) +rxt(k,371)* y(k,94) +rxt(k,402)* y(k,95) &
1226 : +rxt(k,356)* y(k,96) +rxt(k,403)* y(k,97) +rxt(k,404)* y(k,98) &
1227 : +rxt(k,426)* y(k,99) +rxt(k,416)* y(k,100) +rxt(k,417)* y(k,101) &
1228 : +rxt(k,424)* y(k,102) +rxt(k,427)* y(k,103) +rxt(k,379)* y(k,105) &
1229 : +rxt(k,380)* y(k,106) +rxt(k,389)* y(k,107) +rxt(k,390)* y(k,108) &
1230 : +rxt(k,391)* y(k,109) +rxt(k,484)* y(k,110) +rxt(k,393)* y(k,111) &
1231 : +rxt(k,202)* y(k,112) +rxt(k,428)* y(k,115) +rxt(k,429)* y(k,116) &
1232 : +rxt(k,519)* y(k,120) +rxt(k,210)* y(k,125) +rxt(k,201)* y(k,126) &
1233 : +rxt(k,357)* y(k,127) +rxt(k,489)* y(k,128) +rxt(k,182)* y(k,134) &
1234 : +rxt(k,183)* y(k,136) +rxt(k,505)* y(k,139) +rxt(k,342)* y(k,141) &
1235 : +rxt(k,461)* y(k,144) +rxt(k,464)* y(k,145) +rxt(k,360)* y(k,148) &
1236 : +rxt(k,364)* y(k,149) +rxt(k,511)* y(k,150) +rxt(k,516)* y(k,152) &
1237 : +rxt(k,507)* y(k,153) +rxt(k,493)* y(k,177) +rxt(k,494)* y(k,178) &
1238 : +rxt(k,498)* y(k,179) +rxt(k,500)* y(k,180) +rxt(k,501)* y(k,181) &
1239 : +rxt(k,468)* y(k,182) +rxt(k,469)* y(k,183) +rxt(k,435)* y(k,184) &
1240 : +rxt(k,471)* y(k,185) +rxt(k,474)* y(k,186) +rxt(k,477)* y(k,187) &
1241 : +rxt(k,478)* y(k,188) + 2._r8*(rxt(k,184) +rxt(k,185))* y(k,228) &
1242 315224452 : + het_rates(k,228))* y(k,228)
1243 : prod(k,226) = (2.000_r8*rxt(k,173)*y(k,76) +rxt(k,176)*y(k,134) + &
1244 : rxt(k,177)*y(k,136) +rxt(k,198)*y(k,126) +rxt(k,203)*y(k,124) + &
1245 : rxt(k,219)*y(k,56) +.450_r8*rxt(k,332)*y(k,200) + &
1246 : .150_r8*rxt(k,362)*y(k,232) +.450_r8*rxt(k,383)*y(k,215) + &
1247 : .200_r8*rxt(k,387)*y(k,217) +.400_r8*rxt(k,436)*y(k,191) + &
1248 : .400_r8*rxt(k,450)*y(k,202) +.400_r8*rxt(k,456)*y(k,216))*y(k,90) &
1249 : + (rxt(k,178)*y(k,76) +.130_r8*rxt(k,318)*y(k,25) + &
1250 : .360_r8*rxt(k,347)*y(k,29) +.240_r8*rxt(k,378)*y(k,105) + &
1251 : .360_r8*rxt(k,392)*y(k,111) +.320_r8*rxt(k,425)*y(k,99) + &
1252 : .630_r8*rxt(k,480)*y(k,6) +.630_r8*rxt(k,483)*y(k,110))*y(k,136) &
1253 : + (rxt(k,170)*y(k,77) +rxt(k,171)*y(k,79) +rxt(k,234)*y(k,85) + &
1254 : rxt(k,237)*y(k,93) +rxt(k,263)*y(k,81) +rxt(k,265)*y(k,92) + &
1255 : rxt(k,296)*y(k,42))*y(k,134) + (.300_r8*rxt(k,303)*y(k,53) + &
1256 : .650_r8*rxt(k,316)*y(k,24) +.500_r8*rxt(k,324)*y(k,27) + &
1257 : .500_r8*rxt(k,360)*y(k,148) +.100_r8*rxt(k,380)*y(k,106) + &
1258 : .600_r8*rxt(k,427)*y(k,103) +.500_r8*rxt(k,435)*y(k,184))*y(k,228) &
1259 : + (rxt(k,153)*y(k,77) +2.000_r8*rxt(k,154)*y(k,241) + &
1260 : rxt(k,246)*y(k,85) +rxt(k,269)*y(k,81) +rxt(k,310)*y(k,54) + &
1261 : rxt(k,313)*y(k,86))*y(k,224) + (rxt(k,3) +rxt(k,273)*y(k,73)) &
1262 : *y(k,241) +rxt(k,21)*y(k,2) +rxt(k,22)*y(k,8) +rxt(k,29)*y(k,23) &
1263 : +rxt(k,30)*y(k,27) +rxt(k,31)*y(k,30) +rxt(k,32)*y(k,32) +rxt(k,38) &
1264 : *y(k,51) +rxt(k,39)*y(k,53) +.330_r8*rxt(k,40)*y(k,54) +rxt(k,43) &
1265 : *y(k,72) +2.000_r8*rxt(k,4)*y(k,79) +rxt(k,9)*y(k,89) +rxt(k,10) &
1266 : *y(k,91) +rxt(k,106)*y(k,92) +rxt(k,107)*y(k,93) +rxt(k,47)*y(k,95) &
1267 : +rxt(k,50)*y(k,103) +rxt(k,54)*y(k,108) +.500_r8*rxt(k,529)*y(k,125) &
1268 : +rxt(k,59)*y(k,128) +rxt(k,62)*y(k,145) +rxt(k,63)*y(k,148) &
1269 : +rxt(k,64)*y(k,149) +rxt(k,66)*y(k,177) +rxt(k,68)*y(k,179) &
1270 : +rxt(k,71)*y(k,182) +rxt(k,72)*y(k,184) +rxt(k,73)*y(k,186) &
1271 315224452 : +rxt(k,74)*y(k,188)
1272 : loss(k,171) = (rxt(k,591)* y(k,63) +rxt(k,593)* y(k,135) + rxt(k,592) &
1273 315224452 : + het_rates(k,229))* y(k,229)
1274 : prod(k,171) = (rxt(k,120) +rxt(k,121) +rxt(k,122) +rxt(k,123) +rxt(k,124) + &
1275 : rxt(k,125) +rxt(k,584)*y(k,220) +rxt(k,585)*y(k,222))*y(k,134) &
1276 : + (rxt(k,127) +rxt(k,128) +rxt(k,129) +rxt(k,131) +rxt(k,132) + &
1277 315224452 : rxt(k,133))*y(k,135)
1278 : loss(k,127) = (rxt(k,459)* y(k,90) +rxt(k,460)* y(k,124) + het_rates(k,230)) &
1279 315224452 : * y(k,230)
1280 : prod(k,127) = (.200_r8*rxt(k,449)*y(k,66) +.140_r8*rxt(k,461)*y(k,144) + &
1281 315224452 : rxt(k,464)*y(k,145))*y(k,228)
1282 : loss(k,172) = (rxt(k,358)* y(k,90) +rxt(k,359)* y(k,124) + het_rates(k,231)) &
1283 315224452 : * y(k,231)
1284 315224452 : prod(k,172) = (.500_r8*rxt(k,360)*y(k,148) +rxt(k,365)*y(k,29))*y(k,228)
1285 : loss(k,206) = (rxt(k,362)* y(k,90) +rxt(k,363)* y(k,124) +rxt(k,361) &
1286 315224452 : * y(k,201) + het_rates(k,232))* y(k,232)
1287 : prod(k,206) = (rxt(k,364)*y(k,149) +rxt(k,366)*y(k,48) + &
1288 : .150_r8*rxt(k,501)*y(k,181))*y(k,228) + (.060_r8*rxt(k,480)*y(k,6) + &
1289 315224452 : .060_r8*rxt(k,483)*y(k,110))*y(k,136) +.150_r8*rxt(k,70)*y(k,181)
1290 : loss(k,203) = (rxt(k,491)* y(k,90) +rxt(k,492)* y(k,124) +rxt(k,490) &
1291 315224452 : * y(k,201) + het_rates(k,233))* y(k,233)
1292 : prod(k,203) = (.500_r8*rxt(k,499)*y(k,126) +rxt(k,500)*y(k,228))*y(k,180) &
1293 315224452 : +rxt(k,493)*y(k,228)*y(k,177)
1294 : loss(k,193) = (rxt(k,496)* y(k,90) +rxt(k,497)* y(k,124) +rxt(k,495) &
1295 315224452 : * y(k,201) + het_rates(k,234))* y(k,234)
1296 : prod(k,193) = (rxt(k,481)*y(k,6) +rxt(k,484)*y(k,110) +rxt(k,498)*y(k,179)) &
1297 315224452 : *y(k,228)
1298 : loss(k,166) = (rxt(k,466)* y(k,90) +rxt(k,467)* y(k,124) + het_rates(k,235)) &
1299 315224452 : * y(k,235)
1300 315224452 : prod(k,166) = (rxt(k,468)*y(k,182) +.650_r8*rxt(k,469)*y(k,183))*y(k,228)
1301 : loss(k,52) = (rxt(k,557)* y(k,90) +rxt(k,558)* y(k,124) + het_rates(k,236)) &
1302 315224452 : * y(k,236)
1303 315224452 : prod(k,52) =rxt(k,556)*y(k,228)*y(k,183)
1304 : loss(k,207) = (rxt(k,432)* y(k,90) +rxt(k,433)* y(k,124) +rxt(k,434) &
1305 : * y(k,126) +rxt(k,430)* y(k,200) +rxt(k,431)* y(k,201) &
1306 315224452 : + het_rates(k,237))* y(k,237)
1307 : prod(k,207) = (rxt(k,402)*y(k,95) +rxt(k,403)*y(k,97) +rxt(k,404)*y(k,98) + &
1308 315224452 : .400_r8*rxt(k,427)*y(k,103) +.500_r8*rxt(k,435)*y(k,184))*y(k,228)
1309 : loss(k,167) = (rxt(k,472)* y(k,90) +rxt(k,473)* y(k,124) + het_rates(k,238)) &
1310 315224452 : * y(k,238)
1311 315224452 : prod(k,167) = (.560_r8*rxt(k,471)*y(k,185) +rxt(k,474)*y(k,186))*y(k,228)
1312 : loss(k,53) = (rxt(k,560)* y(k,90) +rxt(k,561)* y(k,124) + het_rates(k,239)) &
1313 315224452 : * y(k,239)
1314 315224452 : prod(k,53) =rxt(k,559)*y(k,228)*y(k,185)
1315 : loss(k,137) = (rxt(k,475)* y(k,90) +rxt(k,476)* y(k,124) + het_rates(k,240)) &
1316 315224452 : * y(k,240)
1317 315224452 : prod(k,137) = (.300_r8*rxt(k,477)*y(k,187) +rxt(k,478)*y(k,188))*y(k,228)
1318 : loss(k,239) = (rxt(k,273)* y(k,73) +rxt(k,518)* y(k,154) +rxt(k,154) &
1319 : * y(k,224) + rxt(k,1) + rxt(k,2) + rxt(k,3) + het_rates(k,241)) &
1320 315224452 : * y(k,241)
1321 : prod(k,239) = (rxt(k,179)*y(k,77) +rxt(k,180)*y(k,79) +rxt(k,181)*y(k,90) + &
1322 : rxt(k,184)*y(k,228) +rxt(k,187)*y(k,91) +rxt(k,211)*y(k,89) + &
1323 : rxt(k,235)*y(k,85) +rxt(k,238)*y(k,93) +rxt(k,264)*y(k,81) + &
1324 : rxt(k,278)*y(k,41) +rxt(k,280)*y(k,43) +rxt(k,281)*y(k,44) + &
1325 : rxt(k,283)*y(k,46) +rxt(k,288)*y(k,84) +rxt(k,297)*y(k,42) + &
1326 : rxt(k,303)*y(k,53) +rxt(k,304)*y(k,54) +rxt(k,306)*y(k,87) + &
1327 : rxt(k,326)*y(k,28) +rxt(k,328)*y(k,45) +rxt(k,334)*y(k,50) + &
1328 : rxt(k,335)*y(k,51) +rxt(k,351)*y(k,30) +rxt(k,352)*y(k,31) + &
1329 : rxt(k,354)*y(k,49) +rxt(k,360)*y(k,148) +rxt(k,364)*y(k,149) + &
1330 : rxt(k,366)*y(k,48) +.500_r8*rxt(k,379)*y(k,105) +rxt(k,519)*y(k,120)) &
1331 : *y(k,228) + (rxt(k,563)*y(k,93) +rxt(k,569)*y(k,93) + &
1332 : rxt(k,570)*y(k,92) +rxt(k,574)*y(k,93) +rxt(k,575)*y(k,92))*y(k,85) &
1333 : + (rxt(k,521) +rxt(k,174)*y(k,76))*y(k,90) +.050_r8*rxt(k,40) &
1334 325174377 : *y(k,54) +rxt(k,136)*y(k,80)
1335 : end do
1336 9949925 : end subroutine imp_prod_loss
1337 : end module mo_prod_loss
|