File name: Pandora238s1_Granada_20231201_L2Fit_fout2c14d20230901p1-8.txt File generation date: 20250305T233159.5Z Data description: Level 2 spectral fitting results file Data file version: fout2c14p1-8 Local principal investigator: Inmaculada Foyo Moreno Network principal investigator: Alexander Cede Instrument type: Pandora Instrument number: 238 Spectrometer number: 1 Processing software version used: BlickP v1.8.74 Instrument operation file used: Pandora238_OF_v2d20220614.txt Instrument calibration file used: Pandora238s1_CF_v14d20230901.txt Level 1 file used: Pandora238s1_Granada_20231201_L1_smca1c14d20230901p1-8.txt Full location name: Andalusian Institute for Earth System Research (IISTA) - University of Granada Short location name: Granada Country of location: Spain Location latitude [deg]: 37.1640 Location longitude [deg]: -3.6050 Location altitude [m]: 680 Local noon date: 20231201 Notes on s-number (L1 configuration): Corrections NOT applied although requested by the s-number are blind pixel subtraction, dark map for dark correction (replaced by immediate dark method), latency correction, matrix method stray light correction (replaced by simple method) Data caveats: Update session 10 for SO2 First and last pixel inside fitting window: 202 343 Nominal wavelengths inside fitting window [nm]: 305.0915 305.2330 305.3745 305.5159 305.6574 305.7988 305.9402 306.0815 306.2229 306.3642 306.5055 306.6468 306.7881 306.9293 307.0705 307.2117 307.3529 307.4940 307.6352 307.7763 307.9174 308.0585 308.1995 308.3406 308.4816 308.6226 308.7636 308.9045 309.0455 309.1864 309.3273 309.4682 309.6090 309.7499 309.8907 310.0315 310.1723 310.3131 310.4538 310.5946 310.7353 310.8760 311.0167 311.1573 311.2980 311.4386 311.5792 311.7198 311.8604 312.0009 312.1415 312.2820 312.4225 312.5630 312.7035 312.8439 312.9844 313.1248 313.2652 313.4056 313.5460 313.6863 313.8267 313.9670 314.1073 314.2476 314.3879 314.5281 314.6684 314.8086 314.9488 315.0890 315.2292 315.3694 315.5095 315.6497 315.7898 315.9299 316.0700 316.2101 316.3501 316.4902 316.6302 316.7702 316.9102 317.0502 317.1902 317.3301 317.4701 317.6100 317.7499 317.8898 318.0297 318.1696 318.3094 318.4493 318.5891 318.7289 318.8687 319.0085 319.1483 319.2880 319.4278 319.5675 319.7072 319.8469 319.9866 320.1263 320.2659 320.4056 320.5452 320.6848 320.8244 320.9640 321.1036 321.2432 321.3827 321.5223 321.6618 321.8013 321.9408 322.0803 322.2198 322.3592 322.4987 322.6381 322.7775 322.9169 323.0563 323.1957 323.3351 323.4744 323.6138 323.7531 323.8924 324.0317 324.1710 324.3103 324.4495 324.5888 324.7280 324.8673 --------------------------------------------------------------------------------------- Column 1: Two letter code of measurement routine Column 2: UT date and time for center-time of measurement, yyyymmddThhmmssZ (ISO 8601) Column 3: Fractional days since 1-Jan-2000 midnight for center-time of measurement Column 4: Routine count (1 for the first routine of the day, 2 for the second, etc.) Column 5: Repetition count (1 for the first set in the routine, 2 for the second, etc.) Column 6: Total duration of measurement set in seconds Column 7: Data processing type index Column 8: Solar zenith angle for center-time of measurement in degree Column 9: Solar azimuth for center-time of measurement in degree, 0=north, increases clockwise Column 10: Lunar zenith angle for center-time of measurement in degree Column 11: Lunar azimuth for center-time of measurement in degree, 0=north, increases clockwise Column 12: Pointing zenith angle in degree, absolute or relative (see next column), 999=tracker not used Column 13: Zenith pointing mode: zenith angle is... 0=absolute, 1=relative to sun, 2=relative to moon Column 14: Pointing azimuth in degree, increases clockwise, absolute (0=north) or relative (see next column), 999=tracker not used Column 15: Azimuth pointing mode: like zenith angle mode but also fixed scattering angles relative to sun (3) or moon (4) Column 16: Fitting result index: 0=no error or warning, 1, 2=warning, >2=error Column 17: Number of function evaluations used, 0=linear fitting or fitting not successful or no fitting done Column 18: rms of unweighted fitting residuals, -9=fitting not successful Column 19: Normalized rms of fitting residuals weighted with independent uncertainty, -9=fitting not successful or no uncertainty used Column 20: Expected rms based on independent uncertainty, -9=fitting not successful or no uncertainty given Column 21: Expected normalized weighted rms based on independent uncertainty, -9=fitting not successful or no uncertainty given Column 22: Ozone slant column amount [moles per square meter], -9e99=fitting not successful Column 23: Independent uncertainty of ozone slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -3=spectral fitting was done, but no independent uncertainty could be retrieved, -5=no independent uncertainty input was given, -9=spectral fitting not successful Column 24: Structured uncertainty of ozone slant column amount [moles per square meter], -7=not given since method "MEAS" was chosen, -9=spectral fitting not successful Column 25: Common uncertainty of ozone slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -6=no common uncertainty input was given, -7=not given since method "MEAS" was chosen, -9=spectral fitting not successful Column 26: rms-based uncertainty of ozone slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -3=spectral fitting was done, but no rms-based uncertainty could be retrieved, -9=spectral fitting not successful Column 27: Ozone effective temperature [K] Column 28: Independent uncertainty of ozone effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -2=no temperature fitting was requested and output for effective temperature and structured uncertainty of it is based on f-code, -3=spectral fitting was done, but no independent uncertainty could be retrieved, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -5=no independent uncertainty input was given, -9=spectral fitting not successful Column 29: Structured uncertainty of ozone effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -9=spectral fitting not successful Column 30: Common uncertainty of ozone effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -2=no temperature fitting was requested and output for effective temperature and structured uncertainty of it is based on f-code, -3=spectral fitting was done, but no common uncertainty could be retrieved, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -6=no common uncertainty input was given, -9=spectral fitting not successful Column 31: Effective ozone fitting wavelength [nm], -9=fitting not successful Column 32: Diffuse correction applied before fitting at effective fitting wavelength for ozone [%], 0=no diffuse correction applied or fitting not requested, >0=measured diffuse correction, <0=(negative value of) calculated diffuse correction Column 33: Sulfur dioxide slant column amount [moles per square meter], -9e99=fitting not successful Column 34: Independent uncertainty of sulfur dioxide slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -3=spectral fitting was done, but no independent uncertainty could be retrieved, -5=no independent uncertainty input was given, -9=spectral fitting not successful Column 35: Structured uncertainty of sulfur dioxide slant column amount [moles per square meter], -7=not given since method "MEAS" was chosen, -9=spectral fitting not successful Column 36: Common uncertainty of sulfur dioxide slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -6=no common uncertainty input was given, -7=not given since method "MEAS" was chosen, -9=spectral fitting not successful Column 37: rms-based uncertainty of sulfur dioxide slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -3=spectral fitting was done, but no rms-based uncertainty could be retrieved, -9=spectral fitting not successful Column 38: Sulfur dioxide effective temperature [K] Column 39: Independent uncertainty of sulfur dioxide effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -2=no temperature fitting was requested and output for effective temperature and structured uncertainty of it is based on f-code, -3=spectral fitting was done, but no independent uncertainty could be retrieved, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -5=no independent uncertainty input was given, -9=spectral fitting not successful Column 40: Structured uncertainty of sulfur dioxide effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -9=spectral fitting not successful Column 41: Common uncertainty of sulfur dioxide effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -2=no temperature fitting was requested and output for effective temperature and structured uncertainty of it is based on f-code, -3=spectral fitting was done, but no common uncertainty could be retrieved, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -6=no common uncertainty input was given, -9=spectral fitting not successful Column 42: Effective sulfur dioxide fitting wavelength [nm], -9=fitting not successful Column 43: Diffuse correction applied before fitting at effective fitting wavelength for sulfur dioxide [%], 0=no diffuse correction applied or fitting not requested, >0=measured diffuse correction, <0=(negative value of) calculated diffuse correction Column 44: Molecular scattering air mass factor used for molecular scattering subtraction before the fitting Column 45: Estimated uncertainty of molecular scattering air mass factor, -9=molecular scattering was not subtracted before the fitting Column 46: Lower limit used for wavelength scaling [nm] Column 47: Upper limit used for wavelength scaling [nm] Column 48: Order of smoothing polynomial Column 49: Smoothing polynomial coefficient, order 0 Column 50: Independent uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 51: Structured uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 52: Common uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 53: rms-based uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 54: Smoothing polynomial coefficient, order 1 Column 55: Independent uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 56: Structured uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 57: Common uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 58: rms-based uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 59: Smoothing polynomial coefficient, order 2 Column 60: Independent uncertainty of smoothing polynomial coefficient, order 2, -9=fitting not successful Column 61: Structured uncertainty of smoothing polynomial coefficient, order 2, -9=fitting not successful Column 62: Common uncertainty of smoothing polynomial coefficient, order 2, -9=fitting not successful Column 63: rms-based uncertainty of smoothing polynomial coefficient, order 2, -9=fitting not successful Column 64: Smoothing polynomial coefficient, order 3 Column 65: Independent uncertainty of smoothing polynomial coefficient, order 3, -9=fitting not successful Column 66: Structured uncertainty of smoothing polynomial coefficient, order 3, -9=fitting not successful Column 67: Common uncertainty of smoothing polynomial coefficient, order 3, -9=fitting not successful Column 68: rms-based uncertainty of smoothing polynomial coefficient, order 3, -9=fitting not successful Column 69: Mean value of measured data inside fitting window [same units as measurements] Column 70: Wavelength change polynomial coefficient, order 0 Column 71: Independent uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 72: Structured uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 73: Common uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 74: rms-based uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 75: L1 based wavelength change coefficient, order 0, -9=no L1 wavelength change determination Column 76: L1 based wavelength change coefficient, order 1, -9=no L1 wavelength change determination Column 77: L2Fit data quality flag: 0=assured high quality, 1=assured medium quality, 2=assured low quality, 10=not-assured high quality, 11=not-assured medium quality, 12=not-assured low quality Column 78: Sum over 2^i using those i, for which the corresponding L2Fit data quality parameter exceeds the DQ1 limit, 0=L1 data quality above 0, 1=Spectral fitting was not successful, 2=Wavelength shift too large, 3=Normalized rms of fitting residuals weighted with independent uncertainty too large Column 79: Sum over 2^i using those i, for which the corresponding L2Fit data quality parameter exceeds the DQ2 limit (same parameters as for DQ1) Column 80: L1 data quality flag: 0=assured high quality, 1=assured medium quality, 2=assured low quality, 10=not-assured high quality, 11=not-assured medium quality, 12=not-assured low quality Column 81: Sum over 2^i using those i, for which the corresponding L1 data quality parameter exceeds the DQ1 limit, 0=Saturated data, 1=Too few dark counts measurements, 2=No temperature given or effective temperature too different from the reference temperature, 3=Dark count too high, 4=Unsuccessful dark background fitting, 5=The dark count differs significantly from the dark map for too many pixels, 6=Absolute value of estimated average residual stray light level too high, 7=Although attempted, no wavelength change could be retrieved, 8=Absolute value of retrieved wavelength shift too large, 9=Retrieved wavelength shift differs too much from the shift predicted by the effective temperature Column 82: Sum over 2^i using those i, for which the corresponding L1 data quality parameter exceeds the DQ2 limit (same parameters as for DQ1) Column 83: Atmospheric variability [%], 999=no atmospheric variability was determined Column 84: Wavelength effective temperature [°C], 999=no effective temperature given Column 85: Estimated average residual stray light level [%] (only valid for stray light correction methods 2 and higher) Column 86: Retrieved wavelength shift from L1 data [nm], -9=no wavelength change determination Column 87: Retrieved total wavelength shift [nm], -9=no wavelength change fitting Column 88: Number of bright count cycles Column 89: Number of dark count cycles Column 90: Effective position of filterwheel #1, 0=filterwheel not used, 1-9 are valid positions Column 91: Effective position of filterwheel #2, 0=filterwheel not used, 1-9 are valid positions Column 92: Sum over 2^i, 0=spectra were interpolated in time, 1=spectra are corrected for off-target signal Column 93: Integration time [ms] Column 94: Temperature at electronics board [°C], 999=no temperature signal Column 95: Spectrometer control temperature [°C], 999=no temperature signal Column 96: Auxiliary spectrometer temperature [°C], 999=no temperature signal Column 97: Temperature in head sensor [°C], 999=no temperature signal Column 98: Humidity in head sensor [%], -9=no humidity signal Column 99: Pressure in head sensor [hPa], -9=no pressure signal Columns 100-241: Unweighted slant column residuals for each pixel inside the fitting window multiplied by 1e5, 9e5=pixel was not used for fitting Columns 242-383: Normalized slant column residuals weighted with independent instrumental uncertainty for each pixel inside the fitting window multiplied by 1e5, 9e5=pixel was not used for fitting, 0=no independent instrumental uncertainty was given --------------------------------------------------------------------------------------- SS 20231201T155148.6Z 8735.660979 651 1 54.11 2 79.79 232.37 118.50 1.30 0.00 1 0.00 1 0 17 3.011e-01 5.566e-03 2.272e-01 2.262e-03 6.895571e-01 1.156636e-03 7.629784e-03 -6.000000e+00 1.611553e-01 223.42 -2.00 3.42 -2.00 321.08 0.00e+00 3.380386e-03 5.249363e-04 5.722740e-04 -6.000000e+00 7.067486e-03 277.78 -2.00 20.00 -2.00 320.61 0.00e+00 5.413 0.003 302.00 328.00 3 3.9313e-01 3.7401e-03 6.0885e-02 -6.0000e+00 3.9811e-01 8.9340e-02 5.4942e-03 3.6382e-02 -6.0000e+00 4.3335e-01 -6.5540e-02 6.3092e-03 2.0418e-02 -6.0000e+00 2.5792e-01 1.2434e-02 3.0430e-03 8.2433e-03 -6.0000e+00 1.1427e-01 1.1625e-03 -8.1709e-04 5.1895e-05 5.7465e-05 -6.0000e+00 6.3557e-03 -2.7662e-04 0.0000e+00 10 0 0 10 0 0 95.69 20.43 0.99 -0.00028 -0.00642 410 45 1 5 0 131.608 24.50 15.00 20.43 15.82 4.4 884.55 -80340 -73976 -124131 -109730 -96160 -112825 -128921 -92769 -84300 -47718 15433 -46037 -69754 -57230 -62501 -39781 -16515 -21409 -27763 -43472 22481 -26119 -8449 3710 -29705 -22061 -18578 -15409 -30354 -32908 -11461 -23643 1273 -9514 -12704 -23573 -1806 -3996 -5312 -8491 -9769 -543 -3235 -368 376 -1341 -362 768 -1921 721 1487 3196 1667 386 2152 1539 321 -1212 -991 -799 1922 1940 2870 3747 2866 1931 1150 746 753 2293 1468 -5 1897 623 1050 -678 386 -1303 120 -1202 -410 -1508 -1734 -1099 -1293 -877 -458 175 517 404 699 180 245 -259 163 227 964 49 -204 -408 294 66 -662 -93 -1538 -850 -1159 -601 -98 492 299 690 143 536 473 -35 403 -53 -207 -410 -377 186 286 -285 443 389 290 559 -183 183 -195 68 -517 194 5 -574 377 -103 120 -653 428 278 -177 -176 -548 -391 -261 -415 -616 -358 -341 -149 34 -169 -449 -421 -632 -395 -143 -224 -332 -706 194 -351 -107 46 -548 -407 -339 -256 -578 -746 -273 -730 38 -330 -397 -870 -62 -167 -283 -570 -731 -39 -258 -31 33 -123 -34 75 -204 80 193 450 260 68 389 290 64 -237 -186 -148 367 421 714 1071 834 531 335 221 259 919 634 -2 801 250 372 -215 119 -421 46 -565 -220 -958 -1164 -771 -988 -705 -348 134 394 317 498 124 162 -185 137 201 975 59 -258 -492 327 64 -585 -89 -1626 -910 -1338 -747 -133 766 472 1128 230 896 737 -53 632 -88 -336 -686 -597 307 461 -438 652 560 391 846 -297 327 -391 141 -1137 416 11 -1400 899 -258 292 -1490 917 602 SS 20231201T160422.0Z 8735.669699 655 1 54.15 2 81.80 234.50 118.41 4.40 0.00 1 0.00 1 0 25 6.028e-01 9.669e-03 3.216e-01 4.793e-03 8.061535e-01 2.617410e-03 8.828538e-03 -6.000000e+00 3.240668e-01 223.42 -2.00 3.42 -2.00 321.67 0.00e+00 7.824425e-03 1.337468e-03 1.511639e-03 -6.000000e+00 1.518401e-02 277.78 -2.00 20.00 -2.00 321.24 0.00e+00 6.584 0.005 302.00 328.00 3 3.2443e-01 9.0560e-03 7.4991e-02 -6.0000e+00 8.0236e-01 2.5417e-01 1.4840e-02 4.8903e-02 -6.0000e+00 8.7433e-01 -2.9594e-01 1.7565e-02 2.6200e-02 -6.0000e+00 5.1867e-01 1.1287e-01 8.0976e-03 1.3550e-02 -6.0000e+00 2.3059e-01 3.9050e-04 -1.2479e-03 1.1091e-04 8.7447e-05 -6.0000e+00 1.3457e-02 -5.6824e-04 0.0000e+00 12 8 8 10 0 0 35.46 20.42 1.03 -0.00057 -0.00995 266 29 1 5 0 203.190 24.50 15.00 20.42 16.92 4.4 893.33 -96685 -237843 -271216 900000 -274849 -19087 -276207 900000 -82728 -262052 900000 -100162 -224281 65027 900000 -21031 -32479 74704 -10261 -33530 -50453 -119187 -33249 -67048 12617 -27968 -54253 -16810 31979 2064 -55413 -33555 4159 -8091 -40536 -2527 20576 -1328 2891 13913 -1751 -5476 -1819 -3509 -3205 18399 -5438 8323 20941 10885 5825 6504 13048 6688 -491 6137 5013 -1203 4510 -7568 3510 512 3524 4894 3009 7528 4026 1754 5475 3256 5116 3883 1412 -423 315 -3870 1743 -2176 -1268 -965 -1050 -770 -1124 -1084 387 -2262 118 -71 2051 246 336 -394 1795 1320 -1220 -487 388 -487 -4 -512 426 22 -1133 -557 -2642 -2066 -1133 -1646 -397 515 220 313 562 729 334 287 -263 -347 -347 206 -191 942 610 -141 885 529 1063 396 34 556 235 -78 -423 432 -403 -411 517 -295 -208 -417 109 195 -451 -1623 -2923 900000 -2468 -89 -2778 900000 -386 -3365 900000 -468 -3491 304 900000 -98 -190 349 -60 -261 -471 -2087 -272 -835 80 -285 -737 -150 176 18 -1200 -684 72 -168 -1015 -45 332 -32 82 431 -70 -242 -83 -165 -155 750 -293 462 1105 691 460 568 1170 715 -59 708 606 -149 505 -943 416 72 564 914 589 1336 792 361 1263 909 1498 1190 418 -119 77 -871 370 -505 -359 -333 -423 -365 -568 -574 219 -1387 73 -45 1271 146 181 -207 904 731 -828 -387 358 -540 -5 -576 441 19 -915 -479 -2543 -2015 -1184 -1889 -502 687 330 495 877 1101 475 396 -382 -528 -564 339 -296 1380 875 -194 1163 686 1289 551 57 1028 491 -173 -996 1075 -1003 -1067 1310 -788 -538 -999 266 439 SS 20231201T161402.6Z 8735.676419 658 1 54.01 2 83.39 236.11 118.24 6.78 0.00 1 0.00 1 0 25 6.784e-01 2.132e-02 4.281e-01 1.163e-02 9.276512e-01 6.521381e-03 1.022175e-02 -6.000000e+00 3.694127e-01 223.42 -2.00 3.42 -2.00 322.13 0.00e+00 5.115571e-03 3.687866e-03 7.817115e-03 -6.000000e+00 1.671392e-02 277.78 -2.00 20.00 -2.00 321.74 0.00e+00 7.907 0.008 302.00 328.00 3 2.9734e-01 2.5044e-02 1.0021e-01 -6.0000e+00 9.1442e-01 8.0578e-01 3.9015e-02 8.5323e-02 -6.0000e+00 9.9892e-01 -1.0509e+00 4.1547e-02 3.5396e-02 -6.0000e+00 5.9495e-01 4.1963e-01 1.8754e-02 1.0688e-02 -6.0000e+00 2.6312e-01 1.0557e-04 -1.7453e-03 2.6918e-04 1.0135e-04 -6.0000e+00 1.4558e-02 -4.7715e-04 0.0000e+00 12 8 8 10 0 0 24.99 20.43 1.16 -0.00048 -0.01359 134 15 1 5 0 402.664 24.50 15.00 20.43 17.43 4.5 896.75 -85671 900000 -1193 -225926 -195899 -36919 -267502 -157142 -21459 94377 -268360 -176513 -79961 -163244 81747 900000 80482 900000 -43894 900000 -90855 41744 -105208 4994 900000 -49343 -94606 -41220 -44284 900000 -56541 -7173 43085 39881 294256 900000 -41490 -29643 145931 138941 14289 41527 111806 60337 64666 58126 52773 48365 33998 100675 21222 31493 37831 58427 43047 17015 48010 35799 27546 8744 27370 28532 30202 19969 25277 20310 9904 25602 16474 16111 13749 11374 6629 15746 18991 -260 4247 6336 -5289 -2025 2918 -1794 582 2393 1401 1234 346 2221 4272 2011 -519 -2201 -1091 -3119 -4490 -1243 -3963 -2148 -1430 -1308 -2497 -138 -2084 -3151 -647 -754 -2477 -1912 -2755 -1674 -362 462 153 1207 -1140 -1317 -513 201 -1704 795 1094 2743 1201 99 2306 1258 1638 194 -77 1140 1638 997 55 530 507 -347 924 -1038 102 -1037 -958 -802 -970 900000 -14 -4919 -2559 -418 -7222 -1779 -243 1068 -8413 -2613 -905 -4158 925 900000 911 900000 -663 900000 -2143 473 -2779 57 900000 -838 -2499 -636 -684 900000 -1571 -140 621 616 3331 900000 -1366 -1007 1652 1573 500 1154 1611 1428 1464 1436 1412 1444 1260 2072 1136 1685 2064 3175 2827 1483 3152 2513 1965 723 1952 2293 2804 2442 3033 2437 1424 3246 2648 3064 2848 2472 1411 2799 2838 -41 644 1004 -1168 -534 859 -658 223 943 601 557 164 1041 1949 881 -212 -892 -437 -1419 -2506 -802 -3100 -1984 -1407 -1235 -2211 -102 -1387 -2251 -495 -589 -2137 -1810 -2979 -2132 -520 705 235 1803 -1653 -1863 -754 308 -2831 1305 1669 3891 1688 135 2977 1622 1984 275 -132 2004 3271 2135 128 1307 1357 -978 2525 -3014 283 -2667 -2291 -1754