File name: Pandora171s1_Beijing-RADI_20231012_L2Fit_fwvt1c1d20210725p1-8.txt File generation date: 20231013T004254.8Z Data description: Level 2 spectral fitting results file Data file version: fwvt1c1p1-8 Local principal investigator: Li Kaitao Network principal investigator: Alexander Cede Instrument type: Pandora Instrument number: 171 Spectrometer number: 1 Processing software version used: BlickP v1.8.49 Instrument operation file used: Pandora171_OF_v3d20200518.txt Instrument calibration file used: Pandora171s1_CF_v1d20210725.txt Level 1 file used: Pandora171s1_Beijing-RADI_20231012_L1_smca1c1d20210725p1-8.txt Full location name: Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing Short location name: Beijing-RADI Country of location: China Location latitude [deg]: 40.0048 Location longitude [deg]: 116.3786 Location altitude [m]: 59 Local noon date: 20231012 Notes on s-number (L1 configuration): Corrections NOT applied although requested by the s-number are latency correction, matrix method stray light correction (replaced by simple method) Data caveats: None First and last pixel inside fitting window: 1645 1795 Nominal wavelengths inside fitting window [nm]: 492.0172 492.1381 492.2591 492.3800 492.5009 492.6218 492.7427 492.8635 492.9843 493.1052 493.2259 493.3467 493.4675 493.5882 493.7090 493.8297 493.9504 494.0710 494.1917 494.3123 494.4330 494.5536 494.6742 494.7947 494.9153 495.0358 495.1563 495.2768 495.3973 495.5178 495.6382 495.7587 495.8791 495.9995 496.1199 496.2402 496.3606 496.4809 496.6012 496.7215 496.8417 496.9620 497.0822 497.2024 497.3226 497.4428 497.5630 497.6831 497.8033 497.9234 498.0435 498.1635 498.2836 498.4036 498.5237 498.6437 498.7636 498.8836 499.0036 499.1235 499.2434 499.3633 499.4832 499.6030 499.7229 499.8427 499.9625 500.0823 500.2021 500.3218 500.4416 500.5613 500.6810 500.8006 500.9203 501.0399 501.1596 501.2792 501.3988 501.5183 501.6379 501.7574 501.8769 501.9964 502.1159 502.2354 502.3548 502.4742 502.5936 502.7130 502.8324 502.9518 503.0711 503.1904 503.3097 503.4290 503.5482 503.6675 503.7867 503.9059 504.0251 504.1443 504.2634 504.3825 504.5016 504.6207 504.7398 504.8589 504.9779 505.0969 505.2159 505.3349 505.4539 505.5728 505.6917 505.8106 505.9295 506.0484 506.1673 506.2861 506.4049 506.5237 506.6425 506.7612 506.8800 506.9987 507.1174 507.2361 507.3547 507.4734 507.5920 507.7106 507.8292 507.9478 508.0663 508.1849 508.3034 508.4219 508.5403 508.6588 508.7772 508.8957 509.0141 509.1324 509.2508 509.3691 509.4875 509.6058 509.7241 509.8423 509.9606 --------------------------------------------------------------------------------------- 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: 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: Temperature at electronics board 1 [°C], 999=no temperature signal Column 99: Temperature at electronics board 2 [°C], 999=no temperature signal Column 100: Humidity in head sensor [%], -9=no humidity signal Column 101: Pressure in head sensor [hPa], -9=no pressure signal Columns 102-252: Unweighted slant column residuals for each pixel inside the fitting window multiplied by 1e5, 9e5=pixel was not used for fitting Columns 253-403: 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 20231012T034951.7Z 8685.159626 247 1 48.52 2 47.36 176.20 39.37 217.38 0.00 1 0.00 1 0 15 4.247e-03 4.240e-03 3.349e-04 3.341e-04 2.240244e+03 1.937076e+01 1.607801e+02 -6.000000e+00 2.284012e+02 285.12 -2.00 20.00 -2.00 505.25 0.00e+00 1.174631e-03 9.084611e-06 6.953751e-05 -6.000000e+00 1.179500e-04 285.12 -2.00 20.00 -2.00 500.40 0.00e+00 -3.916035e-02 6.490314e-03 1.856985e-02 -6.000000e+00 8.234894e-02 226.43 -2.00 3.49 -2.00 501.99 0.00e+00 1.474 0.000 489.00 513.00 1 5.1259e+00 4.9652e-04 3.0761e-03 -6.0000e+00 6.3018e-03 -2.7000e-02 1.3064e-04 4.0516e-04 -6.0000e+00 1.6588e-03 9.3029e-03 -3.0854e-03 2.8314e-05 2.7670e-05 -6.0000e+00 3.6917e-04 -5.5221e-02 1.9159e-04 11 8 0 10 0 0 99.96 22.68 3.52 -0.05522 -0.07662 133 31 1 2 0 364.760 20.02 22.68 26.95 25.79 25.68 34.1 1031.49 900000 -964 -539 -177 282 703 797 880 723 607 504 205 -31 -230 -297 -404 -178 29 -16 7 259 113 52 187 69 84 239 226 -367 -420 -780 -957 -735 -214 -188 72 295 114 -41 -101 42 -79 217 449 634 624 610 527 265 -27 -392 -704 -1035 -793 -654 -435 -94 248 335 446 398 384 454 513 511 381 304 73 -150 -453 -430 -639 -477 -260 -481 -549 -457 -396 -563 -851 -683 -638 -585 -305 -136 -179 -11 -2 52 -60 56 267 360 267 511 564 252 144 -216 -630 -759 -685 -426 -294 -436 -108 191 -32 -154 -161 -323 -156 138 291 378 625 727 660 482 501 342 210 290 150 -121 285 335 195 378 264 -27 35 -25 -318 -446 -219 -529 -394 25 170 97 291 307 309 330 264 72 -50 -43 -332 -510 900000 -951 -532 -178 290 734 837 926 768 642 532 215 -32 -238 -303 -411 -182 29 -16 7 274 119 55 198 73 89 251 237 -387 -435 -794 -966 -735 -218 -196 76 312 119 -43 -103 43 -81 223 466 656 652 638 548 272 -27 -399 -695 -1017 -770 -643 -435 -96 254 344 456 406 393 461 523 520 385 307 73 -150 -449 -419 -621 -467 -257 -477 -550 -454 -390 -546 -820 -656 -617 -565 -295 -133 -175 -11 -2 51 -59 56 261 358 267 511 563 249 139 -207 -597 -707 -639 -403 -284 -429 -106 190 -32 -149 -156 -313 -154 136 287 377 624 729 659 482 502 338 208 282 146 -118 276 325 190 365 255 -26 34 -23 -299 -422 -207 -500 -378 24 165 95 287 300 303 327 262 71 -48 -40 -314 -483 SQ 20231012T035907.1Z 8685.166055 251 1 50.43 2 47.29 179.32 40.49 220.43 0.00 1 0.00 1 0 15 4.605e-03 4.594e-03 1.561e-04 1.558e-04 2.149729e+03 8.973526e+00 1.513077e+02 -6.000000e+00 2.476787e+02 285.12 -2.00 20.00 -2.00 505.26 0.00e+00 1.324131e-03 4.237340e-06 7.786190e-05 -6.000000e+00 1.279047e-04 285.12 -2.00 20.00 -2.00 500.45 0.00e+00 -7.158344e-02 3.017793e-03 1.936398e-02 -6.000000e+00 8.929379e-02 226.43 -2.00 3.49 -2.00 502.04 0.00e+00 1.472 0.000 489.00 513.00 1 3.6255e+00 2.3083e-04 3.2098e-03 -6.0000e+00 6.8334e-03 -3.3700e-02 6.0725e-05 4.0367e-04 -6.0000e+00 1.7985e-03 4.1745e-02 -3.0856e-03 1.3250e-05 2.8885e-05 -6.0000e+00 4.0032e-04 -5.5721e-02 1.9968e-04 12 8 8 10 0 0 99.62 22.76 3.68 -0.05572 -0.07712 374 88 1 2 0 129.600 20.02 22.76 28.21 25.79 25.79 33.9 1035.07 900000 -1040 -609 -196 299 730 797 881 847 662 544 229 -8 -210 -383 -444 -276 -10 -1 71 299 195 50 157 105 134 277 250 -366 -497 -833 -1014 -770 -278 -216 71 301 131 -55 -111 8 -22 244 507 656 694 674 580 336 -33 -436 -747 -1071 -912 -724 -435 -139 249 349 453 470 454 490 576 567 431 320 123 -169 -474 -520 -643 -546 -346 -469 -514 -476 -368 -679 -886 -787 -668 -619 -338 -229 -174 4 -4 -26 -82 61 215 337 322 487 548 289 158 -285 -648 -792 -715 -552 -309 -356 -4 198 -77 -139 -206 -324 -212 180 292 484 709 781 700 533 549 378 245 284 224 -71 271 371 225 373 434 7 -33 -29 -403 -480 -303 -556 -462 22 190 120 363 420 347 332 298 125 -119 -152 -445 -561 900000 -1021 -600 -198 308 755 830 923 894 696 570 237 -8 -215 -387 -448 -280 -10 -1 74 315 206 53 164 110 142 292 262 -381 -507 -843 -1019 -768 -283 -223 75 316 136 -57 -114 8 -23 251 522 684 722 702 600 348 -33 -439 -740 -1052 -888 -713 -431 -140 254 357 463 478 463 498 587 578 436 322 123 -168 -469 -513 -631 -533 -342 -465 -515 -472 -362 -658 -852 -755 -645 -605 -328 -223 -171 4 -4 -26 -81 60 212 334 322 486 546 284 154 -274 -616 -750 -669 -526 -297 -349 -4 198 -76 -136 -200 -318 -209 178 289 485 713 780 705 538 548 376 240 277 218 -70 265 362 220 364 422 7 -32 -28 -387 -457 -287 -530 -447 21 187 118 358 416 342 329 296 124 -116 -146 -425 -536 SQ 20231012T040718.5Z 8685.171741 255 1 49.12 2 47.31 182.08 41.55 223.03 0.00 1 0.00 1 0 22 4.873e-03 4.862e-03 2.422e-04 2.418e-04 2.163112e+03 1.390575e+01 1.508069e+02 -6.000000e+00 2.620518e+02 285.12 -2.00 20.00 -2.00 505.25 0.00e+00 1.389777e-03 6.574638e-06 8.135630e-05 -6.000000e+00 1.353242e-04 285.12 -2.00 20.00 -2.00 500.45 0.00e+00 -7.534289e-02 4.683311e-03 1.996860e-02 -6.000000e+00 9.445929e-02 226.43 -2.00 3.49 -2.00 502.04 0.00e+00 1.472 0.000 489.00 513.00 1 4.4881e+00 3.5819e-04 3.2810e-03 -6.0000e+00 7.2290e-03 -3.5232e-02 9.4274e-05 4.1011e-04 -6.0000e+00 1.9022e-03 1.7607e-02 -3.3325e-03 2.0564e-05 2.9745e-05 -6.0000e+00 4.2349e-04 -5.5964e-02 2.1394e-04 12 8 8 10 0 0 99.93 22.86 3.74 -0.05596 -0.07908 220 52 1 2 0 220.220 20.02 22.86 29.26 25.90 25.90 34.2 1038.73 900000 -1135 -666 -270 307 733 847 939 920 735 524 226 -68 -206 -412 -446 -330 -34 -3 109 308 259 64 204 112 151 321 260 -324 -518 -831 -1091 -814 -323 -259 116 282 149 -52 -64 26 -31 263 534 666 740 690 576 303 -52 -471 -763 -1164 -955 -792 -423 -98 265 385 468 494 480 573 581 578 442 349 124 -151 -509 -578 -612 -580 -323 -481 -514 -465 -484 -731 -903 -810 -724 -628 -338 -257 -251 -5 -54 -15 -76 69 251 347 315 523 555 317 151 -262 -690 -837 -780 -602 -329 -332 -9 185 -90 -204 -265 -372 -168 246 338 514 740 822 734 607 605 466 258 297 271 -83 301 419 218 369 401 16 -39 -60 -407 -543 -325 -654 -465 66 225 141 399 407 305 369 317 134 -140 -203 -503 -569 900000 -1110 -662 -270 312 760 883 983 971 772 551 235 -70 -211 -418 -451 -336 -34 -3 114 324 272 67 215 117 159 337 272 -337 -528 -831 -1099 -813 -329 -268 121 295 156 -54 -66 27 -32 270 555 691 775 714 601 313 -53 -474 -756 -1147 -932 -771 -419 -99 272 393 476 500 488 586 597 593 450 350 125 -150 -501 -568 -600 -565 -318 -475 -512 -459 -474 -707 -879 -785 -698 -612 -330 -249 -248 -5 -53 -15 -74 68 249 345 315 523 555 313 148 -254 -653 -791 -738 -571 -319 -327 -9 184 -88 -200 -258 -362 -167 244 335 515 743 829 737 611 610 463 252 290 264 -81 294 409 213 359 390 16 -38 -58 -387 -512 -308 -625 -450 64 221 138 392 402 303 364 314 132 -136 -195 -480 -543