File name: Pandora166s1_PhiladelphiaPA_20231122_L2Fit_fout2c5d20220925p1-8.txt
File generation date: 20231123T052143.6Z
Data description: Level 2 spectral fitting results file
Data file version: fout2c5p1-8
Local principal investigator: Lukas Valin
Network principal investigator: Alexander Cede
Instrument type: Pandora
Instrument number: 166
Spectrometer number: 1
Processing software version used: BlickP v1.8.49
Instrument operation file used: Pandora166_OF_v5d20210217.txt
Instrument calibration file used: Pandora166s1_CF_v5d20220925.txt
Level 1 file used: Pandora166s1_PhiladelphiaPA_20231122_L1_smca1c5d20220925p1-8.txt
Full location name: Philadelphia, Pennsylvania; Philadelphia-Camden-Wilmington; North East Waste (NE
Short location name: PhiladelphiaPA
Country of location: United States
Location latitude [deg]: 39.9919
Location longitude [deg]: -75.0811
Location altitude [m]: 6
Local noon date: 20231122
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: HCHO modification of session 4 "Field cal based on session3"
First and last pixel inside fitting window: 223 365
Nominal wavelengths inside fitting window [nm]: 305.0654 305.2061 305.3468 305.4875 305.6282 305.7689 305.9096 306.0503 306.1910 306.3316 306.4722 306.6129 306.7535 306.8941 307.0347 307.1753 307.3158 307.4564 307.5969 307.7375 307.8780 308.0185 308.1590 308.2995 308.4400 308.5805 308.7209 308.8614 309.0018 309.1422 309.2827 309.4231 309.5635 309.7038 309.8442 309.9846 310.1249 310.2652 310.4056 310.5459 310.6862 310.8265 310.9668 311.1070 311.2473 311.3875 311.5278 311.6680 311.8082 311.9484 312.0886 312.2288 312.3690 312.5091 312.6493 312.7894 312.9295 313.0696 313.2097 313.3498 313.4899 313.6300 313.7700 313.9101 314.0501 314.1901 314.3301 314.4701 314.6101 314.7501 314.8901 315.0300 315.1700 315.3099 315.4498 315.5897 315.7296 315.8695 316.0094 316.1493 316.2891 316.4290 316.5688 316.7086 316.8484 316.9882 317.1280 317.2678 317.4076 317.5473 317.6871 317.8268 317.9665 318.1062 318.2459 318.3856 318.5253 318.6650 318.8046 318.9443 319.0839 319.2235 319.3631 319.5027 319.6423 319.7819 319.9215 320.0610 320.2006 320.3401 320.4796 320.6191 320.7586 320.8981 321.0376 321.1770 321.3165 321.4559 321.5954 321.7348 321.8742 322.0136 322.1530 322.2923 322.4317 322.5711 322.7104 322.8497 322.9891 323.1284 323.2677 323.4069 323.5462 323.6855 323.8247 323.9640 324.1032 324.2424 324.3816 324.5208 324.6600 324.7992 324.9383
---------------------------------------------------------------------------------------
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 20231122T194652.1Z 8726.824215 384 1 41.87 2 73.41 223.92 83.38 98.16 0.00 1 0.00 1 0 17 5.143e-01 1.620e-02 2.202e-01 1.145e-02 4.403946e-01 5.856017e-03 5.102105e-03 -6.000000e+00 2.909759e-01 223.32 -2.00 3.49 -2.00 320.05 0.00e+00 3.392767e-03 1.643587e-03 2.420676e-04 -6.000000e+00 1.219225e-02 280.46 -2.00 20.00 -2.00 319.73 0.00e+00 3.446 0.001 302.00 328.00 3 4.5491e+00 1.6315e-02 4.0896e-02 -6.0000e+00 7.2121e-01 7.3106e-02 1.9319e-02 2.2373e-02 -6.0000e+00 7.8772e-01 -1.7749e-02 1.7880e-02 1.1710e-02 -6.0000e+00 4.6536e-01 5.5267e-03 9.5884e-03 3.6698e-03 -6.0000e+00 1.9828e-01 1.0285e-04 -1.3254e-04 3.4714e-04 4.3532e-05 -6.0000e+00 1.4319e-02 -1.8282e-02 -2.7578e-05 12 9 9 12 9 1 99.99 19.40 0.71 -0.01828 -0.01928 20 5 1 5 0 1993.720 26.21 19.98 19.40 19.06 23.5 971.87 -54046 900000 -82606 900000 14357 116754 450840 -17784 36850 -8436 302055 -18088 48132 -14138 -25815 28778 6631 110200 -26889 26338 12317 10637 69476 1296 65089 -12535 4960 -4749 3407 40220 48070 23963 -17631 17664 -10602 17251 -12434 -7075 18095 -3927 3379 -20654 4557 -1465 1760 -753 4196 -1110 6619 -2020 -8247 7119 4079 2004 4936 2976 2258 -875 9074 -4879 -2525 1860 -3745 2995 -2792 -1798 -4155 2137 1601 1069 1035 1952 5491 -1710 3924 -4920 -169 -2935 -2173 4349 -3924 -1174 -1309 -833 -2769 -798 424 695 1890 1055 378 754 -2505 -356 348 91 1192 1595 -1345 2054 -154 -2182 1198 -1209 600 -2821 -11 -431 106 1515 -441 -47 -227 270 -369 1209 1696 -803 526 -1628 -484 -637 1523 -184 901 711 -1618 1312 257 807 -104 -323 736 -478 398 -1371 -58 -574 1114 -45 -330 -7 292 -1402 900000 -3296 900000 221 1304 5034 -404 572 -207 3373 -525 917 -557 -1337 1058 322 2319 -2002 1276 693 606 2418 87 2618 -1045 359 -364 236 2064 2543 1773 -2122 1633 -1223 1513 -1536 -959 2313 -721 644 -5253 956 -334 393 -175 936 -262 1471 -504 -2340 1907 1213 650 1626 1025 800 -320 2962 -1803 -940 718 -1689 1388 -1416 -869 -1956 942 758 576 615 1207 3285 -1054 2130 -2634 -79 -1375 -1074 2311 -2579 -843 -1020 -681 -2452 -740 410 694 1873 1025 346 637 -2045 -280 291 85 1250 1869 -1749 2606 -189 -2506 1235 -1276 668 -3437 -14 -566 148 2255 -718 -80 -381 431 -555 1695 2345 -1176 807 -2684 -802 -1039 2370 -281 1306 992 -2213 1756 376 1283 -179 -581 1355 -925 809 -2973 -127 -1291 2482 -98 -698 -14 569
SS 20231122T195121.7Z 8726.827335 387 1 50.82 2 74.01 224.82 82.54 98.86 0.00 1 0.00 1 0 17 8.503e-02 5.665e-03 5.698e-02 2.542e-03 4.659353e-01 1.297507e-03 5.324881e-03 -6.000000e+00 4.813453e-02 223.32 -2.00 3.49 -2.00 320.24 0.00e+00 1.398831e-03 3.871016e-04 2.879301e-04 -6.000000e+00 1.964496e-03 280.46 -2.00 20.00 -2.00 319.92 0.00e+00 3.568 0.001 302.00 328.00 3 1.9746e+00 3.6523e-03 4.2617e-02 -6.0000e+00 1.1931e-01 5.5062e-03 4.3996e-03 2.3686e-02 -6.0000e+00 1.3027e-01 5.1924e-02 4.2760e-03 1.2558e-02 -6.0000e+00 7.6843e-02 -3.3302e-02 2.2740e-03 4.0002e-03 -6.0000e+00 3.2401e-02 1.1996e-03 1.2777e-04 7.6917e-05 4.5976e-05 -6.0000e+00 2.3348e-03 -2.1611e-02 5.6010e-05 10 0 0 10 0 0 100.00 19.39 0.13 -0.02161 -0.02065 362 65 1 5 0 140.370 26.21 20.00 19.39 19.22 23.4 973.43 -37981 25247 -26943 9098 -8052 -13165 65402 -3417 25428 -25124 2201 -6723 8273 -1394 2205 15112 4080 15615 2119 12237 -7508 -4282 1617 -3781 6712 -3579 3194 -3012 -528 3189 2123 7140 1487 4667 -1189 -1199 166 1004 356 -2233 -2857 -4177 -883 -1824 -688 36 -1523 -532 -2134 -523 613 1600 940 1164 531 -584 1036 -2200 581 -1700 1000 1668 1580 3137 -841 427 -96 -315 811 325 528 38 839 332 623 -57 148 -599 -944 -213 -1203 -501 -2394 -2095 -823 -653 37 381 774 527 920 953 323 1228 600 713 106 1377 466 -131 -239 -1172 -222 426 -594 -1549 -715 -783 -430 276 182 93 169 120 11 266 427 41 485 -214 -470 -294 299 449 242 91 -534 1039 416 182 -108 -674 -66 -201 81 -455 -12 384 157 150 1 -243 -25 -626 286 -488 123 -116 -194 546 -55 360 -567 41 -144 182 -39 71 498 160 598 97 537 -400 -224 83 -217 377 -231 199 -193 -32 193 143 521 127 414 -108 -107 16 109 46 -351 -506 -799 -172 -368 -139 7 -317 -111 -459 -115 144 408 261 347 167 -192 345 -759 192 -568 334 598 628 1351 -392 190 -41 -135 368 168 301 23 504 193 336 -28 67 -266 -447 -115 -747 -347 -1839 -1693 -699 -591 35 372 758 503 822 791 254 943 496 660 111 1601 591 -168 -293 -1331 -232 441 -663 -1851 -888 -1025 -596 414 294 157 283 194 17 381 606 60 751 -351 -784 -482 474 689 357 129 -735 1417 616 296 -188 -1243 -125 -398 168 -1000 -26 875 362 339 3 -506 -51
SS 20231122T195513.5Z 8726.830017 390 1 1.74 2 74.53 225.59 81.82 99.46 0.00 1 0.00 1 0 17 8.822e-02 8.466e-03 5.206e-02 4.131e-03 4.815986e-01 2.066930e-03 5.618341e-03 -6.000000e+00 4.992760e-02 223.32 -2.00 3.49 -2.00 319.11 0.00e+00 7.759313e-04 4.636506e-04 2.059586e-04 -6.000000e+00 2.037902e-03 280.46 -2.00 20.00 -2.00 318.81 0.00e+00 3.681 0.001 302.00 328.00 3 8.4622e-01 5.5891e-03 4.3893e-02 -6.0000e+00 1.2375e-01 5.3095e-02 6.4185e-03 2.4177e-02 -6.0000e+00 1.3514e-01 -2.8677e-02 5.2351e-03 1.2415e-02 -6.0000e+00 7.9712e-02 3.9101e-03 2.8492e-03 3.7720e-03 -6.0000e+00 3.3600e-02 3.2662e-03 -3.7590e-04 1.1657e-04 4.0636e-05 -6.0000e+00 2.4226e-03 -2.0350e-02 2.5785e-05 12 9 1 12 1 1 97.56 19.40 0.35 -0.02035 -0.02317 8 48 1 5 0 192.560 26.21 20.01 19.40 19.51 23.4 975.46 21871 -24962 15810 -44287 -13601 -5179 -36609 21137 -40991 2074 -42092 7072 -18151 -12483 8058 -13120 10695 -11500 17437 -2328 -121 4535 2458 4421 -8563 7767 -8472 -2253 -2923 -4699 1876 -3253 3294 -1433 -1615 -4611 -2390 1127 -3244 1779 -3841 328 -1071 -669 1362 -533 705 -49 1217 -1201 2708 544 1590 3957 256 1841 -357 395 544 70 1042 1918 3003 3340 1173 -648 -194 -219 1893 469 325 933 25 410 -1233 -398 -1143 -1917 451 -1714 -702 -2301 -1789 -1627 -1473 -1072 -541 865 -45 1029 -735 667 134 212 1086 -515 815 668 761 -710 19 -964 -776 594 -1075 -1341 -685 -493 157 -45 366 484 558 88 296 219 592 739 353 859 -1237 147 252 276 666 106 173 1243 433 -437 -202 -910 -210 -412 -347 -59 -345 762 -27 33 -12 -146 -134 357 -679 322 -1526 -332 -112 -1099 426 -1660 62 -1942 226 -845 -642 410 -964 685 -950 1211 -198 -10 372 213 406 -968 789 -1033 -250 -313 -523 219 -440 473 -231 -253 -726 -398 212 -745 447 -1152 101 -337 -215 432 -174 229 -16 404 -425 996 217 682 1786 122 896 -183 201 271 35 518 1020 1701 2014 756 -411 -117 -132 1197 331 248 722 19 316 -896 -267 -727 -1194 293 -1267 -558 -2039 -1676 -1566 -1467 -1138 -588 958 -51 1111 -775 658 128 199 1061 -562 947 837 1013 -945 25 -1193 -902 678 -1305 -1695 -903 -671 222 -66 576 781 896 139 445 318 843 1091 543 1359 -1982 232 394 417 977 153 245 1753 639 -689 -337 -1548 -368 -742 -651 -113 -677 1509 -52 64 -22 -268 -243