********************************************************************* ;
* A citation to the MCM website and the relevant mechanism          * ;
* construction protocols should be given in any publication using   * ;
* information obtained from this source, using the following or     * ;
* comparable wording:                                               * ;
* The chemical mechanistic information was taken from the Master    * ;
* Chemical Mechanism, MCM v3.2 (ref), via website:                  * ;
* http://mcm.leeds.ac.uk/MCM.                                       * ;
* The reference should be: (Jenkin et al., Atmos. Environ., 31, 81, * ;
* 1997; Saunders et al., Atmos. Chem. Phys., 3, 161, 2003), for     * ;
* non aromatic schemes; (Jenkin et al., Atmos. Chem. Phys., 3,  * ;
* 181, 2003; Bloss et al., Atmos. Chem. Phys., 5, 641, 2005), for   * ;
* aromatic schemes; and (Jenkin et al., Atmos. Chem. Phys.,  12, * ;
* 5275, 2012), for the beta-caryophyllene scheme.                  * ;
********************************************************************* ;
Output file for Fortran Subset Extractor


***********************************************************************
Species present in subset
1 	 C4PAN5 
2 	 C4PAN6 
3 	 ETHENO3O 
4 	 NC3OO 
5 	 CHOPRNO3 
6 	 CHOCOHCO 
7 	 C524CO 
8 	 ISOPDOOH 
9 	 IPRHOCO2H 
10 	 ETHGLY 
11 	 MACROH 
12 	 MACROO 
13 	 ETHO2HNO3 
14 	 HMGLOOA 
15 	 INB1NBCO2H 
16 	 MACRO2 
17 	 HMACO3H 
18 	 ACO3H 
19 	 ACO3B 
20 	 INAHCO3 
21 	 C525OOH 
22 	 INAHCHO 
23 	 ISOPAOOH 
24 	 IPROPOLO 
25 	 ISOPCOOH 
26 	 NC4OOA 
27 	 INB1NACO3H 
28 	 IEACO3H 
29 	 C3DIOLOOH 
30 	 INAHPCO3 
31 	 INDOOH 
32 	 NO2 
33 	 NO3 
34 	 VGLYOX 
35 	 C58OOH 
36 	 H2O2 
37 	 C58NO3 
38 	 C59O 
39 	 INCCO 
40 	 PAN 
41 	 C524O 
42 	 INB1GLYOX 
43 	 PRONO3AO2 
44 	 CONM2CO3H 
45 	 IEB1O2 
46 	 CO2H3CO3H 
47 	 CH3CHO 
48 	 PRONO3BO 
49 	 NOA 
50 	 INCOOH 
51 	 ISOPBOOH 
52 	 INAOOH 
53 	 IEPOXC 
54 	 C524NO3 
55 	 IEPOXA 
56 	 MACROOA 
57 	 C58NO3CO3 
58 	 MACROOH 
59 	 INAHPCO3H 
60 	 CO2N3CHO 
61 	 C510OH 
62 	 INANCOCO2H 
63 	 CO2N3CO3 
64 	 PHAN 
65 	 HOCH2COCHO 
66 	 HC4CCO2H 
67 	 HMVKAO 
68 	 MACO2H 
69 	 MGLOOA 
70 	 ACRPAN 
71 	 HMACROH 
72 	 ISOPCNO3 
73 	 MVK 
74 	 H14CO23C4 
75 	 C510O2 
76 	 INDO2 
77 	 INB1NACHO 
78 	 MACRNBPAN 
79 	 COHM2CO3H 
80 	 NC524OOH 
81 	 HOCH2CHO 
82 	 ISOPANO3 
83 	 INANCOCHO 
84 	 PRNO3CO3H 
85 	 INB1OOH 
86 	 INANCO3H 
87 	 MACRNCO2H 
88 	 INAOH 
89 	 C3DIOLO 
90 	 PRONO3AO 
91 	 IECCO3 
92 	 PRNO3PAN 
93 	 HYPROPO2H 
94 	 ISOPAO 
95 	 INCOH 
96 	 INCO 
97 	 MMALNBCO3H 
98 	 IPRHOCO3 
99 	 INCO2 
100 	 MACRO 
101 	 INANCOCO3 
102 	 IEB2OOH 
103 	 MGLYOOB 
104 	 MGLYOOA 
105 	 INCNO3 
106 	 INB1O 
107 	 CH3CHOO 
108 	 INB2OOH 
109 	 INCNCHO 
110 	 INANPAN 
111 	 INAO2 
112 	 INCNCO3 
113 	 ISOPDNO3 
114 	 IBUTALOH 
115 	 C32OH13CO 
116 	 CH3CO3H 
117 	 HIEB1O2 
118 	 C57O 
119 	 CH3C2H2O2 
120 	 MVKOHAO 
121 	 H1CO23CHO 
122 	 INANO3 
123 	 INB1HPCHO 
124 	 ACO2H 
125 	 CO2N3PAN 
126 	 ISOPBNO3 
127 	 CH3CHOOA 
128 	 MMALNACO3 
129 	 IECCO3H 
130 	 HC4ACHO 
131 	 CH4 
132 	 HMGLYOO 
133 	 HCOCO3 
134 	 NC4CO3 
135 	 HCOCH2O2 
136 	 MVKO 
137 	 MMALNACO3H 
138 	 ACLOO 
139 	 IEB2O2 
140 	 NC524NO3 
141 	 CONM2CO3 
142 	 PROPALO 
143 	 PRONO3BO2 
144 	 HCOOH 
145 	 C58NO3PAN 
146 	 HCOCO2H 
147 	 ETHENO3O2 
148 	 HIEPOXB 
149 	 HCHO 
150 	 HMVKAO2 
151 	 CONM2CO2H 
152 	 MACRNO3 
153 	 HMACRO2 
154 	 IEB1O 
155 	 IPROPOLPER 
156 	 ALLYLOH 
157 	 IECCHO 
158 	 HOCH2CO3 
159 	 MACROHNO3 
160 	 HMVKNO3 
161 	 CH2OOG 
162 	 CH2OOE 
163 	 A2PAN 
164 	 CH2OOC 
165 	 CH2OOB 
166 	 CH2OOA 
167 	 CONM2CHO 
168 	 CH3O 
169 	 HOCH2CO2H 
170 	 MACRNB 
171 	 IPRHOCO3H 
172 	 INB1NBCHO 
173 	 HYPERACET 
174 	 MACRNOO 
175 	 C510O 
176 	 CH3OOH 
177 	 OCCOHCOOH 
178 	 PROPGLY 
179 	 IEC1OOH 
180 	 NC4CHO 
181 	 HMVKNGLYOX 
182 	 C57NO3CO3H 
183 	 CH3CHOHCHO 
184 	 NISOPNO3 
185 	 IEC2OOH 
186 	 MVKOHAOH 
187 	 MVKOHBO2 
188 	 MGLOO 
189 	 INDHPCO3 
190 	 INDOH 
191 	 C524OH 
192 	 MGLYOO 
193 	 C58NO3CO2H 
194 	 C524O2 
195 	 NOAOOA 
196 	 H13CO2CO3H 
197 	 INB1NACO2H 
198 	 CO2H3CHO 
199 	 IEAPAN 
200 	 MGLYOX 
201 	 MVKOHAOOH 
202 	 MVKOHAO2 
203 	 INANCOPAN 
204 	 HC4CCHO 
205 	 INB1NO3 
206 	 HOCHOCOOH 
207 	 CH3O2NO2 
208 	 MACRNBCO3H 
209 	 ISOPDO 
210 	 C57NO3 
211 	 C525O 
212 	 MMALNACO2H 
213 	 MMALNBCO2H 
214 	 GAOO 
215 	 MACRNBCO3 
216 	 IEACO3 
217 	 HMGLOO 
218 	 CH3COCO2H 
219 	 IEPOXB 
220 	 C57O2 
221 	 INCGLYOX 
222 	 DNC524CO 
223 	 OCCOHCOH 
224 	 IECPAN 
225 	 HIEB1O 
226 	 C2OHOCOOH 
227 	 INB1HPCO3H 
228 	 C510OOH 
229 	 C5H8 
230 	 OCCOHCO2 
231 	 C59O2 
232 	 ACETOL 
233 	 INB1HPCO2H 
234 	 IPROPOLO2 
235 	 CONM2PAN 
236 	 CH3CO3 
237 	 INB1HPPAN 
238 	 C57OH 
239 	 MMALNBPAN 
240 	 HIEB2O2 
241 	 HPNC524CO 
242 	 CH3COCH2O2 
243 	 MACO3H 
244 	 HMVKBO 
245 	 GAOOB 
246 	 CH3NO3 
247 	 ISOPAOH 
248 	 HMVKANO3 
249 	 MACROHO2 
250 	 HMVKBO2 
251 	 HYPROPO 
252 	 INB1NBCO3H 
253 	 H13CO2CHO 
254 	 COHM2CO3 
255 	 HCOCO3H 
256 	 HMGLYOOA 
257 	 C57NO3CO3 
258 	 ISOPCO2 
259 	 INANCOCO3H 
260 	 IPROPOLPAN 
261 	 HNO3 
262 	 MACRNCO3 
263 	 HYETHO2H 
264 	 ETHOHNO3 
265 	 GLYOO 
266 	 ISOPAO2 
267 	 A2PANOO 
268 	 HIEB1OOH 
269 	 GLYOX 
270 	 HC4CCO3 
271 	 HMACR 
272 	 INCNCO2H 
273 	 C3H6 
274 	 CO2N3CO3H 
275 	 HC4ACO3 
276 	 MACROHOOH 
277 	 CO23C4NO3 
278 	 C57OOH 
279 	 CH2OO 
280 	 COHM2PAN 
281 	 CH3CHOHCO3 
282 	 COHM2CO2H 
283 	 INB1HPCO3 
284 	 C5PAN19 
285 	 C5PAN18 
286 	 C5PAN17 
287 	 C525O2 
288 	 C2H4 
289 	 HCOCOHCO3 
290 	 INB1NBCO3 
291 	 HCOC5 
292 	 NISOPOOH 
293 	 MACRNOOA 
294 	 HO12CO3C4 
295 	 HCOCH2O 
296 	 CH3COCO3 
297 	 H13CO2CO3 
298 	 INDHCO3H 
299 	 INANCO 
300 	 IEC1O2 
301 	 BIACETO 
302 	 NO3CH2CO3 
303 	 C3DIOLO2 
304 	 NO3CH2CHO 
305 	 INB2O2 
306 	 MVKOHANO3 
307 	 ISOPCO 
308 	 PROPOLNO3 
309 	 CO 
310 	 CL 
311 	 INAHPPAN 
312 	 HOCH2CH2O 
313 	 C42AOH 
314 	 IEB2O 
315 	 IEB1OOH 
316 	 HC4ACO2H 
317 	 HIEB2OOH 
318 	 HMACO2H 
319 	 HNC524CO 
320 	 HCOCH2OOH 
321 	 MVKNO3 
322 	 C2OHOCO2H 
323 	 ISOPBOH 
324 	 MPAN 
325 	 ISOPDOH 
326 	 C58O2 
327 	 INDHPPAN 
328 	 A2PANO 
329 	 NO3CH2CO3H 
330 	 C58OH 
331 	 ISOPDO2 
332 	 SO3 
333 	 SO2 
334 	 IEC1O 
335 	 ISOPBO2 
336 	 HMACRO 
337 	 C57NO3CO2H 
338 	 HOCH2CO3H 
339 	 PR2O2HNO3 
340 	 HMVKAOOH 
341 	 INAHCO3H 
342 	 NC4CO3H 
343 	 HCOCOHCO3H 
344 	 C58NO3CO3H 
345 	 CH3OH 
346 	 INDHCO3 
347 	 INDHCHO 
348 	 HMACO3 
349 	 O3 
350 	 INDHPCHO 
351 	 CH3COCH2O 
352 	 INB1NAPAN 
353 	 ACLOOA 
354 	 IPROPOLO2H 
355 	 INAO 
356 	 INAHPCHO 
357 	 INANCHO 
358 	 INCNPAN 
359 	 NC524O 
360 	 MVKOHBO 
361 	 HCOCOHPAN 
362 	 CH3COCO3H 
363 	 OH 
364 	 GLYOOB 
365 	 GLYOOC 
366 	 CH3O2 
367 	 HMVKBOOH 
368 	 BIACETOH 
369 	 HIEB2O 
370 	 PR1O2HNO3 
371 	 HYPROPO2 
372 	 MVKO2 
373 	 C59OOH 
374 	 NC4OO 
375 	 HOCH2CH2O2 
376 	 IEACHO 
377 	 INCNCO3H 
378 	 MACROHO 
379 	 H13CO2C3 
380 	 INB1NACO3 
381 	 MVKOH 
382 	 NC3OOA 
383 	 ACO3 
384 	 C57NO3PAN 
385 	 PROLNO3 
386 	 MVKOO 
387 	 INB1O2 
388 	 H2 
389 	 CHOMOHCO3 
390 	 C58O 
391 	 NC2OO 
392 	 HC4CCO3H 
393 	 OCCOHCO 
394 	 NMGLYOX 
395 	 CH3COPAN 
396 	 INB2O 
397 	 CHOMOHPAN 
398 	 MMALNBCO3 
399 	 HOCH2COCO2H 
400 	 ACR 
401 	 INDHPAN 
402 	 INB1NBPAN 
403 	 NO 
404 	 HMACROOH 
405 	 HC4ACO3H 
406 	 NC524OH 
407 	 CH3CO2H 
408 	 MVKOOH 
409 	 NC524O2 
410 	 C4M2ALOHNO3 
411 	 MVKOOA 
412 	 MMALNAPAN 
413 	 NISOPO 
414 	 GLYPAN 
415 	 INAHPCO2H 
416 	 CO23C3CHO 
417 	 CH3COCH3 
418 	 PRNO3CO3 
419 	 MACRNBCO2H 
420 	 NOAOO 
421 	 CO2H3CO3 
422 	 INANCO3 
423 	 NO3CH2CO2H 
424 	 INB1OH 
425 	 MVKOHBOOH 
426 	 INAHCO2H 
427 	 CHOMOHCO3H 
428 	 INAHPAN 
429 	 HO2 
430 	 NISOPO2 
431 	 NC4CO2H 
432 	 PRNO3CO2H 
433 	 MACR 
434 	 INANCO2H 
435 	 MACO3 
436 	 MACRNCO3H 
437 	 INB1CO 
438 	 NC2OOA 
439 	 C4PAN10 
440 	 NO3CH2PAN 
441 	 C524OOH 
442 	 INDHPCO3H 
443 	 MACRNPAN 
444 	 HMPAN 
445 	 C3MDIALOH 
446 	 INDO 
447 	 ISOPBO 

