File name: Pandora111s1_Bucharest_20241201_L2Fit_fwvt1c1d20200901p1-8.txt File generation date: 20241201T221015.5Z Data description: Level 2 spectral fitting results file Data file version: fwvt1c1p1-8 Local principal investigator: Anca Nemuc Network principal investigator: Alexander Cede Instrument type: Pandora Instrument number: 111 Spectrometer number: 1 Processing software version used: BlickP v1.8.60 Instrument operation file used: Pandora111_OF_v1d20200312.txt Instrument calibration file used: Pandora111s1_CF_v1d20200901.txt Level 1 file used: Pandora111s1_Bucharest_20241201_L1_smca1c1d20200901p1-8.txt Full location name: MARS Atmospheric Observatory Short location name: Bucharest Country of location: Romania Location latitude [deg]: 44.3440 Location longitude [deg]: 26.0120 Location altitude [m]: 93 Local noon date: 20241201 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: 1632 1780 Nominal wavelengths inside fitting window [nm]: 492.1036 492.2257 492.3477 492.4698 492.5918 492.7138 492.8357 492.9577 493.0796 493.2016 493.3235 493.4454 493.5672 493.6891 493.8109 493.9327 494.0545 494.1763 494.2981 494.4199 494.5416 494.6633 494.7850 494.9067 495.0284 495.1500 495.2716 495.3933 495.5149 495.6364 495.7580 495.8796 496.0011 496.1226 496.2441 496.3656 496.4871 496.6085 496.7299 496.8513 496.9727 497.0941 497.2155 497.3368 497.4582 497.5795 497.7008 497.8220 497.9433 498.0645 498.1858 498.3070 498.4282 498.5493 498.6705 498.7916 498.9127 499.0339 499.1549 499.2760 499.3971 499.5181 499.6391 499.7601 499.8811 500.0021 500.1230 500.2440 500.3649 500.4858 500.6067 500.7275 500.8484 500.9692 501.0900 501.2108 501.3316 501.4524 501.5731 501.6938 501.8145 501.9352 502.0559 502.1766 502.2972 502.4178 502.5384 502.6590 502.7796 502.9002 503.0207 503.1412 503.2617 503.3822 503.5027 503.6231 503.7436 503.8640 503.9844 504.1048 504.2251 504.3455 504.4658 504.5861 504.7064 504.8267 504.9470 505.0672 505.1874 505.3076 505.4278 505.5480 505.6682 505.7883 505.9084 506.0285 506.1486 506.2687 506.3888 506.5088 506.6288 506.7488 506.8688 506.9888 507.1087 507.2287 507.3486 507.4685 507.5884 507.7082 507.8281 507.9479 508.0677 508.1875 508.3073 508.4270 508.5468 508.6665 508.7862 508.9059 509.0256 509.1452 509.2649 509.3845 509.5041 509.6237 509.7432 509.8628 509.9823 --------------------------------------------------------------------------------------- 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-250: Unweighted slant column residuals for each pixel inside the fitting window multiplied by 1e5, 9e5=pixel was not used for fitting Columns 251-399: 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 20241201T072506.8Z 9101.309107 99 1 40.44 2 75.66 142.00 80.56 144.39 0.00 1 0.00 1 0 15 3.991e-03 3.948e-03 1.771e-03 1.758e-03 1.862058e+03 1.014613e+02 1.352135e+02 -6.000000e+00 2.084176e+02 275.34 -2.00 20.00 -2.00 505.33 0.00e+00 8.282801e-04 4.720436e-05 4.701887e-05 -6.000000e+00 1.088302e-04 275.34 -2.00 20.00 -2.00 500.27 0.00e+00 3.306417e-01 3.451156e-02 1.690406e-02 -6.000000e+00 7.759156e-02 222.28 -2.00 3.58 -2.00 501.92 0.00e+00 3.951 0.001 489.00 513.00 1 6.8989e+00 2.6315e-03 5.8705e-03 -6.0000e+00 5.9272e-03 1.1942e-02 6.8364e-04 4.8189e-04 -6.0000e+00 1.5511e-03 1.0915e-03 2.7429e-04 1.5100e-04 2.1401e-05 -6.0000e+00 3.5280e-04 2.1802e-02 -1.4264e-04 11 8 0 10 0 0 99.95 15.97 0.80 0.02180 0.02371 95 31 2 2 0 425.580 15.01 15.97 9.89 23.43 23.11 37.9 1023.16 -1253 -600 -633 194 304 645 681 680 490 220 273 -225 -382 -555 -406 -412 -153 -32 222 274 228 232 278 276 371 323 186 -134 -211 -793 -653 -384 -69 199 213 457 46 -55 -223 -51 87 132 461 279 425 457 342 36 -336 -542 -881 -778 -731 -452 -146 25 164 289 328 339 609 472 430 343 163 148 -268 -175 -411 -261 -249 -109 -187 -204 -168 -219 -264 -754 -523 -556 -552 -292 -271 -47 -284 -125 -129 -189 -123 132 333 230 446 437 409 173 -397 -445 -997 -736 -620 -180 -163 110 33 -160 -319 -343 -351 -105 82 335 540 773 563 782 404 497 322 171 181 45 214 124 176 227 369 63 -41 23 -245 -211 -266 -319 -229 -78 131 299 296 217 323 308 270 145 -219 -253 -661 -581 -587 -1186 -580 -636 200 320 691 742 744 533 236 288 -235 -394 -566 -412 -423 -160 -35 243 301 251 255 307 306 412 358 205 -146 -221 -809 -661 -394 -73 219 238 506 51 -59 -240 -55 93 143 505 309 473 506 375 39 -355 -557 -882 -765 -721 -456 -153 24 163 286 324 334 602 470 428 340 159 142 -254 -164 -381 -239 -229 -101 -177 -195 -159 -204 -241 -683 -472 -505 -505 -269 -252 -44 -270 -119 -123 -179 -117 127 323 227 441 427 392 162 -365 -400 -882 -648 -557 -168 -157 108 32 -155 -303 -324 -335 -102 81 334 541 777 570 795 411 499 318 167 173 43 208 121 171 221 357 61 -40 22 -230 -195 -244 -294 -216 -76 129 296 295 217 324 312 274 147 -219 -246 -630 -546 -556 SQ 20241201T112418.2Z 9101.475211 176 1 2.79 2 68.69 199.70 73.96 196.63 0.00 1 0.00 1 0 15 4.483e-03 4.428e-03 2.249e-03 2.244e-03 1.095627e+03 1.212853e+02 9.003939e+01 -6.000000e+00 2.341219e+02 275.33 -2.00 20.00 -2.00 505.29 0.00e+00 7.916999e-04 5.969007e-05 4.410718e-05 -6.000000e+00 1.222558e-04 275.33 -2.00 20.00 -2.00 500.54 0.00e+00 1.970375e-01 4.323129e-02 9.707301e-03 -6.000000e+00 8.716667e-02 222.28 -2.00 3.58 -2.00 502.12 0.00e+00 2.726 0.000 489.00 513.00 1 6.1802e+00 3.2942e-03 4.0793e-03 -6.0000e+00 6.6586e-03 1.1204e-02 8.6441e-04 2.9226e-04 -6.0000e+00 1.7424e-03 2.6932e-03 3.8249e-04 1.9495e-04 1.8573e-05 -6.0000e+00 3.9634e-04 2.1765e-02 -1.4146e-04 12 9 1 12 1 1 999.00 16.11 0.76 0.02176 0.02442 2 15 1 2 0 930.170 15.00 16.11 11.13 23.53 23.22 37.8 1022.63 -1108 -584 -234 -49 727 752 558 668 500 284 214 -59 -655 -657 -371 -480 -180 -250 495 115 -217 124 194 656 394 546 761 -187 -86 -847 -625 -244 38 138 94 359 228 42 -297 -353 100 55 354 568 467 313 371 174 -506 -603 -1089 -742 -990 -658 -437 -54 411 373 692 413 399 421 199 250 256 -54 -400 -275 -206 -742 -210 -302 97 -40 -120 -666 -22 -693 -267 -736 -514 -287 153 -132 -230 -182 -333 162 19 212 111 480 370 359 443 209 -263 -333 -844 -1175 -573 -391 -63 -362 -114 -51 -341 -65 -100 -462 150 782 328 722 796 511 429 554 329 66 -82 140 134 192 136 332 509 494 148 14 -320 -194 -495 -320 -134 -164 525 302 -32 307 459 510 387 -158 -537 86 -418 -750 -614 -1059 -571 -237 -48 732 773 584 702 522 292 217 -59 -649 -644 -361 -472 -180 -258 519 121 -230 130 205 693 419 578 797 -194 -86 -827 -605 -239 39 145 101 380 239 43 -305 -361 103 57 370 600 496 332 387 179 -512 -592 -1042 -731 -933 -636 -437 -56 423 384 706 422 411 436 207 258 260 -54 -395 -268 -198 -709 -201 -291 95 -39 -118 -646 -21 -653 -250 -696 -490 -275 148 -129 -227 -180 -329 159 19 211 113 493 382 367 444 206 -253 -312 -780 -1086 -538 -380 -64 -355 -116 -51 -336 -64 -99 -467 153 768 326 719 796 515 431 549 321 67 -82 140 135 194 137 334 509 494 148 14 -312 -185 -471 -306 -130 -164 504 294 -32 302 453 506 385 -158 -530 83 -409 -728 -600 SQ 20241201T113215.5Z 9101.480735 180 1 2.29 2 69.19 201.60 74.39 198.36 0.00 1 0.00 1 0 15 3.764e-03 3.751e-03 1.117e-03 1.116e-03 1.172197e+03 5.990098e+01 7.859502e+01 -6.000000e+00 1.965721e+02 275.33 -2.00 20.00 -2.00 505.27 0.00e+00 7.639383e-04 2.988099e-05 4.403779e-05 -6.000000e+00 1.026465e-04 275.33 -2.00 20.00 -2.00 500.55 0.00e+00 1.593291e-01 2.149409e-02 1.326075e-02 -6.000000e+00 7.317457e-02 222.28 -2.00 3.58 -2.00 502.13 0.00e+00 2.788 0.000 489.00 513.00 1 4.8695e+00 1.6439e-03 4.2198e-03 -6.0000e+00 5.5899e-03 1.4118e-02 4.2880e-04 3.3894e-04 -6.0000e+00 1.4625e-03 9.9314e-03 1.4652e-04 9.6272e-05 1.7330e-05 -6.0000e+00 3.3274e-04 1.9513e-02 -1.2173e-04 12 1 1 12 1 1 99.03 16.12 0.77 0.01951 0.02053 8 41 1 2 0 254.030 15.00 16.12 11.39 23.53 23.22 38.1 1024.78 -1085 -653 -313 176 344 740 858 532 520 167 227 -87 -203 -416 -511 -442 -82 -216 238 303 262 103 362 339 257 276 162 -193 -316 -666 -561 -439 -89 164 281 318 180 55 -80 -107 -95 49 455 562 492 550 294 -41 -263 -427 -798 -693 -802 -439 -275 74 243 310 525 415 578 340 274 317 170 181 -84 -236 -233 -290 -442 -489 -309 -79 -254 -233 -334 -538 -509 -622 -444 -384 -181 -155 -148 -116 -116 -151 -73 191 306 245 396 337 186 72 -464 -474 -758 -547 -344 -195 -87 -160 -181 -65 -118 -222 -183 -129 198 327 551 633 840 661 509 422 214 150 63 185 190 111 172 283 338 247 56 -126 -49 -129 -177 -255 -123 -147 73 243 186 213 340 222 183 47 -127 -154 -608 -631 -540 -1065 -623 -310 179 357 781 875 548 531 177 236 -89 -205 -417 -512 -446 -85 -219 244 312 270 106 373 351 268 286 167 -196 -325 -666 -557 -442 -93 169 293 329 184 55 -80 -107 -95 50 464 580 511 568 300 -41 -272 -429 -781 -665 -773 -434 -282 74 245 312 524 414 580 344 278 318 168 185 -85 -235 -229 -283 -433 -485 -311 -81 -256 -231 -325 -518 -489 -602 -433 -377 -180 -155 -149 -117 -117 -153 -70 186 302 246 398 335 182 73 -458 -456 -717 -515 -330 -194 -85 -160 -181 -64 -113 -213 -177 -127 198 330 560 646 859 681 524 429 214 148 61 179 186 109 170 278 331 242 54 -122 -47 -126 -171 -249 -116 -143 72 242 186 215 343 226 187 48 -128 -151 -581 -596 -513