File name: Pandora163s1_Seoul-SNU_20210412_L2Fit_fout2c14d20210114p1-8.txt File generation date: 20250326T105515.7Z Data description: Level 2 spectral fitting results file Data file version: fout2c14p1-8 Local principal investigator: Jinsoo Park Network principal investigator: Alexander Cede Instrument type: Pandora Instrument number: 163 Spectrometer number: 1 Processing software version used: BlickP v1.8.75 Instrument operation file used: Pandora163_OF_v8d20190605.txt Instrument calibration file used: Pandora163s1_CF_v14d20210114.txt Level 1 file used: Pandora163s1_Seoul-SNU_20210412_L1_smca1c14d20210114p1-8.txt Full location name: School of Earth and Environmental Sciences, Seoul National University Short location name: Seoul-SNU Country of location: Korea, Republic of Location latitude [deg]: 37.4580 Location longitude [deg]: 126.9510 Location altitude [m]: 116 Local noon date: 20210412 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, radiometric temperature correction, matrix method stray light correction (replaced by simple method) Data caveats: None First and last pixel inside fitting window: 225 367 Nominal wavelengths inside fitting window [nm]: 305.0608 305.2017 305.3426 305.4835 305.6244 305.7653 305.9062 306.0470 306.1879 306.3287 306.4695 306.6103 306.7511 306.8919 307.0327 307.1735 307.3143 307.4550 307.5958 307.7365 307.8772 308.0180 308.1587 308.2994 308.4400 308.5807 308.7214 308.8620 309.0027 309.1433 309.2839 309.4245 309.5651 309.7057 309.8463 309.9869 310.1274 310.2680 310.4085 310.5490 310.6895 310.8300 310.9705 311.1110 311.2515 311.3919 311.5324 311.6728 311.8133 311.9537 312.0941 312.2345 312.3748 312.5152 312.6556 312.7959 312.9363 313.0766 313.2169 313.3572 313.4975 313.6378 313.7780 313.9183 314.0585 314.1988 314.3390 314.4792 314.6194 314.7596 314.8998 315.0399 315.1801 315.3202 315.4603 315.6005 315.7406 315.8807 316.0207 316.1608 316.3009 316.4409 316.5810 316.7210 316.8610 317.0010 317.1410 317.2810 317.4209 317.5609 317.7008 317.8408 317.9807 318.1206 318.2605 318.4004 318.5402 318.6801 318.8199 318.9598 319.0996 319.2394 319.3792 319.5190 319.6588 319.7985 319.9383 320.0780 320.2177 320.3574 320.4971 320.6368 320.7765 320.9162 321.0558 321.1955 321.3351 321.4747 321.6143 321.7539 321.8935 322.0331 322.1726 322.3121 322.4517 322.5912 322.7307 322.8702 323.0097 323.1491 323.2886 323.4280 323.5675 323.7069 323.8463 323.9857 324.1251 324.2645 324.4038 324.5432 324.6825 324.8218 324.9611 --------------------------------------------------------------------------------------- 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 20210412T015408.8Z 7772.079269 207 1 51.70 2 36.41 135.84 40.74 137.55 0.00 1 0.00 1 0 17 1.089e-02 7.000e-03 8.426e-03 3.563e-03 1.917590e-01 1.761443e-03 2.289415e-03 -6.000000e+00 6.090997e-03 228.65 -2.00 3.44 -2.00 316.78 0.00e+00 -1.499253e-04 2.177025e-04 3.466876e-05 -6.000000e+00 2.675267e-04 281.74 -2.00 20.00 -2.00 316.06 0.00e+00 1.242 0.000 302.00 328.00 3 6.0014e+00 4.6259e-03 1.5734e-02 -6.0000e+00 1.5181e-02 1.8836e-02 5.1935e-03 9.7189e-03 -6.0000e+00 1.6502e-02 3.8545e-03 3.0302e-03 4.7781e-03 -6.0000e+00 9.8146e-03 -5.4140e-03 1.4408e-03 1.3389e-03 -6.0000e+00 4.1708e-03 3.1058e-04 -3.0933e-04 9.8264e-05 1.4708e-05 -6.0000e+00 2.7515e-04 -2.8818e-03 0.0000e+00 11 8 0 10 0 0 99.76 26.00 0.42 -0.00288 -0.00521 17 3 1 5 0 3040.750 26.00 20.01 22.32 25.05 31.4 1051.60 988 1166 -1393 -1169 -4145 -1995 833 1600 2883 -1022 -446 -979 162 -2769 629 -842 -228 -1645 689 741 -1173 491 81 4596 353 664 -1769 878 170 598 -10 2250 129 -1909 -1422 -2678 -1711 584 -980 -1570 -804 236 -657 95 -880 331 -411 -360 -1204 191 1416 703 1792 469 1882 645 286 1142 -939 348 394 825 957 1262 -187 21 -835 -1221 252 -268 1220 546 915 287 321 8 -509 -566 -125 -1525 -510 -716 -2447 -935 -811 -1587 653 -178 807 403 783 333 308 -949 -473 -337 898 1225 795 -50 -425 -821 -962 -864 -848 -669 -426 86 13 -123 809 367 153 -404 2 -494 -432 -801 569 -317 -265 164 144 363 209 446 481 573 1283 747 -101 424 50 -199 -6 534 224 112 102 -315 -443 -556 -709 163 209 -259 -203 -644 -286 122 263 525 -202 -86 -187 34 -665 164 -242 -67 -497 210 234 -383 156 26 1486 126 244 -661 299 53 181 -3 741 46 -717 -492 -904 -616 241 -489 -911 -503 152 -436 63 -583 216 -264 -222 -728 113 862 458 1213 336 1362 481 217 857 -702 256 296 663 825 1149 -171 17 -638 -903 195 -235 1148 539 900 274 288 6 -374 -391 -89 -1230 -453 -695 -2548 -992 -921 -1929 822 -232 1048 506 914 357 299 -898 -469 -362 1045 1559 1065 -66 -532 -955 -1069 -997 -1065 -881 -580 121 19 -190 1287 588 236 -583 3 -617 -541 -1047 786 -470 -404 246 212 511 282 572 587 697 1642 1021 -144 612 73 -300 -9 886 379 195 178 -541 -740 -901 -1132 SS 20210412T020738.6Z 7772.088641 216 1 51.17 2 34.62 140.69 39.00 141.93 0.00 1 0.00 1 0 17 8.171e-03 5.693e-03 4.373e-03 1.977e-03 1.807591e-01 9.775876e-04 2.153323e-03 -6.000000e+00 4.570051e-03 228.65 -2.00 3.44 -2.00 316.60 0.00e+00 4.277136e-04 1.159867e-04 3.047288e-05 -6.000000e+00 2.007290e-04 281.74 -2.00 20.00 -2.00 315.88 0.00e+00 1.214 0.000 302.00 328.00 3 5.0297e+00 2.5613e-03 1.5156e-02 -6.0000e+00 1.1391e-02 6.3901e-03 2.8714e-03 9.1568e-03 -6.0000e+00 1.2382e-02 5.2512e-03 1.6664e-03 4.4924e-03 -6.0000e+00 7.3639e-03 -5.0330e-03 7.8479e-04 1.2855e-03 -6.0000e+00 3.1293e-03 8.6375e-04 -1.3823e-04 5.4404e-05 1.3923e-05 -6.0000e+00 2.0645e-04 -2.7152e-03 0.0000e+00 10 0 0 10 0 0 100.00 26.00 0.31 -0.00272 -0.00375 122 22 1 5 0 415.950 26.00 20.01 22.32 25.33 31.5 1053.32 294 -1339 -2556 -2068 -2017 -990 564 2467 -126 935 -256 677 -741 -322 415 246 515 -157 424 58 874 -306 827 501 901 436 81 547 -1489 564 703 1855 1776 -574 -403 -1951 -1305 -430 -416 -826 -1293 -1102 -1060 -935 -222 -434 -1207 -793 -814 -501 304 1164 1440 663 795 1027 576 755 251 800 1291 1036 1064 1472 473 -403 -197 -522 99 734 541 674 137 395 325 -130 -24 -747 -323 -773 -626 -1184 -2288 -1430 -1189 -781 -25 58 176 485 225 339 343 180 403 331 449 575 368 -112 -379 -1003 -305 -168 -621 -875 -686 -360 -237 384 282 435 263 -199 -114 -197 240 83 62 0 -134 -31 221 302 424 376 326 1001 323 665 59 83 93 -161 182 -98 157 -52 81 -286 -138 -495 -454 54 -271 -530 -395 -339 -154 90 438 -26 197 -53 136 -166 -81 117 75 160 -50 139 20 300 -104 280 178 337 169 31 196 -505 180 232 645 664 -225 -144 -681 -490 -187 -216 -493 -840 -732 -721 -638 -152 -295 -805 -506 -507 -307 193 772 1000 485 593 775 445 578 188 594 982 848 930 1356 434 -342 -153 -389 78 647 516 669 137 380 296 -108 -18 -524 -234 -623 -563 -1161 -2372 -1529 -1360 -946 -32 76 229 612 265 365 336 169 397 354 523 729 496 -148 -473 -1168 -340 -193 -777 -1160 -933 -508 -344 587 446 689 402 -284 -149 -244 296 107 86 0 -203 -47 325 424 571 481 397 1206 418 903 83 119 134 -241 284 -162 262 -90 140 -484 -228 -792 -720 SS 20210412T021309.4Z 7772.092470 219 1 50.88 2 33.94 142.79 38.33 143.82 0.00 1 0.00 1 0 17 8.567e-03 5.894e-03 3.058e-03 1.469e-03 1.808032e-01 7.273231e-04 2.174499e-03 -6.000000e+00 4.791505e-03 228.65 -2.00 3.44 -2.00 316.38 0.00e+00 3.034363e-04 8.236966e-05 2.886222e-05 -6.000000e+00 2.104512e-04 281.74 -2.00 20.00 -2.00 315.68 0.00e+00 1.205 0.000 302.00 328.00 3 4.8224e+00 1.8999e-03 1.5107e-02 -6.0000e+00 1.1942e-02 7.4096e-03 2.1281e-03 9.1940e-03 -6.0000e+00 1.2982e-02 1.7097e-03 1.2302e-03 4.5181e-03 -6.0000e+00 7.7207e-03 -3.4234e-03 5.7434e-04 1.2753e-03 -6.0000e+00 3.2810e-03 1.0741e-03 -2.5769e-04 4.0179e-05 1.4268e-05 -6.0000e+00 2.1645e-04 -2.8818e-03 0.0000e+00 10 0 0 10 0 0 99.99 26.00 0.31 -0.00288 -0.00482 63 11 1 5 0 807.490 26.00 20.00 22.32 26.05 31.5 1057.62 643 -1652 -1594 12 -4330 -1383 486 398 1267 750 -685 -542 -400 -1082 385 -596 302 -174 25 786 716 853 -14 1457 -256 485 -86 798 -389 1377 773 1220 2561 -901 -572 -2713 -1148 -1159 -473 -409 -1124 -414 -807 -869 -863 -316 -1005 -768 -781 136 473 735 1391 523 919 702 670 790 144 616 1020 877 1093 1468 521 -233 -408 -325 567 921 299 108 325 468 259 -24 27 -772 -216 -757 -1174 -916 -1640 -1550 -1266 -850 -294 43 106 357 448 217 442 91 730 129 745 592 527 -300 -582 -730 -695 -144 -631 -814 -922 -290 -314 320 648 277 106 -70 -120 -292 123 271 236 -138 -305 -42 -29 464 527 576 621 757 956 797 15 281 10 -128 -136 54 -18 53 -78 -120 -400 -335 -600 126 -360 -351 2 -798 -233 84 78 275 169 -152 -119 -96 -294 115 -195 100 -59 9 280 261 304 -5 545 -102 198 -35 304 -138 461 269 451 997 -371 -217 -1014 -454 -533 -255 -253 -755 -284 -570 -617 -612 -222 -695 -509 -502 86 309 509 1004 395 707 545 530 621 112 473 803 737 978 1376 487 -201 -325 -249 454 825 291 109 328 458 240 -20 21 -562 -161 -628 -1074 -906 -1718 -1677 -1446 -1029 -372 55 137 451 527 236 437 87 724 139 867 753 703 -395 -729 -847 -775 -167 -787 -1063 -1249 -406 -453 482 1008 433 160 -100 -157 -363 153 347 322 -200 -453 -62 -43 646 703 739 758 919 1219 1071 21 397 14 -189 -211 88 -29 88 -131 -199 -646 -527 -933