***********************************************************************


 Rate expressions fortran routine

subroutine constants(p,t,y)  
implicit none  
!     calculates rate constants from arrhenius informtion 
integer    mnsp,mre
parameter (mnsp = 3487,mre=20000)
double precision:: ro2 , kro2no3, temp, p(mre),t
double precision:: hmix, area, gn2o5, gho2, gch3o2, goh,ghno3, gh2o2, gch3oh, gno3
double precision:: m, o2, n2, h2o, k0, ki, fc, f, k1, k2, k3, k4, kmt01, kmt02 
double precision:: kmt03, kmt04, kmt05, kmt06, kmt07, kmt08, kmt09,kmt10,kmt11 
double precision:: kmt12, kmt13, kmt14, kmt15, kmt16, kmt17, kfpan, kbpan 
double precision:: fcc,krc,fcd,krd,fc2,fc1,fc3,fc4,fc5,kr1,kr2,kr3,kr4
double precision:: fc7,fc8,fc9,fc10,fc13,fc14,kr7,kr8,kr9,kr10,kr13,kr14
double precision:: kc0, kci,kd0,kdi,fd,k10,k1i,f1,k20,k2i,f2,k30,k3i,f3
double precision:: k40,k4i,f4,k70,k7i,f7,ck80,ck8i,ckr8,cfc8,k80,k8i,f8
double precision:: k90,k9i,f9,k100,k10i,f10,k130,k13i,f13,k140,k14i,f14
double precision:: k160,k16i,kr16,fc16,f16, ssign,x
double precision:: y(mnsp)
open(13,file='photolysis.txt', status='old')    
!  calculate zenith angle and temperature for time of day

call zenith(theta, secx, cosx, t)

call temperature(temp, h2o, t)

! ***************************************************************************
!    define ro2 concentration by adding up concentrations of all 
!    contributing species. 
! ***************************************************************************

ro2 = 0.0d+0


	 ro2=ro2+y(430) 
	 ro2=ro2+y(266) 
	 ro2=ro2+y(335) 
	 ro2=ro2+y(258) 
	 ro2=ro2+y(331) 
	 ro2=ro2+y(119) 
	 ro2=ro2+y(435) 
	 ro2=ro2+y(16) 
	 ro2=ro2+y(249) 
	 ro2=ro2+y(236) 
	 ro2=ro2+y(150) 
	 ro2=ro2+y(250) 
	 ro2=ro2+y(366) 
	 ro2=ro2+y(372) 
	 ro2=ro2+y(99) 
	 ro2=ro2+y(134) 
	 ro2=ro2+y(75) 
	 ro2=ro2+y(43) 
	 ro2=ro2+y(143) 
	 ro2=ro2+y(371) 
	 ro2=ro2+y(234) 
	 ro2=ro2+y(111) 
	 ro2=ro2+y(194) 
	 ro2=ro2+y(275) 
	 ro2=ro2+y(326) 
	 ro2=ro2+y(387) 
	 ro2=ro2+y(305) 
	 ro2=ro2+y(270) 
	 ro2=ro2+y(220) 
	 ro2=ro2+y(76) 
	 ro2=ro2+y(158) 
	 ro2=ro2+y(231) 
	 ro2=ro2+y(133) 
	 ro2=ro2+y(389) 
	 ro2=ro2+y(135) 
	 ro2=ro2+y(421) 
	 ro2=ro2+y(383) 
	 ro2=ro2+y(202) 
	 ro2=ro2+y(187) 
	 ro2=ro2+y(302) 
	 ro2=ro2+y(45) 
	 ro2=ro2+y(139) 
	 ro2=ro2+y(262) 
	 ro2=ro2+y(300) 
	 ro2=ro2+y(98) 
	 ro2=ro2+y(418) 
	 ro2=ro2+y(281) 
	 ro2=ro2+y(216) 
	 ro2=ro2+y(409) 
	 ro2=ro2+y(287) 
	 ro2=ro2+y(57) 
	 ro2=ro2+y(141) 
	 ro2=ro2+y(91) 
	 ro2=ro2+y(63) 
	 ro2=ro2+y(112) 
	 ro2=ro2+y(215) 
	 ro2=ro2+y(242) 
	 ro2=ro2+y(297) 
	 ro2=ro2+y(19) 
	 ro2=ro2+y(230) 
	 ro2=ro2+y(30) 
	 ro2=ro2+y(422) 
	 ro2=ro2+y(20) 
	 ro2=ro2+y(117) 
	 ro2=ro2+y(240) 
	 ro2=ro2+y(348) 
	 ro2=ro2+y(153) 
	 ro2=ro2+y(283) 
	 ro2=ro2+y(380) 
	 ro2=ro2+y(290) 
	 ro2=ro2+y(257) 
	 ro2=ro2+y(189) 
	 ro2=ro2+y(346) 
	 ro2=ro2+y(147) 
	 ro2=ro2+y(375) 
	 ro2=ro2+y(101) 
	 ro2=ro2+y(128) 
	 ro2=ro2+y(398) 
	 ro2=ro2+y(267) 
	 ro2=ro2+y(289) 
	 ro2=ro2+y(296) 
	 ro2=ro2+y(254) 
	 ro2=ro2+y(303) 

c 	 WARNING: The following species do not have SMILES strings in the database.
c 	 If any of these are peroxy radicals the RO2 sum will be wrong. 
! ***************************************************************************
m=2.55e+19
o2=0.2095*m 

***********************************************************************
c 	 Generic Rate Coefficients 

***********************************************************************
	KRO2NO = 2.7D-12*EXP(360/TEMP)
	KRO2HO2 = 2.91D-13*EXP(1300/TEMP)
	KAPHO2 = 5.2D-13*EXP(980/TEMP)
	KAPNO = 7.5D-12*EXP(290/TEMP)
	KRO2NO3 = 2.3D-12
	KNO3AL = 1.4D-12*EXP(-1860/TEMP)
	KDEC = 1.00D+06
	KROPRIM = 2.50D-14*EXP(-300/TEMP)
	KROSEC = 2.50D-14*EXP(-300/TEMP)
	KCH3O2 = 1.03D-13*EXP(365/TEMP)
	K298CH3O2 = 3.5D-13

***********************************************************************
c 	  Complex reactions 

