File name: Pandora118s1_Cabauw_20220209_L2Fit_fout2c7d20210820p1-8.txt
File generation date: 20241001T103113.3Z
Data description: Level 2 spectral fitting results file
Data file version: fout2c7p1-8
Local principal investigator: Ankie Piters
Network principal investigator: Alexander Cede
Instrument type: Pandora
Instrument number: 118
Spectrometer number: 1
Processing software version used: BlickP v1.8.60
Instrument operation file used: Pandora118_OF_v3d20200505.txt
Instrument calibration file used: Pandora118s1_CF_v7d20210820.txt
Level 1 file used: Pandora118s1_Cabauw_20220209_L1_smca1c7d20210820p1-8.txt
Full location name: Cabauw
Short location name: Cabauw
Country of location: Netherlands
Location latitude [deg]: 51.9704
Location longitude [deg]: 4.9262
Location altitude [m]: 0
Local noon date: 20220209
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: newDate as OF is out there
First and last pixel inside fitting window: 234 375
Nominal wavelengths inside fitting window [nm]: 305.0994 305.2405 305.3816 305.5226 305.6636 305.8047 305.9457 306.0867 306.2277 306.3686 306.5096 306.6506 306.7915 306.9324 307.0734 307.2143 307.3552 307.4961 307.6369 307.7778 307.9186 308.0595 308.2003 308.3411 308.4819 308.6227 308.7635 308.9043 309.0451 309.1858 309.3266 309.4673 309.6080 309.7487 309.8894 310.0301 310.1708 310.3114 310.4521 310.5927 310.7333 310.8739 311.0145 311.1551 311.2957 311.4363 311.5769 311.7174 311.8579 311.9985 312.1390 312.2795 312.4200 312.5604 312.7009 312.8414 312.9818 313.1222 313.2627 313.4031 313.5435 313.6839 313.8242 313.9646 314.1050 314.2453 314.3856 314.5259 314.6663 314.8065 314.9468 315.0871 315.2274 315.3676 315.5079 315.6481 315.7883 315.9285 316.0687 316.2089 316.3490 316.4892 316.6293 316.7695 316.9096 317.0497 317.1898 317.3299 317.4700 317.6101 317.7501 317.8902 318.0302 318.1702 318.3102 318.4502 318.5902 318.7302 318.8701 319.0101 319.1500 319.2900 319.4299 319.5698 319.7097 319.8496 319.9894 320.1293 320.2691 320.4090 320.5488 320.6886 320.8284 320.9682 321.1080 321.2477 321.3875 321.5272 321.6670 321.8067 321.9464 322.0861 322.2258 322.3654 322.5051 322.6448 322.7844 322.9240 323.0636 323.2032 323.3428 323.4824 323.6220 323.7615 323.9011 324.0406 324.1801 324.3197 324.4592 324.5986 324.7381 324.8776
---------------------------------------------------------------------------------------
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 20220209T110423.3Z 8075.461381 245 1 49.68 2 67.45 166.86 85.10 63.90 0.00 1 0.00 1 0 17 2.874e-01 2.806e-02 1.230e-01 1.295e-02 3.847266e-01 6.569955e-03 5.069475e-03 -6.000000e+00 1.642595e-01 223.30 -2.00 3.84 -2.00 318.87 0.00e+00 -1.050473e-03 1.405608e-03 9.093716e-04 -6.000000e+00 7.224498e-03 266.96 -2.00 20.00 -2.00 318.45 0.00e+00 2.587 0.000 302.00 328.00 3 6.7114e+00 1.7721e-02 3.3471e-02 -6.0000e+00 4.0786e-01 1.4028e-01 2.0353e-02 2.1506e-02 -6.0000e+00 4.4545e-01 -1.8918e-03 1.5283e-02 1.1206e-02 -6.0000e+00 2.6015e-01 -1.0309e-02 8.1819e-03 3.5184e-03 -6.0000e+00 1.1098e-01 2.8383e-05 -4.4449e-04 3.7677e-04 3.7731e-05 -6.0000e+00 7.9057e-03 1.0561e-02 0.0000e+00 12 8 8 10 0 0 99.98 21.17 3.22 0.01056 0.00722 21 5 1 5 0 2365.270 21.17 15.00 17.25 14.88 29.0 1006.18 25560 -14425 18634 -37695 900000 -64452 900000 900000 279844 900000 -61703 -48420 3438 10156 17275 40463 3237 44480 -20037 17180 11429 -2579 16841 -7430 -14349 -16719 32092 19042 34764 -50737 42522 7678 10431 36711 -5091 -8161 -16332 21774 9278 6190 3449 1786 6806 -7438 14919 6340 -3387 2344 -2131 3184 995 27557 -2860 -1607 -4987 -393 11110 2603 4784 -128 1585 353 1593 554 187 -2597 -695 -3082 -5363 -545 -719 1065 -1681 -2129 5082 3232 699 -8091 1201 -2337 1054 -2157 -4722 -2062 -2840 2110 -2883 -1528 2426 -320 -998 -1909 1920 1630 1115 6445 3901 534 -467 -650 2839 812 -2785 -1844 264 3003 4086 -2736 127 -198 -728 -19 48 -1876 -379 -1970 -44 1145 -953 -2580 1751 -773 1199 4855 877 -2123 -1019 -1636 1000 -485 896 -372 2573 -1900 938 504 114 482 784 935 -1855 -1838 720 -645 620 -1963 900000 -3999 900000 900000 3530 900000 -5527 -4170 213 688 1265 2635 315 3110 -2680 1683 1194 -307 1778 -1062 -2316 -2875 3438 2159 3253 -11031 4272 1176 1711 4697 -905 -1532 -3719 4115 2368 1868 1148 628 2345 -2934 4773 2200 -1269 830 -790 1172 405 9391 -1379 -809 -2674 -208 5353 1332 2345 -66 837 205 998 371 122 -1630 -415 -1914 -3762 -411 -578 859 -1357 -1634 3339 1952 404 -5091 766 -1749 854 -1963 -4628 -2068 -3066 2297 -3402 -1806 2733 -349 -1012 -1815 1677 1462 1098 6715 4629 711 -642 -860 3383 907 -3115 -2174 325 3753 5280 -3944 190 -318 -1225 -32 76 -2863 -534 -2730 -61 1655 -1486 -4198 2721 -1197 1770 6635 1193 -2833 -1312 -2233 1445 -755 1443 -619 4288 -3445 1750 970 223 955 1525 1748 -3414 -3273
SS 20220209T111954.4Z 8075.472158 256 1 2.02 2 66.99 170.90 82.96 66.73 0.00 1 0.00 1 0 17 1.769e-02 6.506e-03 1.081e-02 2.408e-03 3.630462e-01 1.202469e-03 4.817031e-03 -6.000000e+00 9.938359e-03 223.30 -2.00 3.84 -2.00 317.76 0.00e+00 4.308219e-04 1.886980e-04 1.117094e-04 -6.000000e+00 4.174576e-04 266.96 -2.00 20.00 -2.00 317.34 0.00e+00 2.538 0.000 302.00 328.00 3 1.4291e+00 3.1635e-03 3.2242e-02 -6.0000e+00 2.4585e-02 4.0757e-02 3.5606e-03 2.0056e-02 -6.0000e+00 2.6796e-02 -1.6158e-02 2.3241e-03 1.0122e-02 -6.0000e+00 1.5862e-02 -3.1132e-04 1.2176e-03 3.0008e-03 -6.0000e+00 6.8174e-03 6.4413e-03 -1.0399e-04 6.7215e-05 3.2801e-05 -6.0000e+00 4.7147e-04 1.1514e-02 0.0000e+00 12 9 1 12 1 1 99.93 21.17 0.26 0.01151 0.01073 15 93 1 5 0 114.830 21.17 15.00 17.26 15.19 29.0 1008.73 -4746 -2781 -5050 -502 -6567 -7737 -8395 -5718 -1709 -4450 -1228 -5191 -3107 -6 -2631 -995 -264 1141 777 434 -308 -664 -54 -854 -121 265 -1117 -3049 -976 -175 946 2227 1248 -1085 -2983 -2227 -1621 -676 -3 -410 -433 -786 -35 -223 1098 -54 523 -321 107 953 766 526 1185 1052 1208 1772 270 -111 163 733 785 1417 1920 838 -403 -412 -699 851 882 23 780 173 91 258 -174 -368 -621 134 -690 -717 -1294 -1960 -1968 -1238 -671 -303 -144 -41 409 -213 124 393 885 931 576 343 470 572 -273 -297 -350 -485 151 -581 -1357 -1219 -935 -344 187 225 -72 -106 278 150 359 404 752 270 -338 -301 -290 280 551 729 677 507 1347 980 461 303 -66 59 53 -183 -58 -494 337 -251 -717 44 -244 -499 -331 -207 -389 -34 -425 -510 -613 -461 -147 -408 -110 -529 -365 -1 -388 -160 -44 200 143 83 -60 -130 -11 -189 -28 64 -266 -691 -208 -39 230 591 339 -300 -849 -648 -490 -242 -1 -189 -215 -409 -17 -111 580 -27 275 -168 56 483 423 314 747 659 778 1165 183 -74 106 475 537 1054 1470 687 -324 -307 -535 654 711 20 721 165 84 225 -146 -278 -466 100 -546 -615 -1258 -1999 -2118 -1328 -771 -351 -173 -50 491 -241 135 413 875 946 616 393 579 775 -366 -398 -432 -574 172 -717 -1732 -1623 -1251 -483 269 349 -112 -165 428 215 496 540 1028 393 -508 -468 -448 418 796 1042 924 689 1771 1373 678 458 -102 93 86 -303 -99 -867 599 -456 -1287 75 -415 -818
SS 20220209T113655.4Z 8075.483974 264 1 48.55 2 66.67 175.37 80.56 69.80 0.00 1 0.00 1 0 17 2.767e-01 1.704e-02 1.279e-01 8.876e-03 3.681583e-01 4.451980e-03 4.817203e-03 -6.000000e+00 1.559912e-01 223.30 -2.00 3.84 -2.00 318.79 0.00e+00 -1.240316e-03 9.311940e-04 3.340353e-04 -6.000000e+00 6.649665e-03 266.96 -2.00 20.00 -2.00 318.36 0.00e+00 2.506 0.000 302.00 328.00 3 6.2342e+00 1.1978e-02 3.2191e-02 -6.0000e+00 3.8559e-01 8.5240e-02 1.3707e-02 2.0368e-02 -6.0000e+00 4.2027e-01 2.1779e-02 1.0241e-02 1.0492e-02 -6.0000e+00 2.4839e-01 -1.9502e-02 5.5103e-03 3.2844e-03 -6.0000e+00 1.0680e-01 5.2018e-05 4.4820e-04 2.5684e-04 3.4845e-05 -6.0000e+00 7.4023e-03 7.2418e-03 0.0000e+00 12 8 8 10 0 0 100.00 21.17 1.62 0.00724 0.01061 32 7 1 5 0 1516.980 21.17 15.00 17.29 15.22 29.0 1008.33 17592 -19048 14310 -3466 900000 -12436 239938 18860 121263 145944 -15945 -9881 16542 7501 -12963 55562 5166 16723 5824 14047 -595 -2694 704 -17289 -12821 -16371 3053 16851 8873 -33500 13689 3898 -1542 8653 -5189 -5212 -5649 4246 4345 1142 610 3272 455 -7382 3570 1532 -4293 1741 -4293 1042 3897 11043 952 -1679 -1011 -1343 4629 3503 5678 -233 1691 2917 8 2357 2270 -4942 658 -4332 465 -2052 -1435 -205 -1092 308 2251 2291 -577 -3779 189 -923 902 -2477 -2122 -1622 -172 351 -1148 -691 1272 1413 21 649 -270 407 1243 4103 2950 307 837 -1769 28 -87 -3266 -766 -1933 1774 705 -1723 -815 597 -178 -302 911 -702 363 -538 137 527 -1342 -1334 581 -362 1117 3241 -212 -2014 1064 -263 1536 -405 863 -639 1283 -832 -252 703 -525 637 -76 421 -557 -1056 568 -937 524 -138 900000 -502 2076 720 1987 1984 -985 -628 956 548 -1314 3279 512 1575 635 1471 -72 -334 90 -2819 -2129 -2923 447 2029 1103 -6430 1867 641 -291 1482 -950 -985 -1207 981 1204 373 215 1174 172 -3002 1312 574 -1673 638 -1672 403 1587 4515 455 -868 -535 -733 2430 1827 2841 -123 920 1703 5 1588 1491 -3237 396 -2778 316 -1599 -1184 -171 -896 236 1552 1430 -346 -2347 125 -697 746 -2304 -2062 -1647 -184 394 -1353 -824 1469 1542 22 614 -245 376 1245 4437 3561 413 1142 -2370 34 -99 -3731 -909 -2469 2277 954 -2490 -1241 958 -300 -506 1453 -1062 510 -742 192 769 -2114 -2161 918 -564 1671 4555 -295 -2725 1372 -360 2234 -635 1401 -1071 2175 -1500 -476 1355 -1039 1264 -149 794 -1023 -1884