File name: Pandora132s1_Berlin_20220116_L2Fit_fout2c13d20210422p1-8.txt File generation date: 20241123T162925.7Z Data description: Level 2 spectral fitting results file Data file version: fout2c13p1-8 Local principal investigator: Thomas Ruhtz Network principal investigator: Alexander Cede Instrument type: Pandora Instrument number: 132 Spectrometer number: 1 Processing software version used: BlickP v1.8.60 Instrument operation file used: Pandora132_OF_v1d20200910.txt Instrument calibration file used: Pandora132s1_CF_v13d20210422.txt Level 1 file used: Pandora132s1_Berlin_20220116_L1_smca1c13d20210422p1-8.txt Full location name: Institute of Space Sciences FU Berlin Short location name: Berlin Country of location: Germany Location latitude [deg]: 52.4600 Location longitude [deg]: 13.3100 Location altitude [m]: 100 Local noon date: 20220116 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: Modification of session 11, but different RefDay (April 2021). First and last pixel inside fitting window: 243 385 Nominal wavelengths inside fitting window [nm]: 305.1071 305.2480 305.3889 305.5298 305.6707 305.8116 305.9524 306.0932 306.2341 306.3749 306.5157 306.6565 306.7973 306.9380 307.0788 307.2195 307.3603 307.5010 307.6417 307.7824 307.9231 308.0638 308.2044 308.3451 308.4857 308.6264 308.7670 308.9076 309.0482 309.1888 309.3293 309.4699 309.6104 309.7510 309.8915 310.0320 310.1725 310.3130 310.4534 310.5939 310.7344 310.8748 311.0152 311.1556 311.2960 311.4364 311.5768 311.7172 311.8575 311.9979 312.1382 312.2785 312.4188 312.5591 312.6994 312.8397 312.9800 313.1202 313.2604 313.4007 313.5409 313.6811 313.8213 313.9615 314.1016 314.2418 314.3819 314.5221 314.6622 314.8023 314.9424 315.0825 315.2225 315.3626 315.5026 315.6427 315.7827 315.9227 316.0627 316.2027 316.3427 316.4827 316.6226 316.7626 316.9025 317.0424 317.1823 317.3222 317.4621 317.6020 317.7418 317.8817 318.0215 318.1613 318.3012 318.4410 318.5807 318.7205 318.8603 319.0000 319.1398 319.2795 319.4192 319.5589 319.6986 319.8383 319.9780 320.1177 320.2573 320.3969 320.5366 320.6762 320.8158 320.9554 321.0949 321.2345 321.3741 321.5136 321.6531 321.7927 321.9322 322.0717 322.2111 322.3506 322.4901 322.6295 322.7689 322.9084 323.0478 323.1872 323.3266 323.4659 323.6053 323.7447 323.8840 324.0233 324.1627 324.3020 324.4413 324.5805 324.7198 324.8591 324.9983 --------------------------------------------------------------------------------------- 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 20220116T092859.7Z 8051.395135 342 1 16.80 2 77.01 154.32 101.80 352.50 0.00 1 0.00 1 0 17 3.371e-01 2.132e-02 2.309e-01 1.274e-02 5.564394e-01 6.659557e-03 7.054162e-03 -6.000000e+00 1.913831e-01 221.58 -2.00 3.85 -2.00 320.54 0.00e+00 5.483378e-03 2.305382e-03 9.201485e-04 -6.000000e+00 8.788064e-03 265.91 -2.00 20.00 -2.00 320.12 0.00e+00 4.336 0.002 302.00 328.00 3 3.9608e+00 1.9339e-02 5.2421e-02 -6.0000e+00 4.7288e-01 2.0313e-01 2.3814e-02 3.1214e-02 -6.0000e+00 5.1922e-01 -1.6557e-01 2.2871e-02 1.7527e-02 -6.0000e+00 3.0478e-01 6.1903e-02 1.1792e-02 8.1338e-03 -6.0000e+00 1.2846e-01 7.3387e-05 -6.7979e-04 3.7555e-04 9.0464e-05 -6.0000e+00 9.4620e-03 -2.4399e-02 0.0000e+00 12 9 9 12 1 1 99.95 15.48 0.07 -0.02440 -0.02951 6 4 1 8 0 2375.280 15.01 15.48 7.71 20.04 20.46 22.8 984.64 -187543 900000 -42961 443 900000 -99280 900000 -87411 900000 -87828 900000 -111049 -61389 150431 -50684 900000 -27278 31408 -43673 66660 23564 3148 5738 -26812 121703 -27180 68874 -931 29967 -902 17402 46276 -24110 21876 -21622 37726 -17690 62407 8864 8763 10201 -2018 24076 -11131 22116 -5814 21165 685 16802 39 8798 17107 2750 18695 -5909 15566 1326 -2527 -244 -1727 5631 -4534 14301 -2178 8681 -3260 -1099 -390 1381 4380 -855 4948 -5293 8764 -3206 3023 -3473 800 453 -3647 3623 -5807 623 -4692 167 858 -1081 4404 -2605 95 -2231 312 -3660 -1013 -948 -818 1215 561 980 -1187 920 -3527 583 -1114 -2507 -486 -2078 -436 -694 1013 486 1509 750 -54 -1481 -1966 1605 -502 1391 -1813 1006 732 1637 1453 961 2090 -1813 2309 699 367 -318 -153 -178 102 -41 -1361 1055 -149 564 5 757 -354 -2073 -8921 900000 -630 5 900000 -2128 900000 -1972 900000 -2279 900000 -4259 -1801 1867 -2116 900000 -1169 815 -2562 1354 744 121 220 -1512 1812 -1686 1632 -43 1014 -45 844 1984 -2366 1407 -2049 1958 -1756 3309 970 1117 1298 -309 3232 -1991 3144 -1022 3152 118 2767 7 1674 3308 665 4118 -1736 3799 378 -739 -68 -484 1532 -1510 4451 -881 3241 -1321 -423 -153 590 2066 -462 2599 -3016 4078 -1532 1219 -1397 313 200 -1965 2098 -4149 451 -3784 137 741 -1000 3939 -2473 84 -1878 237 -2770 -771 -804 -796 1225 641 1191 -1456 1022 -3680 591 -1178 -2659 -534 -2439 -537 -936 1486 777 2476 1222 -83 -2198 -2817 2242 -750 2165 -3015 1613 1134 2423 2079 1326 2729 -2393 3071 1055 617 -582 -292 -354 214 -84 -2994 2313 -334 1247 10 1629 -728 -4261 SS 20220116T094555.2Z 8051.406889 351 1 49.20 2 75.97 158.24 102.05 356.23 0.00 1 0.00 1 0 17 2.612e-01 1.321e-02 2.245e-01 1.137e-02 5.311202e-01 5.873201e-03 6.639655e-03 -6.000000e+00 1.473081e-01 221.59 -2.00 3.85 -2.00 320.52 0.00e+00 3.986775e-03 1.985535e-03 7.047733e-04 -6.000000e+00 6.289087e-03 265.91 -2.00 20.00 -2.00 320.11 0.00e+00 4.035 0.001 302.00 328.00 3 4.4422e+00 1.6964e-02 4.9470e-02 -6.0000e+00 3.6274e-01 2.1221e-01 2.0895e-02 3.0732e-02 -6.0000e+00 3.9759e-01 -1.4082e-01 2.0209e-02 1.5543e-02 -6.0000e+00 2.3431e-01 4.9012e-02 1.0474e-02 6.2439e-03 -6.0000e+00 9.8209e-02 6.1016e-05 -1.0095e-03 3.3924e-04 6.4279e-05 -6.0000e+00 7.2652e-03 -2.4397e-02 0.0000e+00 12 8 8 10 0 0 99.98 15.50 0.14 -0.02440 -0.03198 20 4 1 8 0 2459.370 15.00 15.50 8.15 20.04 20.46 23.0 986.38 -107761 -57946 -53565 -9565 900000 -141249 900000 -122966 50620 -71836 13318 -82226 -57834 -13084 -20526 87616 -29048 67338 -23179 25639 20812 1272 14382 1314 35667 -15833 6552 -10843 -10377 16588 16205 7211 2067 26068 3598 8125 7946 13158 19328 -5544 3056 10501 7927 3050 8407 -1675 12356 -12044 -1312 5070 3511 13669 -3561 3634 619 7193 -4084 849 2466 -6968 4370 920 4329 -306 -640 646 2745 1670 487 3670 -1476 2266 -2062 2124 -2826 -3402 -499 1928 -3038 -536 -2267 -2416 -1920 -1447 -686 1653 181 367 179 1334 400 806 -313 1663 -2371 -927 413 494 765 83 63 -1670 -1119 29 -124 -1049 -1290 -777 -543 1407 -277 -83 -199 -1157 1519 -817 219 84 810 -1228 -340 99 2037 460 1338 569 -182 1150 773 98 -68 -673 -154 436 112 -433 -138 485 657 -190 -21 -397 -1124 -2824 -1051 -972 -106 900000 -5126 900000 -4462 561 -1810 162 -2652 -1807 -330 -683 1263 -1347 1419 -1168 827 713 53 536 60 1222 -942 310 -590 -554 719 817 465 158 1629 268 573 617 1154 1929 -827 462 1567 1262 520 1373 -307 1993 -2525 -254 971 747 2893 -969 982 180 2016 -1301 254 706 -2170 1264 303 1549 -126 -271 259 1040 653 213 1768 -814 1242 -1158 1077 -1373 -1492 -197 752 -1384 -278 -1380 -1651 -1409 -1116 -558 1398 163 338 164 1150 323 600 -224 1199 -1972 -869 436 585 964 104 73 -1773 -1105 29 -133 -1200 -1556 -990 -752 2104 -456 -141 -332 -1850 2213 -1167 312 125 1274 -2044 -558 156 3034 671 1854 759 -237 1544 1157 164 -122 -1274 -301 897 246 -996 -324 1144 1535 -430 -45 -804 -2251 SS 20220116T095627.1Z 8051.414203 357 1 49.10 2 75.41 160.71 102.12 358.55 0.00 1 0.00 1 0 17 1.667e-01 6.591e-03 1.029e-01 4.142e-03 5.142048e-01 2.103297e-03 6.416752e-03 -6.000000e+00 9.312264e-02 221.59 -2.00 3.85 -2.00 320.47 0.00e+00 1.743408e-03 7.129987e-04 6.375105e-04 -6.000000e+00 3.935998e-03 265.91 -2.00 20.00 -2.00 320.05 0.00e+00 3.889 0.001 302.00 328.00 3 3.1583e+00 6.0447e-03 4.7613e-02 -6.0000e+00 2.2938e-01 9.7763e-02 7.4657e-03 2.9183e-02 -6.0000e+00 2.5122e-01 -3.9126e-02 7.3036e-03 1.4976e-02 -6.0000e+00 1.4888e-01 4.3046e-03 3.7919e-03 5.0613e-03 -6.0000e+00 6.2672e-02 2.6998e-04 -2.9820e-04 1.2164e-04 5.6639e-05 -6.0000e+00 4.5471e-03 -2.5280e-02 0.0000e+00 11 8 0 10 0 0 99.97 15.51 0.19 -0.02528 -0.02752 62 13 1 8 0 791.860 15.00 15.51 8.38 20.04 20.46 23.0 987.87 -4311 -50113 29478 -81219 -83572 -45639 -62363 -55561 -62974 -20033 -63772 -7394 -31834 -31069 -8456 12125 -3508 -16467 11266 -8819 1755 -3796 603 693 5927 -4446 -5256 4175 -12900 -2150 -4995 -313 1483 7206 2858 -1329 1387 -4701 7558 -536 -1409 -1945 896 -804 3023 -802 -2857 1939 -3364 2012 1500 2236 1271 1021 2374 1694 2112 796 689 -690 2015 2092 2261 1843 352 991 844 444 308 1977 174 725 1163 -733 197 -1661 -1486 -1161 -2627 -1556 -1961 -1423 -2006 -1668 -1319 -603 78 670 1310 986 1457 244 -334 -1025 -32 250 266 1054 259 -219 -377 -757 -300 206 437 -1130 -378 -824 225 25 286 -46 189 -111 -397 -38 390 193 -45 -28 -551 241 579 -48 601 -204 -62 930 689 -218 -152 -195 -244 -27 56 -16 -37 749 -102 -86 297 -423 -510 -35 -729 202 -1608 -1553 -571 -1008 -965 -1323 -291 -1494 -111 -746 -904 -239 323 -120 -698 387 -388 71 -163 26 32 282 -246 -293 204 -735 -116 -321 -22 117 559 222 -107 121 -524 910 -75 -219 -324 151 -140 513 -143 -527 347 -652 391 321 525 327 281 677 504 638 242 204 -205 607 693 838 751 149 399 329 178 137 985 96 411 645 -391 95 -734 -610 -480 -1223 -836 -1218 -984 -1506 -1320 -1109 -535 73 631 1212 876 1198 188 -246 -782 -26 238 288 1267 334 -280 -447 -811 -301 211 480 -1333 -465 -1083 318 39 480 -79 321 -179 -600 -55 566 294 -73 -47 -927 390 900 -73 862 -283 -83 1289 1062 -376 -280 -376 -490 -57 115 -34 -82 1658 -227 -184 652 -885 -1053