***********************************************************************
	KD0 = 4.90D-3*EXP(-12100/TEMP)*M
	KDI = 5.4D+16*EXP(-13830/TEMP)
	KRD = KD0/KDI
	FCD = 0.30
	NCD = 0.75-1.27*(LOG10(FCD))
	FD = 10**(LOG10(FCD)/(1+(LOG10(KRD)/NCD)**2))
	KBPAN = (KD0*KDI)*FD/(KD0+KDI)
	KC0 = 2.7D-28*M*(TEMP/300)**-7.1
	KCI = 1.2D-11*(TEMP/300)**-0.9
	KRC = KC0/KCI
	FCC = 0.30
	NC = 0.75-1.27*(LOG10(FCC))
	FC = 10**(LOG10(FCC)/(1+(LOG10(KRC)/NC)**2))
	KFPAN = (KC0*KCI)*FC/(KC0+KCI)
	K10 = 1.0D-31*M*(TEMP/300)**-1.6
	K1I = 3.00D-11*(TEMP/300)**0.3
	KR1 = K10/K1I
	FC1 = 0.85
	NC1 = 0.75-1.27*(LOG10(FC1))
	F1 = 10**(LOG10(FC1)/(1+(LOG10(KR1)/NC1)**2))
	KMT01 = (K10*K1I)*F1/(K10+K1I)
	K20 = 1.3D-31*M*(TEMP/300)**-1.5
	K2I = 2.3D-11*(TEMP/300)**0.24
	KR2 = K20/K2I
	FC2 = 0.6
	NC2 = 0.75-1.27*(LOG10(FC2))
	F2 = 10**(LOG10(FC2)/(1+(LOG10(KR2)/NC2)**2))
	KMT02 = (K20*K2I)*F2/(K20+K2I)
	K30 = 3.6D-30*M*(TEMP/300)**-4.1
	K3I = 1.9D-12*(TEMP/300)**0.2
	KR3 = K30/K3I
	FC3 = 0.35
	NC3 = 0.75-1.27*(LOG10(FC3))
	F3 = 10**(LOG10(FC3)/(1+(LOG10(KR3)/NC3)**2))
	KMT03 = (K30*K3I)*F3/(K30+K3I)
	K40 = 1.3D-3*M*(TEMP/300)**-3.5*EXP(-11000/TEMP)
	K4I = 9.7D+14*(TEMP/300)**0.1*EXP(-11080/TEMP)
	KR4 = K40/K4I
	FC4 = 0.35
	NC4 = 0.75-1.27*(LOG10(FC4))
	F4 = 10**(LOG10(FC4)/(1+(LOG10(KR4)/NC4)**2))
	KMT04 = (K40*K4I)*F4/(K40+K4I)
	KMT05 = 1.44D-13*(1+(M/4.2D+19))
	KMT06 = 1 + (1.40D-21*EXP(2200/TEMP)*H2O)
	K70 = 7.4D-31*M*(TEMP/300)**-2.4
	K7I = 3.3D-11*(TEMP/300)**-0.3
	KR7 = K70/K7I
	FC7 = EXP(-TEMP/1420)
	NC7 = 0.75-1.27*(LOG10(FC7))
	F7 = 10**(LOG10(FC7)/(1+(LOG10(KR7)/NC7)**2))
	KMT07 = (K70*K7I)*F7/(K70+K7I)
	K80 = 3.3D-30*M*(TEMP/300)**-3.0
	K8I = 4.1D-11
	KR8 = K80/K8I
	FC8 = 0.4
	NC8 = 0.75-1.27*(LOG10(FC8))
	F8 = 10**(LOG10(FC8)/(1+(LOG10(KR8)/NC8)**2))
	KMT08 = (K80*K8I)*F8/(K80+K8I)
	K90 = 1.8D-31*M*(TEMP/300)**-3.2
	K9I = 4.7D-12
	KR9 = K90/K9I
	FC9 = 0.6
	NC9 = 0.75-1.27*(LOG10(FC9))
	F9 = 10**(LOG10(FC9)/(1+(LOG10(KR9)/NC9)**2))
	KMT09 = (K90*K9I)*F9/(K90+K9I)
	K100 = 4.10D-05*M*EXP(-10650/TEMP)
	K10I = 4.8D+15*EXP(-11170/TEMP)
	KR10 = K100/K10I
	FC10 = 0.6
	NC10 = 0.75-1.27*(LOG10(FC10))
	F10 = 10**(LOG10(FC10)/(1+(LOG10(KR10)/NC10)**2))
	KMT10 = (K100*K10I)*F10/(K100+K10I)
	K1 = 2.40D-14*EXP(460/TEMP)
	K3 = 6.50D-34*EXP(1335/TEMP)
	K4 = 2.70D-17*EXP(2199/TEMP)
	K2 = (K3*M)/(1+(K3*M/K4))
	KMT11 = K1 + K2
	K120 = 4.5D-31*M*(TEMP/300)**-3.9
	K12I = 1.3D-12*(TEMP/300)**-0.7
	KR12 = K120/K12I
	FC12 = 0.525
	NC12 = 0.75-1.27*(LOG10(FC12))
	F12 = 10**(LOG10(FC12)/(1.0+(LOG10(KR12)/NC12)**2))
	KMT12 = (K120*K12I*F12)/(K120+K12I)
	K130 = 2.5D-30*M*(TEMP/300)**-5.5
	K13I = 1.8D-11
	KR13 = K130/K13I
	FC13 = 0.36
	NC13 = 0.75-1.27*(LOG10(FC13))
	F13 = 10**(LOG10(FC13)/(1+(LOG10(KR13)/NC13)**2))
	KMT13 = (K130*K13I)*F13/(K130+K13I)
	K140 = 9.0D-5*EXP(-9690/TEMP)*M
	K14I = 1.1D+16*EXP(-10560/TEMP)
	KR14 = K140/K14I
	FC14 = 0.4
	NC14 = 0.75-1.27*(LOG10(FC14))
	F14 = 10**(LOG10(FC14)/(1+(LOG10(KR14)/NC14)**2))
	KMT14 = (K140*K14I)*F14/(K140+K14I)
	K150 = 8.6D-29*M*(TEMP/300)**-3.1
	K15I = 9.0D-12*(TEMP/300)**-0.85
	KR15 = K150/K15I
	FC15 = 0.48
	NC15 = 0.75-1.27*(LOG10(FC15))
	F15 = 10**(LOG10(FC15)/(1+(LOG10(KR15)/NC15)**2))
	KMT15 = (K150*K15I)*F15/(K150+K15I)
	K160 = 8D-27*M*(TEMP/300)**-3.5
	K16I = 3.0D-11*(TEMP/300)**-1
	KR16 = K160/K16I
	FC16 = 0.5
	NC16 = 0.75-1.27*(LOG10(FC16))
	F16 = 10**(LOG10(FC16)/(1+(LOG10(KR16)/NC16)**2))
	KMT16 = (K160*K16I)*F16/(K160+K16I)
	K170 = 5.0D-30*M*(TEMP/300)**-1.5
	K17I = 1.0D-12
	KR17 = K170/K17I
	FC17 = 0.17*EXP(-51/TEMP)+EXP(-TEMP/204)
	NC17 = 0.75-1.27*(LOG10(FC17))
	F17 = 10**(LOG10(FC17)/(1.0+(LOG10(KR17)/NC17)**2))
	KMT17 = (K170*K17I*F17)/(K170+K17I)
	KMT18 = 9.5D-39*O2*EXP(5270/TEMP)/(1+7.5D-29*O2*EXP(5610/TEMP))

***********************************************************************
! ***************************************************************************

!    define photolysis reaction rates using derwent method from mcm2box.fac

! ***************************************************************************
rewind(13)
read(13,*) 
do i = 1, 35
	read(13,*) k, l(k), mm(k), nn(k) 
	if (cosx.lt.1.00d-10) then
		j(k) = 1.0d-30
	else
		j(k) = l(k)*cosx**( mm(k))*dexp(-nn(k)*secx) 
	endif
