File name: Pandora192s1_Tokyo-Sophia_20250502_L2Fit_fout2c4d20230711p1-8.txt File generation date: 20250502T144955.9Z Data description: Level 2 spectral fitting results file Data file version: fout2c4p1-8 Local principal investigator: Satoshi Inomata Network principal investigator: Alexander Cede Instrument type: Pandora Instrument number: 192 Spectrometer number: 1 Processing software version used: BlickP v1.8.75 Instrument operation file used: Pandora192_OF_v3d20210309.txt Instrument calibration file used: Pandora192s1_CF_v4d20230711.txt Level 1 file used: Pandora192s1_Tokyo-Sophia_20250502_L1_smca1c4d20230711p1-8.txt Full location name: Sophia University Short location name: Tokyo-Sophia Country of location: Japan Location latitude [deg]: 35.6831 Location longitude [deg]: 139.7317 Location altitude [m]: 45 Local noon date: 20250502 Notes on s-number (L1 configuration): Corrections NOT applied although requested by the s-number are dark map for dark correction (replaced by immediate dark method), latency correction, matrix method stray light correction (replaced by simple method) Data caveats: Based on Labanal session 3 for fieldcal First and last pixel inside fitting window: 198 340 Nominal wavelengths inside fitting window [nm]: 305.0295 305.1706 305.3117 305.4528 305.5938 305.7349 305.8759 306.0169 306.1580 306.2990 306.4400 306.5810 306.7219 306.8629 307.0038 307.1448 307.2857 307.4266 307.5676 307.7085 307.8493 307.9902 308.1311 308.2720 308.4128 308.5536 308.6945 308.8353 308.9761 309.1169 309.2577 309.3984 309.5392 309.6800 309.8207 309.9614 310.1022 310.2429 310.3836 310.5243 310.6650 310.8056 310.9463 311.0869 311.2276 311.3682 311.5088 311.6494 311.7900 311.9306 312.0712 312.2117 312.3523 312.4928 312.6334 312.7739 312.9144 313.0549 313.1954 313.3359 313.4763 313.6168 313.7572 313.8977 314.0381 314.1785 314.3189 314.4593 314.5997 314.7401 314.8804 315.0208 315.1611 315.3014 315.4418 315.5821 315.7224 315.8627 316.0029 316.1432 316.2834 316.4237 316.5639 316.7041 316.8444 316.9846 317.1247 317.2649 317.4051 317.5453 317.6854 317.8255 317.9657 318.1058 318.2459 318.3860 318.5260 318.6661 318.8062 318.9462 319.0863 319.2263 319.3663 319.5063 319.6463 319.7863 319.9263 320.0662 320.2062 320.3461 320.4860 320.6260 320.7659 320.9058 321.0456 321.1855 321.3254 321.4652 321.6051 321.7449 321.8847 322.0245 322.1643 322.3041 322.4439 322.5837 322.7234 322.8632 323.0029 323.1426 323.2823 323.4220 323.5617 323.7014 323.8411 323.9807 324.1204 324.2600 324.3996 324.5392 324.6788 324.8184 324.9580 --------------------------------------------------------------------------------------- Column 1: Two letter code of measurement routine Column 2: UT date and time for center-time of measurement, yyyymmddThhmmssZ (ISO 8601) Column 3: Fractional days since 1-Jan-2000 midnight for center-time of measurement Column 4: Routine count (1 for the first routine of the day, 2 for the second, etc.) Column 5: Repetition count (1 for the first set in the routine, 2 for the second, etc.) Column 6: Total duration of measurement set in seconds Column 7: Data processing type index Column 8: Solar zenith angle for center-time of measurement in degree Column 9: Solar azimuth for center-time of measurement in degree, 0=north, increases clockwise Column 10: Lunar zenith angle for center-time of measurement in degree Column 11: Lunar azimuth for center-time of measurement in degree, 0=north, increases clockwise Column 12: Pointing zenith angle in degree, absolute or relative (see next column), 999=tracker not used Column 13: Zenith pointing mode: zenith angle is... 0=absolute, 1=relative to sun, 2=relative to moon Column 14: Pointing azimuth in degree, increases clockwise, absolute (0=north) or relative (see next column), 999=tracker not used Column 15: Azimuth pointing mode: like zenith angle mode but also fixed scattering angles relative to sun (3) or moon (4) Column 16: Fitting result index: 0=no error or warning, 1, 2=warning, >2=error Column 17: Number of function evaluations used, 0=linear fitting or fitting not successful or no fitting done Column 18: rms of unweighted fitting residuals, -9=fitting not successful Column 19: Normalized rms of fitting residuals weighted with independent uncertainty, -9=fitting not successful or no uncertainty used Column 20: Expected rms based on independent uncertainty, -9=fitting not successful or no uncertainty given Column 21: Expected normalized weighted rms based on independent uncertainty, -9=fitting not successful or no uncertainty given Column 22: Ozone slant column amount [moles per square meter], -9e99=fitting not successful Column 23: Independent uncertainty of ozone slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -3=spectral fitting was done, but no independent uncertainty could be retrieved, -5=no independent uncertainty input was given, -9=spectral fitting not successful Column 24: Structured uncertainty of ozone slant column amount [moles per square meter], -7=not given since method "MEAS" was chosen, -9=spectral fitting not successful Column 25: Common uncertainty of ozone slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -6=no common uncertainty input was given, -7=not given since method "MEAS" was chosen, -9=spectral fitting not successful Column 26: rms-based uncertainty of ozone slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -3=spectral fitting was done, but no rms-based uncertainty could be retrieved, -9=spectral fitting not successful Column 27: Ozone effective temperature [K] Column 28: Independent uncertainty of ozone effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -2=no temperature fitting was requested and output for effective temperature and structured uncertainty of it is based on f-code, -3=spectral fitting was done, but no independent uncertainty could be retrieved, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -5=no independent uncertainty input was given, -9=spectral fitting not successful Column 29: Structured uncertainty of ozone effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -9=spectral fitting not successful Column 30: Common uncertainty of ozone effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -2=no temperature fitting was requested and output for effective temperature and structured uncertainty of it is based on f-code, -3=spectral fitting was done, but no common uncertainty could be retrieved, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -6=no common uncertainty input was given, -9=spectral fitting not successful Column 31: Effective ozone fitting wavelength [nm], -9=fitting not successful Column 32: Diffuse correction applied before fitting at effective fitting wavelength for ozone [%], 0=no diffuse correction applied or fitting not requested, >0=measured diffuse correction, <0=(negative value of) calculated diffuse correction Column 33: Sulfur dioxide slant column amount [moles per square meter], -9e99=fitting not successful Column 34: Independent uncertainty of sulfur dioxide slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -3=spectral fitting was done, but no independent uncertainty could be retrieved, -5=no independent uncertainty input was given, -9=spectral fitting not successful Column 35: Structured uncertainty of sulfur dioxide slant column amount [moles per square meter], -7=not given since method "MEAS" was chosen, -9=spectral fitting not successful Column 36: Common uncertainty of sulfur dioxide slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -6=no common uncertainty input was given, -7=not given since method "MEAS" was chosen, -9=spectral fitting not successful Column 37: rms-based uncertainty of sulfur dioxide slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -3=spectral fitting was done, but no rms-based uncertainty could be retrieved, -9=spectral fitting not successful Column 38: Sulfur dioxide effective temperature [K] Column 39: Independent uncertainty of sulfur dioxide effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -2=no temperature fitting was requested and output for effective temperature and structured uncertainty of it is based on f-code, -3=spectral fitting was done, but no independent uncertainty could be retrieved, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -5=no independent uncertainty input was given, -9=spectral fitting not successful Column 40: Structured uncertainty of sulfur dioxide effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -9=spectral fitting not successful Column 41: Common uncertainty of sulfur dioxide effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -2=no temperature fitting was requested and output for effective temperature and structured uncertainty of it is based on f-code, -3=spectral fitting was done, but no common uncertainty could be retrieved, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -6=no common uncertainty input was given, -9=spectral fitting not successful Column 42: Effective sulfur dioxide fitting wavelength [nm], -9=fitting not successful Column 43: Diffuse correction applied before fitting at effective fitting wavelength for sulfur dioxide [%], 0=no diffuse correction applied or fitting not requested, >0=measured diffuse correction, <0=(negative value of) calculated diffuse correction Column 44: Molecular scattering air mass factor used for molecular scattering subtraction before the fitting Column 45: Estimated uncertainty of molecular scattering air mass factor, -9=molecular scattering was not subtracted before the fitting Column 46: Lower limit used for wavelength scaling [nm] Column 47: Upper limit used for wavelength scaling [nm] Column 48: Order of smoothing polynomial Column 49: Smoothing polynomial coefficient, order 0 Column 50: Independent uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 51: Structured uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 52: Common uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 53: rms-based uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 54: Smoothing polynomial coefficient, order 1 Column 55: Independent uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 56: Structured uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 57: Common uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 58: rms-based uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 59: Smoothing polynomial coefficient, order 2 Column 60: Independent uncertainty of smoothing polynomial coefficient, order 2, -9=fitting not successful Column 61: Structured uncertainty of smoothing polynomial coefficient, order 2, -9=fitting not successful Column 62: Common uncertainty of smoothing polynomial coefficient, order 2, -9=fitting not successful Column 63: rms-based uncertainty of smoothing polynomial coefficient, order 2, -9=fitting not successful Column 64: Smoothing polynomial coefficient, order 3 Column 65: Independent uncertainty of smoothing polynomial coefficient, order 3, -9=fitting not successful Column 66: Structured uncertainty of smoothing polynomial coefficient, order 3, -9=fitting not successful Column 67: Common uncertainty of smoothing polynomial coefficient, order 3, -9=fitting not successful Column 68: rms-based uncertainty of smoothing polynomial coefficient, order 3, -9=fitting not successful Column 69: Mean value of measured data inside fitting window [same units as measurements] Column 70: Wavelength change polynomial coefficient, order 0 Column 71: Independent uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 72: Structured uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 73: Common uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 74: rms-based uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 75: L1 based wavelength change coefficient, order 0, -9=no L1 wavelength change determination Column 76: L1 based wavelength change coefficient, order 1, -9=no L1 wavelength change determination Column 77: L2Fit data quality flag: 0=assured high quality, 1=assured medium quality, 2=assured low quality, 10=not-assured high quality, 11=not-assured medium quality, 12=not-assured low quality Column 78: Sum over 2^i using those i, for which the corresponding L2Fit data quality parameter exceeds the DQ1 limit, 0=L1 data quality above 0, 1=Spectral fitting was not successful, 2=Wavelength shift too large, 3=Normalized rms of fitting residuals weighted with independent uncertainty too large Column 79: Sum over 2^i using those i, for which the corresponding L2Fit data quality parameter exceeds the DQ2 limit (same parameters as for DQ1) Column 80: L1 data quality flag: 0=assured high quality, 1=assured medium quality, 2=assured low quality, 10=not-assured high quality, 11=not-assured medium quality, 12=not-assured low quality Column 81: Sum over 2^i using those i, for which the corresponding L1 data quality parameter exceeds the DQ1 limit, 0=Saturated data, 1=Too few dark counts measurements, 2=No temperature given or effective temperature too different from the reference temperature, 3=Dark count too high, 4=Unsuccessful dark background fitting, 5=The dark count differs significantly from the dark map for too many pixels, 6=Absolute value of estimated average residual stray light level too high, 7=Although attempted, no wavelength change could be retrieved, 8=Absolute value of retrieved wavelength shift too large, 9=Retrieved wavelength shift differs too much from the shift predicted by the effective temperature Column 82: Sum over 2^i using those i, for which the corresponding L1 data quality parameter exceeds the DQ2 limit (same parameters as for DQ1) Column 83: Atmospheric variability [%], 999=no atmospheric variability was determined Column 84: Wavelength effective temperature [°C], 999=no effective temperature given Column 85: Estimated average residual stray light level [%] (only valid for stray light correction methods 2 and higher) Column 86: Retrieved wavelength shift from L1 data [nm], -9=no wavelength change determination Column 87: Retrieved total wavelength shift [nm], -9=no wavelength change fitting Column 88: Number of bright count cycles Column 89: Number of dark count cycles Column 90: Effective position of filterwheel #1, 0=filterwheel not used, 1-9 are valid positions Column 91: Effective position of filterwheel #2, 0=filterwheel not used, 1-9 are valid positions Column 92: Sum over 2^i, 0=spectra were interpolated in time, 1=spectra are corrected for off-target signal Column 93: Integration time [ms] Column 94: Temperature at electronics board [°C], 999=no temperature signal Column 95: Spectrometer control temperature [°C], 999=no temperature signal Column 96: Auxiliary spectrometer temperature [°C], 999=no temperature signal Column 97: Temperature in head sensor [°C], 999=no temperature signal Column 98: Humidity in head sensor [%], -9=no humidity signal Column 99: Pressure in head sensor [hPa], -9=no pressure signal Columns 100-242: Unweighted slant column residuals for each pixel inside the fitting window multiplied by 1e5, 9e5=pixel was not used for fitting Columns 243-385: Normalized slant column residuals weighted with independent instrumental uncertainty for each pixel inside the fitting window multiplied by 1e5, 9e5=pixel was not used for fitting, 0=no independent instrumental uncertainty was given --------------------------------------------------------------------------------------- SS 20250502T003957.7Z 9253.027751 194 1 49.72 2 33.26 119.92 76.24 65.80 0.00 1 0.00 1 0 17 2.142e-02 1.348e-02 9.979e-03 3.960e-03 2.024758e-01 2.198816e-03 2.469440e-03 -6.000000e+00 1.368355e-02 230.35 -2.00 3.41 -2.00 316.92 0.00e+00 1.393553e-03 2.866830e-04 5.234157e-04 -6.000000e+00 6.037554e-04 285.27 -2.00 20.00 -2.00 316.16 0.00e+00 1.195 0.000 302.00 328.00 3 6.7238e+00 5.7664e-03 1.5751e-02 -6.0000e+00 3.4118e-02 6.8725e-02 6.4124e-03 1.0100e-02 -6.0000e+00 3.6723e-02 -6.6536e-03 3.7462e-03 4.9302e-03 -6.0000e+00 2.1877e-02 -7.1351e-03 1.7513e-03 1.3781e-03 -6.0000e+00 9.1761e-03 1.4828e-04 -3.4024e-03 1.3840e-04 1.9905e-05 -6.0000e+00 7.3595e-04 -9.6979e-02 0.0000e+00 12 8 8 10 0 0 99.86 18.11 0.69 -0.09698 -0.12255 19 4 1 5 0 2616.440 22.40 14.99 18.11 22.12 30.2 981.02 900000 -3612 5208 -611 -2254 -4740 -1060 -2289 -3393 -5378 -2647 -2447 -7555 -5959 -2869 -408 147 -438 739 2741 2141 4611 2303 1603 1044 1604 -442 -1299 -781 -3060 -3074 910 2666 -1470 -465 -2626 -7362 -3703 -3351 2994 1298 1576 456 902 677 1699 309 133 -840 -933 -943 813 1726 990 1967 750 1781 570 -493 -927 -32 752 2033 2857 1828 1109 -302 -1585 -1457 -976 1761 774 974 1303 492 160 -3264 -4341 -4728 -4406 -3206 -2301 -2230 -1548 -731 -290 1285 1458 1740 2002 1234 838 -1468 -1499 -1461 -1720 -114 -406 499 647 -123 -772 -1263 -2004 -2219 -1143 -976 -485 730 1090 1718 1914 1961 1035 156 -586 -1235 -908 -431 933 967 1273 1312 963 420 -465 -849 -738 -1319 -702 -435 -1014 -686 -925 -438 61 822 772 567 349 -163 -99 -1021 900000 -542 768 -96 -345 -712 -146 -327 -528 -919 -458 -428 -1434 -1205 -626 -97 37 -118 210 784 631 1344 689 489 330 521 -149 -436 -249 -986 -969 282 855 -514 -163 -918 -2697 -1402 -1442 1404 703 926 284 575 437 1083 197 83 -517 -563 -569 498 1100 667 1363 539 1289 421 -366 -686 -24 567 1603 2370 1584 955 -248 -1288 -1158 -809 1546 710 926 1236 454 143 -2767 -3412 -3577 -3423 -2679 -2100 -2213 -1629 -777 -327 1512 1798 2191 2501 1494 986 -1620 -1534 -1455 -1782 -126 -491 622 836 -158 -983 -1522 -2365 -2693 -1451 -1307 -678 1053 1578 2588 2963 3023 1545 228 -800 -1614 -1179 -578 1308 1425 1913 1959 1404 590 -627 -1095 -956 -1730 -954 -596 -1452 -1004 -1392 -681 99 1380 1335 998 615 -283 -168 -1693 SS 20250502T005744.0Z 9253.040093 206 1 13.06 2 30.22 125.67 73.05 67.91 0.00 1 0.00 1 0 25 2.015e-02 1.364e-02 7.688e-03 4.033e-03 1.950316e-01 2.261361e-03 2.415945e-03 -6.000000e+00 1.287916e-02 230.35 -2.00 3.41 -2.00 316.00 0.00e+00 1.193943e-03 2.393133e-04 4.049206e-04 -6.000000e+00 5.682020e-04 285.28 -2.00 20.00 -2.00 315.30 0.00e+00 1.157 0.000 302.00 328.00 3 6.2043e+00 5.8664e-03 1.5183e-02 -6.0000e+00 3.2113e-02 6.7270e-02 6.4922e-03 9.6630e-03 -6.0000e+00 3.4561e-02 -1.7793e-02 3.7258e-03 4.7100e-03 -6.0000e+00 2.0589e-02 5.5069e-05 1.6758e-03 1.2986e-03 -6.0000e+00 8.6386e-03 2.6417e-04 -3.3078e-03 1.3953e-04 2.2803e-05 -6.0000e+00 6.9252e-04 -9.8779e-02 0.0000e+00 12 9 9 12 1 1 99.97 18.07 0.72 -0.09878 -0.12363 4 4 1 5 0 2609.600 22.40 15.00 18.07 23.41 30.2 986.94 900000 -598 887 1072 -5484 -969 -9875 -1340 312 -3674 -4295 -4617 -3863 -1808 415 -910 2364 871 931 622 2869 3134 3839 1598 2086 872 -914 106 -905 -2108 -2958 -95 -408 -848 -2676 -3231 -5678 -4899 -1431 -604 601 1585 1986 1925 603 606 1375 -645 -800 350 565 1746 814 1241 539 1051 1019 460 798 -1094 -360 742 2856 2135 1741 881 711 -1027 -1144 402 651 1836 773 827 990 -1146 -2176 -4215 -4827 -3999 -3380 -2294 -1838 -1446 -1548 -672 586 1294 1066 1447 764 404 -661 -2234 -2235 -1334 90 1106 723 474 -382 -1055 -667 -1846 -1496 -1447 -898 -504 793 1377 1909 1766 1624 1158 607 -561 -970 -788 257 931 1312 1307 1052 1441 657 599 -402 -478 -1061 -1032 -993 -973 -1385 -803 -564 496 241 564 667 315 -44 -239 -1037 900000 -128 200 243 -1264 -197 -2088 -266 66 -864 -1056 -1151 -981 -485 121 -297 805 315 347 238 1100 1213 1480 637 851 372 -404 46 -376 -842 -1214 -39 -175 -364 -1183 -1402 -2520 -2310 -743 -346 385 1095 1435 1369 438 439 1019 -477 -575 243 391 1231 608 967 421 839 826 378 650 -901 -296 622 2478 1977 1608 841 641 -882 -962 349 612 1765 780 835 963 -1107 -1938 -3658 -4034 -3408 -3112 -2209 -1916 -1540 -1747 -770 705 1622 1324 1846 944 467 -735 -2431 -2364 -1449 100 1332 908 619 -497 -1311 -808 -2200 -1824 -1796 -1171 -686 1073 1938 2796 2574 2360 1680 826 -743 -1227 -1029 332 1266 1804 1840 1476 2034 896 778 -521 -599 -1344 -1363 -1327 -1350 -1977 -1167 -819 746 379 879 1056 500 -68 -377 -1605 SS 20250502T010545.1Z 9253.045661 212 1 49.20 2 28.92 128.57 71.60 68.84 0.00 1 0.00 1 0 17 1.703e-02 1.209e-02 3.088e-03 1.403e-03 1.873272e-01 7.821682e-04 2.305399e-03 -6.000000e+00 1.088151e-02 230.35 -2.00 3.41 -2.00 316.45 0.00e+00 1.389757e-03 9.200048e-05 4.426615e-04 -6.000000e+00 4.800630e-04 285.28 -2.00 20.00 -2.00 315.72 0.00e+00 1.142 0.000 302.00 328.00 3 5.1562e+00 2.0409e-03 1.4833e-02 -6.0000e+00 2.7132e-02 5.2784e-02 2.2665e-03 9.3228e-03 -6.0000e+00 2.9200e-02 -1.4242e-02 1.3109e-03 4.5519e-03 -6.0000e+00 1.7395e-02 -3.2309e-03 6.0106e-04 1.2724e-03 -6.0000e+00 7.2989e-03 7.8067e-04 -3.1378e-03 4.8696e-05 1.9922e-05 -6.0000e+00 5.8508e-04 -9.9030e-02 0.0000e+00 12 8 8 10 0 0 99.98 18.07 0.60 -0.09903 -0.12261 97 21 1 5 0 507.130 22.40 15.00 18.07 24.15 30.1 990.26 900000 170 1247 288 -2128 -2666 -4377 -2097 -2272 -3042 -2975 -4013 -2322 -2987 335 78 1036 1632 1947 2508 2271 2584 1792 1705 1256 666 334 -1224 -2075 -2025 -2782 -1042 -1334 -812 -3280 -4722 -4962 -4327 -1741 -215 1431 527 837 1312 680 1017 853 765 -300 -529 -278 212 1345 1073 1018 851 928 187 230 271 477 1295 1602 2624 2164 1141 -182 -1089 -1450 -4 1205 1302 1065 937 213 -796 -2453 -3652 -3992 -3517 -2817 -2653 -2050 -2077 -1094 -79 737 1499 1401 1519 799 226 -748 -1804 -2015 -1337 -479 341 865 847 161 -737 -1519 -1695 -1483 -1178 -851 -483 313 1062 1738 1877 1637 1120 171 -381 -722 -594 367 575 1077 938 1205 918 467 106 -350 -567 -1036 -653 -648 -801 -804 -643 -219 -2 432 424 498 530 118 -302 -1039 900000 30 236 55 -398 -463 -736 -352 -410 -595 -604 -831 -485 -677 82 21 298 487 597 786 729 837 587 564 429 237 133 -487 -796 -737 -995 -373 -503 -318 -1325 -1890 -1987 -1844 -824 -116 862 317 526 841 443 724 602 528 -203 -351 -183 143 943 790 778 668 739 151 185 216 382 1061 1384 2376 2033 1067 -162 -917 -1199 -3 1103 1281 1091 962 212 -751 -2147 -2964 -3112 -2806 -2425 -2510 -2096 -2266 -1250 -86 850 1801 1727 1861 949 251 -848 -1912 -2080 -1420 -545 421 1040 1057 200 -963 -1886 -2052 -1842 -1543 -1067 -633 424 1496 2542 2820 2452 1616 233 -532 -961 -788 501 829 1480 1324 1686 1254 676 146 -462 -727 -1341 -877 -914 -1175 -1207 -988 -349 -3 677 686 818 871 191 -477 -1613