File name: Pandora69s1_NewBrunswickNJ_20230716_L2Fit_fout2c23d20210421p1-8.txt
File generation date: 20230717T055338.5Z
Data description: Level 2 spectral fitting results file
Data file version: fout2c23p1-8
Local principal investigator: Lukas Valin
Network principal investigator: Alexander Cede
Instrument type: Pandora
Instrument number: 69
Spectrometer number: 1
Processing software version used: BlickP v1.8.49
Instrument operation file used: Pandora69_OF_v2d20210223.txt
Instrument calibration file used: Pandora69s1_CF_v23d20210421.txt
Level 1 file used: Pandora69s1_NewBrunswickNJ_20230716_L1_smca1c23d20210421p1-8.txt
Full location name: Rutgers NJDEP (USEPA AQS ID 34-023-0011)
Short location name: NewBrunswickNJ
Country of location: USA
Location latitude [deg]: 40.4622
Location longitude [deg]: -74.4294
Location altitude [m]: 19
Local noon date: 20230716
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 s13, field cal for NewBrunswick, revisited HCHO
First and last pixel inside fitting window: 259 400
Nominal wavelengths inside fitting window [nm]: 305.1248 305.2657 305.4066 305.5475 305.6883 305.8292 305.9700 306.1109 306.2517 306.3925 306.5333 306.6741 306.8148 306.9556 307.0963 307.2371 307.3778 307.5185 307.6592 307.7999 307.9406 308.0812 308.2219 308.3625 308.5032 308.6438 308.7844 308.9250 309.0656 309.2062 309.3467 309.4873 309.6278 309.7683 309.9089 310.0494 310.1898 310.3303 310.4708 310.6113 310.7517 310.8921 311.0326 311.1730 311.3134 311.4537 311.5941 311.7345 311.8748 312.0152 312.1555 312.2958 312.4361 312.5764 312.7167 312.8570 312.9972 313.1375 313.2777 313.4179 313.5581 313.6983 313.8385 313.9787 314.1189 314.2590 314.3991 314.5393 314.6794 314.8195 314.9596 315.0997 315.2397 315.3798 315.5198 315.6599 315.7999 315.9399 316.0799 316.2199 316.3598 316.4998 316.6397 316.7797 316.9196 317.0595 317.1994 317.3393 317.4792 317.6190 317.7589 317.8987 318.0386 318.1784 318.3182 318.4580 318.5978 318.7375 318.8773 319.0170 319.1568 319.2965 319.4362 319.5759 319.7156 319.8552 319.9949 320.1345 320.2742 320.4138 320.5534 320.6930 320.8326 320.9722 321.1117 321.2513 321.3908 321.5303 321.6698 321.8093 321.9488 322.0883 322.2278 322.3672 322.5067 322.6461 322.7855 322.9249 323.0643 323.2037 323.3430 323.4824 323.6217 323.7611 323.9004 324.0397 324.1790 324.3183 324.4575 324.5968 324.7360 324.8753
---------------------------------------------------------------------------------------
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
---------------------------------------------------------------------------------------
SS 20230716T135503.7Z 8597.579904 223 1 5.16 2 44.14 101.09 30.68 104.50 0.00 1 0.00 1 0 17 9.911e-03 7.232e-03 5.860e-03 3.889e-03 1.813093e-01 2.062879e-03 2.193084e-03 -6.000000e+00 5.886038e-03 232.84 -2.00 3.50 -2.00 315.11 0.00e+00 -2.864136e-05 1.610878e-04 1.805522e-05 -6.000000e+00 2.310618e-04 289.82 -2.00 20.00 -2.00 314.82 0.00e+00 1.391 0.000 302.00 328.00 3 4.1333e+00 5.3688e-03 1.6914e-02 -6.0000e+00 1.4914e-02 5.7960e-02 5.9290e-03 9.4737e-03 -6.0000e+00 1.6155e-02 -5.3710e-03 3.3811e-03 4.6040e-03 -6.0000e+00 9.4570e-03 -1.9965e-03 1.4803e-03 1.2747e-03 -6.0000e+00 3.9371e-03 1.7364e-03 4.1001e-04 1.1396e-04 1.7612e-05 -6.0000e+00 2.8066e-04 1.5526e-02 -8.3914e-06 12 9 1 12 1 1 999.00 22.97 0.35 0.01553 0.01861 2 6 1 5 0 1661.110 26.11 20.01 22.97 31.97 28.1 1064.62 -3716 -533 -1698 -127 -363 -2688 1532 -1358 2220 -1835 1313 -1853 -205 666 1011 201 87 1694 417 2083 1042 412 837 881 313 609 467 -734 1611 461 847 272 1556 -2895 -891 -2422 -2351 -673 -916 -997 -1588 122 -712 -677 -132 -350 202 -798 -272 327 1265 -1 1595 332 622 590 128 749 1253 1094 1026 1298 968 541 -659 -17 -115 403 481 256 55 -690 297 532 65 -131 -1358 -251 -916 -1043 -1581 -2302 -1746 -133 -785 -235 814 98 313 86 539 -326 237 508 -86 1157 1057 601 -843 459 -636 -219 430 -397 -1246 -415 -537 -253 228 -4 383 36 175 -484 -148 441 126 810 -351 -286 -511 -80 22 436 743 774 740 637 591 232 -118 65 -181 324 -40 -286 663 59 -108 0 -1046 -695 -1226 -177 -552 -37 -99 -757 447 -441 743 -635 437 -663 -78 275 448 93 41 802 204 1021 513 204 427 463 169 332 248 -368 762 222 421 141 829 -1549 -450 -1285 -1359 -430 -658 -770 -1273 98 -584 -553 -107 -281 158 -614 -205 250 997 0 1336 287 542 519 113 655 1085 959 931 1229 955 540 -637 -15 -101 355 451 256 57 -730 308 535 62 -116 -1157 -211 -813 -992 -1605 -2455 -1911 -150 -926 -286 1009 123 388 103 610 -351 243 522 -93 1313 1273 758 -1075 565 -759 -251 490 -475 -1558 -530 -706 -340 313 -6 548 51 238 -629 -183 534 156 1031 -470 -392 -705 -109 30 563 931 941 887 779 746 302 -157 87 -246 452 -58 -421 986 90 -163 0 -1525 -997
SS 20230716T140238.1Z 8597.585163 229 1 6.23 2 42.73 102.63 29.32 106.34 0.00 1 0.00 1 0 17 8.830e-03 6.833e-03 5.009e-03 3.421e-03 1.756952e-01 1.815155e-03 2.108274e-03 -6.000000e+00 5.244635e-03 232.84 -2.00 3.50 -2.00 314.97 0.00e+00 2.850931e-04 1.381470e-04 1.695881e-05 -6.000000e+00 2.058394e-04 289.82 -2.00 20.00 -2.00 314.70 0.00e+00 1.360 0.000 302.00 328.00 3 4.1331e+00 4.7183e-03 1.6453e-02 -6.0000e+00 1.3289e-02 5.2118e-02 5.2090e-03 9.1302e-03 -6.0000e+00 1.4394e-02 -4.4054e-03 2.9677e-03 4.4320e-03 -6.0000e+00 8.4260e-03 -1.5694e-03 1.2941e-03 1.2262e-03 -6.0000e+00 3.5087e-03 1.8132e-03 4.9369e-04 9.9940e-05 1.6929e-05 -6.0000e+00 2.5002e-04 1.4277e-02 2.8013e-07 12 9 1 12 1 1 999.00 22.98 0.20 0.01428 0.01799 2 5 1 5 0 2018.250 26.11 20.00 22.98 32.62 28.3 1071.17 769 127 -1849 -1787 -478 -158 -69 -493 247 -1708 -976 -451 -1266 2758 -1320 -160 -230 737 -25 292 838 1945 912 -134 2101 -804 944 -896 1508 1191 970 2719 -1185 -1355 -1350 -916 -2379 -794 -64 -883 -130 -1002 -347 -600 151 -1024 -993 -766 -654 862 506 983 1097 998 501 962 -404 754 385 715 953 810 2027 -56 -181 -206 -153 258 696 347 356 478 -584 238 -382 49 -231 -778 -1244 -554 -1019 -2256 -1251 -1178 -1089 -31 310 -286 -19 459 458 362 879 1276 -492 412 695 850 -769 -232 73 -243 933 -179 -789 -631 -420 -303 75 -59 271 -97 -336 -221 361 512 464 380 -253 -409 -290 123 468 -229 731 194 1057 667 618 116 -363 -214 345 -116 73 337 30 -406 190 -75 -507 -732 258 44 -635 -558 -138 -46 -22 -166 88 -617 -345 -166 -508 1167 -617 -77 -112 363 -12 150 427 983 479 -73 1150 -455 512 -463 735 588 496 1426 -660 -735 -702 -492 -1405 -517 -46 -691 -104 -825 -287 -496 124 -837 -794 -597 -502 665 405 810 931 867 440 851 -362 664 338 633 872 774 1997 -57 -175 -188 -136 229 653 348 371 503 -611 241 -367 44 -197 -660 -1111 -527 -1032 -2406 -1365 -1330 -1284 -37 384 -359 -23 547 518 388 899 1306 -531 468 836 1065 -976 -286 87 -277 1058 -213 -980 -804 -549 -404 102 -83 385 -136 -455 -285 443 617 568 481 -336 -557 -396 167 622 -296 912 236 1260 812 776 150 -479 -285 464 -161 104 490 45 -608 282 -110 -731 -1043
SS 20230716T140634.0Z 8597.587894 232 1 13.01 2 42.00 103.45 28.62 107.35 0.00 1 0.00 1 0 17 7.105e-03 5.898e-03 1.290e-03 8.802e-04 1.733347e-01 4.664190e-04 2.098254e-03 -6.000000e+00 4.218793e-03 232.84 -2.00 3.50 -2.00 314.99 0.00e+00 1.458242e-04 3.551567e-05 1.654015e-05 -6.000000e+00 1.656512e-04 289.82 -2.00 20.00 -2.00 314.71 0.00e+00 1.344 0.000 302.00 328.00 3 2.5548e+00 1.2125e-03 1.6278e-02 -6.0000e+00 1.0689e-02 4.9184e-02 1.3392e-03 9.0579e-03 -6.0000e+00 1.1580e-02 -5.9927e-03 7.6277e-04 4.4037e-03 -6.0000e+00 6.7785e-03 -1.9891e-03 3.3223e-04 1.2191e-03 -6.0000e+00 2.8216e-03 8.9943e-03 1.6243e-04 2.5730e-05 1.6791e-05 -6.0000e+00 2.0120e-04 1.5974e-02 -3.2837e-05 12 1 1 12 1 1 99.93 22.97 0.24 0.01597 0.01719 33 27 1 5 0 377.120 26.11 20.00 22.97 33.23 28.3 1075.55 -619 -560 -1279 -2052 -883 -151 461 1360 533 -811 -488 -189 -769 201 344 59 802 364 1027 709 100 1126 946 351 335 386 -199 209 931 1342 736 1384 44 -1707 -973 -1432 -1463 -736 -680 -1117 -1024 -535 -933 -821 -489 -855 -536 -624 -468 708 498 648 1038 462 616 965 501 495 751 921 1170 1062 1277 158 -278 390 -275 434 674 135 211 -42 329 225 87 118 -237 33 -383 -1038 -922 -1818 -1647 -915 -940 -498 -193 146 99 123 449 242 712 726 415 638 959 190 -379 -160 -411 -432 329 -133 -980 -605 -472 -482 141 82 -34 -130 -19 -164 307 447 458 414 60 -490 -97 -67 199 203 630 193 1334 838 322 269 -127 -94 189 289 50 -244 324 -222 -317 -273 -383 -528 -217 -188 -451 -645 -250 -43 142 459 195 -300 -176 -72 -301 87 162 29 400 174 500 354 50 557 487 189 186 217 -108 103 463 673 387 708 24 -903 -518 -747 -861 -471 -481 -887 -851 -428 -762 -665 -392 -717 -432 -484 -354 539 402 547 860 399 539 851 446 434 649 808 1078 992 1265 160 -264 355 -239 376 642 135 225 -43 350 227 85 103 -202 28 -350 -1000 -923 -1978 -1776 -1038 -1091 -609 -245 180 126 145 507 263 715 732 451 720 1131 243 -489 -196 -498 -491 369 -161 -1219 -783 -608 -646 197 114 -48 -180 -26 -215 369 546 574 524 79 -673 -133 -91 261 263 772 237 1576 1028 398 348 -171 -128 250 400 70 -355 481 -337 -475 -398 -544 -764