enddo

	p(1)= 3.15D-12*EXP(-450/TEMP)
	p(2)= 1.03D-14*EXP(-1995/TEMP)*0.3
	p(3)= 1.03D-14*EXP(-1995/TEMP)*0.2
	p(4)= 1.03D-14*EXP(-1995/TEMP)*0.3
	p(5)= 1.03D-14*EXP(-1995/TEMP)*0.2
	p(6)= 2.7D-11*EXP(390/TEMP)*0.148
	p(7)= 2.7D-11*EXP(390/TEMP)*0.444
	p(8)= 2.7D-11*EXP(390/TEMP)*0.102
	p(9)= 2.7D-11*EXP(390/TEMP)*0.306
	p(10)= KRO2HO2*0.706
	p(11)= KRO2NO*0.052
	p(12)= KRO2NO*0.948
	p(13)= KRO2NO3
	p(14)= 1.30D-12*0.2*RO2
	p(15)= 1.30D-12*0.2*RO2
	p(16)= 1.30D-12*0.6*RO2
	p(17)= KDEC*0.22
	p(18)= KDEC*0.51
	p(19)= KDEC*0.27
	p(20)= J(18)
	p(21)= J(19)
	p(22)= 3.4D-15
	p(23)= 1.4D-15*EXP(-2100/TEMP)*0.12
	p(24)= 1.4D-15*EXP(-2100/TEMP)*0.88
	p(25)= 8.0D-12*EXP(380/TEMP)*0.45
	p(26)= 8.0D-12*EXP(380/TEMP)*0.47
	p(27)= 8.0D-12*EXP(380/TEMP)*0.08
	p(28)= J(23)
	p(29)= J(24)
	p(30)= 8.5D-16*EXP(-1520/TEMP)*0.5
	p(31)= 8.5D-16*EXP(-1520/TEMP)*0.5
	p(32)= 2.6D-12*EXP(610/TEMP)*0.3
	p(33)= 2.6D-12*EXP(610/TEMP)*0.7
	p(34)= J(11)
	p(35)= J(12)
	p(36)= 5.5D-16
	p(37)= 5.4D-12*EXP(135/TEMP)
	p(38)= KDEC*0.255
	p(39)= KDEC*0.255
	p(40)= KDEC*0.22
	p(41)= KDEC*0.27
	p(42)= KDEC*0.255
	p(43)= KDEC*0.255
	p(44)= KDEC*0.22
	p(45)= KDEC*0.27
	p(46)= KRO2HO2*0.706
	p(47)= KRO2NO*0.100
	p(48)= KRO2NO*0.900
	p(49)= KRO2NO3
	p(50)= 2.40D-12*0.1*RO2
	p(51)= 2.40D-12*0.8*RO2
	p(52)= 2.40D-12*0.1*RO2
	p(53)= KRO2HO2*0.706
	p(54)= KRO2NO*0.066
	p(55)= KRO2NO*0.934
	p(56)= KRO2NO3
	p(57)= 8.00D-13*0.8*RO2
	p(58)= 8.00D-13*0.2*RO2
	p(59)= KRO2HO2*0.706
	p(60)= KRO2NO*0.100
	p(61)= KRO2NO*0.900
	p(62)= KRO2NO3
	p(63)= 2.00D-12*0.1*RO2
	p(64)= 2.00D-12*0.1*RO2
	p(65)= 2.00D-12*0.8*RO2
	p(66)= KRO2HO2*0.706
	p(67)= KRO2NO*0.134
	p(68)= KRO2NO*0.866
	p(69)= KRO2NO3
	p(70)= 2.90D-12*0.1*RO2
	p(71)= 2.90D-12*0.8*RO2
	p(72)= 2.90D-12*0.1*RO2
	p(73)= J(41)
	p(74)= 1.03D-10
	p(75)= J(53)*2.0
	p(76)= 8.55D-11
	p(77)= KROPRIM*O2
	p(78)= J(53)
	p(79)= 5.30D-17*0.5
	p(80)= 5.30D-17*0.5
	p(81)= 6.93D-11
	p(82)= J(18)
	p(83)= KNO3AL*4.25
	p(84)= 2.40D-17*0.5
	p(85)= 2.40D-17*0.5
	p(86)= 4.16D-11*0.52
	p(87)= 4.16D-11*0.48
	p(88)= 1.20D-15
	p(89)= 1.00D-14
	p(90)= 1.00D-15
	p(91)= 7.00D-14
	p(92)= 6.00D-18*H2O
	p(93)= 1.00D-17*H2O
	p(94)= KDEC*0.35
	p(95)= KDEC*0.65
	p(96)= KAPHO2*0.44
	p(97)= KAPHO2*0.15
	p(98)= KAPHO2*0.41
	p(99)= 8.70D-12*EXP(290/TEMP)
	p(100)= KFPAN
	p(101)= KRO2NO3*1.74
	p(102)= 1.00D-11*0.7*RO2
	p(103)= 1.00D-11*0.3*RO2
	p(104)= KDEC*0.18
	p(105)= KDEC*0.82
	p(106)= J(34)
	p(107)= KNO3AL*2.4
	p(108)= 1.9D-12*EXP(575/TEMP)
	p(109)= KDEC*0.37
	p(110)= KDEC*0.47
	p(111)= KDEC*0.16
	p(112)= KRO2HO2*0.625
	p(113)= KRO2NO
	p(114)= KRO2NO3
	p(115)= 9.20D-14*0.7*RO2
	p(116)= 9.20D-14*0.3*RO2
	p(117)= KRO2HO2*0.625
	p(118)= KRO2NO*0.017
	p(119)= KRO2NO*0.983
	p(120)= KRO2NO3
	p(121)= 1.4D-12*0.2*RO2
	p(122)= 1.4D-12*0.2*RO2
	p(123)= 1.4D-12*0.6*RO2
	p(124)= 4.6D-13*EXP(-1155/TEMP)*0.35
	p(125)= 4.6D-13*EXP(-1155/TEMP)*0.65
	p(126)= 5.5D-15*EXP(-1880/TEMP)*0.5
	p(127)= 5.5D-15*EXP(-1880/TEMP)*0.5
	p(128)= KMT16*0.87
	p(129)= KMT16*0.13
	p(130)= KAPHO2*0.15
	p(131)= KAPHO2*0.41
	p(132)= KAPHO2*0.44
	p(133)= 7.5D-12*EXP(290/TEMP)
	p(134)= KFPAN
	p(135)= 4.0D-12
	p(136)= 2*(K298CH3O2*2.9D-12*EXP(500/TEMP))**0.5*RO2*0.3
	p(137)= 2*(K298CH3O2*2.9D-12*EXP(500/TEMP))**0.5*RO2*0.7
	p(138)= KDEC*0.20
	p(139)= KDEC*0.20
	p(140)= KDEC*0.24
	p(141)= KDEC*0.36
	p(142)= KDEC*0.24
	p(143)= KDEC*0.40
	p(144)= KDEC*0.36
	p(145)= KRO2HO2*0.625
	p(146)= KRO2NO*0.017
	p(147)= KRO2NO*0.983
	p(148)= KRO2NO3
	p(149)= 2.00D-12*0.2*RO2
	p(150)= 2.00D-12*0.6*RO2
	p(151)= 2.00D-12*0.2*RO2
	p(152)= KRO2HO2*0.625
	p(153)= KRO2NO
	p(154)= KRO2NO3
	p(155)= 8.80D-13*0.2*RO2
	p(156)= 8.80D-13*0.6*RO2
	p(157)= 8.80D-13*0.2*RO2
	p(158)= 1.2D-15
	p(159)= 1.0D-14
	p(160)= 1.0D-15
	p(161)= 7.0D-14
	p(162)= 1.0D-17*H2O
	p(163)= 6.0D-18*H2O
	p(164)= 3.8D-13*EXP(780/TEMP)*(1-1/(1+498*EXP(-1160/TEMP)))
	p(165)= 3.8D-13*EXP(780/TEMP)*(1/(1+498*EXP(-1160/TEMP)))
	p(166)= 2.3D-12*EXP(360/TEMP)*0.001
	p(167)= 2.3D-12*EXP(360/TEMP)*0.999
	p(168)= KMT13
	p(169)= 1.2D-12
	p(170)= 2*KCH3O2*RO2*7.18*EXP(-885/TEMP)
	p(171)= 2*KCH3O2*RO2*0.5*(1-7.18*EXP(-885/TEMP))
	p(172)= 2*KCH3O2*RO2*0.5*(1-7.18*EXP(-885/TEMP))
	p(173)= 1.2D-15
	p(174)= 1.0D-14
	p(175)= 1.0D-15
	p(176)= 7.0D-14
	p(177)= 6.0D-18*H2O
	p(178)= KRO2HO2*0.625
	p(179)= KRO2NO
	p(180)= KRO2NO3
	p(181)= 2.00D-12*0.6*RO2
	p(182)= 2.00D-12*0.2*RO2
	p(183)= 2.00D-12*0.2*RO2
	p(184)= J(41)
	p(185)= 1.54D-10*0.05
	p(186)= 1.54D-10*0.93
	p(187)= 1.54D-10*0.02
	p(188)= J(53)
	p(189)= 5.30D-17*0.5
	p(190)= 5.30D-17*0.5
	p(191)= 6.93D-11
	p(192)= KDEC*0.25
	p(193)= KDEC*0.75
	p(194)= J(18)
	p(195)= J(19)
	p(196)= KNO3AL*4.25
	p(197)= 2.40D-17*0.5
	p(198)= 2.40D-17*0.5
	p(199)= 4.52D-11*0.52
	p(200)= 4.52D-11*0.48
	p(201)= 9.30D-11*0.5
	p(202)= 9.30D-11*0.5
	p(203)= J(41)
	p(204)= 5.00D-11*0.92
	p(205)= 5.00D-11*0.08
	p(206)= J(55)
	p(207)= 1.06D-16*0.50
	p(208)= 1.06D-16*0.50
	p(209)= 1.36D-11*0.72
	p(210)= 1.36D-11*0.28
	p(211)= KDEC
	p(212)= 3.85D-11
	p(213)= J(41)
	p(214)= 1.54D-10*0.05
	p(215)= 1.54D-10*0.93
	p(216)= 1.54D-10*0.02
	p(217)= KDEC*0.75
	p(218)= KDEC*0.25
	p(219)= J(18)
	p(220)= J(19)
	p(221)= KNO3AL*4.25
	p(222)= 2.40D-17*0.5
	p(223)= 2.40D-17*0.5
	p(224)= 4.52D-11*0.52
	p(225)= 4.52D-11*0.48
	p(226)= J(41)
	p(227)= 1.15D-10*0.22
	p(228)= 1.15D-10*0.75
	p(229)= 1.15D-10*0.03
	p(230)= J(54)
	p(231)= 1.06D-16*0.50
	p(232)= 1.06D-16*0.50
	p(233)= 2.54D-11
	p(234)= KDEC
	p(235)= J(24)
	p(236)= 3.81D-11
	p(237)= 7.38D-11
	p(238)= KDEC*0.11
	p(239)= KDEC*0.89
	p(240)= 1D-12
	p(241)= J(57)
	p(242)= J(56)
	p(243)= KNO3AL
	p(244)= 1.00D-11*0.200
	p(245)= 1.00D-11*0.800
	p(246)= J(15)
	p(247)= KDEC*0.11
	p(248)= KDEC*0.89
	p(249)= KRO2HO2*0.706
	p(250)= KRO2NO*0.145
	p(251)= KRO2NO*0.855
	p(252)= KRO2NO3
	p(253)= 2.90D-12*0.1*RO2
	p(254)= 2.90D-12*0.8*RO2
	p(255)= 2.90D-12*0.1*RO2
	p(256)= KAPHO2*0.15
	p(257)= KAPHO2*0.41
	p(258)= KAPHO2*0.44
	p(259)= KAPNO
	p(260)= KFPAN
	p(261)= KRO2NO3*1.74
	p(262)= 1.00D-11*0.3*RO2
	p(263)= 1.00D-11*0.7*RO2
	p(264)= KDEC*0.11
	p(265)= KDEC*0.89
	p(266)= KDEC*0.11
	p(267)= KDEC*0.89
	p(268)= J(31)
	p(269)= J(33)
	p(270)= J(32)
	p(271)= KNO3AL*0.6
	p(272)= KNO3AL*0.4
	p(273)= 3.1D-12*EXP(340/TEMP)*0.6
	p(274)= 3.1D-12*EXP(340/TEMP)*0.4
	p(275)= KRO2HO2*0.706
	p(276)= KRO2NO
	p(277)= KRO2NO3
	p(278)= 9.20D-14*0.7*RO2
	p(279)= 9.20D-14*0.3*RO2
	p(280)= 4.5D-13
	p(281)= 1.51D-11
	p(282)= J(41)
	p(283)= 1.661D-11
	p(284)= 1.6D+16*EXP(-13500/TEMP)
	p(285)= 8.2D-18
	p(286)= 2.9D-11
	p(287)= 1.2D-15
	p(288)= 1.0D-14
	p(289)= 1.0D-15
	p(290)= 7.0D-14
	p(291)= 6.0D-18*H2O
	p(292)= J(17)
	p(293)= J(41)
	p(294)= 3.77D-11
	p(295)= KDEC
	p(296)= J(17)
	p(297)= 3.42D-11
	p(298)= 5.55D-11
	p(299)= J(17)
	p(300)= J(41)
	p(301)= 2.18D-11
	p(302)= J(55)
	p(303)= J(17)
	p(304)= KDEC
	p(305)= 1.32D-10
	p(306)= J(17)*2
	p(307)= KRO2HO2*0.520
	p(308)= KRO2NO
	p(309)= KRO2NO3
	p(310)= 6.00D-13*0.2*RO2
	p(311)= 6.00D-13*0.6*RO2
	p(312)= 6.00D-13*0.2*RO2
	p(313)= KRO2HO2*0.520
	p(314)= KRO2NO
	p(315)= KRO2NO3
	p(316)= 4.00D-14*0.2*RO2
	p(317)= 4.00D-14*0.2*RO2
	p(318)= 4.00D-14*0.6*RO2
	p(319)= J(13)
	p(320)= 1.4D-12*EXP(-1860/TEMP)
	p(321)= 4.7D-12*EXP(345/TEMP)*0.95
	p(322)= 4.7D-12*EXP(345/TEMP)*0.05
	p(323)= KDEC*0.24
	p(324)= KDEC*0.36
	p(325)= KDEC*0.20
	p(326)= KDEC*0.20
	p(327)= KRO2HO2*0.520
	p(328)= KRO2NO3
	p(329)= 8.80D-13*0.2*RO2
	p(330)= 8.80D-13*0.6*RO2
	p(331)= 8.80D-13*0.2*RO2
	p(332)= KRO2NO*0.977
	p(333)= KRO2NO*0.023
	p(334)= KRO2HO2*0.520
	p(335)= KRO2NO*0.991
	p(336)= KRO2NO*0.009
	p(337)= KRO2NO3
	p(338)= 2.00D-12*0.2*RO2
	p(339)= 2.00D-12*0.6*RO2
	p(340)= 2.00D-12*0.2*RO2
	p(341)= 8.00D-13
	p(342)= 3.70D-12
	p(343)= J(41)
	p(344)= 3D-14
	p(345)= KBPAN
	p(346)= 1.2D-15
	p(347)= 1.0D-14
	p(348)= 1.0D-15
	p(349)= 7.0D-14
	p(350)= 1.0D-17*H2O
	p(351)= 6.0D-18*H2O
	p(352)= J(41)
	p(353)= 5.77D-11
	p(354)= 2.23D-12
	p(355)= KDEC
	p(356)= J(15)
	p(357)= KNO3AL*4.0
	p(358)= 2.45D-11
	p(359)= J(22)
	p(360)= 1.88D-11
	p(361)= J(41)
	p(362)= 3.95D-11
	p(363)= KDEC
	p(364)= J(35)
	p(365)= 2.69D-12
	p(366)= J(41)
	p(367)= 5.3D-12*EXP(190/TEMP)*0.6
	p(368)= 5.3D-12*EXP(190/TEMP)*0.4
	p(369)= J(51)
	p(370)= 4.0D-13*EXP(-845/TEMP)
	p(371)= 7.2D-14*EXP(-1080/TEMP)*O2
	p(372)= KMT14
	p(373)= 2.85D-12*EXP(-345/TEMP)
	p(374)= J(24)
	p(375)= J(41)
	p(376)= 1.90D-12*EXP(190/TEMP)
	p(377)= 2.55D-11
	p(378)= KDEC
	p(379)= 7.51D-16*EXP(-1521/TEMP)*0.5
	p(380)= 7.51D-16*EXP(-1521/TEMP)*0.5
	p(381)= 4.60D-12*EXP(452/TEMP)*0.3
	p(382)= 4.60D-12*EXP(452/TEMP)*0.7
	p(383)= J(23)
	p(384)= J(24)
	p(385)= KNO3AL*2.0
	p(386)= 2.95D-11
	p(387)= J(34)
	p(388)= 5.23D-12
	p(389)= 1.6D-12*EXP(305/TEMP)
	p(390)= J(22)
	p(391)= KDEC*0.11
	p(392)= KDEC*0.89
	p(393)= KDEC*0.11
	p(394)= KDEC*0.89
	p(395)= KNO3AL
	p(396)= 3.40D-12
	p(397)= J(57)
	p(398)= J(56)
	p(399)= KRO2HO2*0.706
	p(400)= KRO2NO*0.072
	p(401)= KRO2NO*0.928
	p(402)= KRO2NO3
	p(403)= 8.00D-13*0.8*RO2
	p(404)= 8.00D-13*0.2*RO2
	p(405)= KRO2HO2*0.706
	p(406)= KRO2NO*0.134
	p(407)= KRO2NO*0.866
	p(408)= KRO2NO3
	p(409)= 2.90D-12*RO2*0.1
	p(410)= 2.90D-12*RO2*0.8
	p(411)= 2.90D-12*RO2*0.1
	p(412)= KAPHO2*0.44
	p(413)= KAPHO2*0.15
	p(414)= KAPHO2*0.41
	p(415)= KAPNO
	p(416)= KFPAN
	p(417)= KRO2NO3*1.74
	p(418)= 1.00D-11*0.7*RO2
	p(419)= 1.00D-11*0.3*RO2
	p(420)= KRO2HO2*0.706
	p(421)= KRO2NO*0.019
	p(422)= KRO2NO*0.981
	p(423)= KRO2NO3
	p(424)= 9.20D-14*0.7*RO2
	p(425)= 9.20D-14*0.3*RO2
	p(426)= 9.05D-12*0.5
	p(427)= 9.05D-12*0.5
	p(428)= KDEC*0.36
	p(429)= KDEC*0.20
	p(430)= KDEC*0.24
	p(431)= KDEC*0.20
	p(432)= 4.34D-12*0.16
	p(433)= 4.34D-12*0.84
	p(434)= J(17)
	p(435)= KRO2HO2*0.706
	p(436)= KRO2NO*0.145
	p(437)= KRO2NO*0.855
	p(438)= KRO2NO3
	p(439)= 2.90D-12*0.1*RO2
	p(440)= 2.90D-12*0.8*RO2
	p(441)= 2.90D-12*0.1*RO2
	p(442)= KRO2HO2*0.706
	p(443)= KRO2NO*0.064
	p(444)= KRO2NO*0.936
	p(445)= 8.80D-13*0.2*RO2
	p(446)= 8.80D-13*0.2*RO2
	p(447)= 8.80D-13*0.6*RO2
	p(448)= 1.50D-11*0.719
	p(449)= 1.50D-11*0.281
	p(450)= KAPHO2*0.44
	p(451)= KAPHO2*0.15
	p(452)= KAPHO2*0.41
	p(453)= KAPNO
	p(454)= KFPAN
	p(455)= KRO2NO3*1.74
	p(456)= 1.00D-11*0.7*RO2
	p(457)= 1.00D-11*0.3*RO2
	p(458)= KDEC*0.11
	p(459)= KDEC*0.89
	p(460)= KRO2HO2*0.706
	p(461)= KRO2NO
	p(462)= KRO2NO3
	p(463)= 9.20D-14*0.7*RO2
	p(464)= 9.20D-14*0.3*RO2
	p(465)= KDEC*0.18
	p(466)= KDEC*0.82
	p(467)= 1.33D-12*0.33
	p(468)= 1.33D-12*0.67
	p(469)= J(56)
	p(470)= J(57)
	p(471)= KDEC*0.18
	p(472)= KDEC*0.82
	p(473)= KRO2HO2*0.706
	p(474)= KRO2NO*0.072
	p(475)= KRO2NO*0.928
	p(476)= KRO2NO3
	p(477)= 8.00D-13*0.8*RO2
	p(478)= 8.00D-13*0.2*RO2
	p(479)= KAPHO2*0.44
	p(480)= KAPHO2*0.15
	p(481)= KAPHO2*0.41
	p(482)= KAPNO
	p(483)= KFPAN
	p(484)= KRO2NO3*1.74
	p(485)= 1.00D-11*0.7*RO2
	p(486)= 1.00D-11*0.3*RO2
	p(487)= KRO2HO2*0.706
	p(488)= KRO2NO
	p(489)= KRO2NO3
	p(490)= 9.20D-14*RO2
	p(491)= 1.2D-15
	p(492)= 1.0D-14
	p(493)= 1.0D-15
	p(494)= 7.0D-14
	p(495)= 6.0D-18*H2O
	p(496)= 1.0D-17*H2O
	p(497)= KDEC
	p(498)= 1.2D-15
	p(499)= 1.0D-14
	p(500)= 1.0D-15
	p(501)= 7.0D-14
	p(502)= 6.0D-18*H2O
	p(503)= J(41)
	p(504)= 3.31D-11*0.89
	p(505)= 3.31D-11*0.11
	p(506)= 1.98D-12*0.445
	p(507)= 1.98D-12*0.414
	p(508)= 1.98D-12*0.141
	p(509)= KDEC
	p(510)= 3.30D-12
	p(511)= 1.53D-11
	p(512)= 2.16D-11
	p(513)= J(41)
	p(514)= 2.52D-11
	p(515)= KBPAN
	p(516)= 2.16D-11
	p(517)= 1.2D-15
	p(518)= 1.0D-14
	p(519)= 1.0D-15
	p(520)= 7.0D-14
	p(521)= 6.0D-18*H2O
	p(522)= 1.0D-17*H2O
	p(523)= 1.2D-15
	p(524)= 1.0D-14
	p(525)= 1.0D-15
	p(526)= 7.0D-14
	p(527)= 6.0D-18*H2O
	p(528)= KAPHO2*0.15
	p(529)= KAPHO2*0.41
	p(530)= KAPHO2*0.44
	p(531)= KAPNO
	p(532)= KFPAN
	p(533)= KRO2NO3*1.74
	p(534)= 1.00D-11*0.7*RO2
	p(535)= 1.00D-11*0.3*RO2
	p(536)= J(41)
	p(537)= 2.81D-11
	p(538)= KDEC
	p(539)= 2.69D-11
	p(540)= 1.4D-11
	p(541)= J(17)
	p(542)= KAPHO2*0.56
	p(543)= KAPHO2*0.44
	p(544)= KAPNO
	p(545)= KFPAN
	p(546)= KRO2NO3*1.74
	p(547)= 1.00D-11*RO2
	p(548)= 1.69D-12
	p(549)= 1.90D-12*EXP(190/TEMP)
	p(550)= J(41)
	p(551)= KROPRIM*O2
	p(552)= 7.00D+03
	p(553)= KNO3AL*2.4
	p(554)= 3.55D-12
	p(555)= J(57)
	p(556)= J(56)
	p(557)= 9.16D-13
	p(558)= 3.47D-12
	p(559)= 1.90D-12*EXP(190/TEMP)
	p(560)= J(41)
	p(561)= 7.00D+03
	p(562)= KROSEC*O2
	p(563)= 1.71D-12
	p(564)= KRO2HO2*0.387
	p(565)= KRO2NO
	p(566)= KRO2NO3
	p(567)= 2.00D-12*0.2*RO2
	p(568)= 2.00D-12*0.6*RO2
	p(569)= 2.00D-12*0.2*RO2
	p(570)= 1.20D-15
	p(571)= 1.00D-14
	p(572)= 1.00D-15
	p(573)= 7.00D-14
	p(574)= 6.00D-18*H2O
	p(575)= 1.00D-17*H2O
	p(576)= 6.6D-12*EXP(-1240/TEMP)
	p(577)= 1.85D-12*EXP(-1690/TEMP)
	p(578)= 2.44D-11
	p(579)= 1.90D-12*EXP(190/TEMP)
	p(580)= J(41)
	p(581)= 2.00D+14*EXP(-6410/TEMP)
	p(582)= 1.20D-11*0.613
	p(583)= 1.20D-11*0.387
	p(584)= 1.83D-11
	p(585)= 1.90D-12*EXP(190/TEMP)
	p(586)= J(41)
	p(587)= 2.00D+14*EXP(-5505/TEMP)
	p(588)= KNO3AL*2.4
	p(589)= 1.7D-11
	p(590)= J(17)
	p(591)= J(34)
	p(592)= 8.0D-13
	p(593)= KAPHO2*0.56
	p(594)= KAPHO2*0.44
	p(595)= KAPNO
	p(596)= KFPAN
	p(597)= KRO2NO3*1.74
	p(598)= 1.00D-11*RO2
	p(599)= J(34)
	p(600)= J(35)
	p(601)= KNO3AL*4.0
	p(602)= 1.23D-11
	p(603)= KAPHO2*0.15
	p(604)= KAPHO2*0.41
	p(605)= KAPHO2*0.44
	p(606)= KAPNO
	p(607)= KFPAN
	p(608)= KRO2NO3*1.74
	p(609)= 1.00D-11*0.3*RO2
	p(610)= 1.00D-11*0.7*RO2
	p(611)= KDEC*0.24
	p(612)= KDEC*0.20
	p(613)= KDEC*0.20
	p(614)= KDEC*0.36
	p(615)= J(34)
	p(616)= KNO3AL*2.4
	p(617)= 1.44D-11
	p(618)= KRO2HO2*0.625
	p(619)= KRO2NO*0.017
	p(620)= KRO2NO*0.983
	p(621)= KRO2NO3
	p(622)= 2.00D-12*RO2*0.2
	p(623)= 2.00D-12*RO2*0.6
	p(624)= 2.00D-12*RO2*0.2
	p(625)= KRO2HO2*0.625
	p(626)= KRO2NO
	p(627)= KRO2NO3
	p(628)= 8.80D-13*RO2*0.2
	p(629)= 8.80D-13*RO2*0.2
	p(630)= 8.80D-13*RO2*0.6
	p(631)= 2.59D-11
	p(632)= J(17)
	p(633)= KNO3AL*7.5
	p(634)= 2.20D-11
	p(635)= 1.2D-15
	p(636)= 1.0D-14
	p(637)= 1.0D-15
	p(638)= 7.0D-14
	p(639)= 6.0D-18*H2O
	p(640)= 1.0D-17*H2O
	p(641)= 1.2D-15
	p(642)= 1.0D-14
	p(643)= 1.0D-15
	p(644)= 7.0D-14
	p(645)= 6.0D-18*H2O
	p(646)= KAPHO2*0.44
	p(647)= KAPHO2*0.15
	p(648)= KAPHO2*0.41
	p(649)= KAPNO
	p(650)= KFPAN
	p(651)= KRO2NO3*1.74
	p(652)= 1.00D-11*0.7*RO2
	p(653)= 1.00D-11*0.3*RO2
	p(654)= KDEC
	p(655)= J(41)
	p(656)= 1.01D-11*0.65
	p(657)= 1.01D-11*0.35
	p(658)= 2.00D-12*0.07
	p(659)= 2.00D-12*0.39
	p(660)= 2.00D-12*0.07
	p(661)= 2.00D-12*0.33
	p(662)= 2.00D-12*0.14
	p(663)= KDEC
	p(664)= 6.68D-12
	p(665)= 1.18D-10*0.22
	p(666)= 1.18D-10*0.03
	p(667)= 1.18D-10*0.75
	p(668)= J(41)
	p(669)= J(54)
	p(670)= 2.94D-11
	p(671)= KDEC
	p(672)= 4.21D-11
	p(673)= J(24)
	p(674)= 7.78D-11
	p(675)= 2.52D-11
	p(676)= J(41)
	p(677)= 2.88D-11
	p(678)= KBPAN
	p(679)= 2.52D-11
	p(680)= J(41)
	p(681)= 3.16D-11
	p(682)= J(15)
	p(683)= 2.32D-11
	p(684)= KDEC
	p(685)= 3.04D-11
	p(686)= KRO2HO2*0.706
	p(687)= KRO2NO
	p(688)= KRO2NO3
	p(689)= 9.20D-14*RO2
	p(690)= KRO2HO2*0.706
	p(691)= KRO2NO
	p(692)= KRO2NO3
	p(693)= 8.80D-13*RO2
	p(694)= 1.2D-15
	p(695)= 1.0D-14
	p(696)= 1.0D-15
	p(697)= 7.0D-14
	p(698)= 1.0D-17*H2O
	p(699)= 6.0D-18*H2O
	p(700)= 6.78D-12
	p(701)= J(56)
	p(702)= J(57)*2
	p(703)= KAPHO2*0.44
	p(704)= KAPHO2*0.15
	p(705)= KAPHO2*0.41
	p(706)= KAPNO
	p(707)= KFPAN
	p(708)= KRO2NO3*1.74
	p(709)= 1.00D-11*0.7*RO2
	p(710)= 1.00D-11*0.3*RO2
	p(711)= 1.27D-11*0.35
	p(712)= 1.27D-11*0.34
	p(713)= 1.27D-11*0.31
	p(714)= J(41)
	p(715)= 1.63D-12*0.50
	p(716)= 1.63D-12*0.50
	p(717)= KDEC
	p(718)= 3.27D-12
	p(719)= J(56)
	p(720)= J(57)
	p(721)= 6.65D-12*0.71
	p(722)= 6.65D-12*0.29
	p(723)= 1.59D-11*0.73
	p(724)= 1.59D-11*0.27
	p(725)= J(41)
	p(726)= KDEC
	p(727)= KRO2HO2*0.706
	p(728)= KRO2NO
	p(729)= KRO2NO3
	p(730)= 9.20D-14*RO2
	p(731)= J(15)
	p(732)= KNO3AL*7.5
	p(733)= 2.76D-11
	p(734)= 2.52D-11
	p(735)= J(41)
	p(736)= 2.88D-11
	p(737)= KBPAN
	p(738)= 2.52D-11
	p(739)= J(41)
	p(740)= 3.16D-11
	p(741)= KDEC
	p(742)= 3.04D-11
	p(743)= 7.20D-12
	p(744)= J(15)
	p(745)= 1.2D-15
	p(746)= 1.0D-14
	p(747)= 1.0D-15
	p(748)= 7.0D-14
	p(749)= 6.0D-18*H2O
	p(750)= 9.20D-12*0.61
	p(751)= 9.20D-12*0.39
	p(752)= J(41)
	p(753)= 1.80D+13*(TEMP/298)**1.7*EXP(-4733/TEMP)
	p(754)= 1.80D+13*(TEMP/298)**1.7*EXP(-4079/TEMP)
	p(755)= 5.60D-12
	p(756)= 2.73D-12
	p(757)= 6.19D-12
	p(758)= J(41)
	p(759)= 1.12D-12
	p(760)= KBPAN
	p(761)= J(41)
	p(762)= J(22)
	p(763)= 9.70D-12
	p(764)= KDEC
	p(765)= 4.52D-12*0.19
	p(766)= 4.52D-12*0.81
	p(767)= J(15)
	p(768)= J(34)
	p(769)= 1.34D-11
	p(770)= J(34)
	p(771)= 1.23D-11
	p(772)= J(41)
	p(773)= J(15)
	p(774)= 1.58D-11
	p(775)= KBPAN
	p(776)= 1.22D-11
	p(777)= KAPHO2*0.44
	p(778)= KAPHO2*0.15
	p(779)= KAPHO2*0.41
	p(780)= KAPNO
	p(781)= KFPAN
	p(782)= KRO2NO3*1.74
	p(783)= 1.00D-11*RO2*0.7
	p(784)= 1.00D-11*RO2*0.3
	p(785)= 8.8D-12*EXP(-1320/TEMP) + 1.7D-14*EXP(423/TEMP)
	p(786)= J(21)
	p(787)= 6.99D-11
	p(788)= J(41)
	p(789)= J(17)
	p(790)= 6.64D-11
	p(791)= KBPAN
	p(792)= KAPHO2*0.44
	p(793)= KAPHO2*0.15
	p(794)= KAPHO2*0.41
	p(795)= KAPNO
	p(796)= KFPAN
	p(797)= KRO2NO3*1.74
	p(798)= 1.00D-11*0.7*RO2
	p(799)= 1.00D-11*0.3*RO2
	p(800)= KDEC
	p(801)= 2.91D-11
	p(802)= 1.90D-12*EXP(190/TEMP)
	p(803)= J(41)
	p(804)= J(15)
	p(805)= KAPHO2*0.44
	p(806)= KAPHO2*0.56
	p(807)= KAPNO
	p(808)= KFPAN
	p(809)= KRO2NO3*1.74
	p(810)= 1.00D-11*RO2
	p(811)= J(22)
	p(812)= J(41)
	p(813)= 7.34D-12
	p(814)= KBPAN
	p(815)= 3.74D-12
	p(816)= 8.66D-12
	p(817)= J(41)
	p(818)= 1.22D-11
	p(819)= KBPAN
	p(820)= 8.63D-12
	p(821)= 1.20D-15
	p(822)= 1.00D-14
	p(823)= 1.00D-15
	p(824)= 7.00D-14
	p(825)= 6.00D-18*H2O
	p(826)= 1.00D-17*H2O
	p(827)= 5.98D-11
	p(828)= J(41)
	p(829)= J(22)
	p(830)= 4.37D-12
	p(831)= KDEC
	p(832)= J(15)
	p(833)= KNO3AL*4.0
	p(834)= 2.66D-11
	p(835)= 2.10D-11
	p(836)= J(22)
	p(837)= 4.39D-12
	p(838)= J(22)
	p(839)= J(41)
	p(840)= KDEC
	p(841)= 4.44D-12
	p(842)= J(35)
	p(843)= 1.72D-13*EXP(-1190/TEMP)
	p(844)= 0.68*2.00E-11
	p(845)= 0.255*2.00E-11
	p(846)= 0.065*2.00E-11
	p(847)= 0.3*J(18)
	p(848)= 0.4*J(18)
	p(849)= 0.3*J(18)
	p(850)= 0.5*2.9E-19
	p(851)= 0.5*2.9E-19
	p(852)= KAPHO2*0.44
	p(853)= KAPHO2*0.56
	p(854)= KAPNO
	p(855)= KFPAN
	p(856)= KRO2NO3*1.74
	p(857)= 1.00D-11*RO2
	p(858)= 1.68D-13
	p(859)= 3.63D-12
	p(860)= J(41)
	p(861)= 1.12D-14
	p(862)= KBPAN
	p(863)= J(17)
	p(864)= J(41)
	p(865)= 2.67D-11
	p(866)= J(17)
	p(867)= 4.22D-12
	p(868)= J(56)
	p(869)= J(57)
	p(870)= 1.21D-12*0.56
	p(871)= 1.21D-12*0.44
	p(872)= J(17)
	p(873)= 2.29D-11
	p(874)= 1.31D-11*0.667
	p(875)= 1.31D-11*0.333
	p(876)= KRO2HO2*0.706
	p(877)= KRO2NO*0.072
	p(878)= KRO2NO*0.928
	p(879)= KRO2NO3
	p(880)= 8.00D-13*RO2*0.8
	p(881)= 8.00D-13*RO2*0.2
	p(882)= 3.40D-15
	p(883)= 7.80D-19*0.5
	p(884)= 7.80D-19*0.5
	p(885)= 4.83D-11*0.376
	p(886)= 4.83D-11*0.624
	p(887)= KRO2HO2*0.706
	p(888)= KRO2NO
	p(889)= KRO2NO3
	p(890)= 9.20D-14*RO2
	p(891)= KAPHO2*0.15
	p(892)= KAPHO2*0.41
	p(893)= KAPHO2*0.44
	p(894)= KAPNO
	p(895)= KFPAN
	p(896)= 1.00D-11*0.3*RO2
	p(897)= 1.00D-11*0.7*RO2
	p(898)= J(17)
	p(899)= J(41)
	p(900)= 3.90D-11
	p(901)= KDEC
	p(902)= J(41)
	p(903)= J(15)
	p(904)= 5.34D-11
	p(905)= KDEC
	p(906)= 1.34D-12*0.44
	p(907)= 1.34D-12*0.15
	p(908)= KAPHO2*0.15
	p(909)= KAPHO2*0.41
	p(910)= KAPHO2*0.44
	p(911)= KAPNO
	p(912)= KFPAN
	p(913)= KRO2NO3*1.74
	p(914)= 1.00D-11*0.3*RO2
	p(915)= 1.00D-11*0.7*RO2
	p(916)= 4.42D-12*0.15
	p(917)= 4.42D-12*0.85
	p(918)= J(41)
	p(919)= 8.21D-13
	p(920)= KBPAN
	p(921)= J(15)
	p(922)= 2.41D-11
	p(923)= J(15)
	p(924)= 1.85D-11
	p(925)= J(17)
	p(926)= 1.85D-11
	p(927)= 1.35D-11
	p(928)= J(34)
	p(929)= 9.37D-12*0.54
	p(930)= 9.37D-12*0.46
	p(931)= J(17)
	p(932)= J(22)
	p(933)= J(41)
	p(934)= 3.60D-12
	p(935)= 1.57D-11
	p(936)= KDEC
	p(937)= KAPHO2*0.56
	p(938)= KAPHO2*0.44
	p(939)= KAPNO
	p(940)= KFPAN
	p(941)= KRO2NO3*1.74
	p(942)= 1.00D-11*RO2
	p(943)= KAPHO2*0.56
	p(944)= KAPHO2*0.44
	p(945)= KAPNO
	p(946)= KFPAN
	p(947)= KRO2NO3*1.74
	p(948)= 1.00D-11*RO2
	p(949)= 2.58D-11
	p(950)= J(17)
	p(951)= J(41)
	p(952)= 2.27D-11
	p(953)= J(17)
	p(954)= KAPHO2*0.15
	p(955)= KAPHO2*0.41
	p(956)= KAPHO2*0.44
	p(957)= KAPNO
	p(958)= KFPAN
	p(959)= KRO2NO3*1.74
	p(960)= 1.00D-11*0.3*RO2
	p(961)= 1.00D-11*0.7*RO2
	p(962)= J(17)
	p(963)= 2.15D-11
	p(964)= KAPHO2*0.15
	p(965)= KAPHO2*0.41
	p(966)= KAPHO2*0.44
	p(967)= KAPNO
	p(968)= KFPAN
	p(969)= KRO2NO3*1.74
	p(970)= 1.00D-11*0.3*RO2
	p(971)= 1.00D-11*0.7*RO2
	p(972)= 1.72D-12
	p(973)= J(41)
	p(974)= 4.80D-12
	p(975)= KBPAN
	p(976)= 4.75D-13
	p(977)= 1.36D-13*EXP(1250/TEMP)*0.15
	p(978)= 1.36D-13*EXP(1250/TEMP)*0.85
	p(979)= KRO2NO
	p(980)= KRO2NO3
	p(981)= 2*(K298CH3O2*8.0D-12)**0.5*RO2*0.2
	p(982)= 2*(K298CH3O2*8.0D-12)**0.5*RO2*0.6
	p(983)= 2*(K298CH3O2*8.0D-12)**0.5*RO2*0.2
	p(984)= 3.14D-13
	p(985)= 3.77D-12
	p(986)= J(41)
	p(987)= 1.43D-13
	p(988)= KBPAN
	p(989)= 9.34D-12
	p(990)= J(41)
	p(991)= 2.34D-12
	p(992)= KBPAN
	p(993)= 2.89D-12
	p(994)= J(34)
	p(995)= KAPHO2*0.56
	p(996)= KAPHO2*0.44
	p(997)= KAPNO
	p(998)= KFPAN
	p(999)= KRO2NO3*1.74
	p(1000)= 1.00D-11*RO2
	p(1001)= 1.44D-11
	p(1002)= J(34)
	p(1003)= J(35)
	p(1004)= KRO2HO2*0.52
	p(1005)= KRO2NO
	p(1006)= KRO2NO3
	p(1007)= 8.8D-13*RO2*0.6
	p(1008)= 8.8D-13*RO2*0.2
	p(1009)= 8.8D-13*RO2*0.2
	p(1010)= 0.4*0.52*KRO2HO2
	p(1011)= 0.6*0.52*KRO2HO2
	p(1012)= 0.05*KRO2NO
	p(1013)= 0.95*KRO2NO
	p(1014)= KRO2NO3
	p(1015)= 2.0D-12*RO2*0.2
	p(1016)= 2.0D-12*RO2*0.6
	p(1017)= 2.0D-12*RO2*0.2
	p(1018)= 3.3D-12*EXP(-2880/TEMP)
	p(1019)= 9.1D-15*EXP(-2580/TEMP)
	p(1020)= KMT15
	p(1021)= KDEC*0.24
	p(1022)= KDEC*0.20
	p(1023)= KDEC*0.20
	p(1024)= KDEC*0.36
	p(1025)= J(41)
	p(1026)= 4.81D-12
	p(1027)= KBPAN
	p(1028)= 1.21D-12
	p(1029)= KAPHO2*0.44
	p(1030)= KAPHO2*0.15
	p(1031)= KAPHO2*0.41
	p(1032)= KAPNO
	p(1033)= KFPAN
	p(1034)= KRO2NO3*1.74
	p(1035)= 1.00D-11*0.7*RO2
	p(1036)= 1.00D-11*0.3*RO2
	p(1037)= KAPHO2*0.44
	p(1038)= KAPHO2*0.15
	p(1039)= KAPHO2*0.41
	p(1040)= KAPNO
	p(1041)= KFPAN
	p(1042)= KRO2NO3*1.74
	p(1043)= 1.00D-11*0.7*RO2
	p(1044)= 1.00D-11*0.3*RO2
	p(1045)= 3.79D-12
	p(1046)= KAPHO2*0.44
	p(1047)= KAPHO2*0.15
	p(1048)= KAPHO2*0.41
	p(1049)= KAPNO
	p(1050)= KFPAN
	p(1051)= KRO2NO3*1.74
	p(1052)= 1.00D-11*0.7*RO2
	p(1053)= 1.00D-11*0.3*RO2
	p(1054)= KRO2HO2*0.706
	p(1055)= KRO2NO
	p(1056)= KRO2NO3
	p(1057)= 9.20D-14*RO2
	p(1058)= KRO2HO2*0.706
	p(1059)= KRO2NO
	p(1060)= KRO2NO3
	p(1061)= 8.80D-13*RO2
	p(1062)= 1.32D-11*0.728
	p(1063)= 1.32D-11*0.272
	p(1064)= J(41)
	p(1065)= 2.43D-12
	p(1066)= 1.80D+13*(TEMP/298)**1.7*EXP(-4079/TEMP)
	p(1067)= 1.80D+13*(TEMP/298)**1.7*EXP(-4733/TEMP)
	p(1068)= 9.60D-12
	p(1069)= KAPHO2*0.15
	p(1070)= KAPHO2*0.41
	p(1071)= KAPHO2*0.44
	p(1072)= KAPNO
	p(1073)= KFPAN
	p(1074)= KRO2NO3*1.74
	p(1075)= 1.00D-11*RO2*0.3
	p(1076)= 1.00D-11*RO2*0.7
	p(1077)= KDEC*0.37
	p(1078)= KDEC*0.50
	p(1079)= KDEC*0.13
	p(1080)= KDEC*0.18
	p(1081)= KDEC*0.82
	p(1082)= KRO2HO2*0.625
	p(1083)= KRO2NO
	p(1084)= KRO2NO3
	p(1085)= 9.20D-14*RO2*0.7
	p(1086)= 9.20D-14*RO2*0.3
	p(1087)= 1.37D-11
	p(1088)= J(41)
	p(1089)= KDEC
	p(1090)= 2.49D-12
	p(1091)= J(41)
	p(1092)= 5.57D-12*0.68
	p(1093)= 5.57D-12*0.32
	p(1094)= KBPAN
	p(1095)= 1.97D-12
	p(1096)= 3.70D-12
	p(1097)= J(56)
	p(1098)= J(57)
	p(1099)= 6.78D-12
	p(1100)= J(56)
	p(1101)= J(57)
	p(1102)= J(41)
	p(1103)= 3.18D-12
	p(1104)= J(56)
	p(1105)= J(57)
	p(1106)= KBPAN
	p(1107)= KAPHO2*0.15
	p(1108)= KAPHO2*0.41
	p(1109)= KAPHO2*0.44
	p(1110)= KAPNO
	p(1111)= KFPAN
	p(1112)= KRO2NO3*1.74
	p(1113)= 1.00D-11*0.3*RO2
	p(1114)= 1.00D-11*0.7*RO2
	p(1115)= KAPHO2*0.15
	p(1116)= KAPHO2*0.41
	p(1117)= KAPHO2*0.44
	p(1118)= KAPNO
	p(1119)= KFPAN
	p(1120)= KRO2NO3*1.74
	p(1121)= 1.00D-11*0.3*RO2
	p(1122)= 1.00D-11*0.7*RO2
	p(1123)= KAPHO2*0.15
	p(1124)= KAPHO2*0.41
	p(1125)= KAPHO2*0.44
	p(1126)= KAPNO
	p(1127)= KFPAN
	p(1128)= KRO2NO3*1.74
	p(1129)= 1.00D-11*0.3*RO2
	p(1130)= 1.00D-11*0.7*RO2
	p(1131)= 2.53D-11*0.86
	p(1132)= 2.53D-11*0.14
	p(1133)= J(17)
	p(1134)= J(15)
	p(1135)= KAPHO2*0.15
	p(1136)= KAPHO2*0.41
	p(1137)= KAPHO2*0.44
	p(1138)= KAPNO
	p(1139)= KFPAN
	p(1140)= KRO2NO3*1.74
	p(1141)= KRO2NO3*1.74
	p(1142)= 1.00D-11*0.3*RO2
	p(1143)= 1.00D-11*0.7*RO2
	p(1144)= J(17)
	p(1145)= J(22)
	p(1146)= 2.73D-11
	p(1147)= J(41)
	p(1148)= 1.04D-11
	p(1149)= KBPAN
	p(1150)= 6.80D-12
	p(1151)= 4.11D-12
	p(1152)= J(41)
	p(1153)= 5.10D-13
	p(1154)= KBPAN
	p(1155)= J(41)
	p(1156)= KAPHO2*0.56
	p(1157)= KAPHO2*0.44
	p(1158)= KAPNO
	p(1159)= KFPAN
	p(1160)= KRO2NO3*1.74
	p(1161)= 1.00D-11*RO2
	p(1162)= KAPHO2*0.56
	p(1163)= KAPHO2*0.44
	p(1164)= KAPNO
	p(1165)= KFPAN
	p(1166)= KRO2NO3*1.74
	p(1167)= 1.00D-11*RO2
	p(1168)= 1.66D-12
	p(1169)= J(41)
	p(1170)= 4.74D-12
	p(1171)= KBPAN
	p(1172)= 1.14D-12
	p(1173)= 1.23D-12*0.32
	p(1174)= 1.23D-12*0.68
	p(1175)= 4.31D-12*0.09
	p(1176)= 4.31D-12*0.91
	p(1177)= J(41)
	p(1178)= 7.10D-13
	p(1179)= KBPAN
	p(1180)= 1.90D-12*EXP(190/TEMP)
	p(1181)= 8.39D-12
	p(1182)= J(22)
	p(1183)= J(41)
	p(1184)= J(41)
	p(1185)= 9.43D-12
	p(1186)= KBPAN
	p(1187)= 5.83D-12
	p(1188)= 4.77D-11
	p(1189)= J(41)
	p(1190)= KDEC
	p(1191)= 6.22D-11
	p(1192)= 1.36D-10
	p(1193)= J(15)*2
	p(1194)= 9.258E-11
	p(1195)= J(41)
	p(1196)= J(15)
	p(1197)= 2.92D-11
	p(1198)= KDEC
	p(1199)= KRO2HO2*0.387
	p(1200)= KRO2NO
	p(1201)= KRO2NO3
	p(1202)= 6.00D-13*0.6*RO2
	p(1203)= 6.00D-13*0.2*RO2
	p(1204)= 6.00D-13*0.2*RO2
	p(1205)= 1.53D-13*EXP(1300/TEMP)
	p(1206)= KRO2NO*0.005
	p(1207)= KRO2NO*0.995
	p(1208)= KRO2NO3
	p(1209)= 2*(KCH3O2*7.8D-14*EXP(1000/TEMP))**0.5*RO2*0.2
	p(1210)= 2*(KCH3O2*7.8D-14*EXP(1000/TEMP))**0.5*RO2*0.6
	p(1211)= 2*(KCH3O2*7.8D-14*EXP(1000/TEMP))**0.5*RO2*0.2
	p(1212)= J(41)
	p(1213)= 6.50D-12
	p(1214)= J(41)*2
	p(1215)= 9.58D-12
	p(1216)= J(41)
	p(1217)= KBPAN
	p(1218)= 5.98D-12
	p(1219)= 1.36D-12
	p(1220)= J(41)
	p(1221)= 4.08D-12
	p(1222)= KBPAN
	p(1223)= 4.85D-13
	p(1224)= KAPHO2*0.15
	p(1225)= KAPHO2*0.41
	p(1226)= KAPHO2*0.44
	p(1227)= KAPNO
	p(1228)= KFPAN
	p(1229)= KRO2NO3*1.74
	p(1230)= 1.00D-11*0.3*RO2
	p(1231)= 1.00D-11*0.7*RO2
	p(1232)= 3.04D-12
	p(1233)= J(41)
	p(1234)= 6.12D-12
	p(1235)= KBPAN
	p(1236)= 2.52D-12
	p(1237)= 4.30D-11
	p(1238)= J(41)
	p(1239)= J(17)
	p(1240)= KDEC
	p(1241)= 5.74D-11
	p(1242)= J(41)
	p(1243)= J(15)
	p(1244)= KDEC
	p(1245)= 2.98D-11
	p(1246)= J(17)
	p(1247)= 1.93D-11
	p(1248)= J(17)
	p(1249)= 3.85D-12
	p(1250)= J(56)
	p(1251)= J(57)
	p(1252)= 5.25D-12
	p(1253)= J(22)
	p(1254)= 2.67D-11
	p(1255)= J(17)
	p(1256)= 1.84D-11
	p(1257)= 1.99D-11
	p(1258)= J(41)
	p(1259)= 2.90D-11
	p(1260)= KBPAN
	p(1261)= 1.20D-15
	p(1262)= 1.00D-14
	p(1263)= 1.00D-15
	p(1264)= 7.00D-14
	p(1265)= 6.00D-18*H2O
	p(1266)= 4.17D-11
	p(1267)= J(17)
	p(1268)= J(41)
	p(1269)= KDEC
	p(1270)= 3.82D-11
	p(1271)= J(17)
	p(1272)= 4.93D-12
	p(1273)= J(56)
	p(1274)= J(57)
	p(1275)= 8.01D-12
	p(1276)= J(41)
	p(1277)= J(56)
	p(1278)= J(57)
	p(1279)= 4.41D-12
	p(1280)= J(56)
	p(1281)= J(57)
	p(1282)= KBPAN
	p(1283)= 3.69D-12
	p(1284)= J(41)
	p(1285)= J(35)
	p(1286)= 1.02D-13
	p(1287)= KBPAN
	p(1288)= 7.40D-12
	p(1289)= 1.09D-11
	p(1290)= J(41)
	p(1291)= 7.26D-12
	p(1292)= KBPAN
	p(1293)= 1.72D-12
	p(1294)= 5.18D-12
	p(1295)= J(41)
	p(1296)= 1.58D-12
	p(1297)= KBPAN
	p(1298)= 1.73D-12
	p(1299)= 5.19D-12
	p(1300)= J(41)
	p(1301)= 1.59D-12
	p(1302)= KBPAN
	p(1303)= KAPHO2*0.44
	p(1304)= KAPHO2*0.15
	p(1305)= KAPHO2*0.41
	p(1306)= KAPNO
	p(1307)= KFPAN
	p(1308)= KRO2NO3*1.74
	p(1309)= 1.00D-11*0.7*RO2
	p(1310)= 1.00D-11*0.3*RO2
	p(1311)= KAPHO2*0.44
	p(1312)= KAPHO2*0.15
	p(1313)= KAPHO2*0.41
	p(1314)= KAPNO
	p(1315)= KFPAN
	p(1316)= KRO2NO3*1.74
	p(1317)= 1.00D-11*0.7*RO2
	p(1318)= 1.00D-11*0.3*RO2
	p(1319)= 6.52D-12
	p(1320)= 6.52D-12*0.39
	p(1321)= 6.52D-12*0.61
	p(1322)= J(41)
	p(1323)= 6.00D-12
	p(1324)= KBPAN
	p(1325)= 8.64D-12
	p(1326)= J(41)*2
	p(1327)= 5.04D-12
	p(1328)= KBPAN
	p(1329)= J(41)
	p(1330)= 5.66D-12
	p(1331)= J(41)
	p(1332)= 1.96D-12
	p(1333)= KBPAN
	p(1334)= 2.16D-11
	p(1335)= J(17)
	p(1336)= 2.47D-11
	p(1337)= J(41)
	p(1338)= J(17)
	p(1339)= 2.11D-11
	p(1340)= KBPAN
	p(1341)= J(17)
	p(1342)= KAPHO2*0.44
	p(1343)= KAPHO2*0.15
	p(1344)= KAPHO2*0.41
	p(1345)= KAPNO
	p(1346)= KFPAN
	p(1347)= KRO2NO3*1.74
	p(1348)= 1.00E-11*0.7*RO2
	p(1349)= 1.00E-11*0.3*RO2
	p(1350)= KAPHO2*0.44
	p(1351)= KAPHO2*0.56
	p(1352)= KAPNO
	p(1353)= KFPAN
	p(1354)= KRO2NO3*1.74
	p(1355)= 1.00D-11*RO2
	p(1356)= J(53)
	p(1357)= J(34)
	p(1358)= KNO3AL*2.4
	p(1359)= 1.24D-11
	p(1360)= 1.90D-12*EXP(190/TEMP)
	p(1361)= 1.62D-12
	p(1362)= J(41)
	p(1363)= 7.00D+03
	p(1364)= KROPRIM*O2
	p(1365)= 8.40D-13
	p(1366)= 1.90D-12*EXP(190/TEMP)
	p(1367)= 1.38D-11
	p(1368)= J(41)
	p(1369)= 9.50D+13*EXP(-5988/TEMP)
	p(1370)= KROPRIM*O2
	p(1371)= 1.45D-11
	p(1372)= J(56)
	p(1373)= J(57)
	p(1374)= 9.35D-13
	p(1375)= J(41)
	p(1376)= J(56)
	p(1377)= J(57)
	p(1378)= 4.02D-12
	p(1379)= 1.30D-13
	p(1380)= J(56)
	p(1381)= J(57)
	p(1382)= J(35)
	p(1383)= KBPAN
	p(1384)= J(56)
	p(1385)= J(57)
	p(1386)= 4.15D-13
	p(1387)= 1.36D-11
	p(1388)= J(34)
	p(1389)= KAPHO2*0.44
	p(1390)= KAPHO2*0.56
	p(1391)= KAPNO
	p(1392)= KFPAN
	p(1393)= KRO2NO3*1.74
	p(1394)= 1.00D-11*RO2
	p(1395)= 2.23D-11
	p(1396)= J(15)
	p(1397)= 2.59D-11
	p(1398)= J(41)
	p(1399)= J(15)
	p(1400)= 2.18D-11
	p(1401)= J(15)
	p(1402)= KBPAN
	p(1403)= KAPHO2*0.15
	p(1404)= KAPHO2*0.41
	p(1405)= KAPHO2*0.44
	p(1406)= KAPNO
	p(1407)= KFPAN
	p(1408)= KRO2NO3*1.74
	p(1409)= 1.00D-11*0.3*RO2
	p(1410)= 1.00D-11*0.7*RO2
	p(1411)= KDEC
	p(1412)= 1.867E-11
	p(1413)= 1.513E-11
	p(1414)= J(41)
	p(1415)= 1.865E-11
	p(1416)= KBPAN
	p(1417)= 7.33D-11
	p(1418)= J(41)
	p(1419)= 6.97D-11
	p(1420)= KBPAN
	p(1421)= KDEC
	p(1422)= KRO2HO2*0.520
	p(1423)= KRO2NO
	p(1424)= KRO2NO3
	p(1425)= 2.00D-12*RO2
	p(1426)= 2.78D-11
	p(1427)= J(41)
	p(1428)= KDEC
	 return 
 	 end
