File name: Pandora159s1_Wakkerstroom_20201031_L2Fit_fout2c6d20191220p1-8.txt
File generation date: 20240926T165938.9Z
Data description: Level 2 spectral fitting results file
Data file version: fout2c6p1-8
Local principal investigator: B. Scholes
Network principal investigator: Alexander Cede
Instrument type: Pandora
Instrument number: 159
Spectrometer number: 1
Processing software version used: BlickP v1.8.60
Instrument operation file used: Pandora159_OF_v5d20191004.txt
Instrument calibration file used: Pandora159s1_CF_v6d20191220.txt
Level 1 file used: Pandora159s1_Wakkerstroom_20201031_L1_smca1c6d20191220p1-8.txt
Full location name: Global Change Institute, University of Witwatersrand field site in Wakkerstroom
Short location name: Wakkerstroom
Country of location: South Africa
Location latitude [deg]: -27.3493
Location longitude [deg]: 30.1438
Location altitude [m]: 1760
Local noon date: 20201031
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: Update session 4 based on session2 for SO2
First and last pixel inside fitting window: 222 364
Nominal wavelengths inside fitting window [nm]: 305.0490 305.1898 305.3306 305.4713 305.6121 305.7529 305.8936 306.0343 306.1750 306.3157 306.4564 306.5971 306.7378 306.8785 307.0191 307.1597 307.3004 307.4410 307.5816 307.7222 307.8628 308.0033 308.1439 308.2845 308.4250 308.5655 308.7060 308.8465 308.9870 309.1275 309.2680 309.4084 309.5489 309.6893 309.8297 309.9701 310.1105 310.2509 310.3913 310.5317 310.6720 310.8124 310.9527 311.0930 311.2333 311.3736 311.5139 311.6542 311.7944 311.9347 312.0749 312.2152 312.3554 312.4956 312.6358 312.7760 312.9161 313.0563 313.1964 313.3366 313.4767 313.6168 313.7569 313.8970 314.0371 314.1771 314.3172 314.4572 314.5973 314.7373 314.8773 315.0173 315.1573 315.2973 315.4372 315.5772 315.7171 315.8570 315.9969 316.1369 316.2767 316.4166 316.5565 316.6964 316.8362 316.9760 317.1159 317.2557 317.3955 317.5352 317.6750 317.8148 317.9545 318.0943 318.2340 318.3737 318.5134 318.6531 318.7928 318.9325 319.0721 319.2118 319.3514 319.4910 319.6307 319.7703 319.9098 320.0494 320.1890 320.3285 320.4681 320.6076 320.7471 320.8866 321.0261 321.1656 321.3051 321.4445 321.5840 321.7234 321.8628 322.0023 322.1417 322.2810 322.4204 322.5598 322.6991 322.8385 322.9778 323.1171 323.2564 323.3957 323.5350 323.6743 323.8135 323.9528 324.0920 324.2312 324.3704 324.5096 324.6488 324.7880 324.9272
---------------------------------------------------------------------------------------
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 20201031T102927.7Z 7609.437126 5 1 49.32 2 16.92 317.98 159.70 143.42 0.00 1 0.00 1 0 17 1.181e-02 7.725e-03 8.070e-03 4.688e-03 1.476044e-01 2.370040e-03 1.604646e-03 -6.000000e+00 6.734655e-03 228.16 -2.00 3.06 -2.00 315.69 0.00e+00 2.452300e-04 2.295737e-04 1.702201e-05 -6.000000e+00 2.931606e-04 279.24 -2.00 20.00 -2.00 315.08 0.00e+00 1.045 0.000 302.00 328.00 3 6.4502e+00 6.1380e-03 1.0855e-02 -6.0000e+00 1.6772e-02 5.4637e-02 6.8364e-03 6.6886e-03 -6.0000e+00 1.8210e-02 -4.0924e-03 3.9132e-03 3.3013e-03 -6.0000e+00 1.0822e-02 -9.8427e-04 1.7654e-03 9.4172e-04 -6.0000e+00 4.5931e-03 3.1145e-04 -2.3500e-04 1.3343e-04 1.1843e-05 -6.0000e+00 3.1554e-04 3.9544e-03 0.0000e+00 11 8 0 10 0 0 99.98 19.77 1.08 0.00395 0.00219 40 9 1 8 0 1232.640 15.22 19.77 22.54 21.98 21.98 14.3 934.94 3269 -1615 -736 -5471 630 -85 -384 1142 571 -490 -1211 3103 -1421 -463 -1292 -1157 518 -2481 -49 1522 1707 -547 332 -882 1539 20 -480 202 805 3841 3442 934 3506 -1960 1293 1066 -1678 -488 -289 -58 -1564 -863 -607 -592 -1461 -300 -1141 -1080 -565 211 1120 444 2587 -198 808 1230 479 712 592 110 1240 690 838 1160 -728 12 -18 816 1349 54 57 461 -14 -750 192 -675 295 -298 1099 -1071 -129 -926 -1595 -762 -1084 -906 -258 -276 -60 -628 811 417 172 904 1810 660 136 726 -1012 -100 -391 -734 389 185 -1093 -1002 -710 -518 -54 686 125 -282 -56 227 -95 58 328 -284 819 -358 -750 273 482 224 599 1761 862 1653 473 325 99 -552 278 486 68 -658 211 649 -347 -106 -249 -870 -746 862 -480 -218 -1571 151 -19 -89 290 162 -146 -352 858 -436 -157 -489 -468 210 -1049 -20 648 744 -240 145 -403 718 10 -232 90 332 1475 1346 394 1510 -906 544 430 -742 -245 -172 -39 -1173 -668 -477 -467 -1154 -232 -864 -789 -393 143 773 326 1948 -159 657 1008 401 590 487 90 1029 606 790 1145 -724 10 -15 634 1085 49 56 474 -15 -768 185 -589 229 -217 799 -870 -115 -906 -1667 -823 -1239 -1111 -333 -366 -80 -818 976 460 172 844 1724 691 155 905 -1352 -130 -486 -851 422 207 -1347 -1299 -956 -723 -77 1021 195 -444 -84 317 -122 69 385 -349 1068 -504 -1103 397 681 307 774 2158 1012 1889 573 414 132 -751 377 677 100 -1022 334 1047 -566 -172 -391 -1330 -1122
SS 20201031T110028.1Z 7609.458659 21 1 2.25 2 22.23 301.93 155.16 128.40 0.00 1 0.00 1 0 17 9.675e-03 8.166e-03 1.508e-03 1.203e-03 1.398993e-01 6.211106e-04 1.535986e-03 -6.000000e+00 5.515574e-03 228.16 -2.00 3.06 -2.00 314.01 0.00e+00 1.287581e-04 4.415238e-05 2.009349e-05 -6.000000e+00 2.401532e-04 279.24 -2.00 20.00 -2.00 313.56 0.00e+00 1.080 0.000 302.00 328.00 3 2.2397e+00 1.5820e-03 1.0853e-02 -6.0000e+00 1.3736e-02 1.6985e-02 1.7489e-03 6.3453e-03 -6.0000e+00 1.4915e-02 -8.4478e-04 1.0054e-03 3.1320e-03 -6.0000e+00 8.8639e-03 -6.8362e-04 4.3368e-04 8.8724e-04 -6.0000e+00 3.7610e-03 2.1043e-02 -2.9325e-04 3.3265e-05 1.2546e-05 -6.0000e+00 2.5848e-04 5.2716e-03 0.0000e+00 12 9 1 12 1 1 96.56 19.06 0.29 0.00527 0.00307 16 92 1 8 0 115.210 15.04 19.06 24.27 21.38 21.27 12.9 938.05 -1247 -1394 -1892 -3293 -581 950 1561 1730 1050 -504 -476 1654 357 379 378 -936 -89 231 260 755 -124 1239 1482 78 483 224 -591 -550 1171 2120 2265 1949 1847 -1460 -547 1195 -216 166 91 -1448 -1963 -1088 -1356 -1416 -444 -1404 -1019 -617 -647 581 588 464 812 944 493 1079 522 367 614 860 1270 1174 813 446 -845 -116 169 709 1764 298 -488 -263 -401 8 -344 216 558 828 1513 -133 -899 -1141 -2139 -814 -1089 -1129 -489 -660 -467 -543 132 218 685 1933 1650 741 843 319 -805 -681 -238 -225 941 735 -910 -946 -504 -780 -137 -96 -416 -400 -96 163 424 1216 870 635 539 -463 -806 -143 -91 169 903 709 1005 1881 587 507 115 0 78 433 -342 -585 -124 -34 -574 -352 -198 -531 -133 -569 -682 -933 -1525 -246 368 622 768 485 -246 -225 752 176 194 215 -542 -52 136 156 464 -78 758 907 50 324 146 -385 -354 721 1227 1326 1203 1135 -941 -325 711 -131 110 68 -1206 -1700 -970 -1236 -1297 -399 -1264 -887 -537 -539 468 488 407 708 855 458 1009 498 349 571 792 1195 1122 801 464 -878 -114 156 648 1585 292 -505 -285 -436 9 -358 211 507 693 1271 -119 -865 -1204 -2313 -905 -1241 -1329 -584 -816 -579 -654 152 237 723 1896 1671 787 941 377 -982 -854 -278 -253 1035 836 -1097 -1164 -619 -998 -181 -127 -560 -540 -128 208 509 1398 988 751 658 -593 -1078 -189 -117 217 1093 851 1180 2148 681 625 142 0 98 561 -454 -790 -172 -48 -798 -499 -272 -736 -179
SS 20201031T112657.8Z 7609.477057 32 1 17.51 2 27.46 293.00 150.52 119.20 0.00 1 0.00 1 0 25 9.490e-03 7.545e-03 1.339e-03 9.329e-04 1.469230e-01 4.759100e-04 1.600606e-03 -6.000000e+00 5.410697e-03 228.17 -2.00 3.06 -2.00 314.82 0.00e+00 8.036376e-05 3.894235e-05 1.354851e-05 -6.000000e+00 2.355562e-04 279.25 -2.00 20.00 -2.00 314.29 0.00e+00 1.127 0.000 302.00 328.00 3 3.5685e+00 1.2206e-03 1.1341e-02 -6.0000e+00 1.3475e-02 2.3021e-02 1.3557e-03 6.6506e-03 -6.0000e+00 1.4631e-02 -6.2591e-04 7.7548e-04 3.2803e-03 -6.0000e+00 8.6951e-03 -1.4645e-03 3.4101e-04 9.2335e-04 -6.0000e+00 3.6899e-03 5.2739e-03 -2.2878e-04 2.6022e-05 1.2390e-05 -6.0000e+00 2.5354e-04 4.5522e-03 0.0000e+00 12 9 1 12 1 1 100.00 18.99 0.27 0.00455 0.00283 61 40 1 8 0 266.630 15.02 18.99 23.70 21.38 21.27 12.7 934.86 -686 -1943 -1943 -3047 -417 991 1556 1811 1012 -501 -412 1819 157 177 280 -1147 -329 38 195 716 -372 700 1639 377 724 201 -598 -509 1106 2415 2027 2168 1901 -1372 -355 1021 -41 414 188 -1398 -1698 -1131 -1165 -1219 -770 -1302 -969 -572 -733 528 946 523 841 828 420 1145 594 260 737 871 1553 1137 918 478 -943 -84 229 676 1704 483 -371 -208 -457 -57 -415 425 602 936 1485 38 -907 -1135 -2073 -956 -975 -1153 -508 -747 -425 -494 257 363 602 1799 1397 736 816 454 -790 -563 -295 -240 383 931 -896 -753 -730 -588 -119 15 -83 -252 36 91 249 1070 782 584 281 -439 -657 -127 -44 41 879 571 965 1555 651 411 142 3 89 332 -289 -408 -112 76 -503 -296 -149 -520 -250 -248 -762 -765 -1124 -132 287 464 601 375 -196 -156 666 63 78 126 -561 -163 19 100 372 -198 365 848 205 411 118 -350 -278 546 1219 1042 1074 999 -757 -181 536 -21 248 124 -1087 -1438 -919 -966 -1016 -634 -1067 -828 -466 -573 399 734 431 725 742 358 983 520 227 629 793 1337 1046 909 500 -989 -80 210 586 1524 457 -388 -228 -503 -60 -418 382 518 741 1184 34 -847 -1162 -2275 -1081 -1098 -1398 -649 -928 -529 -636 305 417 618 1718 1380 795 904 554 -1037 -729 -358 -268 432 1075 -1070 -955 -965 -756 -158 21 -121 -370 51 126 310 1315 948 689 357 -606 -885 -170 -62 55 1101 681 1156 1818 805 505 182 4 117 447 -388 -579 -164 113 -763 -441 -215 -729 -341