File name: Pandora147s1_SWDetroitMI_20231010_L2Fit_fout2c9d20221206p1-8.txt
File generation date: 20260323T005057.6Z
Data description: Level 2 spectral fitting results file
Data file version: fout2c9p1-8
Local principal investigator: Lukas Valin
Network principal investigator: Alexander Cede
Instrument type: Pandora
Instrument number: 147
Spectrometer number: 1
Processing software version used: BlickP v1.8.88
Instrument operation file used: Pandora147_OF_v4d20210429.txt
Instrument calibration file used: Pandora147s1_CF_v9d20221206.txt
Level 1 file used: Pandora147s1_SWDetroitMI_20231010_L1_smca1c9d20221206p1-8.txt.bz2
Full location name: SW Detroit (SWHS) (AQS-26-163-0015)
Short location name: SWDetroitMI
Country of location: USA
Location latitude [deg]: 42.3026
Location longitude [deg]: -83.1068
Location altitude [m]: 178
Local noon date: 20231010
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: 236 378
Nominal wavelengths inside fitting window [nm]: 305.0989 305.2398 305.3807 305.5215 305.6624 305.8032 305.9440 306.0848 306.2256 306.3664 306.5072 306.6479 306.7887 306.9294 307.0701 307.2108 307.3516 307.4922 307.6329 307.7736 307.9143 308.0549 308.1955 308.3362 308.4768 308.6174 308.7580 308.8985 309.0391 309.1797 309.3202 309.4607 309.6013 309.7418 309.8823 310.0227 310.1632 310.3037 310.4441 310.5846 310.7250 310.8654 311.0058 311.1462 311.2866 311.4270 311.5673 311.7077 311.8480 311.9883 312.1286 312.2689 312.4092 312.5495 312.6898 312.8300 312.9703 313.1105 313.2507 313.3909 313.5311 313.6713 313.8115 313.9517 314.0918 314.2319 314.3721 314.5122 314.6523 314.7924 314.9325 315.0725 315.2126 315.3526 315.4927 315.6327 315.7727 315.9127 316.0527 316.1927 316.3326 316.4726 316.6125 316.7525 316.8924 317.0323 317.1722 317.3121 317.4519 317.5918 317.7317 317.8715 318.0113 318.1511 318.2909 318.4307 318.5705 318.7103 318.8500 318.9898 319.1295 319.2692 319.4089 319.5486 319.6883 319.8280 319.9677 320.1073 320.2469 320.3866 320.5262 320.6658 320.8054 320.9450 321.0845 321.2241 321.3636 321.5032 321.6427 321.7822 321.9217 322.0612 322.2007 322.3401 322.4796 322.6190 322.7584 322.8979 323.0373 323.1767 323.3160 323.4554 323.5948 323.7341 323.8734 324.0128 324.1521 324.2914 324.4307 324.5699 324.7092 324.8485 324.9877
---------------------------------------------------------------------------------------
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 20231010T152930.0Z 8683.645487 304 1 14.88 2 55.04 145.98 30.12 205.95 0.00 1 0.00 1 0 17 1.070e-02 6.536e-03 3.068e-03 1.298e-03 2.448098e-01 6.431411e-04 2.961068e-03 -6.000000e+00 5.981699e-03 226.22 -2.00 3.54 -2.00 316.63 0.00e+00 5.268720e-04 7.926363e-05 9.101030e-05 -6.000000e+00 2.670345e-04 282.89 -2.00 20.00 -2.00 315.90 0.00e+00 1.740 0.000 302.00 328.00 3 3.1926e+00 1.6746e-03 2.1156e-02 -6.0000e+00 1.4810e-02 3.4082e-02 1.8858e-03 1.2503e-02 -6.0000e+00 1.6141e-02 -1.6001e-03 1.1025e-03 6.1614e-03 -6.0000e+00 9.6262e-03 -1.3312e-03 5.2526e-04 1.7470e-03 -6.0000e+00 4.0835e-03 2.9298e-03 5.3920e-05 3.6181e-05 1.9267e-05 -6.0000e+00 2.8093e-04 3.6991e-04 3.0898e-07 12 9 1 12 1 1 99.95 16.13 0.23 0.00037 0.00078 32 25 1 5 0 450.690 22.29 15.00 16.13 19.33 30.2 937.03 -1407 -820 -755 -4456 -2856 -2841 -1760 -820 -1279 -3144 -2703 -2405 -1666 -398 692 -584 65 164 1155 1273 1013 204 891 828 472 964 -220 -532 -759 857 714 1949 1211 -1456 -2760 -2288 -2517 -719 -335 66 -480 -165 -540 -406 54 -230 -28 -273 -559 46 824 791 1208 592 1134 1222 694 278 172 354 1306 1505 2061 1453 430 -431 -1232 -425 652 412 168 283 378 372 -153 -256 -1064 -1030 -756 -830 -1114 -1658 -2034 -1644 -1000 -548 24 373 363 296 434 -69 -266 1 69 116 787 613 369 -292 -667 -1018 -795 -318 -690 -878 -672 -691 165 274 709 637 361 139 -181 58 97 445 263 128 -68 260 317 758 579 269 591 786 487 436 -102 -216 -193 -101 46 -244 231 176 -25 -194 -320 -528 -683 -216 -134 -123 -685 -422 -411 -272 -142 -223 -573 -526 -488 -347 -94 183 -169 19 50 363 413 332 66 297 293 176 368 -83 -187 -250 278 248 723 488 -600 -1063 -879 -1081 -329 -183 38 -299 -106 -357 -268 36 -151 -18 -169 -341 28 534 548 878 453 883 970 561 223 135 278 1058 1223 1806 1347 390 -358 -1036 -350 541 389 160 277 365 368 -140 -208 -855 -800 -629 -728 -1096 -1675 -2174 -1837 -1119 -652 29 473 456 354 507 -73 -279 1 72 134 938 805 475 -389 -826 -1237 -930 -369 -871 -1164 -874 -938 235 392 1052 942 538 191 -243 74 125 602 362 187 -101 376 444 1071 778 344 766 976 653 624 -144 -310 -283 -154 71 -392 378 294 -42 -314 -501 -843 -1079
SS 20231010T163214.0Z 8683.689051 341 1 3.17 2 50.20 164.66 37.37 229.69 0.00 1 0.00 1 0 17 1.048e-02 6.985e-03 3.709e-03 2.081e-03 2.200500e-01 1.037140e-03 2.677258e-03 -6.000000e+00 5.857569e-03 226.22 -2.00 3.54 -2.00 315.77 0.00e+00 5.159181e-04 1.054015e-04 7.080153e-05 -6.000000e+00 2.615109e-04 282.89 -2.00 20.00 -2.00 315.10 0.00e+00 1.559 0.000 302.00 328.00 3 2.8944e+00 2.6734e-03 1.8961e-02 -6.0000e+00 1.4503e-02 2.6709e-02 2.9948e-03 1.1209e-02 -6.0000e+00 1.5806e-02 -2.4256e-03 1.7207e-03 5.5240e-03 -6.0000e+00 9.4265e-03 -1.5649e-03 7.8747e-04 1.5629e-03 -6.0000e+00 3.9986e-03 4.9289e-03 -1.1051e-04 5.7287e-05 1.8230e-05 -6.0000e+00 2.7513e-04 4.5184e-04 -1.1975e-05 12 9 1 12 1 1 98.68 16.13 0.33 0.00045 -0.00038 9 36 1 5 0 316.430 22.29 15.00 16.13 20.30 30.4 944.55 -357 -328 -1723 -2619 -4843 -397 -716 1289 -1127 -1301 -3609 -1684 -773 -514 1195 -530 84 872 709 1415 1479 -54 1811 774 891 267 182 -656 66 1003 1432 1249 1432 -1815 -3160 -2248 -2484 -459 -195 -58 -1141 -134 -325 -42 209 -206 -478 -552 -714 71 720 971 961 951 617 691 1118 441 311 -84 1216 1532 2247 1260 113 -331 -503 -444 1111 156 589 115 436 -199 219 -491 -1204 -1211 -814 -1158 -1107 -1358 -2050 -1463 -1334 -424 -507 341 535 248 509 144 -143 205 -194 218 515 828 236 -318 -670 -1085 -643 -150 -933 -893 -677 -536 334 397 672 550 434 230 -132 327 316 231 218 154 -99 529 415 766 663 195 412 679 479 354 -161 -365 -353 46 38 -161 71 343 -46 -149 -503 -485 -705 -82 -80 -425 -644 -1078 -82 -161 321 -289 -342 -942 -448 -229 -173 446 -199 32 342 288 589 619 -23 765 348 420 130 87 -291 27 406 612 575 707 -919 -1493 -1053 -1206 -258 -123 -39 -829 -100 -248 -32 158 -155 -351 -393 -496 49 528 757 741 766 509 580 945 371 256 -69 1037 1329 2091 1242 109 -290 -433 -377 1014 154 591 119 445 -195 209 -444 -994 -1022 -689 -1050 -1118 -1424 -2165 -1611 -1579 -506 -637 415 675 298 589 160 -145 204 -208 245 607 1021 305 -398 -817 -1282 -731 -180 -1169 -1120 -885 -729 453 572 964 784 613 312 -173 404 395 304 296 212 -139 725 578 1011 872 244 514 864 624 469 -223 -515 -506 65 58 -245 110 543 -73 -230 -751 -728 -1049
SS 20231010T164127.1Z 8683.695452 348 1 0.75 2 49.79 167.60 38.70 232.53 0.00 1 0.00 1 0 17 1.282e-02 8.577e-03 9.494e-03 5.678e-03 2.158137e-01 2.839026e-03 2.663606e-03 -6.000000e+00 7.162180e-03 226.22 -2.00 3.54 -2.00 315.55 0.00e+00 4.112773e-04 2.747297e-04 2.740136e-05 -6.000000e+00 3.200114e-04 282.89 -2.00 20.00 -2.00 314.89 0.00e+00 1.546 0.000 302.00 328.00 3 2.9025e+00 7.2922e-03 1.8780e-02 -6.0000e+00 1.7732e-02 2.8217e-02 8.1650e-03 1.1092e-02 -6.0000e+00 1.9330e-02 4.7373e-03 4.6836e-03 5.4733e-03 -6.0000e+00 1.1527e-02 -2.4471e-03 2.1174e-03 1.5470e-03 -6.0000e+00 4.8875e-03 4.9714e-03 -6.9405e-04 1.5624e-04 1.8954e-05 -6.0000e+00 3.3669e-04 1.4351e-03 -1.1876e-04 12 9 1 12 1 1 999.00 16.15 0.32 0.00144 -0.00378 1 31 1 5 0 362.270 22.29 15.00 16.15 21.26 30.5 953.01 307 2489 -1129 -1850 -3370 2026 -3730 5335 126 -1281 -4999 383 -366 -119 -731 890 -81 901 325 -107 -321 -1104 1209 771 250 1626 961 -2309 -515 806 1381 1500 2968 -442 -1089 -4067 -1982 -960 605 -1109 -1044 -1158 -311 -43 -953 -102 -953 825 -1370 659 -541 2018 468 1605 -227 1888 1152 1082 46 921 104 1173 2257 742 -387 317 -2090 136 1070 -12 708 132 81 806 842 -1070 -263 -1858 -78 -693 -1112 -1015 -1658 -1775 122 -1497 -329 -101 932 721 -163 1114 -913 -49 -370 384 1248 618 229 -979 -443 -699 -1282 332 -1516 -529 30 217 -941 1051 17 262 138 824 664 -599 336 144 -179 -584 -296 27 277 846 621 182 436 1399 1145 346 351 -287 -230 -247 60 -66 -51 510 -565 311 -213 -944 -631 80 674 -315 -470 -806 448 -945 1394 36 -384 -1508 114 -118 -42 -290 358 -34 383 144 -47 -143 -488 557 363 124 817 474 -1088 -230 353 639 717 1500 -229 -523 -2016 -1033 -565 393 -787 -769 -884 -240 -33 -739 -78 -712 604 -998 475 -406 1572 378 1338 -191 1586 987 915 38 766 91 1065 2133 726 -370 288 -1846 116 959 -12 719 136 85 817 806 -956 -223 -1539 -68 -647 -1104 -1073 -1811 -1986 141 -1794 -401 -126 1139 858 -188 1202 -941 -50 -400 434 1477 774 292 -1247 -544 -809 -1479 390 -1887 -669 39 290 -1281 1489 25 379 193 1098 846 -749 421 186 -241 -818 -411 38 380 1123 802 229 533 1725 1482 459 474 -396 -322 -354 88 -99 -78 776 -870 470 -319 -1395 -919