!**************** end of subroutine constants *****************

subroutine zenith(theta, secx, cosx, ttime)
double precision:: lat,pi,radian,dec,lha,ttime,theta,cosx,secx
double precision:: sinld,cosld
! solar declination angle from july 1st - harwell traj model
dec = 23.79
! latitude
lat = 50.0
pi = 4.0*atan(1.0)
! local hour angle - representing time of day
lha = (1.0+ttime/4.32d+4)*pi
radian = 180.0/pi
lat = lat/radian
dec = dec/radian
theta = acos(cos(lha)*cos(dec)*cos(lat)+sin(dec)*sin(lat))
sinld = sin(lat)*sin(dec)
cosld = cos(lat)*cos(dec)
cosx = (cos(lha)*cosld)+sinld
cosx = cos(theta)
secx = 1.0d+0/(cosx+1.0d-30)
return

end

!**************** end of subroutine zenith *****************

subroutine temperature(temp, h2o, ttime)
!c subroutine to calculate diurnal variations in temperature
double precision:: temp,ttime,sin, rh, h2o
temp = 289.86 + 8.3*sin((7.2722d-5*ttime)-1.9635)
rh=23.0*sin((7.2722d-5*ttime)+1.1781)+66.5 
h2o=6.1078*dexp(-1.0d+0*(597.3-0.57*(temp-273.16))*18.0/1.986*(1.0/temp-1.0/273.16))*10./(1.38d-16*temp)*rh 
return
end
!**************** end of subroutine temperature *************


