File name: Pandora187s1_PittsburghPA_20231223_L2Fit_fout2c3d20220705p1-8.txt File generation date: 20240123T221735.5Z Data description: Level 2 spectral fitting results file Data file version: fout2c3p1-8 Local principal investigator: Lukas Valin Network principal investigator: Alexander Cede Instrument type: Pandora Instrument number: 187 Spectrometer number: 1 Processing software version used: BlickP v1.8.59 Instrument operation file used: Pandora187_OF_v4d20210517.txt Instrument calibration file used: Pandora187s1_CF_v3d20220705.txt Level 1 file used: Pandora187s1_PittsburghPA_20231223_L1_smca1c3d20220705p1-8.txt Full location name: Pittsburgh, PA - Lawrenceville ACHD (USEPA AQS ID 42-003-0008) Short location name: PittsburghPA Country of location: United States Location latitude [deg]: 40.4655 Location longitude [deg]: -79.9608 Location altitude [m]: 265 Local noon date: 20231223 Notes on s-number (L1 configuration): Corrections NOT applied although requested by the s-number are latency correction, matrix method stray light correction (replaced by simple method) Data caveats: None First and last pixel inside fitting window: 194 335 Nominal wavelengths inside fitting window [nm]: 305.0947 305.2360 305.3774 305.5187 305.6601 305.8014 305.9427 306.0840 306.2253 306.3666 306.5079 306.6491 306.7904 306.9316 307.0728 307.2141 307.3553 307.4965 307.6376 307.7788 307.9200 308.0611 308.2023 308.3434 308.4845 308.6256 308.7667 308.9078 309.0489 309.1899 309.3310 309.4720 309.6131 309.7541 309.8951 310.0361 310.1771 310.3181 310.4591 310.6000 310.7410 310.8819 311.0228 311.1637 311.3046 311.4455 311.5864 311.7273 311.8681 312.0090 312.1498 312.2907 312.4315 312.5723 312.7131 312.8539 312.9946 313.1354 313.2761 313.4169 313.5576 313.6983 313.8390 313.9797 314.1204 314.2611 314.4018 314.5424 314.6831 314.8237 314.9643 315.1049 315.2455 315.3861 315.5267 315.6673 315.8078 315.9484 316.0889 316.2294 316.3699 316.5104 316.6509 316.7914 316.9319 317.0723 317.2128 317.3532 317.4936 317.6340 317.7744 317.9148 318.0552 318.1956 318.3359 318.4763 318.6166 318.7569 318.8972 319.0375 319.1778 319.3181 319.4584 319.5986 319.7389 319.8791 320.0193 320.1596 320.2998 320.4400 320.5801 320.7203 320.8605 321.0006 321.1407 321.2809 321.4210 321.5611 321.7012 321.8412 321.9813 322.1214 322.2614 322.4014 322.5415 322.6815 322.8215 322.9615 323.1014 323.2414 323.3814 323.5213 323.6612 323.8012 323.9411 324.0810 324.2209 324.3607 324.5006 324.6405 324.7803 324.9201 --------------------------------------------------------------------------------------- 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: Ozone slant column amount [moles per square meter], -9e99=fitting not successful Column 23: 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 24: 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 25: 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 26: 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 27: Ozone effective temperature [K] Column 28: 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 29: 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 30: 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 31: Effective ozone fitting wavelength [nm], -9=fitting not successful Column 32: 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 33: Sulfur dioxide slant column amount [moles per square meter], -9e99=fitting not successful Column 34: Independent uncertainty of sulfur 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 sulfur 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 sulfur 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 sulfur 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: Sulfur dioxide effective temperature [K] Column 39: Independent uncertainty of sulfur 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 sulfur 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 sulfur 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 sulfur dioxide fitting wavelength [nm], -9=fitting not successful Column 43: Diffuse correction applied before fitting at effective fitting wavelength for sulfur dioxide [%], 0=no diffuse correction applied or fitting not requested, >0=measured diffuse correction, <0=(negative value of) calculated diffuse correction Column 44: Molecular scattering air mass factor used for molecular scattering subtraction before the fitting Column 45: Estimated uncertainty of molecular scattering air mass factor, -9=molecular scattering was not subtracted before the fitting Column 46: Lower limit used for wavelength scaling [nm] Column 47: Upper limit used for wavelength scaling [nm] Column 48: Order of smoothing polynomial Column 49: Smoothing polynomial coefficient, order 0 Column 50: Independent uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 51: Structured uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 52: Common uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 53: rms-based uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 54: Smoothing polynomial coefficient, order 1 Column 55: Independent uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 56: Structured uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 57: Common uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 58: rms-based uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 59: Smoothing polynomial coefficient, order 2 Column 60: Independent uncertainty of smoothing polynomial coefficient, order 2, -9=fitting not successful Column 61: Structured uncertainty of smoothing polynomial coefficient, order 2, -9=fitting not successful Column 62: Common uncertainty of smoothing polynomial coefficient, order 2, -9=fitting not successful Column 63: rms-based uncertainty of smoothing polynomial coefficient, order 2, -9=fitting not successful Column 64: Smoothing polynomial coefficient, order 3 Column 65: Independent uncertainty of smoothing polynomial coefficient, order 3, -9=fitting not successful Column 66: Structured uncertainty of smoothing polynomial coefficient, order 3, -9=fitting not successful Column 67: Common uncertainty of smoothing polynomial coefficient, order 3, -9=fitting not successful Column 68: rms-based uncertainty of smoothing polynomial coefficient, order 3, -9=fitting not successful Column 69: Mean value of measured data inside fitting window [same units as measurements] Column 70: Wavelength change polynomial coefficient, order 0 Column 71: Independent uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 72: Structured uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 73: Common uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 74: rms-based uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 75: L1 based wavelength change coefficient, order 0, -9=no L1 wavelength change determination Column 76: L1 based wavelength change coefficient, order 1, -9=no L1 wavelength change determination Column 77: 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 78: 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 79: 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 80: 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 81: 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 82: 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 83: Atmospheric variability [%], 999=no atmospheric variability was determined Column 84: Wavelength effective temperature [°C], 999=no effective temperature given Column 85: Estimated average residual stray light level [%] (only valid for stray light correction methods 2 and higher) Column 86: Retrieved wavelength shift from L1 data [nm], -9=no wavelength change determination Column 87: Retrieved total wavelength shift [nm], -9=no wavelength change fitting Column 88: Number of bright count cycles Column 89: Number of dark count cycles Column 90: Effective position of filterwheel #1, 0=filterwheel not used, 1-9 are valid positions Column 91: Effective position of filterwheel #2, 0=filterwheel not used, 1-9 are valid positions Column 92: Sum over 2^i, 0=spectra were interpolated in time, 1=spectra are corrected for off-target signal Column 93: Integration time [ms] Column 94: Temperature at electronics board [°C], 999=no temperature signal Column 95: Spectrometer control temperature [°C], 999=no temperature signal Column 96: Auxiliary spectrometer temperature [°C], 999=no temperature signal Column 97: Temperature in head sensor [°C], 999=no temperature signal Column 98: Humidity in head sensor [%], -9=no humidity signal Column 99: Pressure in head sensor [hPa], -9=no pressure signal Columns 100-241: Unweighted slant column residuals for each pixel inside the fitting window multiplied by 1e5, 9e5=pixel was not used for fitting Columns 242-383: 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 --------------------------------------------------------------------------------------- SS 20231223T144137.2Z 8757.612236 246 1 6.02 2 73.62 142.71 120.11 5.62 0.00 1 0.00 1 0 25 1.919e-01 1.362e-02 3.224e-01 1.948e-02 4.870153e-01 1.005803e-02 5.855527e-03 -6.000000e+00 1.094646e-01 222.99 -2.00 3.49 -2.00 319.88 0.00e+00 -4.312052e-04 3.110318e-03 8.390543e-04 -6.000000e+00 4.785075e-03 275.53 -2.00 20.00 -2.00 319.35 0.00e+00 3.489 0.001 302.00 328.00 3 3.8548e+00 2.8703e-02 4.2006e-02 -6.0000e+00 2.6969e-01 1.5461e-02 3.4428e-02 2.5397e-02 -6.0000e+00 2.9353e-01 -1.8953e-02 3.0019e-02 1.2938e-02 -6.0000e+00 1.7471e-01 6.4528e-03 1.5671e-02 3.8367e-03 -6.0000e+00 7.5565e-02 2.1623e-04 -1.1775e-03 5.6737e-04 4.9272e-05 -6.0000e+00 5.3670e-03 -2.4977e-02 0.0000e+00 12 9 9 12 1 1 999.00 18.82 0.47 -0.02498 -0.03383 1 5 1 5 0 3003.770 22.19 14.99 18.82 13.15 28.0 953.37 900000 -46730 -31716 84617 -54270 -64085 -84064 12203 58734 66753 -39329 46309 -54379 22592 -45748 23029 -4583 11202 -40610 117 -40351 -15165 -34264 1942 9686 -18580 21742 -11699 2006 -3490 4532 -13732 8706 -1245 -767 1670 2663 -7323 -6761 5243 1622 6000 7445 5701 -3898 10387 1714 -2672 -4292 5079 2181 297 -549 3743 -1485 4338 -278 -3209 3030 -5648 1974 7237 8512 -697 7451 -2231 2434 -471 215 -1317 433 -42 2169 136 1315 -279 -505 892 -3469 1074 -3120 97 -2551 -2238 -3333 -789 873 -761 1008 3072 965 740 -171 -638 -243 -2455 704 127 -1058 -230 569 -852 -185 -687 -597 668 -1077 524 534 582 -1162 2174 193 372 -124 1126 -838 350 -267 -723 -501 -448 -470 -744 421 -1284 943 1466 1068 183 1679 736 441 105 621 -914 -946 -14 -50 160 -1040 199 900000 -887 -602 1607 -1030 -1217 -2128 232 1115 1267 -933 879 -1979 465 -2244 656 -189 425 -2763 6 -2937 -864 -2494 117 585 -1440 1238 -888 128 -247 336 -1398 796 -131 -77 164 290 -1049 -1179 959 313 1177 1496 1182 -894 2277 377 -621 -1039 1289 573 87 -176 1238 -543 1551 -107 -1254 1084 -2189 740 2828 3663 -357 3549 -1116 1132 -232 117 -825 289 -29 1438 87 755 -147 -249 440 -1998 687 -2379 80 -2138 -1981 -3198 -785 893 -806 1035 2930 869 673 -151 -586 -228 -2686 847 159 -1396 -291 711 -981 -202 -793 -740 859 -1383 705 783 932 -1889 3464 325 589 -184 1605 -1246 542 -442 -1235 -843 -730 -741 -1124 602 -1786 1265 2090 1692 319 2866 1317 826 209 1223 -1892 -1998 -29 -103 337 -2103 381 SS 20231223T145229.6Z 8757.619787 249 1 45.91 2 72.40 144.98 119.84 8.42 0.00 1 0.00 1 0 25 9.588e-02 9.030e-03 1.235e-01 5.120e-03 4.672784e-01 2.615973e-03 5.518616e-03 -6.000000e+00 5.469235e-02 222.99 -2.00 3.49 -2.00 319.65 0.00e+00 4.271327e-04 7.576610e-04 4.327009e-04 -6.000000e+00 2.378903e-03 275.53 -2.00 20.00 -2.00 319.12 0.00e+00 3.261 0.001 302.00 328.00 3 4.1417e+00 7.3753e-03 3.9471e-02 -6.0000e+00 1.3475e-01 3.0106e-02 8.7981e-03 2.4022e-02 -6.0000e+00 1.4663e-01 1.4182e-02 7.5776e-03 1.2243e-02 -6.0000e+00 8.7270e-02 -1.6985e-02 3.9951e-03 3.6643e-03 -6.0000e+00 3.7468e-02 2.0076e-04 -1.0835e-03 1.4925e-04 4.2145e-05 -6.0000e+00 2.6668e-03 -2.4043e-02 0.0000e+00 11 9 0 11 40 0 99.38 18.81 0.44 -0.02404 -0.03218 13 4 1 5 0 3531.040 22.19 15.00 18.81 14.10 28.1 960.76 -17768 38743 18727 -2974 -43940 34529 -19543 -1545 21253 -5098 -25392 31635 -25214 12367 -6340 5205 -16304 10769 -8570 9767 8510 -15532 -26047 4161 -9916 -15417 7982 -2685 -6563 -10700 -10402 -4960 34088 3962 -6114 735 3967 -5931 -1927 -224 -4685 -26 2022 363 -104 4992 -2003 -700 -2732 226 2191 -331 -2786 1798 4356 2357 2864 1292 2079 -2596 644 1731 3622 1693 3039 -1109 197 -1 702 861 874 54 170 932 -1589 -598 -2495 -1268 -2420 -1083 -2071 -776 -2890 -2359 -644 -104 881 1159 1440 1804 1170 -1124 702 316 251 -383 265 -564 79 -72 -640 -187 -87 271 -1334 -1227 -1119 -251 992 345 -85 831 339 -29 59 -188 106 -167 -96 -117 -609 533 -604 145 374 -1012 904 878 623 227 579 157 215 623 213 -894 51 -461 -196 -120 387 -396 -234 332 207 -39 -799 283 -300 -23 280 -95 -598 498 -782 347 -240 218 -861 505 -520 518 467 -1092 -2141 275 -811 -1403 579 -207 -509 -908 -999 -514 2684 438 -714 81 471 -925 -352 -47 -1150 -6 512 95 -28 1267 -547 -191 -775 66 681 -108 -1001 662 1634 933 1151 523 813 -1060 266 777 1776 917 1630 -589 100 0 404 561 617 40 122 625 -996 -334 -1330 -675 -1452 -743 -1637 -629 -2562 -2200 -633 -108 946 1264 1545 1826 1102 -1004 632 276 244 -424 309 -737 107 -94 -776 -219 -98 318 -1602 -1546 -1486 -351 1488 534 -142 1379 546 -47 88 -274 159 -263 -151 -190 -978 820 -966 220 541 -1420 1226 1270 997 373 1007 286 384 1171 426 -1873 108 -991 -412 -240 735 -770 SS 20231223T145949.2Z 8757.624874 253 1 46.51 2 71.61 146.54 119.60 10.30 0.00 1 0.00 1 0 25 1.445e-01 1.578e-02 7.153e-02 7.804e-03 4.419413e-01 3.950226e-03 5.258941e-03 -6.000000e+00 8.240136e-02 222.99 -2.00 3.49 -2.00 319.49 0.00e+00 2.277995e-03 1.046673e-03 3.278773e-04 -6.000000e+00 3.584601e-03 275.53 -2.00 20.00 -2.00 318.93 0.00e+00 3.130 0.001 302.00 328.00 3 4.8412e+00 1.0941e-02 3.7586e-02 -6.0000e+00 2.0302e-01 1.4340e-01 1.2814e-02 2.2458e-02 -6.0000e+00 2.2092e-01 -1.4760e-01 9.9175e-03 1.1586e-02 -6.0000e+00 1.3149e-01 6.2058e-02 5.2111e-03 3.9820e-03 -6.0000e+00 5.6443e-02 1.1412e-04 -1.0238e-03 2.2885e-04 4.4201e-05 -6.0000e+00 4.0188e-03 -2.2761e-02 0.0000e+00 12 9 8 11 32 0 97.48 18.81 -0.62 -0.02276 -0.03045 12 4 1 5 0 3874.990 22.19 14.99 18.81 14.85 28.4 966.38 -41885 -782 -17361 -26607 -66661 -61347 -38429 -55230 -18203 -44897 -35266 -26238 -47083 -19268 -13482 -8676 -14787 20034 2817 13671 20074 -18980 -20203 -1762 -13240 -13234 -9647 2385 -5740 -11107 -12029 -6029 45236 5964 -4118 -3481 3092 2300 2257 5521 -1698 2903 6667 1851 4066 3679 5039 173 3187 -2845 9368 -1489 -769 4845 5911 4725 3279 3068 2451 -3069 565 2799 5666 3749 4430 -1369 -1487 -313 1780 523 1344 410 -183 3080 -2646 -3523 -4183 -4101 -3944 -3215 -3966 -1898 -3500 -3185 364 -653 1978 2059 1474 2113 387 246 144 165 -1956 -1142 -1376 -1048 -119 610 -883 -1537 -699 -697 -2069 -2119 -882 -92 1613 806 221 1935 -69 -167 -619 -483 442 -338 841 445 -759 1444 239 242 -684 -127 1261 2177 1303 364 774 1117 -3 793 -502 -556 -144 -608 42 -469 -53 -1492 -1502 -21 -566 -897 -3007 -2534 -1337 -2521 -633 -2147 -1672 -1229 -2994 -1207 -947 -660 -1290 1313 231 1040 1445 -2000 -2223 -171 -1519 -1572 -1101 229 -574 -1222 -1500 -797 3971 796 -574 -481 442 389 460 1267 -458 787 1799 530 1142 1044 1407 51 920 -909 2873 -554 -304 1923 2405 1994 1422 1314 1029 -1356 252 1335 2909 2111 2471 -763 -802 -170 1044 351 965 305 -134 2053 -1701 -2052 -2289 -2261 -2413 -2254 -3185 -1537 -3072 -2943 347 -669 2061 2182 1549 2089 359 228 127 151 -1872 -1228 -1678 -1411 -164 809 -1097 -1750 -760 -797 -2574 -2769 -1201 -131 2296 1259 367 3191 -112 -253 -943 -712 661 -501 1311 718 -1225 2212 384 371 -1012 -179 1717 3124 2077 594 1327 1975 -5 1545 -985 -1127 -296 -1269 85 -976 -104 -2858