File name: Pandora122s1_Downsview_20230106_L2Fit_fwvt1c1d20220801p1-8.txt File generation date: 20260318T010618.8Z Data description: Level 2 spectral fitting results file Data file version: fwvt1c1p1-8 Local principal investigator: Vitali Fioletov Network principal investigator: Alexander Cede Instrument type: Pandora Instrument number: 122 Spectrometer number: 1 Processing software version used: BlickP v1.8.88 Instrument operation file used: Pandora122_OF_v6d20161030.txt Instrument calibration file used: Pandora122s1_CF_v1d20220801.txt Level 1 file used: Pandora122s1_Downsview_20230106_L1_smca1c1d20220801p1-8.txt Full location name: Environment Canada HQ - Downsview Short location name: Downsview Country of location: Canada Location latitude [deg]: 43.7810 Location longitude [deg]: -79.4680 Location altitude [m]: 187 Local noon date: 20230106 Notes on s-number (L1 configuration): Corrections NOT applied although requested by the s-number are latency correction, radiometric temperature correction, matrix method stray light correction (replaced by simple method) Data caveats: None First and last pixel inside fitting window: 1652 1801 Nominal wavelengths inside fitting window [nm]: 492.1089 492.2301 492.3513 492.4725 492.5936 492.7148 492.8359 492.9570 493.0781 493.1992 493.3202 493.4412 493.5623 493.6833 493.8043 493.9252 494.0462 494.1671 494.2880 494.4090 494.5298 494.6507 494.7716 494.8924 495.0132 495.1340 495.2548 495.3756 495.4963 495.6171 495.7378 495.8585 495.9792 496.0998 496.2205 496.3411 496.4617 496.5823 496.7029 496.8235 496.9440 497.0646 497.1851 497.3056 497.4261 497.5465 497.6670 497.7874 497.9078 498.0282 498.1486 498.2690 498.3893 498.5096 498.6300 498.7502 498.8705 498.9908 499.1110 499.2313 499.3515 499.4717 499.5918 499.7120 499.8321 499.9523 500.0724 500.1925 500.3125 500.4326 500.5526 500.6726 500.7927 500.9126 501.0326 501.1526 501.2725 501.3924 501.5123 501.6322 501.7521 501.8719 501.9918 502.1116 502.2314 502.3512 502.4709 502.5907 502.7104 502.8301 502.9498 503.0695 503.1892 503.3088 503.4284 503.5481 503.6677 503.7872 503.9068 504.0263 504.1459 504.2654 504.3849 504.5043 504.6238 504.7432 504.8627 504.9821 505.1014 505.2208 505.3402 505.4595 505.5788 505.6981 505.8174 505.9367 506.0559 506.1752 506.2944 506.4136 506.5328 506.6519 506.7711 506.8902 507.0093 507.1284 507.2475 507.3665 507.4856 507.6046 507.7236 507.8426 507.9616 508.0805 508.1995 508.3184 508.4373 508.5562 508.6751 508.7939 508.9128 509.0316 509.1504 509.2692 509.3879 509.5067 509.6254 509.7441 509.8628 509.9815 --------------------------------------------------------------------------------------- Column 1: Two letter code of measurement routine Column 2: UT date and time for center-time of measurement, yyyymmddThhmmssZ (ISO 8601) Column 3: Fractional days since 1-Jan-2000 midnight for center-time of measurement Column 4: Routine count (1 for the first routine of the day, 2 for the second, etc.) Column 5: Repetition count (1 for the first set in the routine, 2 for the second, etc.) Column 6: Total duration of measurement set in seconds Column 7: Data processing type index Column 8: Solar zenith angle for center-time of measurement in degree Column 9: Solar azimuth for center-time of measurement in degree, 0=north, increases clockwise Column 10: Lunar zenith angle for center-time of measurement in degree Column 11: Lunar azimuth for center-time of measurement in degree, 0=north, increases clockwise Column 12: Pointing zenith angle in degree, absolute or relative (see next column), 999=tracker not used Column 13: Zenith pointing mode: zenith angle is... 0=absolute, 1=relative to sun, 2=relative to moon Column 14: Pointing azimuth in degree, increases clockwise, absolute (0=north) or relative (see next column), 999=tracker not used Column 15: Azimuth pointing mode: like zenith angle mode but also fixed scattering angles relative to sun (3) or moon (4) Column 16: Fitting result index: 0=no error or warning, 1, 2=warning, >2=error Column 17: Number of function evaluations used, 0=linear fitting or fitting not successful or no fitting done Column 18: rms of unweighted fitting residuals, -9=fitting not successful Column 19: Normalized rms of fitting residuals weighted with independent uncertainty, -9=fitting not successful or no uncertainty used Column 20: Expected rms based on independent uncertainty, -9=fitting not successful or no uncertainty given Column 21: Expected normalized weighted rms based on independent uncertainty, -9=fitting not successful or no uncertainty given Column 22: Water vapor slant column amount [moles per square meter], -9e99=fitting not successful Column 23: Independent uncertainty of water vapor slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -3=spectral fitting was done, but no independent uncertainty could be retrieved, -5=no independent uncertainty input was given, -9=spectral fitting not successful Column 24: Structured uncertainty of water vapor slant column amount [moles per square meter], -7=not given since method "MEAS" was chosen, -9=spectral fitting not successful Column 25: Common uncertainty of water vapor slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -6=no common uncertainty input was given, -7=not given since method "MEAS" was chosen, -9=spectral fitting not successful Column 26: rms-based uncertainty of water vapor slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -3=spectral fitting was done, but no rms-based uncertainty could be retrieved, -9=spectral fitting not successful Column 27: Water vapor effective temperature [K] Column 28: Independent uncertainty of water vapor effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -2=no temperature fitting was requested and output for effective temperature and structured uncertainty of it is based on f-code, -3=spectral fitting was done, but no independent uncertainty could be retrieved, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -5=no independent uncertainty input was given, -9=spectral fitting not successful Column 29: Structured uncertainty of water vapor effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -9=spectral fitting not successful Column 30: Common uncertainty of water vapor effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -2=no temperature fitting was requested and output for effective temperature and structured uncertainty of it is based on f-code, -3=spectral fitting was done, but no common uncertainty could be retrieved, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -6=no common uncertainty input was given, -9=spectral fitting not successful Column 31: Effective water vapor fitting wavelength [nm], -9=fitting not successful Column 32: Diffuse correction applied before fitting at effective fitting wavelength for water vapor [%], 0=no diffuse correction applied or fitting not requested, >0=measured diffuse correction, <0=(negative value of) calculated diffuse correction Column 33: Nitrogen dioxide slant column amount [moles per square meter], -9e99=fitting not successful Column 34: Independent uncertainty of nitrogen dioxide slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -3=spectral fitting was done, but no independent uncertainty could be retrieved, -5=no independent uncertainty input was given, -9=spectral fitting not successful Column 35: Structured uncertainty of nitrogen dioxide slant column amount [moles per square meter], -7=not given since method "MEAS" was chosen, -9=spectral fitting not successful Column 36: Common uncertainty of nitrogen dioxide slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -6=no common uncertainty input was given, -7=not given since method "MEAS" was chosen, -9=spectral fitting not successful Column 37: rms-based uncertainty of nitrogen dioxide slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -3=spectral fitting was done, but no rms-based uncertainty could be retrieved, -9=spectral fitting not successful Column 38: Nitrogen dioxide effective temperature [K] Column 39: Independent uncertainty of nitrogen dioxide effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -2=no temperature fitting was requested and output for effective temperature and structured uncertainty of it is based on f-code, -3=spectral fitting was done, but no independent uncertainty could be retrieved, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -5=no independent uncertainty input was given, -9=spectral fitting not successful Column 40: Structured uncertainty of nitrogen dioxide effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -9=spectral fitting not successful Column 41: Common uncertainty of nitrogen dioxide effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -2=no temperature fitting was requested and output for effective temperature and structured uncertainty of it is based on f-code, -3=spectral fitting was done, but no common uncertainty could be retrieved, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -6=no common uncertainty input was given, -9=spectral fitting not successful Column 42: Effective nitrogen dioxide fitting wavelength [nm], -9=fitting not successful Column 43: Diffuse correction applied before fitting at effective fitting wavelength for nitrogen dioxide [%], 0=no diffuse correction applied or fitting not requested, >0=measured diffuse correction, <0=(negative value of) calculated diffuse correction Column 44: Ozone slant column amount [moles per square meter], -9e99=fitting not successful Column 45: Independent uncertainty of ozone slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -3=spectral fitting was done, but no independent uncertainty could be retrieved, -5=no independent uncertainty input was given, -9=spectral fitting not successful Column 46: Structured uncertainty of ozone slant column amount [moles per square meter], -7=not given since method "MEAS" was chosen, -9=spectral fitting not successful Column 47: Common uncertainty of ozone slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -6=no common uncertainty input was given, -7=not given since method "MEAS" was chosen, -9=spectral fitting not successful Column 48: rms-based uncertainty of ozone slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -3=spectral fitting was done, but no rms-based uncertainty could be retrieved, -9=spectral fitting not successful Column 49: Ozone effective temperature [K] Column 50: Independent uncertainty of ozone effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -2=no temperature fitting was requested and output for effective temperature and structured uncertainty of it is based on f-code, -3=spectral fitting was done, but no independent uncertainty could be retrieved, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -5=no independent uncertainty input was given, -9=spectral fitting not successful Column 51: Structured uncertainty of ozone effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -9=spectral fitting not successful Column 52: Common uncertainty of ozone effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -2=no temperature fitting was requested and output for effective temperature and structured uncertainty of it is based on f-code, -3=spectral fitting was done, but no common uncertainty could be retrieved, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -6=no common uncertainty input was given, -9=spectral fitting not successful Column 53: Effective ozone fitting wavelength [nm], -9=fitting not successful Column 54: Diffuse correction applied before fitting at effective fitting wavelength for ozone [%], 0=no diffuse correction applied or fitting not requested, >0=measured diffuse correction, <0=(negative value of) calculated diffuse correction Column 55: Molecular scattering air mass factor used for molecular scattering subtraction before the fitting Column 56: Estimated uncertainty of molecular scattering air mass factor, -9=molecular scattering was not subtracted before the fitting Column 57: Lower limit used for wavelength scaling [nm] Column 58: Upper limit used for wavelength scaling [nm] Column 59: Order of smoothing polynomial Column 60: Smoothing polynomial coefficient, order 0 Column 61: Independent uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 62: Structured uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 63: Common uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 64: rms-based uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 65: Smoothing polynomial coefficient, order 1 Column 66: Independent uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 67: Structured uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 68: Common uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 69: rms-based uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 70: Mean value of measured data inside fitting window [same units as measurements] Column 71: Wavelength change polynomial coefficient, order 0 Column 72: Independent uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 73: Structured uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 74: Common uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 75: rms-based uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 76: L1 based wavelength change coefficient, order 0, -9=no L1 wavelength change determination Column 77: L1 based wavelength change coefficient, order 1, -9=no L1 wavelength change determination Column 78: L2Fit data quality flag: 0=assured high quality, 1=assured medium quality, 2=assured low quality, 10=not-assured high quality, 11=not-assured medium quality, 12=not-assured low quality Column 79: Sum over 2^i using those i, for which the corresponding L2Fit data quality parameter exceeds the DQ1 limit, 0=L1 data quality above 0, 1=Spectral fitting was not successful, 2=Wavelength shift too large, 3=Normalized rms of fitting residuals weighted with independent uncertainty too large Column 80: Sum over 2^i using those i, for which the corresponding L2Fit data quality parameter exceeds the DQ2 limit (same parameters as for DQ1) Column 81: L1 data quality flag: 0=assured high quality, 1=assured medium quality, 2=assured low quality, 10=not-assured high quality, 11=not-assured medium quality, 12=not-assured low quality Column 82: Sum over 2^i using those i, for which the corresponding L1 data quality parameter exceeds the DQ1 limit, 0=Saturated data, 1=Too few dark counts measurements, 2=No temperature given or effective temperature too different from the reference temperature, 3=Dark count too high, 4=Unsuccessful dark background fitting, 5=The dark count differs significantly from the dark map for too many pixels, 6=Absolute value of estimated average residual stray light level too high, 7=Although attempted, no wavelength change could be retrieved, 8=Absolute value of retrieved wavelength shift too large, 9=Retrieved wavelength shift differs too much from the shift predicted by the effective temperature Column 83: Sum over 2^i using those i, for which the corresponding L1 data quality parameter exceeds the DQ2 limit (same parameters as for DQ1) Column 84: Atmospheric variability [%], 999=no atmospheric variability was determined Column 85: Wavelength effective temperature [°C], 999=no effective temperature given Column 86: Estimated average residual stray light level [%] (only valid for stray light correction methods 2 and higher) Column 87: Retrieved wavelength shift from L1 data [nm], -9=no wavelength change determination Column 88: Retrieved total wavelength shift [nm], -9=no wavelength change fitting Column 89: Number of bright count cycles Column 90: Number of dark count cycles Column 91: Effective position of filterwheel #1, 0=filterwheel not used, 1-9 are valid positions Column 92: Effective position of filterwheel #2, 0=filterwheel not used, 1-9 are valid positions Column 93: Sum over 2^i, 0=spectra were interpolated in time, 1=spectra are corrected for off-target signal Column 94: Integration time [ms] Column 95: Temperature at electronics board [°C], 999=no temperature signal Column 96: Spectrometer control temperature [°C], 999=no temperature signal Column 97: Auxiliary spectrometer temperature [°C], 999=no temperature signal Columns 98-247: Unweighted slant column residuals for each pixel inside the fitting window multiplied by 1e5, 9e5=pixel was not used for fitting Columns 248-397: Normalized slant column residuals weighted with independent instrumental uncertainty for each pixel inside the fitting window multiplied by 1e5, 9e5=pixel was not used for fitting, 0=no independent instrumental uncertainty was given --------------------------------------------------------------------------------------- SQ 20230106T213822.8Z 8406.901652 516 1 50.13 2 88.16 235.96 89.15 53.10 0.00 1 0.00 1 0 15 9.641e-03 9.545e-03 6.789e-04 6.768e-04 7.884247e+03 4.310920e+01 5.164729e+02 -6.000000e+00 5.194916e+02 272.33 -2.00 20.00 -2.00 505.41 0.00e+00 5.004213e-03 1.840167e-05 2.749967e-04 -6.000000e+00 2.647995e-04 272.33 -2.00 20.00 -2.00 501.07 0.00e+00 1.142609e+00 1.313430e-02 8.348345e-02 -6.000000e+00 1.915578e-01 223.05 -2.00 3.57 -2.00 502.56 0.00e+00 16.621 0.089 489.00 513.00 1 -4.1158e+00 1.0025e-03 2.8011e-02 -6.0000e+00 1.4655e-02 -4.0834e-02 2.5160e-04 2.2337e-03 -6.0000e+00 3.7815e-03 1.0146e+01 -3.7371e-05 6.7637e-05 1.2216e-04 -6.0000e+00 1.0061e-03 -1.2590e-03 -1.6978e-06 12 8 8 10 0 0 99.96 22.23 6.73 -0.00126 -0.00152 49 9 1 4 0 1022.836 27.45 20.00 22.23 -2577 -1639 -651 -72 553 1312 1439 1348 1411 961 381 -58 -748 -975 -1476 -975 -722 358 476 464 466 563 830 1110 909 913 735 233 -361 -1642 -1455 -1128 -288 434 449 521 251 312 57 -237 194 201 743 1237 1237 1067 1130 441 -113 -730 -1468 -2414 -2138 -1310 -315 294 689 743 628 803 1003 1458 1373 1077 222 403 -455 -475 -613 -1106 -1267 -1095 -877 -509 -281 -500 -554 -1076 -1121 -1301 -1519 -911 -884 -460 -329 -520 -666 -828 -692 -208 -53 338 820 1094 627 467 -248 -801 -1439 -1693 -2114 -1455 -418 219 -176 -347 -520 -737 -346 -293 284 799 1182 1717 1779 1951 1768 1614 1014 307 593 672 875 809 940 688 553 676 436 387 84 -292 -996 -1130 -569 -507 -137 149 472 387 389 467 672 704 184 -258 -857 -1340 -1399 -1380 -2351 -1471 -599 -69 518 1248 1388 1309 1362 920 362 -54 -700 -907 -1379 -915 -689 336 455 449 452 548 809 1089 900 905 726 228 -346 -1558 -1367 -1070 -279 419 442 515 247 305 55 -232 191 200 749 1225 1237 1072 1128 437 -115 -725 -1433 -2343 -2072 -1285 -316 293 697 758 643 823 1032 1503 1382 1084 223 411 -461 -476 -609 -1095 -1254 -1089 -881 -515 -285 -502 -547 -1056 -1094 -1272 -1496 -898 -879 -460 -322 -513 -658 -816 -681 -205 -53 341 830 1100 622 452 -243 -769 -1365 -1592 -2011 -1412 -406 219 -175 -346 -530 -749 -352 -295 290 824 1231 1795 1827 2020 1851 1697 1065 320 608 683 888 829 970 716 578 706 456 405 88 -300 -1021 -1160 -588 -534 -143 158 506 419 424 501 722 758 198 -276 -928 -1432 -1488 -1480 SQ 20230106T214149.2Z 8406.904042 520 1 50.64 2 88.68 236.56 88.67 53.66 0.00 1 0.00 1 0 22 1.220e-02 1.205e-02 5.254e-04 5.234e-04 7.998846e+03 3.335388e+01 5.113435e+02 -6.000000e+00 6.573543e+02 272.33 -2.00 20.00 -2.00 505.43 0.00e+00 5.408318e-03 1.426920e-05 2.963805e-04 -6.000000e+00 3.350606e-04 272.33 -2.00 20.00 -2.00 501.13 0.00e+00 9.977840e-01 1.015655e-02 8.914543e-02 -6.000000e+00 2.422072e-01 223.05 -2.00 3.57 -2.00 502.62 0.00e+00 18.030 0.114 489.00 513.00 1 -4.7866e+00 7.7563e-04 3.1490e-02 -6.0000e+00 1.8531e-02 -4.1599e-02 1.9402e-04 2.4157e-03 -6.0000e+00 4.7776e-03 1.6419e+01 -8.3873e-04 5.2379e-05 1.2781e-04 -6.0000e+00 1.2723e-03 4.1901e-03 1.9320e-05 12 8 8 10 0 0 99.60 22.23 8.67 0.00419 -0.00163 57 11 1 4 0 888.323 27.45 20.00 22.23 -3287 -2294 -1076 -41 684 1545 1777 1600 1615 997 256 -237 -884 -1300 -1636 -1303 -704 305 745 769 626 793 1038 1459 1323 1196 816 228 -700 -1998 -2036 -1590 -563 510 865 791 564 261 -2 -95 247 472 1133 1519 1600 1453 1339 623 -152 -945 -2070 -2774 -2583 -1717 -350 344 843 1048 996 1094 1377 1791 1622 1299 620 334 -59 -595 -887 -1265 -1587 -1291 -850 -691 -413 -654 -996 -1486 -1636 -1756 -1803 -1318 -1065 -663 -418 -791 -906 -940 -897 -231 126 677 1134 1450 1002 432 -263 -1027 -1764 -2257 -2434 -1724 -550 175 55 -318 -710 -870 -426 -215 374 1012 1558 2120 2165 2323 2207 1965 1140 495 610 748 1070 1018 1134 697 796 766 550 477 0 -583 -1238 -1368 -966 -564 -127 244 600 598 551 594 821 726 259 -166 -996 -1772 -2046 -1804 -2910 -2064 -991 -38 640 1465 1706 1503 1509 953 243 -223 -827 -1212 -1529 -1226 -653 287 712 741 607 770 983 1389 1272 1153 785 218 -676 -1909 -1931 -1523 -536 496 855 785 558 257 -2 -93 245 461 1120 1519 1615 1474 1352 627 -151 -926 -2005 -2667 -2484 -1672 -347 340 845 1057 1008 1109 1401 1828 1665 1335 636 339 -60 -594 -879 -1249 -1568 -1280 -851 -699 -418 -655 -986 -1460 -1602 -1722 -1778 -1302 -1059 -663 -421 -782 -897 -928 -886 -228 126 682 1148 1456 993 419 -258 -990 -1681 -2139 -2327 -1682 -538 176 55 -319 -710 -867 -425 -218 384 1051 1592 2175 2241 2427 2330 2036 1183 522 633 769 1100 1057 1185 736 842 811 570 494 -1 -611 -1293 -1431 -1019 -590 -135 263 639 644 598 649 900 797 285 -180 -1074 -1893 -2179 -1936 SQ 20230106T214527.0Z 8406.906563 524 1 49.25 2 89.23 237.20 88.16 54.25 0.00 1 0.00 1 0 29 1.660e-02 1.636e-02 1.579e-03 1.572e-03 7.632435e+03 9.954124e+01 4.675940e+02 -6.000000e+00 8.944698e+02 272.33 -2.00 20.00 -2.00 505.43 0.00e+00 5.591003e-03 4.288334e-05 3.030957e-04 -6.000000e+00 4.559183e-04 272.33 -2.00 20.00 -2.00 501.17 0.00e+00 8.697417e-01 3.038890e-02 9.128949e-02 -6.000000e+00 3.295499e-01 223.05 -2.00 3.57 -2.00 502.65 0.00e+00 19.392 0.142 489.00 513.00 1 -2.9960e+00 2.3184e-03 3.5144e-02 -6.0000e+00 2.5214e-02 -2.6633e-02 5.8060e-04 2.5495e-03 -6.0000e+00 6.4996e-03 2.2853e+00 8.4846e-04 1.5722e-04 1.2737e-04 -6.0000e+00 1.7312e-03 -1.8070e-02 1.6162e-04 12 9 8 11 32 0 97.05 22.23 10.37 -0.01807 -0.01218 16 3 1 4 0 3077.331 27.45 20.00 22.23 -5224 -3105 -1395 -351 119 1856 1836 2223 2357 1386 480 -398 -1565 -2115 -3198 -1959 -245 1134 419 557 736 880 879 2115 1270 869 1168 269 -738 -3150 -2418 -2178 -366 766 667 1701 704 1097 1064 108 650 -683 1881 2451 2534 2238 2442 531 -64 -259 -1790 -4371 -3371 -2221 224 -6 1592 1128 1141 1017 1507 2662 2411 1329 1098 1258 758 -547 -1436 -2036 -1555 -1466 -1363 -308 -804 -1442 -1747 -2086 -1598 -2252 -1951 -1448 -1617 -351 233 -796 -1626 -1072 -1766 44 243 775 1893 2186 1524 -238 -114 -1671 -2623 -2067 -3790 -1105 -242 878 -8 45 -816 -1403 270 -376 909 1051 1587 3666 2167 3148 2384 2394 1327 -862 61 580 1834 569 1386 1048 846 1161 952 398 -966 -104 -2255 -2615 -540 -1014 -75 451 1069 -95 161 501 1204 1443 54 -290 -1805 -1962 -1948 -2273 -4619 -2762 -1265 -326 113 1770 1699 2057 2166 1269 457 -377 -1483 -2000 -2924 -1788 -225 1050 398 533 706 846 850 2039 1244 855 1140 260 -704 -3009 -2280 -2082 -355 759 645 1634 679 1045 1009 103 627 -678 1854 2440 2546 2261 2448 535 -64 -253 -1733 -4286 -3285 -2189 222 -6 1635 1126 1142 1023 1518 2669 2432 1351 1109 1251 747 -541 -1418 -2006 -1522 -1445 -1361 -308 -808 -1440 -1723 -2042 -1550 -2201 -1910 -1421 -1605 -347 231 -800 -1648 -1078 -1788 44 247 795 1936 2216 1528 -237 -111 -1607 -2501 -1942 -3654 -1067 -241 893 -8 45 -818 -1410 268 -383 932 1096 1669 3812 2309 3357 2489 2502 1390 -909 63 614 1925 588 1485 1091 884 1212 997 418 -1017 -107 -2350 -2739 -561 -1080 -80 491 1173 -102 175 547 1311 1569 59 -315 -1952 -2090 -2065 -2444