***********************************************************************


 photolysis.txt, for use when calculating J rate coefficients in constants.f 

      J          L           M          N
      1      6.073D-05     1.743       0.474
      2      4.775D-04     0.298       0.080
      3      1.041D-05     0.723       0.279
      4      1.165D-02     0.244       0.267
      5      2.485D-02     0.168       0.108
      6      1.747D-01     0.155      0.125
      7      2.644D-03     0.261      0.288
      8      9.312D-07     1.230      0.307
     11      4.642D-05     0.762       0.353
     12      6.853D-05     0.477       0.323
     13      7.344D-06     1.202       0.417
     14      2.879D-05     1.067       0.358
     15      2.792D-05     0.805       0.338
     16      1.675D-05     0.805       0.338
     17      7.914D-05     0.764      0.364
     18      1.140D-05     0.396       0.298
     19      1.140D-05     0.396      0.298
     21      7.992D-07     1.578       0.271
     22      5.804D-06     1.092      0.377
     23      1.836D-05     0.395      0.296
     24      1.836D-05     0.395      0.296
     31      6.845D-05     0.130      0.201
     32      1.032D-05     0.130      0.201
     33      3.802D-05     0.644       0.312
     34      1.537D-04     0.170       0.208
     35      3.326D-04     0.148      0.215
     41      7.649D-06     0.682       0.279
     51      1.588D-06     1.154      0.318
     52      1.907D-06     1.244      0.335
     53      2.485D-06     1.196      0.328
     54      4.095D-06     1.111       0.316
     55      1.135D-05     0.974      0.309
     56      7.549D-06     1.015      0.324
     57      3.363D-06     1.296      0.322
     61      7.537D-04     0.499       0.266

