File name: Pandora180s1_BronxNY_20221005_L2Fit_fwvt1c2d20210302p1-8.txt
File generation date: 20260323T084001.6Z
Data description: Level 2 spectral fitting results file
Data file version: fwvt1c2p1-8
Local principal investigator: Lukas Valin
Network principal investigator: Alexander Cede
Instrument type: Pandora
Instrument number: 180
Spectrometer number: 1
Processing software version used: BlickP v1.8.88
Instrument operation file used: Pandora180_OF_v3d20201002.txt
Instrument calibration file used: Pandora180s1_CF_v2d20210302.txt
Level 1 file used: Pandora180s1_BronxNY_20221005_L1_smca1c2d20210302p1-8.txt
Full location name: New York State DEC Bronx-Pfizer PAMS station (36-005-0133)
Short location name: BronxNY
Country of location: USA
Location latitude [deg]: 40.8679
Location longitude [deg]: -73.8781
Location altitude [m]: 31
Local noon date: 20221005
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: 1639 1787
Nominal wavelengths inside fitting window [nm]: 492.1095 492.2310 492.3525 492.4740 492.5954 492.7169 492.8383 492.9597 493.0811 493.2025 493.3238 493.4452 493.5665 493.6878 493.8091 493.9304 494.0516 494.1729 494.2941 494.4153 494.5365 494.6577 494.7788 494.9000 495.0211 495.1422 495.2633 495.3844 495.5054 495.6265 495.7475 495.8685 495.9895 496.1104 496.2314 496.3523 496.4733 496.5942 496.7150 496.8359 496.9568 497.0776 497.1984 497.3192 497.4400 497.5608 497.6815 497.8023 497.9230 498.0437 498.1644 498.2850 498.4057 498.5263 498.6469 498.7675 498.8881 499.0087 499.1292 499.2497 499.3702 499.4907 499.6112 499.7317 499.8521 499.9726 500.0930 500.2134 500.3337 500.4541 500.5744 500.6948 500.8151 500.9353 501.0556 501.1759 501.2961 501.4163 501.5365 501.6567 501.7769 501.8970 502.0172 502.1373 502.2574 502.3775 502.4975 502.6176 502.7376 502.8576 502.9776 503.0976 503.2176 503.3375 503.4574 503.5774 503.6972 503.8171 503.9370 504.0568 504.1766 504.2964 504.4162 504.5360 504.6558 504.7755 504.8952 505.0149 505.1346 505.2543 505.3739 505.4935 505.6132 505.7328 505.8523 505.9719 506.0914 506.2110 506.3305 506.4500 506.5694 506.6889 506.8083 506.9278 507.0472 507.1666 507.2859 507.4053 507.5246 507.6439 507.7632 507.8825 508.0018 508.1210 508.2403 508.3595 508.4787 508.5978 508.7170 508.8361 508.9553 509.0744 509.1935 509.3125 509.4316 509.5506 509.6696 509.7886 509.9076
---------------------------------------------------------------------------------------
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-249: Unweighted slant column residuals for each pixel inside the fitting window multiplied by 1e5, 9e5=pixel was not used for fitting
Columns 250-398: 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 20221005T154809.1Z 8313.658439 227 1 1.05 2 47.52 161.01 144.05 67.08 0.00 1 0.00 1 0 15 4.458e-03 4.469e-03 2.866e-03 2.861e-03 1.675575e+03 1.733333e+02 1.086789e+02 -6.000000e+00 2.565161e+02 285.41 -2.00 20.00 -2.00 505.26 0.00e+00 5.689097e-04 8.002836e-05 3.293813e-05 -6.000000e+00 1.273611e-04 285.41 -2.00 20.00 -2.00 500.46 0.00e+00 1.545652e-01 5.718710e-02 1.540309e-02 -6.000000e+00 8.947328e-02 227.04 -2.00 3.51 -2.00 502.05 0.00e+00 1.478 0.000 489.00 513.00 1 5.3478e+00 4.3924e-03 2.5829e-03 -6.0000e+00 6.8727e-03 3.4217e-03 1.1315e-03 3.4361e-04 -6.0000e+00 1.7752e-03 7.4053e-03 -2.4896e-04 3.0549e-04 1.6264e-05 -6.0000e+00 4.8962e-04 5.6351e-02 -2.2872e-04 12 9 1 12 1 1 999.00 23.97 0.30 0.05635 0.05462 1 28 1 4 0 416.970 27.88 20.00 23.97 19.70 30.8 971.48 -1284 -288 144 175 29 280 655 1105 211 -122 224 -853 -337 -512 -718 -92 -43 62 -472 -110 343 179 668 836 455 377 606 -15 -548 -770 -153 -1025 69 370 324 -217 12 42 -90 597 -372 531 163 154 406 530 211 364 -228 -453 -355 -827 -612 -629 100 294 -443 -235 -159 255 473 805 283 431 618 508 34 -161 -164 -249 -503 -1061 3 291 266 -167 -377 -435 -631 -240 -438 -596 -183 516 -134 -140 -217 17 -208 200 45 100 745 338 228 151 -299 -589 -619 -384 -280 34 234 -393 -407 -60 -107 -757 156 45 -12 497 609 563 636 882 -30 452 114 495 -324 625 203 33 -245 -343 329 -79 370 -75 190 -314 -654 -452 -317 56 -259 39 92 -95 81 453 620 525 500 -297 -596 -685 -575 -1288 -283 145 177 30 295 696 1170 224 -128 231 -876 -339 -512 -721 -93 -44 65 -492 -115 358 187 696 872 477 394 625 -15 -561 -774 -152 -1044 72 381 338 -228 13 43 -92 619 -380 542 169 160 424 551 217 368 -234 -452 -346 -823 -591 -621 100 294 -451 -240 -162 258 481 816 288 435 616 514 34 -160 -160 -243 -491 -1051 3 290 264 -164 -366 -417 -605 -230 -424 -582 -179 508 -130 -136 -211 17 -204 196 45 101 748 338 225 146 -293 -565 -584 -360 -266 33 228 -393 -409 -59 -105 -743 152 45 -12 489 605 562 638 886 -30 451 112 494 -323 613 200 32 -237 -332 325 -78 363 -73 183 -299 -619 -427 -303 55 -249 38 91 -94 81 450 617 522 492 -289 -570 -664 -556
SQ 20221005T162954.1Z 8313.687432 255 1 49.29 2 45.89 175.14 136.85 76.34 0.00 1 0.00 1 0 15 4.106e-03 4.135e-03 1.359e-03 1.353e-03 3.071734e+03 8.605429e+01 2.165757e+02 -6.000000e+00 2.362967e+02 285.40 -2.00 20.00 -2.00 505.20 0.00e+00 4.339226e-04 3.786157e-05 2.616981e-05 -6.000000e+00 1.172975e-04 285.40 -2.00 20.00 -2.00 500.28 0.00e+00 1.146907e-01 2.746055e-02 2.127569e-02 -6.000000e+00 8.233428e-02 227.04 -2.00 3.51 -2.00 501.87 0.00e+00 1.435 0.000 489.00 513.00 1 7.9350e+00 2.1135e-03 2.7985e-03 -6.0000e+00 6.3256e-03 6.9335e-03 5.4238e-04 5.1924e-04 -6.0000e+00 1.6317e-03 5.6181e-04 -1.3473e-03 1.4329e-04 2.4434e-05 -6.0000e+00 4.5064e-04 4.4525e-02 -8.4388e-05 11 8 0 10 0 0 99.87 23.98 0.20 0.04452 0.03518 14 3 1 4 0 3519.620 27.88 20.00 23.98 20.30 30.9 975.47 -819 -704 -673 -227 210 301 764 715 774 196 -209 -482 -696 -741 -762 -447 76 185 145 -94 213 -203 451 575 456 746 462 -90 -461 -880 -819 -383 -145 557 127 40 268 -291 -194 -289 -37 174 260 527 962 532 376 575 -15 -312 -870 -935 -446 -558 172 29 327 189 260 753 267 260 297 169 297 -77 -149 13 -396 -120 -651 -374 265 4 48 -593 -455 -400 -174 -214 -179 -216 93 -110 -41 133 -162 66 162 61 648 428 64 567 250 262 56 -262 -592 -579 -610 84 258 58 -272 -296 -116 -127 288 -66 81 247 383 182 554 537 376 323 -291 -107 -202 -145 -42 5 2 100 102 549 73 -212 -111 -444 -477 -593 -242 -561 -87 133 58 -231 138 545 590 426 289 -502 -320 -521 -394 -832 -731 -669 -231 218 318 815 764 821 206 -217 -496 -707 -747 -769 -455 78 195 155 -102 231 -221 488 623 495 804 493 -94 -473 -888 -825 -390 -151 593 138 43 286 -309 -204 -304 -39 184 278 566 1033 571 400 599 -16 -312 -853 -905 -431 -551 174 29 341 198 272 784 280 273 312 176 305 -79 -149 12 -389 -117 -637 -368 262 4 47 -585 -441 -384 -165 -205 -172 -210 91 -109 -41 132 -162 66 163 62 659 411 61 541 251 257 54 -247 -552 -538 -576 82 257 55 -261 -282 -109 -119 269 -62 78 240 374 179 547 532 373 318 -285 -103 -192 -137 -40 5 2 95 97 515 69 -198 -102 -404 -430 -536 -221 -524 -82 127 56 -225 135 533 578 416 280 -482 -300 -482 -364
SQ 20221005T205306.4Z 8313.870213 412 1 48.85 2 72.95 247.33 90.06 117.14 0.00 1 0.00 1 0 15 3.813e-03 3.802e-03 2.697e-04 2.689e-04 4.754776e+03 1.787803e+01 3.338484e+02 -6.000000e+00 2.193741e+02 285.38 -2.00 20.00 -2.00 505.30 0.00e+00 8.205115e-04 7.543321e-06 4.746408e-05 -6.000000e+00 1.089051e-04 285.38 -2.00 20.00 -2.00 500.42 0.00e+00 4.119577e-01 5.413704e-03 3.247180e-02 -6.000000e+00 7.646470e-02 227.02 -2.00 3.51 -2.00 502.03 0.00e+00 3.359 0.001 489.00 513.00 1 4.6807e+00 4.1451e-04 5.6493e-03 -6.0000e+00 5.8743e-03 -5.7955e-04 1.0636e-04 8.2100e-04 -6.0000e+00 1.5158e-03 1.0730e-02 -1.0785e-03 2.8171e-05 3.7106e-05 -6.0000e+00 4.1851e-04 4.5697e-02 -1.0999e-04 11 8 0 10 0 0 99.93 23.97 0.29 0.04570 0.03822 353 87 1 4 0 133.060 27.88 20.00 23.97 19.62 30.9 970.40 -979 -589 -436 -249 75 411 510 495 377 267 62 -322 -499 -498 -601 -479 -139 103 33 96 219 346 457 593 611 556 482 -102 -520 -660 -572 -343 173 369 71 112 92 -75 -154 -95 31 90 296 451 520 537 389 235 -141 -512 -803 -798 -754 -329 26 65 144 83 253 346 495 459 451 430 224 160 115 -209 -390 -444 -343 -295 -35 76 -125 -412 -431 -375 -343 -293 -244 -232 -47 49 61 47 19 -43 -18 166 557 660 443 459 457 275 -126 -378 -719 -775 -677 -234 141 34 -252 -245 -415 -298 -197 -13 74 218 286 456 598 651 508 253 39 -175 -210 -57 81 77 81 215 153 224 318 63 -192 -314 -502 -564 -213 -130 -149 33 212 121 258 404 498 410 66 -253 -629 -647 -564 -988 -572 -437 -256 79 439 551 536 405 283 65 -334 -513 -508 -614 -494 -146 110 36 99 228 360 475 619 638 578 496 -109 -544 -677 -585 -356 184 380 75 117 96 -77 -157 -97 31 92 307 471 544 560 402 238 -140 -526 -808 -792 -751 -334 26 65 146 85 258 353 506 470 462 438 226 159 113 -203 -375 -423 -327 -284 -34 74 -122 -397 -409 -351 -319 -273 -230 -220 -45 47 59 46 18 -42 -18 163 553 660 445 455 444 261 -117 -345 -649 -698 -620 -220 137 34 -253 -245 -411 -293 -194 -13 75 222 293 470 588 643 503 250 38 -177 -210 -57 80 77 81 215 153 223 315 62 -187 -302 -479 -541 -206 -128 -142 32 207 120 255 400 495 406 65 -246 -600 -609 -530