File name: Pandora145s1_Toronto-Scarborough_20231204_L2Fit_fout2c8d20221111p1-8.txt File generation date: 20231205T054934.2Z Data description: Level 2 spectral fitting results file Data file version: fout2c8p1-8 Local principal investigator: Vitali Fioletov Network principal investigator: Alexander Cede Instrument type: Pandora Instrument number: 145 Spectrometer number: 1 Processing software version used: BlickP v1.8.49 Instrument operation file used: Pandora145_OF_v5d20190718.txt Instrument calibration file used: Pandora145s1_CF_v8d20221111.txt Level 1 file used: Pandora145s1_Toronto-Scarborough_20231204_L1_smca1c8d20221111p1-8.txt Full location name: University of Toronto - Scarborough Short location name: Toronto-Scarborough Country of location: Canada Location latitude [deg]: 43.7843 Location longitude [deg]: -79.1874 Location altitude [m]: 137 Local noon date: 20231204 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: based on session 4 for field cal VP 20221111 First and last pixel inside fitting window: 226 368 Nominal wavelengths inside fitting window [nm]: 305.0546 305.1954 305.3363 305.4771 305.6179 305.7586 305.8994 306.0402 306.1809 306.3217 306.4624 306.6031 306.7438 306.8845 307.0252 307.1659 307.3065 307.4472 307.5878 307.7284 307.8691 308.0097 308.1503 308.2908 308.4314 308.5720 308.7125 308.8531 308.9936 309.1341 309.2746 309.4151 309.5556 309.6961 309.8365 309.9770 310.1174 310.2578 310.3982 310.5386 310.6790 310.8194 310.9598 311.1001 311.2405 311.3808 311.5212 311.6615 311.8018 311.9421 312.0823 312.2226 312.3629 312.5031 312.6433 312.7836 312.9238 313.0640 313.2042 313.3444 313.4845 313.6247 313.7648 313.9049 314.0451 314.1852 314.3253 314.4654 314.6054 314.7455 314.8856 315.0256 315.1656 315.3057 315.4457 315.5857 315.7256 315.8656 316.0056 316.1455 316.2855 316.4254 316.5653 316.7052 316.8451 316.9850 317.1249 317.2647 317.4046 317.5444 317.6842 317.8240 317.9639 318.1036 318.2434 318.3832 318.5229 318.6627 318.8024 318.9422 319.0819 319.2216 319.3612 319.5009 319.6406 319.7802 319.9199 320.0595 320.1991 320.3387 320.4783 320.6179 320.7575 320.8971 321.0366 321.1761 321.3157 321.4552 321.5947 321.7342 321.8737 322.0131 322.1526 322.2920 322.4315 322.5709 322.7103 322.8497 322.9891 323.1285 323.2678 323.4072 323.5465 323.6859 323.8252 323.9645 324.1038 324.2431 324.3824 324.5216 324.6609 324.8001 324.9393 --------------------------------------------------------------------------------------- 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 20231204T181555.3Z 8738.761056 499 1 30.15 2 67.93 197.22 94.79 291.20 0.00 1 0.00 1 0 17 2.812e-02 5.302e-03 1.702e-02 2.138e-03 4.183732e-01 1.042688e-03 5.084859e-03 -6.000000e+00 1.534227e-02 222.34 -2.00 3.57 -2.00 318.62 0.00e+00 5.646354e-04 2.187679e-04 1.498569e-04 -6.000000e+00 6.683448e-04 275.26 -2.00 20.00 -2.00 318.05 0.00e+00 2.638 0.000 302.00 328.00 3 3.7919e+00 2.8055e-03 3.3606e-02 -6.0000e+00 3.7788e-02 5.5105e-02 3.2201e-03 2.1417e-02 -6.0000e+00 4.1259e-02 -1.7724e-02 2.3387e-03 1.0813e-02 -6.0000e+00 2.4563e-02 1.1082e-03 1.2563e-03 3.2052e-03 -6.0000e+00 1.0590e-02 5.1538e-04 -5.2577e-04 5.6603e-05 3.1207e-05 -6.0000e+00 6.9813e-04 -9.5417e-03 -4.1819e-05 10 0 0 10 0 0 66.31 17.94 0.30 -0.00954 -0.01349 20 20 1 5 0 1507.050 22.61 15.00 17.94 11.30 21.6 1029.88 195 -14543 -13390 -9470 -16115 -11940 -4710 -1774 -3407 -695 -2547 617 -3654 -682 -1040 -1332 -286 -268 1063 1862 -4515 1938 -1089 -319 -626 -433 -363 -2737 -2526 -2010 31 12 2619 912 -2598 -3354 -2701 -2232 -904 -632 -746 -114 -395 -1080 908 -617 -645 -409 -674 -501 235 664 1824 1044 1048 582 1636 470 214 -431 957 1257 899 1816 702 674 109 -360 432 845 548 239 558 502 575 -120 -626 -1170 -874 -410 -427 -1032 -2094 -1242 -845 -729 -513 -135 170 622 592 247 251 111 220 7 470 554 382 -62 -81 -1010 -303 -136 -682 -541 -500 -1007 -248 63 114 128 414 -28 -405 228 156 267 166 -220 -320 24 275 512 801 213 597 760 463 488 278 2 -119 -136 -100 36 -28 201 -6 -118 -351 -403 -401 6 -631 -609 -399 -639 -431 -169 -71 -158 -34 -127 31 -226 -50 -83 -123 -28 -28 117 213 -568 225 -133 -42 -89 -65 -55 -406 -350 -278 5 2 491 187 -517 -641 -557 -537 -263 -202 -267 -43 -155 -434 357 -249 -260 -163 -272 -205 103 319 931 531 559 320 918 266 118 -234 523 740 590 1304 491 477 73 -238 302 631 455 209 486 418 443 -90 -424 -763 -603 -303 -362 -924 -2059 -1265 -913 -792 -588 -159 201 706 665 255 240 104 220 8 561 701 527 -85 -105 -1269 -345 -156 -860 -679 -654 -1397 -363 95 185 205 650 -44 -590 327 221 393 247 -349 -516 38 415 745 1169 300 791 1004 672 737 452 3 -196 -234 -174 66 -53 384 -12 -222 -629 -722 -692 SS 20231204T182204.4Z 8738.765329 503 1 29.68 2 68.27 198.70 95.81 292.25 0.00 1 0.00 1 0 17 5.143e-02 6.429e-03 4.247e-02 4.351e-03 4.255911e-01 2.126202e-03 5.142917e-03 -6.000000e+00 2.806233e-02 222.34 -2.00 3.57 -2.00 318.98 0.00e+00 -1.432063e-04 4.959022e-04 1.799872e-04 -6.000000e+00 1.222201e-03 275.26 -2.00 20.00 -2.00 318.43 0.00e+00 2.677 0.000 302.00 328.00 3 4.8049e+00 5.8058e-03 3.4157e-02 -6.0000e+00 6.9120e-02 9.2485e-02 6.7291e-03 2.1840e-02 -6.0000e+00 7.5462e-02 -1.0187e-02 5.1394e-03 1.1082e-02 -6.0000e+00 4.4926e-02 -7.0767e-03 2.7604e-03 3.2644e-03 -6.0000e+00 1.9372e-02 1.7500e-04 -1.0993e-04 1.1688e-04 3.2168e-05 -6.0000e+00 1.2767e-03 -9.6286e-03 5.7131e-06 11 8 0 10 0 0 99.96 17.93 0.17 -0.00963 -0.01045 14 14 1 5 0 2119.770 22.61 15.00 17.93 11.81 21.7 1033.41 10625 -16586 -1440 -27979 -15382 -29755 -13029 -14131 5221 -18320 -10950 4454 489 -8554 -2540 -1971 -2066 2175 -1068 3854 -5805 7073 -161 4880 -832 4297 1872 961 1413 -478 4864 -135 2102 -4143 -1786 -4266 -5217 -1610 1289 -1780 -1020 -586 -1741 -1043 3571 2078 -2402 -1224 -525 -541 25 1176 1486 550 -199 1643 809 1005 459 -622 1019 939 1916 1438 1038 -471 429 -1313 332 520 594 1222 973 747 595 -922 -1257 -2747 -548 -98 -1552 -1214 -1055 -1147 -1121 -883 -812 486 930 458 348 -573 -156 -135 272 223 33 1016 903 635 396 -1207 -581 -222 -925 -973 -397 -599 292 371 22 32 13 193 -551 208 -150 -588 110 -340 -156 610 426 308 1268 68 29 336 562 -24 -186 73 39 29 208 -177 330 -7 -99 -214 -188 -513 23 246 -576 -46 -1112 -475 -1010 -402 -507 177 -855 -477 174 23 -541 -181 -148 -168 181 -97 350 -599 633 -16 499 -97 506 230 113 154 -53 571 -19 325 -733 -289 -674 -906 -316 302 -517 -329 -200 -624 -376 1226 731 -880 -441 -189 -199 10 505 631 253 -96 805 412 510 227 -304 498 497 1115 927 697 -301 255 -784 208 376 476 946 811 597 440 -617 -763 -1623 -336 -69 -1277 -1040 -980 -1114 -1160 -978 -946 578 1102 527 371 -566 -141 -119 259 239 40 1284 1243 872 510 -1428 -667 -257 -1175 -1223 -517 -824 422 549 35 54 21 299 -846 296 -212 -869 172 -533 -247 931 674 469 1819 94 38 440 805 -38 -298 122 66 52 394 -335 634 -14 -196 -412 -364 -948 40 SS 20231204T183129.3Z 8738.771867 509 1 30.16 2 68.85 200.96 97.35 293.88 0.00 1 0.00 1 0 17 7.959e-02 6.487e-03 3.785e-02 3.317e-03 4.349688e-01 1.623581e-03 5.292819e-03 -6.000000e+00 4.342064e-02 222.34 -2.00 3.57 -2.00 318.92 0.00e+00 3.123524e-04 3.708457e-04 2.325871e-04 -6.000000e+00 1.891521e-03 275.26 -2.00 20.00 -2.00 318.37 0.00e+00 2.745 0.000 302.00 328.00 3 4.3901e+00 4.4034e-03 3.4999e-02 -6.0000e+00 1.0695e-01 7.4810e-02 5.0962e-03 2.2356e-02 -6.0000e+00 1.1677e-01 -1.2843e-02 3.8837e-03 1.1344e-02 -6.0000e+00 6.9518e-02 -2.3113e-03 2.0934e-03 3.3485e-03 -6.0000e+00 2.9969e-02 2.4738e-04 -4.3466e-04 8.8038e-05 3.3341e-05 -6.0000e+00 1.9758e-03 -9.0955e-03 -3.0945e-05 11 9 0 11 8 0 99.94 17.95 0.28 -0.00910 -0.01236 13 13 1 5 0 2319.370 22.61 14.99 17.95 12.17 21.7 1034.57 63289 -18275 9240 -28895 -24771 -21108 -21688 -22426 -3031 13083 -15308 6185 -16212 10573 -9305 5721 -4546 6162 -3233 1388 -3154 1601 376 -2060 3452 -1135 1987 -4468 1144 -5228 6998 -2585 6221 -155 583 -6802 -3252 -694 -1014 -523 -1624 -234 -1487 770 -1433 1327 -1958 1904 -647 -241 -1092 1156 438 425 755 188 1431 -478 1439 2 1058 1098 3121 1849 984 404 -552 -281 762 235 195 667 282 336 92 576 -1125 -744 -416 -244 -818 -196 -1982 -1086 -1055 -958 -787 -266 489 227 1014 -145 674 -244 954 -24 626 469 611 -47 -469 -962 -539 -164 -1201 -192 -1279 -550 -33 617 -254 429 471 334 -367 421 207 -138 169 -32 -547 -342 645 296 870 -88 656 402 526 307 167 3 -201 -235 -318 196 -27 260 -161 295 -219 -600 -265 819 -613 250 -1105 -801 -627 -695 -831 -106 423 -662 226 -877 532 -687 422 -402 524 -321 138 -339 161 39 -241 412 -150 260 -591 134 -657 868 -402 1003 -26 92 -1102 -553 -135 -243 -150 -527 -79 -533 274 -520 474 -721 673 -236 -90 -443 506 191 199 371 96 740 -252 721 1 530 594 1843 1223 682 265 -344 -173 494 164 152 541 229 260 71 395 -710 -450 -267 -181 -651 -174 -1955 -1036 -1071 -1043 -901 -313 574 258 1057 -151 642 -229 955 -26 749 596 844 -65 -610 -1214 -616 -188 -1523 -239 -1682 -757 -48 918 -414 722 780 524 -536 605 295 -205 252 -51 -882 -535 968 456 1268 -124 867 533 765 464 272 5 -353 -406 -585 379 -50 497 -311 581 -411 -1073 -455