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