File name: Pandora186s1_MadisonCT_20240401_L2Fit_fout2c2d20210922p1-8.txt File generation date: 20240514T232915.6Z Data description: Level 2 spectral fitting results file Data file version: fout2c2p1-8 Local principal investigator: Lukas Valin Network principal investigator: Alexander Cede Instrument type: Pandora Instrument number: 186 Spectrometer number: 1 Processing software version used: BlickP v1.8.59 Instrument operation file used: Pandora186_OF_v4d20210803.txt Instrument calibration file used: Pandora186s1_CF_v2d20210922.txt Level 1 file used: Pandora186s1_MadisonCT_20240401_L1_smca1c2d20210922p1-8.txt Full location name: Madison - Hammonasset State Park CTDEEP (USEPA AQS ID 09-009-9002) Short location name: MadisonCT Country of location: USA Location latitude [deg]: 41.2568 Location longitude [deg]: -72.5533 Location altitude [m]: 3 Local noon date: 20240401 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: None First and last pixel inside fitting window: 197 339 Nominal wavelengths inside fitting window [nm]: 305.0092 305.1504 305.2917 305.4329 305.5741 305.7153 305.8565 305.9977 306.1389 306.2800 306.4212 306.5623 306.7034 306.8445 306.9856 307.1267 307.2678 307.4089 307.5499 307.6909 307.8320 307.9730 308.1140 308.2550 308.3960 308.5369 308.6779 308.8189 308.9598 309.1007 309.2416 309.3825 309.5234 309.6643 309.8052 309.9460 310.0869 310.2277 310.3685 310.5093 310.6501 310.7909 310.9317 311.0724 311.2132 311.3539 311.4946 311.6354 311.7761 311.9167 312.0574 312.1981 312.3387 312.4794 312.6200 312.7606 312.9012 313.0418 313.1824 313.3230 313.4636 313.6041 313.7447 313.8852 314.0257 314.1662 314.3067 314.4472 314.5876 314.7281 314.8685 315.0090 315.1494 315.2898 315.4302 315.5706 315.7109 315.8513 315.9917 316.1320 316.2723 316.4126 316.5529 316.6932 316.8335 316.9738 317.1140 317.2543 317.3945 317.5347 317.6749 317.8151 317.9553 318.0955 318.2356 318.3758 318.5159 318.6560 318.7961 318.9362 319.0763 319.2164 319.3565 319.4965 319.6366 319.7766 319.9166 320.0566 320.1966 320.3366 320.4765 320.6165 320.7565 320.8964 321.0363 321.1762 321.3161 321.4560 321.5959 321.7357 321.8756 322.0154 322.1552 322.2950 322.4348 322.5746 322.7144 322.8542 322.9939 323.1337 323.2734 323.4131 323.5528 323.6925 323.8322 323.9719 324.1115 324.2512 324.3908 324.5304 324.6700 324.8096 324.9492 --------------------------------------------------------------------------------------- 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 20240401T190843.8Z 8857.797730 470 1 34.67 2 47.16 228.96 135.20 268.54 0.00 1 0.00 1 0 17 1.112e-02 5.261e-03 5.927e-03 2.013e-03 2.681769e-01 1.038797e-03 3.233721e-03 -6.000000e+00 6.562306e-03 227.55 -2.00 3.52 -2.00 317.07 0.00e+00 3.209620e-04 1.413014e-04 9.768002e-05 -6.000000e+00 2.860875e-04 278.77 -2.00 20.00 -2.00 316.42 0.00e+00 1.468 0.000 302.00 328.00 3 5.2876e+00 2.7283e-03 2.0294e-02 -6.0000e+00 1.6302e-02 4.2860e-02 3.0623e-03 1.3749e-02 -6.0000e+00 1.7696e-02 -7.0579e-03 1.8184e-03 6.7902e-03 -6.0000e+00 1.0529e-02 -9.1930e-04 8.7569e-04 1.9098e-03 -6.0000e+00 4.3973e-03 4.2604e-04 -5.4510e-05 6.2019e-05 2.4060e-05 -6.0000e+00 3.2981e-04 1.0207e-03 5.2136e-05 12 1 1 12 1 1 99.53 19.83 0.23 0.00102 0.00061 12 4 1 5 0 2666.040 22.61 15.00 19.83 14.62 24.9 1014.21 -5154 -1495 -4532 -6417 -546 -2420 -1255 -853 859 -586 -1670 -89 1618 1164 -255 -747 -701 710 603 1283 1250 730 322 1269 -498 177 -222 490 603 339 1841 1786 1236 -256 -1570 -1266 -1358 -2018 -380 -815 -520 -1301 -470 -735 -352 -655 -723 -430 553 -296 1015 689 611 910 436 934 1036 752 309 637 699 796 1075 1213 673 43 -175 522 817 1 306 299 -51 -100 73 -411 -148 -672 -224 -796 -704 -1545 -1522 -1619 -861 -701 -28 -16 69 342 119 48 189 26 804 461 219 748 -27 -354 -359 -202 -822 -74 -403 -625 -559 -402 -81 477 583 218 164 -303 218 141 526 604 258 2 -60 -295 10 25 349 443 335 938 522 504 -128 -303 -231 -42 -312 -170 -169 395 78 -205 -51 -315 -351 -600 -187 -610 -857 -63 -268 -140 -104 114 -84 -246 -13 255 208 -52 -171 -170 176 155 340 339 199 89 363 -152 56 -71 150 174 96 528 549 409 -89 -538 -421 -476 -800 -172 -423 -294 -777 -286 -454 -217 -404 -442 -258 324 -175 613 441 410 634 315 683 767 559 228 466 520 621 888 1051 591 36 -140 407 656 1 287 291 -50 -96 66 -351 -117 -512 -174 -670 -646 -1534 -1591 -1761 -975 -831 -35 -21 87 420 139 53 196 26 835 512 262 954 -35 -473 -461 -246 -970 -88 -503 -814 -752 -558 -116 707 893 338 250 -443 302 188 692 814 363 4 -89 -439 15 35 479 588 431 1207 700 705 -186 -450 -348 -65 -499 -281 -284 668 132 -345 -84 -507 -558 SS 20240401T193115.6Z 8857.813375 479 1 4.92 2 50.49 234.96 139.21 272.15 0.00 1 0.00 1 0 17 1.429e-02 7.431e-03 1.441e-02 5.816e-03 2.864577e-01 3.017259e-03 3.438847e-03 -6.000000e+00 8.434621e-03 227.55 -2.00 3.52 -2.00 316.60 0.00e+00 -1.190213e-04 3.662821e-04 9.383178e-05 -6.000000e+00 3.676004e-04 278.77 -2.00 20.00 -2.00 315.96 0.00e+00 1.569 0.000 302.00 328.00 3 5.0428e+00 7.8689e-03 2.1626e-02 -6.0000e+00 2.0954e-02 3.7626e-02 8.8093e-03 1.4596e-02 -6.0000e+00 2.2741e-02 -4.1958e-03 5.1434e-03 7.2026e-03 -6.0000e+00 1.3531e-02 -2.3925e-03 2.4207e-03 2.0149e-03 -6.0000e+00 5.6546e-03 4.7957e-04 5.9925e-04 1.7686e-04 2.5930e-05 -6.0000e+00 4.2368e-04 -1.6495e-03 9.2918e-05 12 9 1 12 1 1 999.00 19.85 0.49 -0.00165 0.00285 1 4 1 5 0 2457.700 22.61 15.00 19.85 15.18 24.8 1017.94 -2881 -1194 -2613 -4971 -4938 -4202 -3628 2189 548 -2830 358 -2353 189 2491 -2597 -1342 1760 2059 -1585 2266 2575 1096 -619 827 -612 1652 1632 -2326 3384 1622 665 3222 1822 -2153 -1932 632 -2026 -1296 -263 -1175 -930 -1069 211 -1812 281 -615 -2039 -1750 -2302 988 943 2378 1036 641 1608 449 1082 1016 434 -355 465 1704 894 678 740 456 -573 482 1038 1530 138 -187 -499 98 363 336 -234 -900 -385 -329 133 -1143 -2131 -626 -1407 -1083 -140 -702 139 634 148 258 -45 687 -201 536 565 -247 204 -804 582 -1204 418 596 -831 -1021 -621 -812 -760 417 325 293 35 89 -97 -350 508 546 213 388 -186 400 212 372 514 752 161 1163 1050 620 83 33 44 -714 -110 169 -170 -378 288 -571 -541 -223 -197 -396 -179 -414 -787 -725 -578 -507 317 89 -505 63 -435 37 534 -651 -366 495 598 -490 700 811 351 -203 279 -219 600 599 -854 1133 539 229 1160 705 -886 -779 245 -836 -589 -135 -676 -577 -690 138 -1216 186 -409 -1359 -1149 -1502 633 622 1627 745 478 1224 350 851 800 339 -278 368 1390 770 609 671 400 -486 399 881 1384 134 -187 -502 96 342 299 -197 -735 -319 -291 126 -1162 -2277 -685 -1606 -1282 -170 -872 172 769 174 287 -48 713 -215 604 677 -314 268 -1061 739 -1477 492 708 -1034 -1317 -822 -1104 -1065 601 482 439 52 127 -132 -463 665 731 296 556 -271 578 302 517 694 989 208 1493 1393 856 118 48 65 -1079 -171 269 -276 -618 470 -924 -859 -346 -303 SS 20240401T193720.4Z 8857.817597 483 1 51.13 2 51.44 236.49 140.29 273.18 0.00 1 0.00 1 0 17 1.695e-02 5.064e-03 1.179e-02 2.654e-03 2.950756e-01 1.364590e-03 3.509855e-03 -6.000000e+00 1.000444e-02 227.55 -2.00 3.52 -2.00 317.99 0.00e+00 -1.581704e-04 2.352812e-04 1.220993e-04 -6.000000e+00 4.361246e-04 278.78 -2.00 20.00 -2.00 317.32 0.00e+00 1.601 0.000 302.00 328.00 3 5.5982e+00 3.6438e-03 2.2220e-02 -6.0000e+00 2.4854e-02 4.5368e-02 4.1218e-03 1.5135e-02 -6.0000e+00 2.6977e-02 -3.9780e-03 2.5962e-03 7.4752e-03 -6.0000e+00 1.6052e-02 -1.8091e-03 1.3070e-03 2.1017e-03 -6.0000e+00 6.7045e-03 2.6160e-04 6.9416e-06 8.3023e-05 2.5654e-05 -6.0000e+00 5.0274e-04 5.9429e-04 6.5677e-05 10 0 0 10 0 0 99.96 19.85 0.44 0.00059 0.00065 26 5 1 5 0 1966.000 22.61 15.00 19.85 15.23 25.0 1019.93 1201 3795 -3306 -1028 -4699 -1569 -10299 -252 -7596 1180 -5024 -1441 -5132 4662 2903 -1586 -2165 943 3049 114 2379 -104 -51 -1197 299 -1239 895 -911 1053 524 1926 -1718 921 -1965 -883 -972 -2464 -1398 -324 -535 -1115 29 -1721 -259 -328 -434 22 -340 -41 547 336 1207 999 870 358 381 1535 722 983 183 1655 690 1689 1008 179 -328 -347 -303 191 644 371 389 2 131 -183 -188 -1031 -101 -1271 -483 -1534 -642 -1582 -1264 -927 -319 -117 473 625 280 327 277 158 -146 468 -317 692 436 536 -477 -158 -806 -404 -190 -647 -681 -598 -434 -226 653 298 366 -111 -91 -69 59 25 289 12 133 61 29 215 72 494 292 256 147 427 304 443 -303 -280 47 -100 -97 151 170 91 68 -159 -397 -488 80 278 -272 -81 -353 -107 -784 -19 -693 107 -493 -140 -564 532 386 -247 -362 160 530 21 443 -20 -10 -247 64 -283 203 -203 217 106 399 -399 226 -516 -224 -240 -653 -418 -113 -218 -507 14 -859 -130 -164 -217 11 -165 -19 263 168 639 561 514 221 240 976 463 621 116 1052 468 1220 777 141 -248 -247 -211 137 508 320 351 2 116 -153 -144 -733 -68 -888 -368 -1310 -594 -1573 -1314 -1014 -367 -140 581 770 337 371 293 156 -141 467 -345 816 559 722 -646 -205 -983 -469 -223 -805 -891 -813 -612 -332 1002 478 595 -177 -138 -99 80 34 402 17 203 95 45 325 106 694 395 334 193 587 442 672 -474 -447 77 -171 -172 274 312 167 122 -281 -681 -826