File name: Pandora166s1_PhiladelphiaPA_20220103_L2Fit_fwvt1c1d20210412p1-8.txt File generation date: 20251125T094100.2Z Data description: Level 2 spectral fitting results file Data file version: fwvt1c1p1-8 Local principal investigator: Lukas Valin Network principal investigator: Alexander Cede Instrument type: Pandora Instrument number: 166 Spectrometer number: 1 Processing software version used: BlickP v1.8.75 Instrument operation file used: Pandora166_OF_v3d20210217.txt Instrument calibration file used: Pandora166s1_CF_v1d20210412.txt Level 1 file used: Pandora166s1_PhiladelphiaPA_20220103_L1_smca1c1d20210412p1-8.txt Full location name: Philadelphia, Pennsylvania; Philadelphia-Camden-Wilmington; North East Waste (NE Short location name: PhiladelphiaPA Country of location: United States Location latitude [deg]: 39.9919 Location longitude [deg]: -75.0811 Location altitude [m]: 6 Local noon date: 20220103 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: 1644 1794 Nominal wavelengths inside fitting window [nm]: 492.0331 492.1540 492.2749 492.3958 492.5166 492.6375 492.7583 492.8791 492.9999 493.1206 493.2414 493.3621 493.4829 493.6036 493.7242 493.8449 493.9655 494.0862 494.2068 494.3274 494.4480 494.5685 494.6891 494.8096 494.9301 495.0506 495.1711 495.2916 495.4120 495.5324 495.6529 495.7732 495.8936 496.0140 496.1343 496.2547 496.3750 496.4953 496.6155 496.7358 496.8560 496.9763 497.0965 497.2166 497.3368 497.4570 497.5771 497.6972 497.8173 497.9374 498.0575 498.1775 498.2976 498.4176 498.5376 498.6576 498.7776 498.8975 499.0174 499.1374 499.2572 499.3771 499.4970 499.6168 499.7367 499.8565 499.9763 500.0961 500.2158 500.3356 500.4553 500.5750 500.6947 500.8144 500.9340 501.0537 501.1733 501.2929 501.4125 501.5321 501.6516 501.7711 501.8907 502.0102 502.1297 502.2491 502.3686 502.4880 502.6074 502.7268 502.8462 502.9656 503.0849 503.2043 503.3236 503.4429 503.5621 503.6814 503.8006 503.9199 504.0391 504.1583 504.2774 504.3966 504.5157 504.6349 504.7540 504.8731 504.9921 505.1112 505.2302 505.3492 505.4682 505.5872 505.7062 505.8251 505.9441 506.0630 506.1819 506.3007 506.4196 506.5385 506.6573 506.7761 506.8949 507.0136 507.1324 507.2511 507.3699 507.4886 507.6073 507.7259 507.8446 507.9632 508.0818 508.2004 508.3190 508.4376 508.5561 508.6746 508.7932 508.9116 509.0301 509.1486 509.2670 509.3854 509.5038 509.6222 509.7406 509.8590 509.9773 --------------------------------------------------------------------------------------- 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: Water vapor slant column amount [moles per square meter], -9e99=fitting not successful Column 23: Independent uncertainty of water vapor 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 water vapor 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 water vapor 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 water vapor 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: Water vapor effective temperature [K] Column 28: Independent uncertainty of water vapor 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 water vapor 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 water vapor 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 water vapor fitting wavelength [nm], -9=fitting not successful Column 32: Diffuse correction applied before fitting at effective fitting wavelength for water vapor [%], 0=no diffuse correction applied or fitting not requested, >0=measured diffuse correction, <0=(negative value of) calculated diffuse correction Column 33: Nitrogen dioxide slant column amount [moles per square meter], -9e99=fitting not successful Column 34: Independent uncertainty of nitrogen 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 nitrogen 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 nitrogen 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 nitrogen 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: Nitrogen dioxide effective temperature [K] Column 39: Independent uncertainty of nitrogen 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 nitrogen 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 nitrogen 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 nitrogen dioxide fitting wavelength [nm], -9=fitting not successful Column 43: Diffuse correction applied before fitting at effective fitting wavelength for nitrogen dioxide [%], 0=no diffuse correction applied or fitting not requested, >0=measured diffuse correction, <0=(negative value of) calculated diffuse correction Column 44: Ozone slant column amount [moles per square meter], -9e99=fitting not successful Column 45: 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 46: 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 47: 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 48: 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 49: Ozone effective temperature [K] Column 50: 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 51: 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 52: 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 53: Effective ozone fitting wavelength [nm], -9=fitting not successful Column 54: 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 55: Molecular scattering air mass factor used for molecular scattering subtraction before the fitting Column 56: Estimated uncertainty of molecular scattering air mass factor, -9=molecular scattering was not subtracted before the fitting Column 57: Lower limit used for wavelength scaling [nm] Column 58: Upper limit used for wavelength scaling [nm] Column 59: Order of smoothing polynomial Column 60: Smoothing polynomial coefficient, order 0 Column 61: Independent uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 62: Structured uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 63: Common uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 64: rms-based uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 65: Smoothing polynomial coefficient, order 1 Column 66: Independent uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 67: Structured uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 68: Common uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 69: rms-based uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 70: Mean value of measured data inside fitting window [same units as measurements] Column 71: Wavelength change polynomial coefficient, order 0 Column 72: Independent uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 73: Structured uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 74: Common uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 75: rms-based uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 76: L1 based wavelength change coefficient, order 0, -9=no L1 wavelength change determination Column 77: L1 based wavelength change coefficient, order 1, -9=no L1 wavelength change determination Column 78: 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 79: 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 80: 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 81: 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 82: 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 83: 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 84: Atmospheric variability [%], 999=no atmospheric variability was determined Column 85: Wavelength effective temperature [°C], 999=no effective temperature given Column 86: Estimated average residual stray light level [%] (only valid for stray light correction methods 2 and higher) Column 87: Retrieved wavelength shift from L1 data [nm], -9=no wavelength change determination Column 88: Retrieved total wavelength shift [nm], -9=no wavelength change fitting Column 89: Number of bright count cycles Column 90: Number of dark count cycles Column 91: Effective position of filterwheel #1, 0=filterwheel not used, 1-9 are valid positions Column 92: Effective position of filterwheel #2, 0=filterwheel not used, 1-9 are valid positions Column 93: Sum over 2^i, 0=spectra were interpolated in time, 1=spectra are corrected for off-target signal Column 94: Integration time [ms] Column 95: Temperature at electronics board [°C], 999=no temperature signal Column 96: Spectrometer control temperature [°C], 999=no temperature signal Column 97: Auxiliary spectrometer temperature [°C], 999=no temperature signal Column 98: Temperature in head sensor [°C], 999=no temperature signal Column 99: Humidity in head sensor [%], -9=no humidity signal Column 100: Pressure in head sensor [hPa], -9=no pressure signal Columns 101-251: Unweighted slant column residuals for each pixel inside the fitting window multiplied by 1e5, 9e5=pixel was not used for fitting Columns 252-402: 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 --------------------------------------------------------------------------------------- SQ 20220103T213319.4Z 8038.898141 318 1 47.34 2 88.47 238.17 80.21 224.56 0.00 1 0.00 1 0 15 5.184e-03 5.134e-03 1.559e-04 1.554e-04 4.317368e+03 9.210899e+00 3.017202e+02 -6.000000e+00 2.779081e+02 276.80 -2.00 20.00 -2.00 505.39 0.00e+00 4.124256e-03 4.224631e-06 2.226078e-04 -6.000000e+00 1.425306e-04 276.80 -2.00 20.00 -2.00 501.03 0.00e+00 1.238615e+00 2.963146e-03 4.981020e-02 -6.000000e+00 1.016317e-01 223.51 -2.00 3.49 -2.00 502.54 0.00e+00 17.410 0.096 489.00 513.00 1 1.0145e+00 2.2642e-04 2.9235e-02 -6.0000e+00 7.7777e-03 -1.6453e-02 5.7520e-05 1.8018e-03 -6.0000e+00 2.0170e-03 5.1088e-02 7.2069e-04 1.4782e-05 8.4230e-05 -6.0000e+00 5.0389e-04 5.7359e-03 -1.5110e-04 12 8 8 10 0 0 98.27 20.15 2.32 0.00574 0.01073 534 106 1 4 0 88.630 26.51 19.99 20.15 5.32 10.2 907.14 -1265 -1082 -704 -190 156 533 905 883 855 664 443 213 -179 -380 -472 -490 -340 -21 237 266 268 415 377 437 580 639 593 448 -15 -504 -702 -729 -604 -28 229 168 185 134 -94 -261 -171 -47 121 363 541 550 563 459 226 -88 -607 -958 -1160 -1078 -643 -343 7 204 312 331 419 574 521 552 481 324 164 3 -223 -429 -554 -479 -377 -190 -15 -180 -333 -460 -625 -690 -724 -670 -648 -523 -287 -307 -313 -276 -299 -255 -84 179 427 495 491 497 302 -63 -438 -784 -1088 -1008 -717 -297 -153 -209 -241 -387 -479 -324 -218 -28 308 571 779 939 1103 1072 906 662 525 347 282 412 356 309 410 382 356 449 224 -78 -221 -412 -683 -493 -277 -247 -22 194 233 276 374 447 384 293 -28 -387 -669 -764 -789 -1133 -978 -627 -175 148 496 856 847 823 635 418 199 -166 -362 -447 -465 -313 -20 228 261 264 410 374 435 559 620 575 433 -15 -487 -663 -688 -578 -28 224 167 185 132 -92 -255 -167 -47 120 367 554 547 561 472 230 -88 -593 -916 -1096 -1022 -621 -342 7 204 314 333 422 580 530 564 490 327 163 3 -218 -430 -551 -476 -378 -185 -15 -178 -325 -460 -618 -679 -714 -665 -625 -508 -281 -304 -311 -274 -296 -253 -84 180 437 513 508 506 302 -62 -420 -767 -1055 -984 -690 -295 -151 -210 -242 -384 -488 -331 -218 -28 319 600 824 1000 1143 1120 949 688 558 364 294 430 376 316 422 394 367 464 232 -80 -232 -428 -708 -514 -284 -259 -23 209 253 293 399 479 413 315 -30 -419 -711 -805 -837 SQ 20220103T213635.6Z 8038.900413 322 1 49.09 2 89.01 238.71 80.63 225.15 0.00 1 0.00 1 0 15 6.282e-03 6.184e-03 3.741e-04 3.721e-04 4.899338e+03 2.223588e+01 3.384436e+02 -6.000000e+00 3.367569e+02 276.80 -2.00 20.00 -2.00 505.42 0.00e+00 5.049509e-03 1.009293e-05 2.709556e-04 -6.000000e+00 1.727264e-04 276.80 -2.00 20.00 -2.00 501.15 0.00e+00 1.219186e+00 7.054650e-03 5.932150e-02 -6.000000e+00 1.231894e-01 223.51 -2.00 3.49 -2.00 502.63 0.00e+00 18.791 0.120 489.00 513.00 1 2.4431e+00 5.3880e-04 3.2981e-02 -6.0000e+00 9.4270e-03 -2.4312e-02 1.3646e-04 2.0428e-03 -6.0000e+00 2.4455e-03 9.9436e-03 1.5272e-03 3.5518e-05 1.0146e-04 -6.0000e+00 6.1058e-04 -9.5428e-03 -8.0465e-05 12 8 8 10 0 0 88.66 20.15 2.97 -0.00954 0.00105 175 35 1 4 0 280.290 26.51 20.00 20.15 5.45 10.3 908.18 -1585 -1376 -883 -231 116 552 1048 1116 1072 801 542 286 -170 -464 -562 -589 -415 -3 267 351 253 405 397 507 680 756 683 554 -35 -577 -938 -996 -753 -152 323 199 286 174 -71 -270 -210 -60 169 489 644 748 768 658 319 -82 -742 -1075 -1348 -1234 -794 -374 31 293 346 430 552 676 673 724 611 487 219 -58 -285 -561 -708 -623 -477 -360 -154 -252 -349 -517 -778 -840 -886 -852 -791 -643 -363 -426 -381 -410 -318 -317 -97 234 499 663 590 577 383 -108 -540 -959 -1300 -1143 -867 -350 -14 -132 -219 -453 -529 -413 -271 36 419 811 939 1131 1343 1262 967 827 610 430 360 429 407 409 408 405 389 467 231 -135 -303 -556 -814 -601 -317 -222 -18 234 286 280 345 496 483 378 6 -407 -797 -963 -869 -1333 -1169 -770 -208 108 524 960 1035 997 741 495 258 -161 -433 -522 -548 -371 -3 248 333 242 388 381 490 661 740 669 541 -34 -542 -864 -917 -704 -146 320 192 275 175 -70 -268 -208 -58 163 481 641 751 773 658 316 -80 -710 -1056 -1311 -1203 -753 -366 31 298 354 440 565 697 699 754 635 502 223 -59 -284 -555 -695 -611 -471 -360 -156 -255 -349 -510 -760 -817 -864 -836 -782 -640 -364 -432 -388 -417 -323 -322 -99 230 500 671 596 574 374 -109 -531 -927 -1247 -1102 -855 -341 -14 -135 -223 -456 -529 -415 -276 37 442 867 970 1176 1407 1335 1026 871 634 442 384 439 421 426 427 426 409 492 243 -141 -312 -569 -832 -619 -333 -238 -19 257 319 314 372 538 526 412 6 -450 -867 -1040 -945 SQ 20220103T213948.3Z 8038.902642 326 1 49.98 2 89.53 239.23 81.05 225.72 0.00 1 0.00 1 0 15 7.791e-03 7.658e-03 7.708e-04 7.638e-04 5.495996e+03 4.591019e+01 3.710903e+02 -6.000000e+00 4.176353e+02 276.80 -2.00 20.00 -2.00 505.45 0.00e+00 6.240169e-03 2.081184e-05 3.336829e-04 -6.000000e+00 2.141948e-04 276.80 -2.00 20.00 -2.00 501.33 0.00e+00 1.188647e+00 1.441633e-02 7.136051e-02 -6.000000e+00 1.527380e-01 223.51 -2.00 3.49 -2.00 502.77 0.00e+00 19.919 0.144 489.00 513.00 1 3.5720e+00 1.1023e-03 3.6534e-02 -6.0000e+00 1.1689e-02 -3.3677e-02 2.7726e-04 2.2878e-03 -6.0000e+00 3.0315e-03 2.7001e-03 1.0195e-03 7.2879e-05 1.2205e-04 -6.0000e+00 7.5724e-04 5.9885e-03 -1.5627e-04 12 8 8 10 0 0 98.63 20.15 3.95 0.00599 0.01306 54 11 1 4 0 925.430 26.51 20.00 20.15 5.58 10.2 908.41 -2170 -1805 -1237 -227 313 771 1148 1230 1320 902 502 31 -167 -600 -703 -820 -659 -82 254 492 567 643 600 609 816 923 978 457 162 -719 -1141 -1051 -987 -126 146 521 433 290 -98 -275 -180 -120 279 545 701 832 1000 661 332 -129 -781 -1344 -1643 -1512 -917 -421 94 402 606 501 708 756 942 846 682 601 283 62 -251 -593 -819 -846 -573 -376 -67 -222 -498 -788 -981 -1194 -1175 -1050 -995 -821 -625 -436 -337 -454 -372 -402 -90 226 614 790 648 846 442 -56 -611 -1155 -1512 -1552 -1138 -446 100 -123 -169 -638 -655 -504 -356 -4 491 1042 1142 1473 1613 1406 1406 835 904 446 284 482 477 610 455 495 406 482 201 -64 -410 -763 -1066 -748 -400 -356 55 309 417 347 376 605 434 327 -8 -591 -962 -1198 -1063 -1825 -1531 -1078 -203 272 684 1037 1125 1209 821 452 28 -148 -525 -612 -717 -582 -74 234 459 532 604 567 579 783 891 946 442 154 -669 -1043 -958 -915 -120 144 492 408 272 -92 -256 -168 -113 267 530 692 829 997 657 327 -125 -744 -1257 -1523 -1407 -870 -412 94 409 620 513 726 781 979 884 712 621 290 63 -252 -591 -810 -836 -570 -378 -68 -227 -503 -785 -967 -1174 -1157 -1039 -992 -823 -632 -444 -344 -464 -378 -410 -92 222 612 797 653 838 431 -56 -605 -1125 -1460 -1509 -1130 -432 99 -125 -171 -643 -655 -506 -364 -4 520 1116 1236 1605 1773 1485 1486 881 989 483 306 520 520 637 479 523 430 511 213 -67 -427 -789 -1101 -777 -423 -385 60 342 467 392 427 690 497 374 -9 -660 -1057 -1307 -1168