File name: Pandora198s1_Kobe_20231008_L2Fit_fout2c1d20220209p1-8.txt File generation date: 20231008T235126.2Z Data description: Level 2 spectral fitting results file Data file version: fout2c1p1-8 Local principal investigator: Kazuyo Yamaji Network principal investigator: Alexander Cede Instrument type: Pandora Instrument number: 198 Spectrometer number: 1 Processing software version used: BlickP v1.8.49 Instrument operation file used: Pandora198_OF_v4d20210227.txt Instrument calibration file used: Pandora198s1_CF_v1d20220209.txt Level 1 file used: Pandora198s1_Kobe_20231008_L1_smca1c1d20220209p1-8.txt Full location name: Kobe University Short location name: Kobe Country of location: Japan Location latitude [deg]: 34.7190 Location longitude [deg]: 135.2900 Location altitude [m]: 23 Local noon date: 20231008 Notes on s-number (L1 configuration): Corrections NOT applied although requested by the s-number are dark map for dark correction (replaced by immediate dark method), latency correction, matrix method stray light correction (replaced by simple method) Data caveats: None First and last pixel inside fitting window: 200 342 Nominal wavelengths inside fitting window [nm]: 305.0180 305.1589 305.2998 305.4407 305.5815 305.7224 305.8632 306.0040 306.1448 306.2856 306.4264 306.5672 306.7080 306.8488 306.9895 307.1302 307.2710 307.4117 307.5524 307.6931 307.8338 307.9745 308.1151 308.2558 308.3964 308.5371 308.6777 308.8183 308.9589 309.0995 309.2401 309.3806 309.5212 309.6617 309.8023 309.9428 310.0833 310.2238 310.3643 310.5048 310.6453 310.7858 310.9262 311.0666 311.2071 311.3475 311.4879 311.6283 311.7687 311.9091 312.0494 312.1898 312.3301 312.4705 312.6108 312.7511 312.8914 313.0317 313.1720 313.3123 313.4525 313.5928 313.7330 313.8732 314.0135 314.1537 314.2939 314.4340 314.5742 314.7144 314.8545 314.9947 315.1348 315.2749 315.4150 315.5551 315.6952 315.8353 315.9754 316.1154 316.2555 316.3955 316.5355 316.6755 316.8155 316.9555 317.0955 317.2355 317.3754 317.5154 317.6553 317.7952 317.9351 318.0751 318.2149 318.3548 318.4947 318.6346 318.7744 318.9142 319.0541 319.1939 319.3337 319.4735 319.6133 319.7530 319.8928 320.0326 320.1723 320.3120 320.4517 320.5914 320.7311 320.8708 321.0105 321.1502 321.2898 321.4295 321.5691 321.7087 321.8483 321.9879 322.1275 322.2671 322.4066 322.5462 322.6857 322.8252 322.9648 323.1043 323.2438 323.3833 323.5227 323.6622 323.8016 323.9411 324.0805 324.2199 324.3593 324.4987 324.6381 324.7775 324.9169 --------------------------------------------------------------------------------------- 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-242: Unweighted slant column residuals for each pixel inside the fitting window multiplied by 1e5, 9e5=pixel was not used for fitting Columns 243-385: 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 20231008T010809.6Z 8681.047334 393 1 2.44 2 46.62 145.24 39.95 268.20 0.00 1 0.00 1 0 33 9.464e-02 8.339e-02 2.611e-03 1.751e-03 1.325759e-01 1.038784e-03 8.138040e-04 -6.000000e+00 6.444337e-02 227.49 -2.00 3.38 -2.00 315.17 0.00e+00 6.978489e-03 8.760446e-05 3.675829e-05 -6.000000e+00 2.656098e-03 289.31 -2.00 20.00 -2.00 314.87 0.00e+00 1.454 0.000 302.00 328.00 3 2.1253e+00 2.6571e-03 1.5243e-02 -6.0000e+00 1.6083e-01 -6.4164e-02 2.9250e-03 4.8878e-03 -6.0000e+00 1.7426e-01 1.0172e-02 1.6805e-03 2.1738e-03 -6.0000e+00 1.0245e-01 -2.8619e-02 7.0872e-04 5.9087e-04 -6.0000e+00 4.1028e-02 1.3787e-02 -4.9274e-03 7.1339e-05 2.2197e-06 -6.0000e+00 3.9260e-03 -2.5773e-02 0.0000e+00 12 9 9 12 1 1 51.31 14.61 0.40 -0.02577 -0.06280 10 52 1 5 0 221.360 20.46 15.00 14.61 29.39 20.3 969.71 900000 14504 16277 12714 7344 -454 -6760 -10582 -11268 -11249 -13080 -14027 -11089 -4832 1961 7280 10742 12744 13990 14134 13518 12098 9104 7190 3734 2041 -1634 -6222 -10780 -14049 -15470 -15759 -14243 -15108 -17286 -19497 -18304 -13555 -6249 1248 7772 10557 12189 13405 12955 11901 9829 6929 3590 784 -502 383 1337 2802 2920 2725 1788 982 -143 -520 825 3969 6785 9218 8559 5277 418 -4244 -5028 -2441 1839 4931 5301 3601 -801 -7182 -14312 -21011 -23628 -21175 -16479 -10267 -5640 -935 3469 7791 11496 14245 14873 12849 8642 2217 -4730 -10663 -12445 -10383 -5583 -680 2798 3506 1457 -2244 -5644 -6401 -4336 -785 2728 6201 9568 12195 14491 14310 12204 7839 2423 -3096 -6169 -5634 -3112 1060 4136 5091 4180 1865 -1662 -5852 -9639 -12105 -12099 -9957 -7544 -5710 -4034 -2643 -173 2399 4562 5864 6370 5274 3110 658 -1041 900000 4101 4754 3851 2257 -143 -2183 -3552 -3919 -4029 -4872 -5414 -4030 -1808 763 2954 4526 5533 6215 6448 6325 5783 4477 3596 1937 996 -822 -3223 -5693 -7533 -8442 -8863 -8186 -8993 -10573 -12194 -11643 -8238 -3942 819 5269 7467 8892 9965 9197 8601 7225 5167 2714 601 -391 303 1085 2320 2478 2353 1472 816 -120 -443 712 3471 6079 8412 7974 4968 393 -4008 -4753 -2338 1798 4649 5125 3530 -793 -7118 -14135 -20830 -23532 -21263 -16953 -10779 -6097 -972 3669 8394 12632 15888 16824 14710 9917 2542 -5392 -12152 -14202 -12016 -6597 -782 3299 4208 1762 -2719 -6830 -7744 -5289 -973 3440 7933 12381 15219 18263 18330 15731 10106 3097 -3937 -7810 -7133 -3995 1379 5467 6834 5663 2532 -2252 -7920 -13049 -16467 -16607 -13819 -10622 -8175 -5570 -3710 -247 3481 6749 8838 9298 7815 4662 995 -1586 SS 20231008T011255.1Z 8681.050638 396 1 28.01 2 46.08 146.69 40.91 268.87 0.00 1 0.00 1 0 33 9.697e-02 8.126e-02 2.721e-03 1.471e-03 1.239776e-01 8.589145e-04 7.189671e-04 -6.000000e+00 6.606919e-02 227.49 -2.00 3.38 -2.00 316.10 0.00e+00 9.369350e-03 8.749056e-05 4.185121e-05 -6.000000e+00 2.721948e-03 289.31 -2.00 20.00 -2.00 315.77 0.00e+00 1.439 0.000 302.00 328.00 3 3.8162e+00 2.2148e-03 1.4979e-02 -6.0000e+00 1.6489e-01 -5.5744e-02 2.4492e-03 4.5172e-03 -6.0000e+00 1.7863e-01 1.4037e-02 1.4130e-03 2.0159e-03 -6.0000e+00 1.0502e-01 -3.1540e-02 6.1838e-04 5.6777e-04 -6.0000e+00 4.2081e-02 2.5844e-03 -4.5337e-03 6.0187e-05 2.7780e-06 -6.0000e+00 4.0231e-03 -2.9377e-02 0.0000e+00 12 9 9 12 1 1 99.99 14.62 0.49 -0.02938 -0.06345 33 14 1 5 0 823.630 20.56 15.01 14.62 29.69 20.4 972.94 900000 16941 19480 14363 9847 -176 -4810 -10238 -11185 -13197 -14891 -16847 -12825 -5931 1321 7806 11513 15166 14777 16682 14734 14040 9715 7537 3072 143 -4040 -7762 -12064 -15119 -15289 -15420 -13610 -14398 -16001 -18657 -17276 -13575 -5987 599 6592 9318 11499 12851 12493 11793 10148 8022 3785 1818 342 1135 1777 3420 2949 2717 1343 96 -1064 -1486 373 3813 7558 9274 9590 5537 833 -4136 -4724 -2459 1928 4773 5569 3690 -475 -6732 -14431 -20643 -23542 -21153 -16087 -10453 -5337 -685 3106 7833 11529 13953 14558 12777 8904 2463 -4682 -10696 -12442 -10903 -5689 -938 2869 3461 1486 -2219 -5494 -6305 -4075 -627 2756 6524 9440 12156 14049 14333 11830 7717 1924 -3384 -6450 -6171 -3475 633 3933 5151 4515 2087 -1398 -5808 -9413 -12045 -11688 -9902 -7486 -5906 -4230 -2653 -408 2230 4508 5813 6023 5445 3100 1014 -826 900000 3418 4051 3155 2176 -40 -1110 -2466 -2767 -3399 -3959 -4682 -3636 -1733 403 2469 3792 5080 5160 5882 5396 5201 3731 2927 1238 59 -1721 -3375 -5358 -6845 -7003 -7303 -6589 -7222 -8198 -9803 -9186 -7473 -3395 355 4044 5958 7545 7820 7799 7500 6559 5234 2529 1222 234 787 1266 2474 2193 2049 1031 74 -838 -1186 300 3113 6277 7910 8275 4863 730 -3641 -4147 -2190 1748 4442 5289 3559 -462 -6544 -14061 -20097 -23092 -20933 -16245 -10826 -5658 -740 3428 8768 13161 14875 15752 13963 9722 2686 -5082 -11599 -14717 -12042 -6389 -1084 3388 4168 1802 -2695 -6655 -7636 -4967 -776 3477 8325 12233 15965 18695 19276 16060 10442 2590 -4517 -8570 -8220 -4686 866 5459 7246 6395 2966 -1982 -8245 -13339 -17178 -16754 -14394 -11034 -8868 -6429 -4087 -638 3554 7317 8876 9361 8549 4937 1624 -1334 SS 20231008T011725.1Z 8681.053762 399 1 48.52 2 45.58 148.09 41.82 269.49 0.00 1 0.00 1 0 33 9.994e-02 7.959e-02 6.585e-03 3.022e-03 1.283861e-01 1.749519e-03 7.691433e-04 -6.000000e+00 6.808849e-02 227.49 -2.00 3.38 -2.00 316.73 0.00e+00 9.610375e-03 2.036190e-04 4.676497e-05 -6.000000e+00 2.805275e-03 289.31 -2.00 20.00 -2.00 316.38 0.00e+00 1.427 0.000 302.00 328.00 3 4.9468e+00 4.5383e-03 1.4924e-02 -6.0000e+00 1.6993e-01 -3.7880e-02 5.0332e-03 4.7003e-03 -6.0000e+00 1.8409e-01 6.6861e-03 2.9350e-03 2.1397e-03 -6.0000e+00 1.0823e-01 -2.5983e-02 1.3098e-03 6.2117e-04 -6.0000e+00 4.3365e-02 8.3212e-04 -4.5734e-03 1.2510e-04 2.4392e-06 -6.0000e+00 4.1463e-03 -2.8893e-02 0.0000e+00 12 8 8 10 0 0 99.99 14.62 0.49 -0.02889 -0.06326 47 11 1 5 0 1031.990 20.56 15.00 14.62 30.11 20.4 976.11 900000 19886 22664 19761 10345 2937 -5973 -7577 -11932 -13172 -16086 -17069 -13740 -7750 2111 7872 13695 14717 17212 17792 16007 13418 11543 5871 2679 649 -4631 -11096 -13119 -14471 -15676 -15648 -13256 -13222 -16883 -17649 -18684 -13450 -5972 188 5814 8592 10559 12113 12066 11508 10454 7916 5388 1884 929 821 1866 2786 2267 1752 815 -378 -1899 -1612 383 4144 6580 9981 8866 5906 587 -3895 -4935 -2888 1356 4027 5171 3151 -756 -7370 -15204 -21600 -23352 -21642 -15846 -10556 -5330 -1025 3296 7593 10892 14024 14163 13218 8422 2602 -5134 -10875 -12458 -10158 -5650 -519 2513 3388 1140 -2155 -5684 -6749 -4166 -716 2796 6097 9609 12089 14247 14334 11922 7774 1772 -3434 -6520 -6045 -3600 817 3837 5282 4166 1872 -2135 -6075 -9587 -11708 -11786 -10042 -7863 -6017 -4423 -3126 -589 1855 4300 5582 6313 5010 3242 768 -1231 900000 3138 3677 3431 1896 546 -1160 -1476 -2465 -2807 -3578 -3930 -3251 -1909 529 2065 3665 4113 4872 5152 4804 4147 3738 1989 929 228 -1697 -4263 -5033 -5568 -6172 -6368 -5485 -5624 -7487 -7873 -8637 -6334 -2899 96 3072 4728 5975 6977 7108 6896 6330 4860 3331 1190 594 538 1249 1907 1596 1254 591 -277 -1414 -1209 290 3162 5175 7929 7229 4850 484 -3211 -4075 -2426 1161 3543 4624 2864 -689 -6750 -13974 -19877 -21413 -20168 -14964 -10261 -5297 -1041 3401 7980 11717 15254 15702 14673 9425 2890 -5709 -12050 -13783 -11367 -6476 -610 3042 4168 1417 -2669 -7034 -8368 -5177 -904 3588 7950 12647 16166 19254 19618 16456 10699 2431 -4668 -8822 -8182 -4943 1133 5414 7531 6005 2703 -3091 -8759 -13784 -16853 -17127 -14807 -11783 -9163 -6822 -4901 -935 3004 7074 9354 10706 8639 5634 1347 -2179