File name: Pandora239s1_Kosetice_20240108_L2Fit_fout2c2d20220901p1-8.txt
File generation date: 20240113T035943.1Z
Data description: Level 2 spectral fitting results file
Data file version: fout2c2p1-8
Local principal investigator: Nadezda Zikova
Network principal investigator: Alexander Cede
Instrument type: Pandora
Instrument number: 239
Spectrometer number: 1
Processing software version used: BlickP v1.8.59
Instrument operation file used: Pandora239_OF_v3d20220426.txt
Instrument calibration file used: Pandora239s1_CF_v2d20220901.txt
Level 1 file used: Pandora239s1_Kosetice_20240108_L1_smca1c2d20220901p1-8.txt
Full location name: National Atmospheric Observatory Kosetice
Short location name: Kosetice
Country of location: Czech Republic
Location latitude [deg]: 49.5733
Location longitude [deg]: 15.0803
Location altitude [m]: 534
Local noon date: 20240108
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: Initial instrument OoB session based on session 1
First and last pixel inside fitting window: 197 339
Nominal wavelengths inside fitting window [nm]: 305.0356 305.1770 305.3183 305.4597 305.6010 305.7424 305.8837 306.0250 306.1663 306.3076 306.4488 306.5901 306.7314 306.8726 307.0138 307.1551 307.2963 307.4374 307.5786 307.7198 307.8610 308.0021 308.1432 308.2844 308.4255 308.5666 308.7077 308.8487 308.9898 309.1309 309.2719 309.4129 309.5540 309.6950 309.8360 309.9769 310.1179 310.2589 310.3998 310.5408 310.6817 310.8226 310.9635 311.1044 311.2453 311.3861 311.5270 311.6678 311.8087 311.9495 312.0903 312.2311 312.3719 312.5126 312.6534 312.7942 312.9349 313.0756 313.2163 313.3570 313.4977 313.6384 313.7791 313.9197 314.0604 314.2010 314.3416 314.4822 314.6228 314.7634 314.9040 315.0445 315.1851 315.3256 315.4661 315.6066 315.7471 315.8876 316.0281 316.1686 316.3090 316.4495 316.5899 316.7303 316.8707 317.0111 317.1515 317.2918 317.4322 317.5725 317.7129 317.8532 317.9935 318.1338 318.2741 318.4144 318.5546 318.6949 318.8351 318.9753 319.1155 319.2557 319.3959 319.5361 319.6763 319.8164 319.9566 320.0967 320.2368 320.3769 320.5170 320.6571 320.7971 320.9372 321.0772 321.2173 321.3573 321.4973 321.6373 321.7773 321.9172 322.0572 322.1971 322.3371 322.4770 322.6169 322.7568 322.8967 323.0366 323.1764 323.3163 323.4561 323.5959 323.7358 323.8756 324.0153 324.1551 324.2949 324.4346 324.5744 324.7141 324.8538 324.9935
---------------------------------------------------------------------------------------
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: Spectrometer control temperature [°C], 999=no temperature signal
Column 95: Auxiliary spectrometer temperature [°C], 999=no temperature signal
Column 96: Temperature in head sensor [°C], 999=no temperature signal
Column 97: Temperature at electronics board 1 [°C], 999=no temperature signal
Column 98: Temperature at electronics board 2 [°C], 999=no temperature signal
Column 99: Humidity in head sensor [%], -9=no humidity signal
Column 100: Pressure in head sensor [hPa], -9=no pressure signal
Columns 101-243: Unweighted slant column residuals for each pixel inside the fitting window multiplied by 1e5, 9e5=pixel was not used for fitting
Columns 244-386: 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 20240108T115312.4Z 8773.495283 218 1 26.68 2 72.60 191.40 90.42 227.84 0.00 1 0.00 1 0 49 1.088e-01 8.623e-02 4.289e-03 3.599e-03 3.026240e-01 1.851649e-03 4.096100e-03 -6.000000e+00 6.110366e-02 221.74 -2.00 3.74 -2.00 313.59 0.00e+00 -1.650069e-03 1.179514e-04 5.213529e-04 -6.000000e+00 2.560292e-03 265.88 -2.00 20.00 -2.00 313.33 0.00e+00 3.298 0.001 302.00 328.00 3 3.8397e+00 4.6126e-03 3.4020e-02 -6.0000e+00 1.5070e-01 7.9559e-01 5.1133e-03 1.5338e-02 -6.0000e+00 1.6424e-01 -6.4161e-01 2.9591e-03 7.8708e-03 -6.0000e+00 9.7267e-02 2.0423e-01 1.2832e-03 3.3066e-03 -6.0000e+00 4.2022e-02 3.8247e-04 7.8914e-05 9.7553e-05 5.0849e-05 -6.0000e+00 2.9091e-03 -3.4189e-02 -1.7760e-03 12 9 8 11 4 0 84.40 22.15 0.50 -0.03419 -0.03360 55 55 1 8 0 484.950 20.00 22.15 -0.16 25.68 25.17 15.8 956.61 21825 23654 18742 5317 -13418 -20942 -18891 -11460 -5061 -6468 -11744 -13825 -6063 2314 11254 14808 13112 12629 11079 11789 7441 2206 1829 4913 6741 5452 412 -10563 -20546 -23957 -21451 -16934 -11953 -14874 -25460 -28780 -22174 -8863 6518 15924 20863 20342 20609 16972 13374 10562 5018 479 -2665 -3312 -47 2747 5881 8025 7904 7852 7755 4585 75 -2838 -53 5338 9759 12287 6901 -2205 -8587 -11957 -4496 3175 7272 7418 6179 1105 -6756 -17134 -22978 -27759 -22173 -14101 -8281 -2413 -1143 1654 5852 7912 10680 12272 9856 6880 -165 -6570 -10176 -13128 -9296 -4561 934 5366 4781 414 -3471 -8108 -8024 -3316 841 2571 5525 6318 8776 10266 10743 8604 4081 -659 -5250 -7122 -5623 -3096 1322 2974 2946 3654 2633 -1069 -3089 -5552 -8259 -6227 -3894 -2963 -3559 -3168 -3210 -2073 743 1450 2407 2991 2949 1422 1136 545 900000 13025 14117 11209 3151 -8034 -12585 -11455 -6934 -3050 -3937 -7073 -8454 -3704 1417 6849 9074 8114 7815 6916 7224 4600 1366 1139 3062 4195 3436 260 -7264 -13972 -16349 -14835 -11711 -8279 -10349 -17643 -20065 -15662 -6263 4638 11650 15271 15132 15251 12661 10000 7872 3767 359 -1992 -2488 -36 2138 4600 6364 6340 6325 6239 3691 61 -2302 -43 4340 8136 10416 5925 -1877 -7279 -10317 -3887 2848 6678 6924 5682 1005 -6019 -14970 -19302 -23669 -19201 -12696 -7828 -2350 -1155 1696 6150 8618 11749 13508 10949 7402 -173 -6720 -9987 -13230 -9740 -5051 1089 5988 5530 480 -3876 -9542 -9040 -3774 992 3109 6797 7352 10549 13017 14208 11588 5474 -857 -6578 -8690 -6862 -3881 1703 3949 3917 4721 3336 -1334 -3739 -7208 -10597 -8135 -4885 -3931 -4984 -4526 -4695 -3146 1166 2345 3935 4958 4837 2275 1752 814 900000
SS 20240108T120025.9Z 8773.500300 222 1 26.88 2 72.85 193.13 91.28 229.15 0.00 1 0.00 1 0 49 1.121e-01 9.365e-02 1.989e-03 1.744e-03 2.835574e-01 8.963040e-04 3.732731e-03 -6.000000e+00 6.291549e-02 221.74 -2.00 3.74 -2.00 313.14 0.00e+00 -1.205007e-03 5.372245e-05 4.270028e-04 -6.000000e+00 2.637702e-03 265.88 -2.00 20.00 -2.00 312.92 0.00e+00 3.343 0.001 302.00 328.00 3 2.3797e+00 2.2260e-03 3.3870e-02 -6.0000e+00 1.5517e-01 8.2947e-01 2.4665e-03 1.4171e-02 -6.0000e+00 1.6915e-01 -6.0908e-01 1.4328e-03 7.2330e-03 -6.0000e+00 1.0017e-01 1.7853e-01 6.1768e-04 2.8918e-03 -6.0000e+00 4.3258e-02 1.6285e-03 -3.4806e-04 4.7015e-05 4.7355e-05 -6.0000e+00 2.9983e-03 -3.4697e-02 -1.9593e-03 12 9 8 11 4 0 89.29 22.16 0.54 -0.03470 -0.03731 290 290 1 8 0 92.680 20.00 22.16 -0.22 25.68 25.17 15.8 958.51 21769 23742 18490 4187 -13337 -20731 -17781 -10867 -4991 -5712 -12354 -13231 -5655 2964 10515 14804 13126 12816 11378 11137 7345 2637 2254 5174 6017 6029 694 -11174 -21989 -24719 -20945 -17445 -12668 -16396 -27007 -29474 -21790 -9284 6248 16140 20325 20839 20303 17056 13495 10636 5465 623 -3258 -3877 -269 3183 5935 8030 7593 7662 7364 4140 163 -2122 546 5151 10235 12746 7225 -2782 -8957 -11475 -4862 3631 7669 8147 5992 1108 -6798 -16757 -25131 -28491 -23450 -14940 -8235 -2522 -1186 1980 6906 10231 12116 13602 11958 7252 397 -6678 -12707 -13970 -9598 -4291 1233 5729 5637 1753 -2961 -8142 -8975 -3507 1119 3010 5517 7386 9762 11733 11798 9416 4785 -445 -5675 -8132 -6708 -3623 698 3487 3732 3998 2385 -1111 -4174 -6896 -9474 -7426 -5233 -3984 -3110 -3780 -3781 -1802 711 1776 3204 4208 3531 2487 1710 815 900000 14155 15397 12046 2723 -8680 -13576 -11655 -7140 -3275 -3748 -8117 -8738 -3740 1960 6998 9821 8778 8555 7655 7433 4919 1761 1514 3473 4061 4063 468 -7551 -14845 -18319 -15498 -13005 -9457 -12343 -20244 -22114 -16407 -7049 4761 12540 15872 16350 15975 13481 10699 8407 4329 494 -2585 -3086 -217 2588 4858 6652 6365 6432 6177 3474 137 -1780 456 4372 8830 11138 6374 -2449 -7812 -10082 -4316 3311 7150 7686 5618 1025 -6164 -14829 -21897 -24796 -20808 -13688 -7831 -2468 -1200 2025 7183 10884 13126 14751 12976 7744 414 -6788 -12662 -14023 -9941 -4659 1407 6857 6310 2148 -3525 -9345 -9975 -3913 1286 3541 6624 9111 11290 14204 14911 12109 6119 -553 -6831 -10544 -8707 -4368 866 4414 4722 4939 2907 -1464 -5384 -8689 -11776 -9382 -6976 -5077 -4114 -5148 -5263 -2581 1056 2709 4920 6528 5459 3738 2493 1155 900000
SS 20240108T124838.1Z 8773.533774 241 1 26.80 2 75.36 204.42 97.39 237.67 0.00 1 0.00 1 0 81 1.176e-01 1.018e-01 3.950e-03 3.549e-03 2.611587e-01 1.785715e-03 3.382829e-03 -6.000000e+00 6.501577e-02 221.74 -2.00 3.74 -2.00 312.84 0.00e+00 -1.137173e-03 1.035035e-04 4.218291e-04 -6.000000e+00 2.754465e-03 265.87 -2.00 20.00 -2.00 312.60 0.00e+00 3.878 0.001 302.00 328.00 3 3.1962e+00 4.4312e-03 3.7862e-02 -6.0000e+00 1.6024e-01 8.9952e-01 4.9769e-03 1.2909e-02 -6.0000e+00 1.7566e-01 -5.7971e-01 2.8764e-03 6.5316e-03 -6.0000e+00 1.0410e-01 1.5437e-01 1.1780e-03 2.6813e-03 -6.0000e+00 4.3748e-02 4.6371e-04 3.0609e-03 9.5435e-05 4.6738e-05 -6.0000e+00 3.1388e-03 -4.2123e-02 -1.6679e-03 12 9 8 11 4 0 95.50 22.16 0.54 -0.04212 -0.01912 96 96 1 8 0 279.060 20.01 22.16 -0.56 25.68 25.17 15.9 962.60 23588 24079 18957 3773 -13936 -20259 -16336 -10293 -4077 -5298 -11765 -11716 -5894 3479 10523 14328 12637 11964 10755 10663 6624 2092 2253 5348 5967 5159 909 -12457 -21749 -25558 -22176 -18296 -14069 -17881 -28503 -29873 -22006 -9081 6404 16323 20215 20455 20725 16747 12734 9721 4526 -384 -4967 -4127 -1021 2345 5276 7768 7180 7300 7439 3828 -201 -2671 1079 5500 10831 13583 6875 -2726 -10535 -13532 -4684 4174 8801 9271 5595 1168 -7953 -18819 -27381 -30601 -25023 -15464 -8741 -3000 -325 2412 9057 12032 14442 15605 13925 8648 874 -7884 -13230 -15746 -10496 -4286 2045 6285 6450 1521 -3822 -9390 -9758 -3404 1188 4148 6706 9240 11269 14043 13881 10420 4945 -967 -7072 -9718 -8498 -5111 736 3377 3982 3833 2697 -2173 -5520 -8075 -11713 -8886 -6507 -5387 -4758 -4591 -5338 -2612 193 1625 14084 4238 3184 2005 1298 720 470 14986 15437 13264 2652 -9852 -14334 -11496 -7335 -2888 -3758 -8328 -8267 -4229 2469 7528 10214 9057 8582 7755 7621 4762 1499 1619 3826 4297 3735 650 -9037 -15539 -18429 -15969 -13176 -10167 -12950 -20553 -21465 -17485 -7256 5130 13269 16495 16791 16905 13811 10585 8057 3751 -318 -4134 -3395 -852 1973 4453 6591 6160 6262 6338 3272 -172 -2297 917 4755 9460 11911 6107 -2397 -9323 -12133 -4163 3787 8090 8601 5230 1072 -7225 -16822 -24144 -26951 -22385 -14096 -8240 -2897 -319 2406 9037 12304 14944 16190 14423 8835 878 -7836 -12792 -15563 -10610 -4491 2228 7171 7498 1770 -4345 -10353 -10452 -3648 1309 4620 7616 10689 13470 17369 17860 13703 6489 -1238 -8837 -11955 -10519 -6460 943 4419 5182 4895 3377 -2706 -6751 -9616 -13943 -10702 -8228 -7148 -6545 -6428 -7698 -3501 268 2300 19789 6093 4572 2805 1941 1045 680