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