File name: Pandora289s1_DeBilt_20260219_L2Fit_fwvt1c2d20250401p1-8.txt File generation date: 20260219T235501.0Z Data description: Level 2 spectral fitting results file Data file version: fwvt1c2p1-8 Local principal investigator: Ankie Piters Network principal investigator: Alexander Cede Instrument type: Pandora Instrument number: 289 Spectrometer number: 1 Processing software version used: BlickP v1.8.88 Instrument operation file used: Pandora289_OF_v3d20240910.txt Instrument calibration file used: Pandora289s1_CF_v2d20250401.txt Level 1 file used: Pandora289s1_DeBilt_20260219_L1_smca1c2d20250401p1-8.txt Full location name: Royal Netherlands Meteorological Institute (KNMI), De Bilt Short location name: DeBilt Country of location: Netherlands Location latitude [deg]: 52.1012 Location longitude [deg]: 5.1783 Location altitude [m]: 20 Local noon date: 20260219 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: OOB session based on clearance session 1 First and last pixel inside fitting window: 1626 1776 Nominal wavelengths inside fitting window [nm]: 492.1097 492.2301 492.3504 492.4707 492.5910 492.7113 492.8315 492.9518 493.0720 493.1922 493.3124 493.4326 493.5527 493.6729 493.7930 493.9131 494.0332 494.1533 494.2733 494.3934 494.5134 494.6334 494.7534 494.8733 494.9933 495.1132 495.2332 495.3531 495.4730 495.5928 495.7127 495.8325 495.9523 496.0721 496.1919 496.3117 496.4315 496.5512 496.6709 496.7906 496.9103 497.0300 497.1496 497.2693 497.3889 497.5085 497.6281 497.7476 497.8672 497.9867 498.1062 498.2257 498.3452 498.4647 498.5841 498.7036 498.8230 498.9424 499.0617 499.1811 499.3005 499.4198 499.5391 499.6584 499.7777 499.8969 500.0162 500.1354 500.2546 500.3738 500.4930 500.6121 500.7313 500.8504 500.9695 501.0886 501.2077 501.3267 501.4458 501.5648 501.6838 501.8028 501.9218 502.0407 502.1597 502.2786 502.3975 502.5164 502.6353 502.7541 502.8729 502.9918 503.1106 503.2293 503.3481 503.4669 503.5856 503.7043 503.8230 503.9417 504.0604 504.1790 504.2976 504.4162 504.5348 504.6534 504.7720 504.8905 505.0090 505.1276 505.2461 505.3645 505.4830 505.6014 505.7198 505.8383 505.9566 506.0750 506.1934 506.3117 506.4300 506.5483 506.6666 506.7849 506.9031 507.0214 507.1396 507.2578 507.3760 507.4941 507.6123 507.7304 507.8485 507.9666 508.0847 508.2028 508.3208 508.4388 508.5568 508.6748 508.7928 508.9108 509.0287 509.1466 509.2645 509.3824 509.5003 509.6182 509.7360 509.8538 509.9716 --------------------------------------------------------------------------------------- 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 --------------------------------------------------------------------------------------- SO 20260219T111731.2Z 9546.470499 1692 1 36.50 2 63.75 170.26 59.54 142.85 0.00 1 0.00 1 0 15 3.048e-03 3.025e-03 8.364e-04 8.347e-04 1.931851e+03 4.357968e+01 1.545015e+02 -6.000000e+00 1.487720e+02 267.18 -2.00 20.00 -2.00 505.27 0.00e+00 4.552930e-04 2.147816e-05 2.527470e-05 -6.000000e+00 8.015323e-05 267.18 -2.00 20.00 -2.00 500.48 0.00e+00 2.098504e-01 1.594956e-02 1.499182e-02 -6.000000e+00 5.849917e-02 223.86 -2.00 3.84 -2.00 502.07 0.00e+00 2.248 0.000 489.00 513.00 1 7.3371e+00 1.2127e-03 3.4950e-03 -6.0000e+00 4.4496e-03 7.2609e-03 3.2028e-04 4.1483e-04 -6.0000e+00 1.1779e-03 9.0663e-04 8.0200e-04 6.4075e-05 2.4452e-05 -6.0000e+00 2.4025e-04 6.8829e-04 3.9914e-05 10 0 0 10 0 0 99.99 24.00 1.53 0.00069 0.00625 71 7 1 4 0 514.040 33.50 19.99 24.00 7.82 0.1 919.39 -630 -281 -270 -25 149 393 592 475 151 90 67 -108 -176 -335 -51 -476 -125 235 219 236 372 354 246 67 511 388 398 7 -368 -227 -339 -289 -23 206 284 22 -169 -61 -299 -359 -277 -57 -39 125 416 253 140 254 -100 -243 -442 -477 -701 -639 -383 -15 45 -37 60 283 169 433 504 309 288 272 6 -183 -383 -400 -264 -186 -44 259 -44 20 18 -543 -197 -161 -394 -393 -451 -27 -177 19 -226 -126 51 -322 -53 234 156 367 251 367 301 -180 -483 -363 -976 -425 -175 -160 131 69 -53 -174 -334 -423 -165 176 179 63 307 626 286 510 287 247 129 133 112 200 134 399 225 373 102 372 51 65 -48 -224 -152 -289 -311 -84 90 -22 288 114 215 278 282 301 -22 -414 -654 -343 -258 -621 -276 -268 -26 155 405 620 501 158 92 69 -111 -179 -337 -51 -479 -128 242 231 250 392 372 260 71 543 412 420 7 -373 -227 -333 -287 -24 216 296 22 -179 -63 -307 -367 -283 -58 -41 132 429 262 148 263 -102 -242 -438 -459 -668 -618 -384 -15 46 -39 62 288 173 445 520 320 295 273 6 -179 -372 -394 -256 -181 -43 254 -44 20 18 -526 -188 -154 -380 -383 -443 -26 -172 18 -223 -124 50 -317 -52 235 156 367 249 365 292 -172 -464 -340 -908 -400 -166 -160 130 69 -52 -172 -325 -415 -164 173 179 64 310 633 292 523 293 248 127 129 110 193 131 388 219 362 99 360 49 64 -46 -211 -143 -275 -304 -82 88 -22 287 114 216 274 279 297 -22 -403 -619 -328 -243 SO 20260219T113508.6Z 9546.482738 1820 1 37.75 2 63.40 175.07 57.96 147.54 0.00 1 0.00 1 0 15 3.968e-03 3.939e-03 1.309e-03 1.306e-03 2.338280e+03 6.921190e+01 1.926754e+02 -6.000000e+00 1.936836e+02 267.18 -2.00 20.00 -2.00 505.27 0.00e+00 4.598861e-04 3.360236e-05 2.485705e-05 -6.000000e+00 1.043582e-04 267.18 -2.00 20.00 -2.00 500.48 0.00e+00 2.701127e-01 2.493995e-02 1.601876e-02 -6.000000e+00 7.619088e-02 223.86 -2.00 3.84 -2.00 502.07 0.00e+00 2.221 0.000 489.00 513.00 1 8.2237e+00 1.8962e-03 3.5043e-03 -6.0000e+00 5.7950e-03 3.0503e-03 4.9932e-04 4.6194e-04 -6.0000e+00 1.5347e-03 3.7301e-04 1.8260e-03 1.0015e-04 2.9492e-05 -6.0000e+00 3.1290e-04 3.2047e-03 1.7337e-07 11 8 0 10 0 0 99.95 24.01 1.56 0.00320 0.01587 11 1 1 4 0 3431.140 33.50 19.89 24.01 8.48 0.1 921.79 -771 49 -747 -328 -69 423 714 404 227 -21 256 -510 314 -607 -153 -291 -193 190 637 203 354 158 172 129 662 160 365 -14 -230 -81 -312 -77 -96 183 326 0 -299 -161 -250 -563 -405 338 107 101 1000 103 -82 456 -191 -446 -436 -684 -636 -1145 -326 -120 279 80 72 310 384 570 505 252 170 558 -1 -647 -570 -159 -422 -108 110 678 -261 367 -30 -896 363 -83 -163 -461 -443 -46 -53 -49 -445 -541 71 -332 7 395 377 613 293 490 522 -293 -542 -217 -942 -442 -438 -102 -3 120 -77 159 -551 -371 -179 308 393 -154 331 697 168 710 -12 138 -65 207 118 82 84 283 -12 361 -26 594 -340 106 266 -109 -11 -473 -308 -415 63 -304 40 255 501 354 418 213 -88 -478 -1090 -52 60 -739 48 -745 -335 -72 448 743 423 236 -22 266 -526 316 -610 -152 -292 -198 195 665 213 371 166 181 136 699 170 383 -14 -240 -81 -307 -76 -100 191 348 -1 -315 -166 -256 -574 -413 346 111 106 1051 110 -86 470 -195 -443 -432 -660 -606 -1114 -327 -121 286 82 74 314 391 583 519 260 174 557 -1 -634 -571 -156 -411 -105 108 659 -260 357 -29 -871 347 -80 -157 -450 -435 -44 -51 -48 -438 -533 69 -326 7 395 384 607 297 484 503 -278 -519 -202 -875 -416 -415 -101 -3 119 -78 155 -536 -363 -179 301 390 -154 333 699 170 721 -12 138 -64 205 116 82 82 274 -11 348 -25 571 -329 103 254 -102 -11 -451 -301 -402 62 -302 40 254 501 357 423 216 -88 -464 -1033 -50 58 SO 20260219T113833.9Z 9546.485115 1846 1 34.69 2 63.36 176.01 57.67 148.47 0.00 1 0.00 1 0 15 3.994e-03 3.960e-03 1.287e-03 1.284e-03 2.253166e+03 6.763171e+01 1.836519e+02 -6.000000e+00 1.949301e+02 267.18 -2.00 20.00 -2.00 505.27 0.00e+00 4.842663e-04 3.304500e-05 2.669231e-05 -6.000000e+00 1.050245e-04 267.18 -2.00 20.00 -2.00 500.47 0.00e+00 2.783171e-01 2.457676e-02 1.629986e-02 -6.000000e+00 7.666298e-02 223.86 -2.00 3.84 -2.00 502.06 0.00e+00 2.218 0.000 489.00 513.00 1 8.1203e+00 1.8679e-03 3.5122e-03 -6.0000e+00 5.8310e-03 2.8150e-03 4.9345e-04 4.5719e-04 -6.0000e+00 1.5439e-03 4.1350e-04 1.2183e-03 9.8515e-05 2.8838e-05 -6.0000e+00 3.1484e-04 1.6270e-03 2.3086e-05 11 8 0 10 0 0 99.89 24.01 1.44 0.00163 0.01008 10 1 1 4 0 3468.270 33.50 19.91 24.01 8.54 0.1 922.05 -575 235 -559 -225 22 416 595 480 12 -78 178 -309 267 -729 -187 -479 -170 151 435 189 485 -3 32 -27 661 217 440 260 -452 154 -351 -85 182 390 479 -354 -295 16 -29 -195 -404 -33 148 241 916 -8 68 146 -367 -164 -579 -947 -914 -1090 -575 -142 323 130 -40 201 546 565 473 189 229 175 392 -396 -334 -263 -291 -7 -8 865 -115 415 -163 -794 511 -138 -426 -411 -653 -164 -115 -18 -425 -468 -22 -416 6 134 67 629 259 576 502 -139 -539 -245 -1198 -519 -500 30 136 564 -239 -91 -220 -492 -212 290 98 -54 128 628 290 778 -79 115 24 143 409 273 136 349 343 447 55 730 -158 52 78 -209 -163 -644 -462 -427 24 -85 231 159 241 290 90 372 -141 -240 -715 -113 -74 -553 232 -559 -232 23 430 624 506 13 -80 186 -320 272 -739 -187 -484 -175 156 459 200 511 -3 34 -29 703 231 465 271 -460 154 -347 -85 182 409 498 -369 -313 16 -30 -199 -415 -34 154 247 942 -8 70 152 -377 -163 -578 -921 -878 -1066 -564 -144 333 135 -41 206 558 582 491 191 236 176 387 -390 -325 -252 -285 -7 -7 846 -115 407 -157 -777 490 -133 -414 -391 -627 -158 -113 -17 -422 -464 -22 -412 6 135 67 628 256 573 488 -132 -520 -230 -1123 -491 -477 30 136 559 -234 -90 -215 -485 -206 285 98 -55 130 636 297 797 -79 116 24 138 392 264 133 341 335 434 53 707 -154 50 75 -198 -154 -619 -441 -417 23 -85 232 160 237 287 89 367 -138 -235 -681 -106 -70