File name: Pandora63s1_LaPorteTX_20211208_L2Fit_fout2c2d20210924p1-8.txt File generation date: 20250514T162750.4Z Data description: Level 2 spectral fitting results file Data file version: fout2c2p1-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.75 Instrument operation file used: Pandora63_OF_v2d20210219.txt Instrument calibration file used: Pandora63s1_CF_v2d20210924.txt Level 1 file used: Pandora63s1_LaPorteTX_20211208_L1_smca1c2d20210924p1-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: 20211208 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: 141 283 Nominal wavelengths inside fitting window [nm]: 305.0348 305.1762 305.3175 305.4588 305.6002 305.7415 305.8828 306.0240 306.1653 306.3066 306.4478 306.5891 306.7303 306.8715 307.0127 307.1539 307.2951 307.4362 307.5774 307.7185 307.8597 308.0008 308.1419 308.2830 308.4241 308.5652 308.7063 308.8473 308.9884 309.1294 309.2704 309.4114 309.5524 309.6934 309.8344 309.9753 310.1163 310.2572 310.3982 310.5391 310.6800 310.8209 310.9618 311.1027 311.2435 311.3844 311.5252 311.6661 311.8069 311.9477 312.0885 312.2293 312.3700 312.5108 312.6515 312.7923 312.9330 313.0737 313.2144 313.3551 313.4958 313.6365 313.7771 313.9178 314.0584 314.1991 314.3397 314.4803 314.6209 314.7615 314.9020 315.0426 315.1831 315.3237 315.4642 315.6047 315.7452 315.8857 316.0262 316.1667 316.3071 316.4476 316.5880 316.7284 316.8688 317.0092 317.1496 317.2900 317.4304 317.5707 317.7111 317.8514 317.9917 318.1321 318.2724 318.4126 318.5529 318.6932 318.8334 318.9737 319.1139 319.2541 319.3944 319.5346 319.6747 319.8149 319.9551 320.0952 320.2354 320.3755 320.5156 320.6557 320.7958 320.9359 321.0760 321.2161 321.3561 321.4962 321.6362 321.7762 321.9162 322.0562 322.1962 322.3362 322.4761 322.6161 322.7560 322.8959 323.0359 323.1758 323.3157 323.4555 323.5954 323.7353 323.8751 324.0150 324.1548 324.2946 324.4344 324.5742 324.7140 324.8538 324.9935 --------------------------------------------------------------------------------------- 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 20211208T173158.8Z 8012.730542 180 1 48.99 2 53.35 168.38 89.73 114.41 0.00 1 0.00 1 0 25 1.787e-02 1.396e-02 3.483e-03 1.372e-03 1.819246e-01 6.674819e-04 2.030857e-03 -6.000000e+00 9.955740e-03 224.45 -2.00 3.24 -2.00 317.08 0.00e+00 1.192537e-03 9.303440e-05 4.785492e-04 -6.000000e+00 4.574839e-04 286.84 -2.00 20.00 -2.00 316.24 0.00e+00 1.671 0.000 302.00 328.00 3 3.9102e+00 1.7484e-03 1.8512e-02 -6.0000e+00 2.4512e-02 3.6306e-02 1.9792e-03 8.2693e-03 -6.0000e+00 2.6725e-02 2.1177e-04 1.1597e-03 4.0507e-03 -6.0000e+00 1.6004e-02 -1.1649e-02 5.5843e-04 1.1783e-03 -6.0000e+00 6.8126e-03 1.6800e-03 4.3503e-04 3.8904e-05 1.1058e-05 -6.0000e+00 4.7135e-04 -1.3656e-03 1.0390e-04 12 8 8 10 0 0 99.98 18.22 0.12 -0.00137 0.00190 76 17 1 5 0 644.490 21.88 14.98 18.22 28.16 8.3 963.66 -1062 1911 -623 -466 -4881 -2044 -3890 -1232 -1273 -3363 -5062 -4241 -1799 -932 614 -39 973 1597 1967 1737 1858 1207 814 1193 277 509 -275 -1527 -2647 -1653 -1119 -668 -371 -3098 -5229 -4749 -3508 -2136 195 340 741 267 1120 671 1390 882 347 188 -1182 -85 372 930 1202 1647 903 1115 618 -188 141 289 1620 1968 2852 2693 868 -676 -1204 -2151 130 1100 563 1169 744 737 -785 -1582 -3401 -4342 -3442 -3112 -2176 -1939 -2424 -1225 -74 262 1300 1384 1881 1023 501 -862 -2432 -2267 -2078 -928 225 1298 741 76 -700 -1968 -2248 -1095 -1080 -1070 232 196 1281 1899 2383 2055 1400 122 -1192 -1861 -1554 -991 363 988 1193 1608 1548 572 177 -1234 -1999 -1351 -1479 -1158 -625 -1242 -743 -227 580 767 1370 1454 944 260 -202 -1204 -1381 -158 304 -103 -73 -722 -272 -543 -187 -214 -607 -913 -761 -340 -199 144 -10 259 431 544 501 546 354 241 365 90 171 -93 -491 -809 -483 -336 -213 -126 -1102 -1803 -1576 -1224 -854 90 180 421 159 678 411 851 541 210 111 -684 -49 220 577 785 1113 628 787 444 -136 100 204 1176 1514 2330 2299 751 -547 -902 -1580 100 931 517 1099 707 682 -695 -1288 -2570 -3152 -2558 -2561 -1966 -1874 -2466 -1288 -81 306 1579 1726 2325 1237 572 -915 -2416 -2204 -2114 -1017 264 1639 979 98 -863 -2293 -2569 -1293 -1374 -1425 316 276 1856 2871 3696 3182 2106 174 -1594 -2392 -2026 -1343 516 1473 1817 2400 2260 815 245 -1655 -2588 -1763 -2051 -1674 -930 -1859 -1125 -354 947 1295 2339 2528 1653 452 -344 -2013 -2300 SS 20211208T174442.8Z 8012.739384 189 1 49.23 2 52.88 172.00 87.30 115.95 0.00 1 0.00 1 0 25 1.738e-02 1.383e-02 2.484e-03 1.000e-03 1.798261e-01 4.866863e-04 1.999447e-03 -6.000000e+00 9.677673e-03 224.45 -2.00 3.24 -2.00 317.00 0.00e+00 1.216438e-03 6.672795e-05 4.526991e-04 -6.000000e+00 4.447439e-04 286.84 -2.00 20.00 -2.00 316.17 0.00e+00 1.653 0.000 302.00 328.00 3 3.3625e+00 1.2733e-03 1.8302e-02 -6.0000e+00 2.3827e-02 3.2123e-02 1.4404e-03 8.1477e-03 -6.0000e+00 2.5979e-02 -2.7355e-04 8.4299e-04 3.9946e-03 -6.0000e+00 1.5558e-02 -1.1376e-02 4.0556e-04 1.1707e-03 -6.0000e+00 6.6218e-03 2.9740e-03 3.1120e-04 2.8391e-05 1.0783e-05 -6.0000e+00 4.5827e-04 -1.2284e-03 7.6278e-05 12 8 8 10 0 0 99.91 18.24 0.07 -0.00123 0.00111 131 29 1 5 0 375.730 21.98 15.00 18.24 28.47 8.3 965.36 649 -861 681 -1051 -3157 -3945 -4141 -2849 -1214 -2348 -3992 -2938 -2719 -270 979 1442 1122 1786 1862 1528 2162 665 1945 951 791 307 -583 -1853 -1704 -1696 -1291 -971 -416 -2562 -5412 -4293 -3995 -1672 -75 423 639 179 675 684 1195 883 -47 -213 -965 -806 579 744 1209 1595 878 1358 414 -5 -155 439 1098 2281 2802 2823 872 -598 -1055 -1955 251 686 823 1098 838 564 -614 -2039 -3146 -4397 -3588 -2757 -2311 -1786 -2352 -1236 -14 407 1109 1528 1983 1134 309 -802 -2323 -2244 -1802 -1087 17 1289 713 -3 -529 -1803 -2074 -1010 -954 -935 81 149 1181 1961 2349 2019 1468 137 -1225 -1914 -1541 -910 524 923 1281 1501 1491 357 304 -1386 -2074 -1383 -1432 -1039 -787 -1154 -1002 -262 451 691 1358 1519 1019 384 -125 -1268 -1372 98 -145 114 -169 -479 -554 -594 -451 -211 -430 -729 -528 -527 -58 235 373 301 487 522 443 657 203 586 300 264 106 -201 -609 -520 -513 -399 -321 -146 -937 -1914 -1439 -1422 -681 -35 225 369 109 416 425 738 544 -28 -129 -569 -470 346 473 806 1101 623 977 305 -3 -111 313 804 1770 2339 2393 744 -489 -783 -1461 197 584 762 1054 803 524 -539 -1684 -2349 -3252 -2746 -2260 -2110 -1708 -2393 -1301 -16 467 1341 1890 2432 1397 354 -867 -2304 -2158 -1827 -1203 20 1628 951 -4 -645 -2115 -2364 -1203 -1191 -1233 110 211 1730 2915 3618 3105 2227 193 -1653 -2458 -1987 -1231 729 1364 1932 2223 2140 514 417 -1825 -2710 -1821 -1945 -1499 -1173 -1741 -1535 -401 731 1168 2341 2664 1801 668 -211 -2085 -2245 SS 20211208T180601.9Z 8012.754189 201 1 50.83 2 52.49 178.16 83.30 118.64 0.00 1 0.00 1 0 25 1.791e-02 1.474e-02 6.400e-03 2.593e-03 1.792837e-01 1.262173e-03 1.993503e-03 -6.000000e+00 9.974078e-03 224.45 -2.00 3.24 -2.00 316.98 0.00e+00 1.343039e-03 1.722047e-04 4.555036e-04 -6.000000e+00 4.583284e-04 286.84 -2.00 20.00 -2.00 316.16 0.00e+00 1.638 0.000 302.00 328.00 3 5.2264e+00 3.3007e-03 1.8151e-02 -6.0000e+00 2.4557e-02 4.7690e-02 3.7338e-03 8.0951e-03 -6.0000e+00 2.6774e-02 3.6303e-04 2.1816e-03 3.9777e-03 -6.0000e+00 1.6034e-02 -1.2055e-02 1.0475e-03 1.1910e-03 -6.0000e+00 6.8251e-03 4.6499e-04 4.0797e-04 7.3374e-05 1.0731e-05 -6.0000e+00 4.7223e-04 -1.4103e-03 9.9656e-05 12 8 8 10 0 0 99.87 18.27 0.16 -0.00141 0.00166 14 3 1 5 0 3629.940 21.98 15.01 18.27 28.82 8.2 966.87 216 -1928 1421 -3468 -3221 -3166 -2676 -1074 -2861 -699 -5864 -2615 -2344 -687 380 617 2992 447 1904 661 1478 1340 1342 25 1048 377 150 -2346 -1982 -113 -2150 -1126 -1135 -1797 -4618 -5006 -3077 -1459 92 264 720 222 749 657 777 816 590 -824 -65 171 349 727 1284 1950 1273 1261 967 324 -225 438 815 2372 2871 2623 836 -562 -1218 -2455 165 680 902 885 448 566 -536 -2119 -3387 -3961 -4129 -3014 -2367 -1853 -2140 -1400 -70 123 1647 1557 2259 1083 382 -1559 -2922 -2766 -2216 -633 261 1285 701 76 -165 -1989 -2227 -997 -868 -925 44 224 1206 2072 2756 2069 1490 -59 -1078 -2303 -1668 -844 547 637 1457 1496 1653 684 213 -930 -1988 -1716 -1757 -1268 -932 -1250 -935 -358 804 525 1637 1359 1274 377 -246 -1309 -1408 33 -329 240 -583 -487 -443 -380 -168 -507 -127 -1099 -476 -460 -150 93 163 802 126 545 197 448 405 407 8 348 130 51 -784 -620 -33 -667 -371 -399 -648 -1621 -1701 -1091 -592 44 142 418 135 464 412 486 510 362 -498 -38 99 209 461 850 1333 895 902 702 236 -162 315 600 1840 2373 2273 730 -461 -925 -1826 128 583 832 843 429 527 -476 -1746 -2568 -2881 -3125 -2496 -2145 -1806 -2191 -1483 -78 145 2004 1942 2803 1319 439 -1679 -2936 -2710 -2270 -693 308 1624 924 99 -203 -2326 -2544 -1180 -1102 -1229 60 314 1746 3119 4259 3202 2229 -83 -1438 -2976 -2159 -1136 773 948 2201 2229 2411 973 293 -1235 -2565 -2232 -2416 -1825 -1374 -1866 -1409 -557 1300 882 2782 2347 2207 648 -414 -2159 -2318