File name: Pandora138s1_Rome-SAP_20201012_L2Fit_fout2c1d20190827p1-8.txt File generation date: 20250203T181932.3Z Data description: Level 2 spectral fitting results file Data file version: fout2c1p1-8 Local principal investigator: Stefano Casadio Network principal investigator: Alexander Cede Instrument type: Pandora Instrument number: 138 Spectrometer number: 1 Processing software version used: BlickP v1.8.70 Instrument operation file used: Pandora138_OF_v6d20190508.txt Instrument calibration file used: Pandora138s1_CF_v1d20190827.txt Level 1 file used: Pandora138s1_Rome-SAP_20201012_L1_smca1c1d20190827p1-8.txt Full location name: Department of Physics, University of Rome La Sapienza Short location name: Rome-SAP Country of location: Italy Location latitude [deg]: 41.9017 Location longitude [deg]: 12.5158 Location altitude [m]: 75 Local noon date: 20201012 Notes on s-number (L1 configuration): Corrections NOT applied although requested by the s-number are latency correction, matrix method stray light correction (replaced by simple method) Data caveats: None First and last pixel inside fitting window: 210 352 Nominal wavelengths inside fitting window [nm]: 305.0730 305.2140 305.3550 305.4960 305.6370 305.7780 305.9189 306.0599 306.2009 306.3418 306.4827 306.6236 306.7645 306.9054 307.0463 307.1872 307.3280 307.4689 307.6097 307.7505 307.8914 308.0322 308.1729 308.3137 308.4545 308.5952 308.7360 308.8767 309.0175 309.1582 309.2989 309.4396 309.5802 309.7209 309.8615 310.0022 310.1428 310.2834 310.4241 310.5647 310.7052 310.8458 310.9864 311.1269 311.2675 311.4080 311.5485 311.6890 311.8295 311.9700 312.1105 312.2509 312.3914 312.5318 312.6723 312.8127 312.9531 313.0935 313.2338 313.3742 313.5146 313.6549 313.7953 313.9356 314.0759 314.2162 314.3565 314.4968 314.6370 314.7773 314.9175 315.0577 315.1980 315.3382 315.4784 315.6186 315.7587 315.8989 316.0390 316.1792 316.3193 316.4594 316.5995 316.7396 316.8797 317.0198 317.1598 317.2999 317.4399 317.5799 317.7199 317.8599 317.9999 318.1399 318.2799 318.4198 318.5598 318.6997 318.8396 318.9795 319.1194 319.2593 319.3992 319.5391 319.6789 319.8187 319.9586 320.0984 320.2382 320.3780 320.5178 320.6575 320.7973 320.9370 321.0768 321.2165 321.3562 321.4959 321.6356 321.7753 321.9149 322.0546 322.1942 322.3338 322.4735 322.6131 322.7527 322.8922 323.0318 323.1714 323.3109 323.4505 323.5900 323.7295 323.8690 324.0085 324.1480 324.2874 324.4269 324.5663 324.7057 324.8452 324.9846 --------------------------------------------------------------------------------------- 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 20201012T074016.1Z 7590.319631 266 1 50.17 2 66.72 125.44 23.63 198.28 0.00 1 0.00 1 0 17 2.617e-02 6.645e-03 8.399e-03 1.647e-03 3.498001e-01 8.142169e-04 4.166486e-03 -6.000000e+00 1.474422e-02 226.03 -2.00 3.53 -2.00 318.09 0.00e+00 -1.855946e-05 1.343765e-04 1.005238e-04 -6.000000e+00 6.003402e-04 283.44 -2.00 20.00 -2.00 317.75 0.00e+00 2.510 0.000 302.00 328.00 3 4.5092e+00 2.1659e-03 3.0637e-02 -6.0000e+00 3.6752e-02 1.6845e-02 2.4432e-03 1.7811e-02 -6.0000e+00 4.0213e-02 1.1437e-02 1.6216e-03 8.8910e-03 -6.0000e+00 2.3678e-02 -1.0398e-02 8.4574e-04 2.5573e-03 -6.0000e+00 9.8016e-03 3.1338e-04 -1.4271e-04 4.6885e-05 2.7607e-05 -6.0000e+00 7.0930e-04 1.1474e-02 0.0000e+00 11 8 0 10 0 0 99.28 12.77 0.16 0.01147 0.01040 38 8 1 8 0 1319.980 15.00 12.77 20.84 20.77 21.38 28.0 952.04 -8054 -10995 -14042 -1205 -13444 3010 -4556 -2501 1733 -5946 -3015 1879 -5134 -2466 -1130 -3827 -988 484 7041 1545 99 1158 -1196 2440 439 337 -224 3132 -2882 1391 4736 1012 -295 -3548 329 -1738 1667 -1941 -741 -474 -1383 -907 83 -700 30 -279 -2162 -970 -1277 -427 7616 444 1295 -15 648 507 984 358 185 285 1401 768 1311 1209 -574 243 -142 -380 1266 84 -403 -370 313 -73 3 -87 195 -524 -221 -137 -1233 -1256 -1823 -1337 -1142 -523 58 -29 555 364 982 191 434 739 323 241 753 703 227 -512 -756 -368 -136 -303 -1173 -923 -631 -433 70 390 336 129 342 76 212 438 120 71 -21 -250 -503 48 262 350 413 -30 688 997 444 330 -299 -108 273 78 -168 -399 161 -84 -172 26 -203 -192 -201 -486 -751 -1012 -73 -827 158 -279 -168 125 -485 -241 154 -515 -281 -146 -560 -149 76 1089 261 17 200 -220 464 90 72 -48 614 -569 271 976 234 -75 -946 80 -437 447 -618 -274 -198 -624 -424 40 -338 15 -135 -1052 -463 -613 -209 3715 244 741 -9 401 321 630 229 116 180 905 529 963 934 -445 178 -101 -271 954 70 -356 -334 279 -63 3 -66 138 -372 -166 -113 -1118 -1220 -1852 -1405 -1257 -598 69 -35 659 419 1067 198 429 732 339 274 919 915 304 -677 -957 -439 -158 -363 -1490 -1208 -852 -602 101 588 524 202 524 111 297 597 166 101 -31 -384 -776 73 391 507 580 -41 920 1361 642 498 -468 -172 441 131 -293 -711 290 -152 -313 45 -353 -326 -340 SS 20201012T081518.5Z 7590.343964 286 1 42.99 2 61.68 133.31 26.64 216.02 0.00 1 0.00 1 0 17 2.782e-02 9.606e-03 7.387e-03 2.398e-03 2.932808e-01 1.187760e-03 3.521970e-03 -6.000000e+00 1.567882e-02 226.03 -2.00 3.53 -2.00 317.24 0.00e+00 1.045273e-03 1.559712e-04 6.341912e-05 -6.000000e+00 6.383063e-04 283.44 -2.00 20.00 -2.00 316.91 0.00e+00 2.098 0.000 302.00 328.00 3 5.4763e+00 3.1183e-03 2.5615e-02 -6.0000e+00 3.9081e-02 5.0255e-02 3.5035e-03 1.4893e-02 -6.0000e+00 4.2759e-02 -2.6769e-03 2.1327e-03 7.4069e-03 -6.0000e+00 2.5177e-02 -4.4686e-03 1.0526e-03 2.1140e-03 -6.0000e+00 1.0425e-02 1.8951e-04 -3.7743e-05 6.8240e-05 2.3559e-05 -6.0000e+00 7.5406e-04 9.7784e-03 0.0000e+00 12 9 9 12 1 1 70.24 12.82 0.18 0.00978 0.00949 10 3 1 8 0 3892.440 15.01 12.82 21.41 20.77 21.38 28.3 956.01 -5507 -12984 -6704 -2080 -9195 -5944 -6866 -4491 1032 -9231 -7948 1431 -3967 -207 1963 -3349 15 1476 11347 2724 4659 2359 2337 5185 715 -1057 -223 3945 -3631 2917 4089 2587 -848 -2604 2281 483 1247 -1174 432 -389 -1138 -541 567 -213 1157 -442 -1535 -645 385 -134 -10183 1127 1273 923 1225 513 363 544 -788 510 1540 1334 1885 1601 152 478 -203 -49 1491 415 257 23 400 545 -650 -898 -875 -1188 -378 -131 -1510 -1457 -1706 -957 -780 -664 133 195 1196 546 1478 -350 -669 -119 10 471 600 1140 771 -115 -640 -626 -291 -642 -839 -1211 -840 -323 431 475 546 123 275 -363 -232 68 -177 109 -153 -11 -81 351 313 122 685 -195 489 958 70 647 -183 -69 -44 148 -4 -406 -20 345 61 -143 -330 -317 -317 -608 -1671 -829 -230 -986 -609 -757 -538 129 -1321 -1114 194 -641 -37 390 -761 3 350 2567 687 1170 601 614 1395 208 -323 -67 1090 -1019 784 1165 801 -290 -914 727 156 438 -479 201 -204 -636 -311 331 -126 678 -262 -908 -374 219 -78 -6712 720 846 640 866 370 265 394 -569 365 1130 1035 1547 1371 129 384 -158 -38 1214 371 242 22 381 505 -573 -737 -676 -904 -302 -114 -1433 -1473 -1792 -1034 -883 -784 161 240 1456 647 1652 -373 -682 -120 11 537 733 1473 1023 -150 -805 -748 -337 -773 -1064 -1587 -1130 -446 614 706 833 188 410 -519 -317 90 -237 150 -221 -16 -121 517 455 172 935 -260 634 1262 98 933 -273 -105 -67 233 -6 -682 -33 593 104 -240 -546 -513 -510 SS 20201012T082107.2Z 7590.348000 289 1 26.47 2 60.91 134.69 27.29 218.60 0.00 1 0.00 1 0 17 2.059e-02 7.637e-03 9.603e-03 3.104e-03 2.883440e-01 1.538760e-03 3.459234e-03 -6.000000e+00 1.160360e-02 226.03 -2.00 3.53 -2.00 317.20 0.00e+00 -1.675728e-04 1.997432e-04 6.162696e-05 -6.000000e+00 4.723216e-04 283.44 -2.00 20.00 -2.00 316.86 0.00e+00 2.047 0.000 302.00 328.00 3 5.5119e+00 4.0394e-03 2.5053e-02 -6.0000e+00 2.8924e-02 3.4492e-02 4.5382e-03 1.4625e-02 -6.0000e+00 3.1643e-02 1.1767e-02 2.7579e-03 7.2765e-03 -6.0000e+00 1.8632e-02 -1.0420e-02 1.3604e-03 2.0780e-03 -6.0000e+00 7.7168e-03 1.9489e-04 1.1546e-04 8.8361e-05 2.3241e-05 -6.0000e+00 5.5793e-04 8.0493e-03 0.0000e+00 12 9 1 12 41 1 49.61 12.84 0.16 0.00805 0.00892 6 3 1 8 0 3757.070 15.00 12.84 21.62 20.77 21.38 28.4 957.40 -2077 -11190 -6285 -2471 -4100 -1541 366 -5098 8148 -5923 -711 -1964 -3181 2527 317 -1725 -1553 -139 7354 1144 497 3382 1333 3059 470 392 -709 -1485 -2403 2494 4641 1574 2856 -4491 965 -1320 1374 -1012 -1008 -203 -1180 -486 478 -204 -208 -1494 -2841 -1847 -817 -1766 2399 719 1554 1201 1078 952 1032 1033 -673 37 1221 1338 1103 1812 51 352 -1046 -723 1251 -659 229 -306 118 430 -741 127 -1749 -1137 -691 -432 -1670 -1533 -924 -1376 -1263 -440 624 126 1155 727 991 335 -14 209 -221 439 620 1003 517 -146 -341 -601 -691 -459 -714 -1151 -864 -334 674 799 533 239 306 -45 -56 -121 -407 364 -236 -121 -695 328 227 472 297 61 228 1043 94 212 -164 -371 139 164 -226 -602 205 510 247 -294 -414 -181 -110 -224 -1418 -778 -276 -421 -153 38 -627 981 -848 -97 -283 -522 448 65 -390 -365 -34 1730 294 130 864 359 855 140 121 -220 -440 -680 683 1328 494 953 -1610 313 -438 487 -416 -480 -108 -670 -284 283 -122 -124 -901 -1707 -1087 -472 -1045 1450 462 1036 836 767 689 759 754 -490 27 904 1046 914 1558 44 285 -822 -566 1022 -595 217 -297 114 400 -659 104 -1368 -869 -554 -378 -1587 -1552 -968 -1495 -1438 -519 755 155 1410 862 1113 357 -14 213 -236 500 756 1296 688 -191 -428 -717 -804 -551 -905 -1508 -1164 -461 957 1183 814 366 455 -65 -77 -160 -543 500 -339 -179 -1040 482 329 663 405 81 296 1369 130 304 -243 -561 213 257 -369 -1007 347 873 421 -494 -682 -292 -176