***********************************************************************


 Fortran routine for calculation of rhs of rate equations

subroutine rate_routine(neqn,t,y,isa,nisa,isb,nisb)
!     routine to calculate right hand sides of rate equations
!
!	neqn		- [i] number of species/ rate equations
!	t		- [i] current timepoint
!	y		- [i] array of species concentrations
!	isa(m,n)	- [i] array holding reactant data m=reaction number, n=species number
!	nisa		- [i] size of isa
!	isb(m,n)	- [i] array holding product data m=reaction number, n=species number
!	nisb		- [i] size of isb
	
parameter mnsp=4000, mre=12000                   
double precision:: y(mnsp),dy(mnsp) 
double precision:: p(mre),r(mre),dy(mnsp)                                  
integer:: i,nsp,nr,isa(2,20000),isb(2,20000)
!	constants is the routine in which the values of the rate coefficients are calculated
external constants

do i = 1,nspecies
	dy(i) = 0.0d+0   
enddo             
 
call constants(p,t,y)  
                                       
do k = 1,nr                                                    
	r(k) = p(k) 
enddo	                                                      

! evaluate the RHS of each rate equation stored in dy()
do is = 1,nisa     
	r(isa(1,is)) = r(isa(1,is))*y(isa(2,is))
enddo
 
do is = 1,nisa
	dy(isa(2,is)) = dy(isa(2,is))- r(isa(1,is)
enddo

***********************************************************************


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

***********************************************************************


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

***********************************************************************

