File name: Pandora239s1_Kosetice_20250425_L2Fit_fout2c2d20220901p1-8.txt File generation date: 20250425T230312.0Z 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.75 Instrument operation file used: Pandora239_OF_v3d20220426.txt Instrument calibration file used: Pandora239s1_CF_v2d20220901.txt Level 1 file used: Pandora239s1_Kosetice_20250425_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: 20250425 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 20250425T100840.4Z 9246.422690 193 1 24.01 2 37.60 160.24 52.02 205.37 0.00 1 0.00 1 0 17 2.575e-01 1.025e-01 1.496e-01 4.903e-02 2.685027e-01 2.492924e-02 3.758058e-03 -6.000000e+00 1.463015e-01 228.62 -2.00 3.74 -2.00 316.63 0.00e+00 8.556711e-03 2.930352e-03 1.229203e-03 -6.000000e+00 6.204657e-03 272.85 -2.00 20.00 -2.00 316.12 0.00e+00 1.261 0.000 302.00 328.00 3 9.0702e+00 6.4933e-02 1.9894e-02 -6.0000e+00 3.6540e-01 4.1752e-01 7.2628e-02 1.9173e-02 -6.0000e+00 3.9793e-01 -1.7004e-01 4.2448e-02 7.5385e-03 -6.0000e+00 2.3556e-01 4.3365e-02 1.9602e-02 7.0680e-03 -6.0000e+00 9.8164e-02 1.1084e-05 -4.7654e-03 1.4071e-03 8.4861e-05 -6.0000e+00 6.8770e-03 -2.6338e-02 -7.2995e-04 12 9 8 11 4 0 67.99 22.06 -0.17 -0.02634 -0.06215 6 6 2 8 0 4000.000 19.97 22.06 15.16 25.58 25.07 23.1 952.13 -46863 45478 -24620 210185 -2069 3204 39869 24180 41805 3035 32120 -32963 -20261 74104 18772 14157 10531 1363 35526 4807 25621 -11966 -19826 47495 50808 47236 2386 17000 -12211 5911 -4249 18927 15111 55848 18395 -3145 -28895 10034 2593 -18795 633 4826 43027 -20366 16125 18046 11934 -10211 10854 -8605 -7209 -11496 14313 -4829 -2377 7139 -377 -7638 8766 -11194 8220 8330 -1712 38403 8147 6815 -14836 -26645 2197 16210 -132 2101 -854 2743 6363 -9825 17909 -17279 -17303 -3254 3001 17320 -2656 -7543 14219 -6470 -3736 15120 9109 4608 20376 -11069 -2953 -2654 2069 -5211 -5585 -3958 4021 8895 11184 -7288 7066 -9367 5027 2232 1162 550 4005 -2050 7946 -18955 -2479 5478 21487 -2349 -2370 -8556 22268 -2714 -13959 18362 12550 400 -12835 1119 -2528 -397 -1234 17890 -7717 5401 -6847 7666 -1 2849 7835 -672 4313 -3624 -4793 -2846 567 -13788 6523 -6882 10049 -420 549 4607 3372 5163 568 4351 -8405 -5086 8562 4413 3948 3231 472 9059 1705 7452 -4910 -8531 10408 10729 10351 817 4809 -4330 1875 -1588 6088 5452 14079 6316 -1244 -15448 4056 1308 -13561 410 3105 19450 -17704 9742 10746 7470 -7767 6511 -6172 -5264 -9294 9331 -3967 -1943 5275 -300 -6440 6096 -9390 5823 6300 -1555 25038 7260 5688 -14057 -27100 1719 12330 -129 2100 -876 2599 5393 -8754 10820 -13894 -14290 -2616 2544 14002 -2793 -8622 13843 -8268 -4872 16824 10808 5508 19130 -12894 -2875 -2446 2088 -5662 -6773 -5247 5244 10864 12493 -8964 7669 -12137 6105 2935 1606 727 5344 -3082 11397 -36004 -3836 7989 24348 -3144 -3138 -12517 25552 -4256 -22691 23352 16441 564 -19453 1402 -3106 -476 -1592 20357 -11917 7514 -10901 10886 -1 4699 12549 -1196 7320 -6446 -8272 -4627 865 SS 20250425T111945.4Z 9246.472053 224 1 24.01 2 36.50 189.07 58.20 225.20 0.00 1 0.00 1 0 145 4.639e-01 4.756e-01 1.686e-01 3.836e-02 1.213371e+00 2.092818e-02 2.128639e-02 -6.000000e+00 2.651201e-01 228.63 -2.00 3.74 -2.00 317.31 0.00e+00 -1.733271e-02 2.654271e-03 8.744400e-03 -6.000000e+00 1.121236e-02 272.85 -2.00 20.00 -2.00 316.77 0.00e+00 1.243 0.000 302.00 328.00 3 6.6521e+00 5.4624e-02 7.9425e-02 -6.0000e+00 6.6222e-01 3.2063e+00 6.0305e-02 8.9355e-02 -6.0000e+00 7.2159e-01 -1.6575e+00 3.6296e-02 3.4756e-02 -6.0000e+00 4.2570e-01 -4.3078e-02 1.8293e-02 2.1588e-02 -6.0000e+00 1.7664e-01 1.3307e-05 -1.2802e-02 7.7337e-04 3.3846e-04 -6.0000e+00 1.2493e-02 -2.4096e-02 -8.1173e-04 12 9 8 11 4 0 26.95 21.99 9.80 -0.02410 -0.12029 6 6 2 8 0 4000.000 19.87 21.99 14.77 25.47 25.07 23.0 952.29 -99380 32819 -52234 329584 -71970 -98053 900000 121678 -85000 6642 6133 -6466 -27030 22240 -23447 22657 -12945 -16290 807 16247 18666 -39081 -49087 -10139 57915 20628 37234 -1524 11124 17324 -10409 57425 56610 74047 47139 24026 -7491 4185 26512 2893 22591 44546 41244 2512 19609 -3837 -8301 -3041 -3162 -11135 -580 7014 20074 13480 22202 33970 22645 15602 13812 -4420 -16478 -24681 -31808 -5049 -11533 2489 1224 -12019 20278 27370 19767 4431 7369 9008 -14466 -26110 -1324 -28733 -23241 -18750 2521 17509 10610 -1775 11204 -5308 -8976 -6513 -20295 -24691 -11808 -42322 11368 -25911 -10470 -7958 606 5311 25407 21065 18909 5741 3543 -24232 -34981 -22842 -29957 -29896 -17095 -8193 -1152 10798 18119 31950 42983 32674 30188 18557 48253 24841 6577 28546 12027 -2324 -6784 9571 7923 8857 14085 48340 41392 67132 63766 73964 71895 72246 -71155 78074 79200 76192 74561 68360 73074 -18895 1944 -7652 12327 -10767 -16503 900000 4551 -17220 621 631 -806 -4802 2911 -5773 4166 -3833 -5331 231 3950 4247 -14739 -19736 -3065 11225 6312 9748 -539 2950 3947 -3468 11715 13474 16365 13303 7597 -3133 1693 10187 1722 12828 21659 20615 1743 10335 -2322 -4939 -1675 -1785 -7155 -365 4555 12422 9534 15325 21599 16400 11777 9580 -3171 -11710 -18295 -27038 -3691 -9882 1862 907 -10462 13514 18798 15996 4264 6564 7014 -12051 -20153 -753 -18855 -15568 -14026 1791 12324 8947 -1774 9981 -5649 -9858 -6821 -23256 -27368 -11350 -45329 7483 -25488 -9869 -8543 694 6067 29173 25842 22136 6700 3525 -28005 -45007 -26826 -38543 -40589 -22442 -11308 -1571 14723 26768 41466 48551 38984 36018 25264 50172 32982 10036 32030 14709 -3129 -9159 10631 8683 9905 17543 51709 55579 79092 85621 97930 106216 114300 -374840 122061 123229 117408 107924 96902 92926 SS 20250425T122147.1Z 9246.515129 251 1 24.01 2 40.10 212.84 65.72 240.08 0.00 1 0.00 1 0 33 3.937e-01 1.411e-01 2.383e-01 6.321e-02 2.383780e-01 3.245098e-02 3.463918e-03 -6.000000e+00 2.242295e-01 228.63 -2.00 3.74 -2.00 316.52 0.00e+00 7.634660e-03 3.804318e-03 5.315476e-04 -6.000000e+00 9.495011e-03 272.86 -2.00 20.00 -2.00 316.02 0.00e+00 1.306 0.000 302.00 328.00 3 9.3550e+00 8.4113e-02 1.9581e-02 -6.0000e+00 5.6009e-01 4.0211e-01 9.3534e-02 1.7873e-02 -6.0000e+00 6.0939e-01 -1.1574e-01 5.5011e-02 7.2662e-03 -6.0000e+00 3.6078e-01 -2.0019e-02 2.6176e-02 5.9305e-03 -6.0000e+00 1.5071e-01 8.4642e-06 -1.9357e-03 1.8578e-03 6.3756e-05 -6.0000e+00 1.0506e-02 -3.1643e-02 -1.9503e-04 12 9 8 11 4 0 94.97 21.98 6.53 -0.03164 -0.04619 6 6 2 8 0 4000.000 20.03 21.98 14.61 25.39 25.07 23.0 952.42 18544 44816 188643 103242 188510 -24549 215208 86982 -14122 -14617 82418 -25755 26824 118622 15475 29943 -30213 19902 -22712 68306 574 -559 -44270 52713 70749 36511 34859 8280 -16836 19300 -1886 6 28561 63012 17927 -11608 -37922 25800 29046 -3308 3170 21148 39619 -5764 32099 21841 7725 10578 13686 14423 -17908 -6774 16532 -6629 -1347 7205 8638 -8188 20648 -48983 8952 -3246 5389 23414 24947 2860 -3421 -35941 18107 19025 12656 -2155 15415 21908 -2523 -16463 9670 -15588 -13722 2089 -9258 17112 -6594 -2188 11553 -1904 690 15050 12142 5148 16806 -8423 -6089 -4939 10739 -5688 566 -10044 10840 12926 13371 -13187 8941 -10416 -1377 9259 -7615 -2834 6878 1063 1890 -18126 -615 19017 17869 10164 -2346 -3042 36500 -295 -16443 12183 10514 -11727 -14522 3109 16772 -1055 5259 15210 -5236 16635 -13155 6978 -1864 2670 -16945 8549 -10287 7193 7813 12895 2307 2814 5951 11629 7833 11621 -5239 13267 6703 -3192 -3454 7198 -6483 4410 9141 3895 7384 -14435 5930 -10735 13335 224 -215 -25700 11645 13085 9566 9312 2765 -6747 5255 -668 2 8657 15538 6530 -4830 -21430 9398 11883 -2169 2132 12326 19761 -4620 17629 13342 5368 6936 8437 8649 -15255 -5466 11026 -5796 -1140 5573 6584 -7296 13423 -63139 6623 -2892 4772 18502 19573 2580 -2991 -41293 12381 14394 11207 -2307 13823 17680 -2418 -16239 6557 -12667 -11151 1616 -8976 14002 -7280 -2392 11655 -2347 870 16954 14154 6203 16576 -9677 -6176 -4677 9027 -6184 641 -14016 13062 14996 14463 -17072 8619 -12375 -1621 10326 -10456 -3875 8904 1542 2859 -33828 -924 23916 20709 11843 -3065 -4163 36002 -449 -27370 16523 14131 -18797 -22561 3853 17163 -1288 6419 17910 -7941 20817 -22464 10066 -3114 4478 -35412 14178 -20737 11840 12331 18682 3626