File name: Pandora63s1_LaPorteTX_20231201_L2Fit_fout2c7d20230718p1-8.txt File generation date: 20231202T063709.6Z Data description: Level 2 spectral fitting results file Data file version: fout2c7p1-8 Local principal investigator: Tom Hanisco Network principal investigator: Alexander Cede Instrument type: Pandora Instrument number: 63 Spectrometer number: 1 Processing software version used: BlickP v1.8.49 Instrument operation file used: Pandora63_OF_v2d20210219.txt Instrument calibration file used: Pandora63s1_CF_v7d20230718.txt Level 1 file used: Pandora63s1_LaPorteTX_20231201_L1_smca1c7d20230718p1-8.txt Full location name: University Of Houston Trailer Short location name: LaPorteTX Country of location: United States Location latitude [deg]: 29.6721 Location longitude [deg]: -95.0647 Location altitude [m]: 22 Local noon date: 20231201 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: Redoing OoB session using Lorentzian slit function due to bad straylight First and last pixel inside fitting window: 142 283 Nominal wavelengths inside fitting window [nm]: 305.1336 305.2750 305.4165 305.5579 305.6993 305.8407 305.9821 306.1235 306.2648 306.4062 306.5475 306.6888 306.8302 306.9715 307.1128 307.2541 307.3953 307.5366 307.6778 307.8191 307.9603 308.1015 308.2427 308.3839 308.5251 308.6663 308.8075 308.9486 309.0898 309.2309 309.3720 309.5131 309.6542 309.7953 309.9364 310.0774 310.2185 310.3595 310.5005 310.6415 310.7825 310.9235 311.0645 311.2055 311.3465 311.4874 311.6283 311.7693 311.9102 312.0511 312.1920 312.3328 312.4737 312.6146 312.7554 312.8963 313.0371 313.1779 313.3187 313.4595 313.6003 313.7410 313.8818 314.0225 314.1632 314.3040 314.4447 314.5854 314.7261 314.8667 315.0074 315.1481 315.2887 315.4293 315.5699 315.7105 315.8511 315.9917 316.1323 316.2729 316.4134 316.5539 316.6945 316.8350 316.9755 317.1160 317.2565 317.3969 317.5374 317.6778 317.8183 317.9587 318.0991 318.2395 318.3799 318.5203 318.6606 318.8010 318.9413 319.0816 319.2220 319.3623 319.5026 319.6429 319.7831 319.9234 320.0636 320.2039 320.3441 320.4843 320.6245 320.7647 320.9049 321.0451 321.1852 321.3254 321.4655 321.6056 321.7457 321.8858 322.0259 322.1660 322.3061 322.4461 322.5861 322.7262 322.8662 323.0062 323.1462 323.2862 323.4262 323.5661 323.7061 323.8460 323.9859 324.1258 324.2657 324.4056 324.5455 324.6854 324.8252 324.9651 --------------------------------------------------------------------------------------- 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-241: Unweighted slant column residuals for each pixel inside the fitting window multiplied by 1e5, 9e5=pixel was not used for fitting Columns 242-383: 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 --------------------------------------------------------------------------------------- SU 20231201T165420.0Z 8735.704399 1425 1 49.09 2 54.57 158.53 90.54 298.87 0.00 1 0.00 1 0 17 1.618e-02 1.157e-02 5.460e-03 2.129e-03 2.111657e-01 1.075403e-03 2.347910e-03 -6.000000e+00 9.303401e-03 224.64 -2.00 3.24 -2.00 316.96 0.00e+00 1.091383e-03 1.488739e-04 4.601902e-04 -6.000000e+00 4.322486e-04 287.37 -2.00 20.00 -2.00 316.15 0.00e+00 1.720 0.000 302.00 328.00 3 5.5412e+00 2.8084e-03 1.9618e-02 -6.0000e+00 2.2918e-02 2.8651e-02 3.1684e-03 9.6377e-03 -6.0000e+00 2.4938e-02 2.3343e-02 1.8604e-03 4.7439e-03 -6.0000e+00 1.4926e-02 -2.2695e-02 8.9621e-04 1.3737e-03 -6.0000e+00 6.3606e-03 2.9439e-04 8.0753e-04 6.1764e-05 1.3264e-05 -6.0000e+00 4.4569e-04 -2.3430e-03 0.0000e+00 12 9 8 11 8 0 99.99 22.19 1.15 -0.00234 0.00373 42 9 1 5 0 1168.620 22.19 15.01 18.51 30.63 20.7 1024.22 1793 689 1092 -3600 3366 -2055 5656 -2473 2929 -4447 -4395 237 -1047 1062 -1760 332 261 3149 -345 2765 -170 1220 1109 604 -152 -1847 448 -1663 -546 -1899 2672 -2677 -2234 -3407 -3251 -3565 -2049 -660 -1188 -525 -422 467 259 799 682 -441 32 -621 -806 88 1358 671 985 311 772 1176 476 333 1089 1541 2343 2278 1596 1548 -267 303 -459 897 1425 257 -28 389 258 -340 -1228 -2521 -2958 -2787 -1933 -2575 -2164 -2126 -1806 -44 112 857 950 1323 1296 676 -271 -1307 -1425 -723 -717 -354 125 478 -209 -328 -1429 -1701 -549 -394 -1217 -174 278 1133 1956 2223 1686 1230 463 -260 -940 -1093 -714 136 687 657 776 780 501 -112 -292 -2412 -1348 -1382 -742 -1036 -1243 -993 -175 713 790 1160 1036 1135 302 47 -1072 -975 262 105 163 -512 429 -274 770 -393 471 -773 -752 42 -212 241 -453 88 71 870 -101 821 -52 363 341 198 -52 -648 147 -525 -163 -579 840 -949 -835 -1231 -1117 -1279 -834 -303 -639 -295 -251 284 160 491 421 -271 20 -377 -484 54 879 443 678 221 556 864 352 242 780 1123 1800 1891 1358 1340 -229 234 -355 716 1232 234 -27 372 240 -314 -1034 -2017 -2220 -2130 -1638 -2349 -2160 -2180 -1923 -49 129 1039 1187 1639 1566 798 -294 -1345 -1409 -747 -811 -419 155 631 -277 -407 -1696 -1935 -641 -499 -1578 -232 385 1626 2854 3383 2600 1826 668 -344 -1229 -1426 -970 188 1013 995 1159 1127 699 -156 -391 -3211 -1780 -1876 -1075 -1508 -1835 -1483 -271 1118 1303 1949 1776 1890 516 78 -1725 -1546 SU 20231201T171109.1Z 8735.716078 1554 1 49.04 2 53.37 163.13 93.62 300.89 0.00 1 0.00 1 0 25 3.395e-02 1.241e-02 1.454e-02 4.499e-03 2.137985e-01 2.264236e-03 2.359953e-03 -6.000000e+00 1.952839e-02 224.64 -2.00 3.24 -2.00 317.45 0.00e+00 1.664560e-04 3.559218e-04 5.400963e-04 -6.000000e+00 9.073219e-04 287.37 -2.00 20.00 -2.00 316.61 0.00e+00 1.672 0.000 302.00 328.00 3 6.5745e+00 5.9773e-03 1.9339e-02 -6.0000e+00 4.8106e-02 2.3671e-02 6.7781e-03 1.0059e-02 -6.0000e+00 5.2347e-02 6.1601e-02 4.0741e-03 4.8177e-03 -6.0000e+00 3.1331e-02 -4.1127e-02 2.0134e-03 1.3784e-03 -6.0000e+00 1.3351e-02 1.0821e-04 9.3568e-04 1.3129e-04 1.3775e-05 -6.0000e+00 9.3554e-04 -2.2718e-03 0.0000e+00 12 8 8 10 0 0 99.99 22.19 2.97 -0.00227 0.00476 43 10 1 5 0 1140.300 22.19 15.00 18.52 30.36 20.8 1022.95 5940 9302 2466 -3292 21492 1550 22431 1077 3533 -5333 909 -2165 3119 756 -2381 -1008 1282 9472 -454 4068 -3553 1710 4049 -756 807 350 -691 2872 -862 -626 -975 -2419 -3112 -3093 -1934 -4952 -2807 -3082 -2342 -1088 -1339 175 -437 -606 -1029 -491 -1708 -391 -2000 -151 121 1913 -131 1211 539 1259 729 -373 464 1448 2063 2813 665 1133 -995 235 135 1524 1465 1125 -183 1192 1019 -895 -503 -2069 -2950 -4111 -767 -1800 -1124 -2094 -1945 491 566 684 1377 1625 605 1531 -104 -966 -2464 -441 -942 -109 440 580 221 -222 -1812 -1934 -53 -230 -593 -224 185 247 2057 1930 1421 1646 171 -438 -1172 -1406 -1345 66 1012 671 138 1135 -223 -393 -557 -2569 -1224 -1492 -965 -1268 -1396 -602 -304 269 602 1123 1177 1487 330 299 -907 -576 638 1037 280 -373 1824 154 2093 136 468 -780 124 -327 507 138 -497 -219 280 1990 -108 970 -913 416 997 -209 230 101 -192 751 -224 -163 -275 -742 -1018 -968 -571 -1575 -1009 -1287 -1133 -574 -723 99 -245 -351 -601 -280 -959 -211 -1083 -83 71 1166 -85 800 366 869 506 -257 315 993 1494 2181 558 963 -805 171 97 1130 1192 999 -172 1109 923 -782 -411 -1540 -2149 -2977 -601 -1589 -1083 -2102 -2030 533 659 817 1683 1969 747 1748 -114 -963 -2393 -441 -1037 -126 550 775 294 -277 -2175 -2221 -62 -294 -771 -303 259 360 3135 2974 2216 2542 248 -607 -1538 -1844 -1841 95 1490 1018 208 1695 -326 -547 -744 -3400 -1602 -2019 -1395 -1915 -2137 -927 -489 439 1031 1954 2086 2648 583 510 -1505 -940 SU 20231201T175323.3Z 8735.745409 1866 1 1.99 2 51.65 175.29 101.03 306.43 0.00 1 0.00 1 0 17 1.491e-02 1.126e-02 6.881e-03 3.944e-03 1.977698e-01 2.005824e-03 2.181768e-03 -6.000000e+00 8.575814e-03 224.64 -2.00 3.24 -2.00 315.81 0.00e+00 1.116659e-03 2.144655e-04 2.295119e-04 -6.000000e+00 3.983655e-04 287.37 -2.00 20.00 -2.00 315.06 0.00e+00 1.608 0.000 302.00 328.00 3 3.6706e+00 5.1494e-03 1.8323e-02 -6.0000e+00 2.1125e-02 1.6696e-02 5.7706e-03 8.9881e-03 -6.0000e+00 2.2986e-02 -5.5878e-03 3.3114e-03 4.3921e-03 -6.0000e+00 1.3758e-02 -7.8698e-03 1.5075e-03 1.3003e-03 -6.0000e+00 5.8640e-03 2.2191e-03 9.1239e-04 1.1257e-04 1.3805e-05 -6.0000e+00 4.1066e-04 -3.6941e-03 0.0000e+00 12 9 9 12 1 1 999.00 22.19 0.34 -0.00369 0.00316 2 16 1 5 0 661.660 22.19 15.00 18.53 30.72 20.7 1024.82 -1020 1626 -2723 -1053 -2677 -2791 -676 -2225 -949 -3299 -2689 -2625 -16 -133 681 297 2596 3437 2296 3086 -884 2706 1182 1061 391 76 -3232 292 -2561 -953 -416 -2402 -837 -4234 -3285 -3625 -572 -407 513 -341 678 650 594 646 923 -332 -393 -150 -113 873 144 1079 758 942 775 1145 1125 489 1556 936 3126 1489 2434 962 -582 -474 398 829 529 655 -164 -817 -457 57 -1116 -2142 -3679 -3045 -2083 -2964 -1925 -2367 -1032 -233 293 388 698 1030 416 433 561 -1033 -1348 -103 -903 35 175 446 -265 -141 -1521 -1242 -768 -895 -1087 278 328 495 1571 2083 1889 825 997 -205 -713 -292 -410 475 566 288 752 1969 464 75 -563 -1096 -987 -728 -617 -865 -959 -1063 -151 316 776 320 722 1119 -2 276 -1038 -1000 -258 412 -705 -242 -606 -624 -162 -576 -258 -918 -729 -763 -5 -46 255 117 1021 1394 915 1270 -376 1099 500 476 184 36 -1501 130 -1115 -408 -191 -1136 -411 -2050 -1497 -1723 -303 -234 325 -239 476 470 436 471 674 -241 -276 -106 -79 621 106 828 609 743 622 937 920 396 1287 765 2651 1331 2281 894 -531 -407 336 696 487 630 -166 -838 -454 56 -1034 -1875 -3048 -2551 -1853 -2851 -1938 -2468 -1106 -264 343 466 863 1264 500 503 615 -1074 -1397 -110 -993 40 213 560 -344 -171 -1831 -1427 -913 -1119 -1390 366 434 682 2206 2949 2699 1176 1359 -267 -908 -368 -524 643 796 405 1073 2687 641 100 -743 -1407 -1268 -962 -845 -1197 -1343 -1510 -216 476 1188 504 1119 1742 -2 425 -1604 -1524