File name: Pandora55s1_QueensNY_20230928_L2Fit_fout2c17d20220526p1-8.txt File generation date: 20230929T222629.6Z Data description: Level 2 spectral fitting results file Data file version: fout2c17p1-8 Local principal investigator: Lukas Valin Network principal investigator: Alexander Cede Instrument type: Pandora Instrument number: 55 Spectrometer number: 1 Processing software version used: BlickP v1.8.49 Instrument operation file used: Pandora55_OF_v5d20190723.txt Instrument calibration file used: Pandora55s1_CF_v17d20220526.txt Level 1 file used: Pandora55s1_QueensNY_20230928_L1_smca1c17d20220526p1-8.txt Full location name: Queens College NYSDEC (USEPA AQS ID 36-081-0124) Short location name: QueensNY Country of location: USA Location latitude [deg]: 40.7361 Location longitude [deg]: -73.8215 Location altitude [m]: 25 Local noon date: 20230928 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: 262 404 Nominal wavelengths inside fitting window [nm]: 305.1320 305.2720 305.4120 305.5520 305.6920 305.8320 305.9719 306.1119 306.2518 306.3917 306.5317 306.6716 306.8115 306.9514 307.0912 307.2311 307.3709 307.5108 307.6506 307.7905 307.9303 308.0701 308.2099 308.3496 308.4894 308.6292 308.7689 308.9087 309.0484 309.1881 309.3278 309.4675 309.6072 309.7469 309.8866 310.0262 310.1659 310.3055 310.4451 310.5847 310.7243 310.8639 311.0035 311.1431 311.2827 311.4222 311.5618 311.7013 311.8408 311.9803 312.1198 312.2593 312.3988 312.5383 312.6777 312.8172 312.9566 313.0960 313.2355 313.3749 313.5143 313.6536 313.7930 313.9324 314.0717 314.2111 314.3504 314.4897 314.6291 314.7684 314.9077 315.0469 315.1862 315.3255 315.4647 315.6040 315.7432 315.8824 316.0216 316.1608 316.3000 316.4392 316.5783 316.7175 316.8566 316.9958 317.1349 317.2740 317.4131 317.5522 317.6913 317.8303 317.9694 318.1085 318.2475 318.3865 318.5255 318.6645 318.8035 318.9425 319.0815 319.2205 319.3594 319.4984 319.6373 319.7762 319.9151 320.0540 320.1929 320.3318 320.4707 320.6095 320.7484 320.8872 321.0260 321.1649 321.3037 321.4425 321.5812 321.7200 321.8588 321.9975 322.1363 322.2750 322.4137 322.5524 322.6911 322.8298 322.9685 323.1072 323.2458 323.3845 323.5231 323.6617 323.8003 323.9389 324.0775 324.2161 324.3547 324.4932 324.6318 324.7703 324.9089 --------------------------------------------------------------------------------------- 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 20230928T175740.7Z 8671.748388 466 1 50.29 2 45.88 205.38 137.36 38.71 0.00 1 0.00 1 0 17 1.005e-02 6.533e-03 1.819e-03 1.196e-03 1.703174e-01 6.323220e-04 2.079018e-03 -6.000000e+00 5.968934e-03 227.81 -2.00 3.51 -2.00 315.17 0.00e+00 -7.744225e-05 5.721143e-05 2.096296e-05 -6.000000e+00 2.686081e-04 286.38 -2.00 20.00 -2.00 314.50 0.00e+00 1.434 0.000 302.00 328.00 3 4.5471e+00 1.6216e-03 1.6751e-02 -6.0000e+00 1.4799e-02 2.8716e-02 1.8007e-03 8.7620e-03 -6.0000e+00 1.6026e-02 5.2084e-03 1.0307e-03 4.2893e-03 -6.0000e+00 9.5696e-03 -4.3094e-03 4.5560e-04 1.1966e-03 -6.0000e+00 4.0626e-03 1.1466e-03 2.1698e-04 3.6802e-05 1.7479e-05 -6.0000e+00 3.0183e-04 1.2749e-02 0.0000e+00 11 9 0 11 4 0 99.99 20.45 0.18 0.01275 0.01438 53 10 1 5 0 948.700 26.00 20.01 20.45 24.26 36.0 1038.82 -3366 256 -5077 -1093 840 630 3429 1237 -555 -1288 128 -51 -2369 1976 -454 -1340 683 1610 1387 672 982 1627 -298 595 191 1272 -832 4 467 844 1971 1323 -2066 -1149 1147 -485 -530 482 -123 -1318 -845 -1573 -1047 -1155 -42 -1432 -513 -359 762 279 273 205 1117 865 693 1297 689 920 885 1363 1569 1169 144 -198 -769 -29 -171 933 788 -258 21 -92 -600 -177 -318 268 666 811 -558 -542 -194 -1332 -1133 -1187 -856 -726 -553 -600 -418 -234 85 93 792 1052 1115 635 594 -443 -172 -174 -547 23 769 160 -596 -333 -717 -172 318 247 -316 35 -222 -344 202 635 726 502 61 -490 -199 81 135 251 186 292 1032 730 167 598 213 257 231 416 -337 -148 99 -378 -508 -241 -348 -323 -58 -1092 82 -1657 -322 228 171 975 390 -187 -440 43 -18 -911 798 -199 -615 314 749 661 327 480 792 -150 309 102 680 -440 2 223 398 951 671 -1102 -601 573 -248 -292 297 -85 -991 -662 -1269 -851 -940 -34 -1152 -402 -274 570 210 212 165 920 733 595 1128 603 802 765 1183 1391 1084 140 -196 -752 -27 -151 809 715 -252 21 -96 -626 -180 -310 245 570 675 -480 -497 -189 -1377 -1215 -1315 -985 -868 -684 -754 -523 -286 100 102 823 1073 1165 700 694 -548 -218 -219 -666 27 874 187 -726 -420 -931 -229 431 345 -452 50 -312 -461 257 778 884 626 80 -664 -275 112 183 334 241 366 1251 881 207 766 280 343 312 571 -478 -216 147 -569 -770 -363 -516 -470 -84 SS 20230928T181539.2Z 8671.760871 479 1 32.69 2 47.50 211.14 135.14 43.99 0.00 1 0.00 1 0 25 1.149e-02 7.116e-03 2.643e-03 1.692e-03 1.756854e-01 8.936515e-04 2.138667e-03 -6.000000e+00 6.827780e-03 227.81 -2.00 3.51 -2.00 315.29 0.00e+00 -1.051985e-05 8.286847e-05 2.940837e-05 -6.000000e+00 3.072440e-04 286.38 -2.00 20.00 -2.00 314.61 0.00e+00 1.478 0.000 302.00 328.00 3 4.7075e+00 2.2944e-03 1.7259e-02 -6.0000e+00 1.6928e-02 3.1124e-02 2.5496e-03 9.0287e-03 -6.0000e+00 1.8332e-02 5.1976e-03 1.4593e-03 4.4188e-03 -6.0000e+00 1.0946e-02 -3.7983e-03 6.4748e-04 1.2339e-03 -6.0000e+00 4.6474e-03 9.2462e-04 3.5807e-04 5.2118e-05 1.7827e-05 -6.0000e+00 3.4524e-04 1.2414e-02 0.0000e+00 12 9 1 12 5 1 99.99 20.27 0.26 0.01241 0.01510 26 8 1 5 0 1210.290 26.00 20.01 20.27 25.29 36.2 1044.08 -2930 564 -4157 -2286 2525 -1429 4189 1228 -984 -1918 0 -1054 -3111 2406 -1151 -2131 -959 2869 1779 1733 983 1603 -6 972 231 2003 -555 -881 620 330 2281 607 -2192 -478 2275 -1081 -1114 958 -375 -1331 -655 -1398 -1125 -991 -99 -849 -284 -1002 614 -394 151 105 1693 831 945 1435 855 796 653 1560 1595 1156 103 -92 -1084 457 -229 779 738 -222 -181 -219 -493 -9 -204 -83 880 556 -398 -646 -55 -1659 -1226 -1006 -783 -1038 -612 -519 -455 -49 -142 203 946 989 1008 501 813 -300 -276 -350 -507 -152 1143 129 -676 -358 -654 -195 519 437 -415 152 -352 -243 76 385 638 620 83 -760 -217 -31 395 114 286 285 1074 615 417 435 264 377 150 473 -78 67 130 -371 -756 -223 -476 -141 -303 -905 173 -1287 -648 644 -377 1132 370 -318 -628 0 -357 -1157 935 -493 -954 -433 1283 820 813 466 759 -3 489 120 1036 -285 -432 289 153 1074 302 -1147 -244 1103 -544 -607 578 -255 -987 -506 -1109 -903 -794 -79 -671 -219 -760 454 -294 116 84 1380 698 805 1238 741 688 561 1341 1404 1065 99 -90 -1057 419 -201 673 667 -215 -184 -228 -513 -9 -198 -76 748 461 -341 -592 -54 -1722 -1313 -1113 -900 -1245 -757 -651 -570 -60 -167 223 983 1009 1056 554 951 -372 -352 -442 -619 -177 1298 151 -826 -453 -851 -260 706 613 -598 219 -496 -328 97 475 781 777 109 -1038 -301 -43 539 153 371 358 1307 747 519 560 349 506 203 653 -111 98 195 -563 -1157 -337 -711 -207 -440 SS 20230928T182713.4Z 8671.768905 488 1 14.42 2 48.69 214.68 133.60 47.19 0.00 1 0.00 1 0 25 8.951e-03 6.635e-03 2.325e-03 1.527e-03 1.776662e-01 8.078350e-04 2.169346e-03 -6.000000e+00 5.318210e-03 227.81 -2.00 3.51 -2.00 315.18 0.00e+00 1.476824e-04 7.329615e-05 4.665899e-05 -6.000000e+00 2.393081e-04 286.38 -2.00 20.00 -2.00 314.51 0.00e+00 1.512 0.000 302.00 328.00 3 3.7270e+00 2.0716e-03 1.7605e-02 -6.0000e+00 1.3186e-02 1.7240e-02 2.3004e-03 9.1434e-03 -6.0000e+00 1.4278e-02 4.8555e-03 1.3168e-03 4.4730e-03 -6.0000e+00 8.5261e-03 -3.0879e-03 5.8214e-04 1.2636e-03 -6.0000e+00 3.6200e-03 2.3856e-03 3.9262e-04 4.7055e-05 1.8372e-05 -6.0000e+00 2.6891e-04 1.2193e-02 0.0000e+00 12 9 1 12 5 1 99.98 20.16 0.23 0.01219 0.01514 25 17 1 5 0 554.170 26.00 19.98 20.16 25.69 36.1 1046.00 -1616 -1856 -3379 -1806 1367 476 2799 1072 -819 -1176 594 368 -903 1162 -1033 -911 9 1424 972 221 1180 1140 616 377 -62 -46 -591 -26 1564 1079 1688 1232 -1575 -602 1122 420 -310 859 -481 -1148 -645 -1484 -1317 -1088 -398 -1443 -991 -608 677 37 547 185 1177 847 729 1080 522 798 938 1397 1332 1068 543 -501 -481 68 -122 1353 818 -48 -134 -633 -473 -677 -68 -7 640 885 -163 -397 -267 -1595 -1260 -1031 -757 -587 -447 -496 -457 -16 36 121 1230 1364 1285 946 680 -354 -486 -650 -409 -97 897 103 -624 -216 -606 -80 -81 -66 -456 -223 -239 -41 535 602 751 723 14 -622 -381 -316 73 250 238 377 1567 845 307 464 5 161 219 375 -371 -26 -28 -228 -492 -120 -321 -326 -144 -514 -605 -1085 -534 370 130 802 338 -274 -398 198 128 -343 471 -454 -415 4 664 463 108 575 556 307 195 -33 -25 -310 -13 742 510 822 630 -843 -317 565 216 -172 529 -332 -855 -500 -1182 -1060 -875 -318 -1149 -774 -463 505 28 423 148 967 714 625 935 455 692 808 1206 1179 985 522 -495 -468 63 -108 1171 742 -46 -137 -659 -491 -688 -66 -6 550 742 -141 -366 -262 -1653 -1351 -1140 -867 -699 -549 -619 -569 -19 42 133 1279 1391 1348 1044 796 -437 -617 -820 -498 -114 1022 120 -760 -272 -784 -106 -110 -92 -650 -319 -334 -56 682 742 918 904 18 -844 -526 -437 99 333 309 474 1904 1027 382 597 6 215 296 517 -527 -38 -41 -344 -747 -181 -477 -474 -208