File name: Pandora239s1_Kosetice_20230413_L2Fit_fout2c2d20220901p1-8.txt File generation date: 20260325T155141.9Z 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.88 Instrument operation file used: Pandora239_OF_v3d20220426.txt Instrument calibration file used: Pandora239s1_CF_v2d20220901.txt Level 1 file used: Pandora239s1_Kosetice_20230413_L1_smca1c2d20220901p1-8.txt.bz2 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: 20230413 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 20230413T141300.4Z 8503.592365 365 1 31.14 2 56.83 241.55 136.19 283.26 0.00 1 0.00 1 0 25 7.638e-02 5.152e-02 2.297e-03 1.738e-03 2.516373e-01 8.840231e-04 3.200672e-03 -6.000000e+00 4.315798e-02 227.43 -2.00 3.74 -2.00 314.63 0.00e+00 9.705742e-04 6.733963e-05 2.219059e-04 -6.000000e+00 1.829368e-03 271.39 -2.00 20.00 -2.00 314.21 0.00e+00 1.822 0.000 302.00 328.00 3 4.2485e+00 2.2484e-03 2.1666e-02 -6.0000e+00 1.0754e-01 3.4422e-01 2.5039e-03 1.3215e-02 -6.0000e+00 1.1727e-01 -3.8412e-01 1.4318e-03 6.4311e-03 -6.0000e+00 6.9420e-02 1.6832e-01 6.1403e-04 1.7662e-03 -6.0000e+00 2.8975e-02 1.0333e-03 -1.5341e-03 4.8089e-05 2.4664e-05 -6.0000e+00 2.0407e-03 -4.4978e-02 -1.5929e-03 12 9 8 11 4 0 98.46 22.22 0.48 -0.04498 -0.05651 87 87 1 8 0 357.900 20.01 22.22 19.02 25.79 25.28 13.5 952.77 11533 14898 12767 3692 -12524 -21647 -20508 -14637 -8616 -8242 -12878 -15480 -8557 -716 6896 11801 11673 11316 9911 10994 8933 5432 3733 4749 5583 5842 2935 -5185 -13295 -15978 -13479 -10857 -6906 -8538 -16171 -20003 -15365 -6907 3344 10119 13018 12731 13264 12301 10485 8810 5754 2340 -1038 -1593 158 1674 3675 4269 3943 3855 3215 1135 -206 -1123 -77 2553 5677 6461 3953 -1575 -6254 -8760 -5006 -371 1787 2498 1652 182 -3718 -8929 -13448 -16681 -14864 -9870 -6847 -3939 -3570 -2058 778 2545 3947 5373 5088 3477 756 -3068 -5695 -7823 -6144 -3472 -369 1935 2390 435 -1029 -3473 -4089 -1927 -115 829 2389 3197 4298 5375 5834 5331 3523 1088 -1291 -2593 -2637 -1399 998 2135 2831 3089 2731 1120 -61 -1698 -3184 -2911 -2112 -1754 -1465 -1515 -2060 -1130 -260 506 966 1183 995 578 -423 -1231 -2065 5171 6814 5893 1685 -5688 -9834 -9341 -6729 -4002 -3886 -6043 -7324 -4101 -352 3480 6072 6081 5947 5333 5880 4825 2915 2020 2612 3103 3272 1654 -2901 -7295 -8684 -7419 -6121 -3972 -5039 -9370 -11505 -9059 -4255 2188 7110 9522 8694 9135 9344 7954 6691 4350 1757 -778 -1187 120 1331 2711 3260 3074 3023 2552 905 -161 -874 -61 2093 4858 5789 3563 -1365 -5239 -7290 -4246 -339 1733 2487 1643 175 -3434 -7737 -10982 -14684 -12248 -8701 -6540 -4015 -3839 -2273 820 2815 4512 6223 5901 3909 875 -3352 -5842 -7963 -6535 -3970 -416 2362 3044 548 -1231 -3919 -4798 -2312 -135 1013 3005 4177 5833 7679 8039 7404 4757 1509 -1679 -3222 -3292 -1808 1354 3043 4063 4350 3760 1492 -78 -2092 -3786 -3506 -2700 -2375 -2072 -2193 -3042 -1586 -384 774 1502 1873 1575 886 -631 -1778 -2944 SS 20230413T141857.2Z 8503.596496 368 1 31.65 2 57.68 242.97 137.07 284.59 0.00 1 0.00 1 0 25 7.759e-02 5.243e-02 2.649e-03 1.982e-03 2.499430e-01 1.003427e-03 3.146406e-03 -6.000000e+00 4.386718e-02 227.43 -2.00 3.74 -2.00 314.71 0.00e+00 1.463856e-03 7.751155e-05 1.932745e-04 -6.000000e+00 1.858693e-03 271.39 -2.00 20.00 -2.00 314.27 0.00e+00 1.864 0.000 302.00 328.00 3 4.4640e+00 2.5534e-03 2.1894e-02 -6.0000e+00 1.0931e-01 3.3997e-01 2.8506e-03 1.3170e-02 -6.0000e+00 1.1917e-01 -3.9000e-01 1.6285e-03 6.3261e-03 -6.0000e+00 7.0545e-02 1.7298e-01 6.9596e-04 1.7891e-03 -6.0000e+00 2.9465e-02 8.0586e-04 -7.1329e-04 5.4869e-05 2.2789e-05 -6.0000e+00 2.0725e-03 -4.7784e-02 -1.5209e-03 12 9 8 11 4 0 91.70 22.22 0.48 -0.04778 -0.05314 72 72 1 8 0 439.550 20.00 22.22 19.39 25.79 25.28 13.6 953.46 11965 15887 13961 4041 -12049 -21010 -20003 -14443 -8837 -8395 -13733 -15933 -8528 -691 7290 12179 12107 11966 10196 11892 9436 5823 3733 4920 6032 4753 2832 -5662 -13009 -15827 -14009 -11003 -7010 -8755 -16463 -20427 -15496 -6964 3747 10259 13290 13132 13083 12033 10289 9196 6165 2436 -938 -1441 300 1932 3363 4575 4063 3462 3129 1002 -596 -1184 83 2789 5658 6994 4033 -1516 -5991 -8813 -4786 -555 1385 2517 1306 98 -3476 -9229 -13446 -16927 -14878 -10194 -7081 -4231 -3937 -2164 916 2703 4221 5263 5174 3737 642 -2863 -5689 -8168 -6121 -3686 -366 1965 2494 798 -947 -3476 -4286 -1804 22 1275 2475 3184 4559 5672 5922 5327 3747 1034 -1406 -3168 -2678 -1482 1027 2411 3259 3098 2952 1225 145 -1638 -3086 -2919 -2127 -1750 -1284 -1549 -1950 -1077 -293 328 909 1354 907 416 -371 -1180 -1906 4877 6571 5816 1846 -5465 -9501 -9054 -6596 -4082 -3923 -6423 -7485 -4046 -337 3640 6204 6243 6212 5438 6265 5037 3089 2001 2668 3286 2646 1570 -3132 -7015 -8479 -7662 -6143 -3991 -5117 -9450 -11653 -9032 -4234 2406 7082 9529 9654 9762 9018 7703 6858 4579 1802 -693 -1058 226 1511 2699 3432 3115 2681 2444 786 -461 -908 65 2250 4778 6150 3581 -1294 -4937 -7237 -3996 -501 1330 2471 1285 93 -3157 -7909 -11914 -14747 -13330 -8914 -6703 -4276 -4200 -2365 1040 2951 4757 6034 5931 4145 734 -3087 -5772 -8258 -6445 -4185 -410 2377 3144 994 -1122 -4239 -4996 -2146 25 1539 3088 4132 6129 8025 8096 7344 5427 1425 -1815 -3932 -3323 -1907 1386 3408 4632 4338 4032 1619 183 -2005 -3978 -3500 -2709 -2360 -1804 -2231 -2864 -1630 -431 500 1408 2131 1428 636 -550 -1836 -2927 SS 20230413T142634.4Z 8503.601787 371 1 30.52 2 58.79 244.76 138.19 286.34 0.00 1 0.00 1 0 25 7.847e-02 5.179e-02 3.017e-03 2.208e-03 2.579810e-01 1.117483e-03 3.272123e-03 -6.000000e+00 4.436968e-02 227.43 -2.00 3.74 -2.00 314.87 0.00e+00 1.500086e-03 8.869007e-05 2.181433e-04 -6.000000e+00 1.879756e-03 271.39 -2.00 20.00 -2.00 314.41 0.00e+00 1.922 0.000 302.00 328.00 3 4.5336e+00 2.8477e-03 2.2616e-02 -6.0000e+00 1.1056e-01 3.5414e-01 3.1805e-03 1.3622e-02 -6.0000e+00 1.2054e-01 -4.0491e-01 1.8148e-03 6.5690e-03 -6.0000e+00 7.1352e-02 1.7864e-01 7.7741e-04 1.8681e-03 -6.0000e+00 2.9800e-02 7.0102e-04 -4.8937e-04 6.1280e-05 2.3433e-05 -6.0000e+00 2.0958e-03 -4.5047e-02 -1.3966e-03 12 9 8 11 4 0 91.93 22.22 0.45 -0.04505 -0.04872 69 69 1 8 0 442.210 20.00 22.22 19.77 25.79 25.28 13.5 952.29 12356 15995 13719 5009 -12543 -20967 -21024 -15157 -9342 -8732 -13147 -16482 -9303 -415 7203 12533 12898 11967 10058 12142 10359 5308 4043 5537 4814 5521 2619 -5436 -13609 -16184 -13827 -10742 -6638 -8350 -15820 -20550 -15458 -7557 3781 10416 13305 13010 13843 12119 10220 8828 5744 2464 -374 -890 349 2319 3393 4009 3937 3723 3000 1648 -486 -1288 83 2479 5506 6948 3730 -1793 -6261 -8821 -5009 -344 1696 2257 4091 -21 -3905 -9123 -13616 -17239 -14651 -10574 -7245 -4125 -3922 -2105 820 2321 3996 5617 5026 3604 212 -3274 -5955 -8292 -6056 -3706 -642 1539 1976 555 -1297 -3871 -4023 -1636 262 1235 2279 3364 4613 5508 5937 5245 3287 904 -1505 -2714 -2768 -1495 1019 2585 2711 3227 2880 1384 4 -1621 -2924 -2653 -2008 -1319 -1155 -1385 -1967 -1309 -299 271 770 1215 1109 554 -533 -1123 -1868 4860 6408 5557 2200 -5546 -9180 -9301 -6738 -4189 -3955 -5907 -7526 -4304 -196 3492 6179 6414 6041 5225 6203 5323 2748 2097 2895 2575 2958 1411 -2909 -7161 -8446 -7338 -5823 -3668 -4739 -8797 -11436 -8769 -4499 2362 6988 9287 9325 9982 8835 7453 6432 4171 1776 -268 -635 256 1765 2661 2958 2956 2813 2295 1256 -403 -1062 70 1959 4550 5977 3252 -1504 -5072 -7110 -4117 -305 1595 2181 3839 -19 -3490 -7652 -11848 -14788 -12874 -9132 -6764 -4102 -4122 -2265 918 2734 4446 6319 5678 4297 240 -3490 -5970 -8283 -6299 -4162 -778 1852 2481 686 -1526 -4691 -4627 -1921 329 1475 2818 4311 6137 7738 8722 7782 4749 1238 -1934 -3335 -3422 -1915 1368 3628 3848 4480 3906 1814 4 -2150 -3739 -3444 -2544 -1765 -1617 -1989 -2885 -1984 -476 413 1193 1912 1740 842 -855 -1740 -2857