File name: Pandora101s1_Izana_20210107_L2Fit_fwvt1c15d20200721p1-8.txt File generation date: 20240504T132039.6Z Data description: Level 2 spectral fitting results file Data file version: fwvt1c15p1-8 Local principal investigator: Alberto Redondas Network principal investigator: Alexander Cede Instrument type: Pandora Instrument number: 101 Spectrometer number: 1 Processing software version used: BlickP v1.8.59 Instrument operation file used: Pandora101_OF_v6d20170928.txt Instrument calibration file used: Pandora101s1_CF_v15d20200721.txt Level 1 file used: Pandora101s1_Izana_20210107_L1_smca1c15d20200721p1-8.txt Full location name: Izana Atmospheric Research Center Short location name: Izana Country of location: Spain Location latitude [deg]: 28.3090 Location longitude [deg]: -16.4994 Location altitude [m]: 2360 Local noon date: 20210107 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: based on session 11 for vp 20200721 First and last pixel inside fitting window: 1643 1792 Nominal wavelengths inside fitting window [nm]: 492.0538 492.1747 492.2957 492.4166 492.5375 492.6584 492.7793 492.9002 493.0210 493.1419 493.2627 493.3835 493.5043 493.6250 493.7458 493.8665 493.9873 494.1080 494.2286 494.3493 494.4700 494.5906 494.7112 494.8318 494.9524 495.0730 495.1936 495.3141 495.4346 495.5551 495.6756 495.7961 495.9166 496.0370 496.1574 496.2778 496.3982 496.5186 496.6390 496.7593 496.8796 496.9999 497.1202 497.2405 497.3608 497.4810 497.6012 497.7215 497.8416 497.9618 498.0820 498.2021 498.3223 498.4424 498.5625 498.6826 498.8026 498.9227 499.0427 499.1627 499.2827 499.4027 499.5227 499.6426 499.7625 499.8824 500.0023 500.1222 500.2421 500.3619 500.4818 500.6016 500.7214 500.8412 500.9609 501.0807 501.2004 501.3201 501.4398 501.5595 501.6792 501.7988 501.9185 502.0381 502.1577 502.2773 502.3968 502.5164 502.6359 502.7554 502.8749 502.9944 503.1139 503.2333 503.3528 503.4722 503.5916 503.7110 503.8304 503.9497 504.0690 504.1884 504.3077 504.4269 504.5462 504.6655 504.7847 504.9039 505.0231 505.1423 505.2615 505.3806 505.4998 505.6189 505.7380 505.8571 505.9762 506.0952 506.2142 506.3333 506.4523 506.5712 506.6902 506.8092 506.9281 507.0470 507.1659 507.2848 507.4037 507.5225 507.6414 507.7602 507.8790 507.9978 508.1166 508.2353 508.3540 508.4728 508.5915 508.7102 508.8288 508.9475 509.0661 509.1847 509.3033 509.4219 509.5405 509.6590 509.7776 509.8961 --------------------------------------------------------------------------------------- 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 Columns 98-247: Unweighted slant column residuals for each pixel inside the fitting window multiplied by 1e5, 9e5=pixel was not used for fitting Columns 248-397: 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 20210107T080455.8Z 7677.336757 218 1 23.57 2 89.71 115.74 36.56 185.05 0.00 1 0.00 1 0 22 5.574e-03 5.506e-03 1.677e-04 1.673e-04 1.052434e+04 1.154751e+01 5.555590e+02 -6.000000e+00 3.022134e+02 272.96 -2.00 20.00 -2.00 505.38 0.00e+00 1.426239e-03 4.546652e-06 9.109827e-05 -6.000000e+00 1.528672e-04 272.96 -2.00 20.00 -2.00 500.91 0.00e+00 1.185033e+00 3.265564e-03 7.187968e-02 -6.000000e+00 1.107840e-01 223.09 -2.00 3.19 -2.00 502.43 0.00e+00 20.524 0.089 489.00 513.00 1 1.7604e+00 2.4878e-04 2.5083e-02 -6.0000e+00 8.4628e-03 -3.4891e-02 6.3098e-05 2.2166e-03 -6.0000e+00 2.2109e-03 3.2713e-02 1.0638e-03 1.5877e-05 1.0241e-04 -6.0000e+00 5.6878e-04 -3.0883e-03 4.7579e-05 12 9 9 12 1 1 99.76 16.85 1.59 -0.00309 0.00429 405 199 1 4 0 57.494 21.48 15.00 16.85 -1183 -841 -410 180 494 806 1086 1070 1031 738 330 176 -151 -352 -477 -400 -198 81 406 411 360 436 373 427 585 516 245 81 -133 -807 -1035 -802 -575 -66 241 169 21 -15 -201 -263 -201 -2 218 419 555 637 528 284 -24 -276 -822 -1065 -1166 -1016 -437 -34 154 299 379 428 461 560 625 561 493 290 15 -278 -369 -515 -681 -602 -505 -299 -149 -329 -623 -776 -793 -823 -870 -692 -668 -529 -406 -403 -324 -331 -307 -192 -6 227 495 604 566 419 58 -360 -572 -1029 -1222 -1037 -680 -246 -4 -157 -299 -389 -439 -220 33 218 503 814 1004 1130 1123 987 752 487 356 278 284 451 465 353 476 404 266 342 282 -168 -337 -429 -548 -279 27 86 268 458 349 334 347 354 284 93 -189 -615 -949 -936 -1060 -765 -374 169 463 773 1059 1028 991 724 320 169 -144 -332 -449 -379 -191 78 398 409 352 427 367 422 581 515 244 80 -132 -785 -993 -772 -566 -66 244 169 21 -15 -202 -263 -202 -2 217 423 565 651 537 287 -24 -276 -805 -1022 -1112 -977 -430 -34 156 304 387 437 472 562 631 567 494 295 15 -277 -364 -503 -663 -588 -499 -299 -150 -330 -616 -757 -767 -794 -842 -674 -657 -525 -408 -399 -321 -328 -304 -190 -6 228 502 611 563 416 56 -349 -544 -969 -1154 -999 -664 -245 -4 -160 -304 -392 -441 -222 34 226 513 833 1031 1170 1176 1021 779 512 368 283 287 458 476 364 491 418 275 354 290 -171 -340 -432 -556 -288 28 90 285 483 372 360 377 378 303 99 -203 -651 -994 -979 SQ 20210107T080825.3Z 7677.339181 223 1 18.52 2 89.02 116.15 36.65 186.46 0.00 1 0.00 1 0 22 4.593e-03 4.550e-03 1.368e-04 1.366e-04 8.822684e+03 9.249482e+00 4.741124e+02 -6.000000e+00 2.490753e+02 272.96 -2.00 20.00 -2.00 505.35 0.00e+00 1.234646e-03 3.703795e-06 7.846875e-05 -6.000000e+00 1.259916e-04 272.96 -2.00 20.00 -2.00 500.80 0.00e+00 1.218225e+00 2.676008e-03 6.048959e-02 -6.000000e+00 9.132215e-02 223.09 -2.00 3.19 -2.00 502.35 0.00e+00 19.132 0.072 489.00 513.00 1 1.2819e+00 2.0398e-04 2.2807e-02 -6.0000e+00 6.9759e-03 -1.9755e-02 5.1912e-05 1.9234e-03 -6.0000e+00 1.8229e-03 6.1451e-02 1.3537e-03 1.3095e-05 8.6475e-05 -6.0000e+00 4.6877e-04 -5.2764e-05 1.2860e-05 11 8 0 10 0 0 98.43 16.86 1.37 -0.00005 0.00934 520 123 1 4 0 35.618 21.58 15.00 16.86 -944 -696 -317 173 413 672 916 860 874 606 245 140 -80 -230 -366 -291 -127 75 353 302 241 359 322 343 440 461 155 34 -101 -663 -880 -661 -475 -82 188 85 -3 -32 -137 -223 -118 -26 169 327 433 494 448 228 -5 -220 -666 -895 -1003 -846 -424 -88 105 238 320 370 382 483 501 475 429 268 20 -222 -278 -450 -541 -472 -399 -256 -120 -235 -488 -636 -609 -690 -805 -624 -573 -405 -346 -334 -295 -215 -253 -147 15 170 400 502 451 316 46 -291 -479 -816 -978 -853 -594 -193 -19 -169 -251 -287 -346 -209 -2 185 432 695 872 943 904 820 606 381 280 279 215 344 362 262 352 341 229 268 236 -126 -272 -362 -419 -234 26 92 231 393 319 291 327 322 215 57 -163 -515 -796 -755 -864 -647 -296 166 398 662 895 849 863 594 238 135 -78 -222 -352 -282 -122 74 347 301 241 360 324 347 448 470 157 35 -100 -656 -857 -646 -463 -82 189 87 -3 -32 -136 -222 -117 -26 171 335 438 500 452 233 -5 -217 -657 -865 -962 -819 -420 -88 107 239 322 372 385 489 510 483 433 268 20 -222 -275 -441 -527 -461 -394 -256 -121 -236 -482 -620 -588 -665 -779 -607 -563 -401 -339 -331 -293 -214 -250 -146 15 172 407 510 451 308 44 -277 -457 -769 -924 -821 -580 -192 -20 -172 -254 -287 -346 -210 -2 186 439 710 894 974 943 845 625 397 288 282 217 347 368 268 361 349 235 275 241 -127 -271 -359 -419 -237 26 95 242 409 336 309 349 338 226 61 -172 -535 -817 -774 SQ 20210107T080841.5Z 7677.339369 223 2 12.95 2 88.96 116.19 36.66 186.57 0.00 1 0.00 1 0 22 4.533e-03 4.492e-03 1.564e-04 1.562e-04 8.772595e+03 1.057506e+01 4.683498e+02 -6.000000e+00 2.458024e+02 272.96 -2.00 20.00 -2.00 505.35 0.00e+00 1.227325e-03 4.238484e-06 7.839767e-05 -6.000000e+00 1.243331e-04 272.96 -2.00 20.00 -2.00 500.81 0.00e+00 1.216671e+00 3.061289e-03 6.064471e-02 -6.000000e+00 9.010610e-02 223.09 -2.00 3.19 -2.00 502.35 0.00e+00 19.000 0.071 489.00 513.00 1 1.1933e+00 2.3343e-04 2.2623e-02 -6.0000e+00 6.8832e-03 -1.7833e-02 5.9351e-05 1.9144e-03 -6.0000e+00 1.7982e-03 6.8118e-02 1.0404e-03 1.4983e-05 8.5776e-05 -6.0000e+00 4.6261e-04 -5.3078e-03 5.9934e-05 12 9 1 12 1 1 99.12 16.86 1.33 -0.00531 0.00191 358 123 1 4 0 35.618 21.58 15.00 16.86 -896 -667 -358 143 371 686 883 848 840 606 267 155 -102 -239 -342 -284 -152 53 354 331 254 357 298 306 472 431 170 39 -92 -646 -869 -636 -501 -85 176 96 -12 -17 -139 -203 -139 2 152 357 439 492 435 232 -7 -203 -661 -890 -977 -854 -418 -78 107 239 305 348 389 483 506 474 400 259 44 -235 -266 -406 -543 -465 -383 -246 -112 -285 -517 -606 -582 -688 -757 -611 -552 -463 -314 -325 -274 -223 -240 -147 51 156 368 456 435 335 33 -262 -462 -822 -1005 -814 -600 -211 -37 -175 -248 -280 -347 -195 41 191 443 705 845 911 923 785 604 354 311 242 195 369 366 279 345 360 218 281 217 -162 -278 -305 -446 -207 28 121 247 393 321 312 322 326 198 50 -153 -519 -799 -751 -815 -616 -333 137 357 673 861 835 827 593 259 152 -99 -229 -328 -273 -149 52 354 329 254 357 300 309 478 438 172 39 -90 -635 -843 -618 -498 -86 177 98 -12 -17 -138 -205 -141 2 153 365 452 497 448 237 -7 -204 -650 -857 -955 -823 -412 -78 109 244 306 349 391 487 513 481 403 264 45 -234 -262 -396 -539 -452 -377 -245 -113 -284 -509 -588 -572 -674 -745 -592 -539 -458 -314 -321 -271 -220 -237 -145 51 157 373 461 433 332 32 -253 -437 -771 -943 -779 -583 -209 -37 -177 -250 -279 -345 -194 41 196 449 718 864 938 960 808 621 368 318 244 195 371 371 284 352 368 222 286 220 -163 -281 -307 -444 -209 29 125 259 408 337 330 343 342 208 53 -161 -537 -834 -783