File name: Pandora111s1_Bucharest_20221219_L2Fit_fwvt1c1d20200901p1-8.txt File generation date: 20240917T180109.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.59 Instrument operation file used: Pandora111_OF_v1d20200312.txt Instrument calibration file used: Pandora111s1_CF_v1d20200901.txt Level 1 file used: Pandora111s1_Bucharest_20221219_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: 20221219 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 20221219T134921.5Z 8388.575943 341 1 49.73 2 83.84 228.36 112.56 272.98 0.00 1 0.00 1 0 15 5.667e-03 5.561e-03 2.299e-04 2.291e-04 1.900708e+03 1.232274e+01 1.327687e+02 -6.000000e+00 2.939221e+02 273.57 -2.00 20.00 -2.00 505.35 0.00e+00 3.048198e-03 6.114559e-06 1.642073e-04 -6.000000e+00 1.538821e-04 273.57 -2.00 20.00 -2.00 500.99 0.00e+00 6.696183e-01 4.350635e-03 3.150485e-02 -6.000000e+00 1.102999e-01 222.26 -2.00 3.58 -2.00 502.48 0.00e+00 8.377 0.011 489.00 513.00 1 3.2834e+00 3.3227e-04 1.2844e-02 -6.0000e+00 8.4206e-03 -3.2449e-02 8.5735e-05 8.2119e-04 -6.0000e+00 2.2044e-03 2.0628e-02 2.4176e-04 2.0053e-05 5.1912e-05 -6.0000e+00 5.0097e-04 1.3174e-02 -9.4160e-05 12 8 8 10 0 0 58.12 14.87 1.65 0.01317 0.01485 498 104 1 2 0 99.650 15.01 14.87 6.55 22.40 22.09 42.1 1031.79 -1601 -1109 -441 102 505 820 966 848 719 344 132 -191 -532 -685 -728 -550 -219 1 278 348 365 372 398 526 532 616 389 5 -493 -918 -927 -714 -129 292 494 426 178 -56 -210 -155 -57 233 467 631 629 632 492 192 -171 -714 -1093 -1267 -1159 -642 -220 204 383 451 454 476 601 627 669 528 278 111 -165 -368 -566 -691 -551 -432 -210 -99 -250 -380 -671 -785 -829 -793 -662 -601 -436 -329 -293 -279 -252 -289 -162 36 308 499 598 562 398 107 -324 -743 -1085 -1157 -854 -443 -52 76 -48 -206 -405 -431 -249 4 363 628 848 1011 1083 1011 892 683 429 245 138 212 273 345 367 380 420 369 275 54 -217 -451 -665 -521 -354 -131 135 328 338 327 397 406 381 221 -58 -439 -815 -1033 -930 -1341 -952 -392 93 471 781 933 825 696 328 125 -178 -490 -625 -664 -505 -205 1 273 344 360 367 394 525 534 618 388 4 -471 -854 -858 -671 -125 294 504 434 179 -56 -208 -154 -57 236 478 654 657 659 509 196 -171 -698 -1043 -1192 -1095 -622 -221 209 397 468 470 493 628 659 704 553 288 113 -166 -368 -560 -677 -540 -429 -212 -101 -254 -378 -657 -764 -806 -775 -652 -596 -436 -333 -299 -285 -257 -295 -166 37 322 493 594 553 416 109 -325 -730 -1049 -1115 -837 -449 -51 75 -47 -201 -389 -414 -242 4 365 642 874 1048 1130 1063 939 713 441 248 138 214 277 352 375 389 430 378 282 55 -217 -445 -652 -514 -357 -135 142 347 360 350 427 439 415 240 -62 -462 -838 -1052 -955 SQ 20221219T135527.9Z 8388.580184 344 1 49.67 2 84.66 229.50 113.62 274.02 0.00 1 0.00 1 0 15 6.324e-03 6.205e-03 3.029e-04 3.014e-04 2.154112e+03 1.619171e+01 1.480491e+02 -6.000000e+00 3.279730e+02 273.57 -2.00 20.00 -2.00 505.34 0.00e+00 3.454232e-03 7.992939e-06 1.838235e-04 -6.000000e+00 1.717092e-04 273.57 -2.00 20.00 -2.00 501.14 0.00e+00 7.148819e-01 5.625024e-03 3.538672e-02 -6.000000e+00 1.230817e-01 222.26 -2.00 3.58 -2.00 502.59 0.00e+00 9.374 0.016 489.00 513.00 1 3.5724e+00 4.2743e-04 1.4445e-02 -6.0000e+00 9.3963e-03 -3.8968e-02 1.1130e-04 9.2463e-04 -6.0000e+00 2.4599e-03 1.3309e-02 3.8830e-04 2.6381e-05 5.9875e-05 -6.0000e+00 5.5902e-04 9.4208e-03 -7.6633e-05 12 8 8 10 0 0 1.49 14.89 2.19 0.00942 0.01211 365 76 1 2 0 136.050 15.01 14.89 6.55 22.40 22.09 41.9 1031.44 -1836 -1257 -489 100 507 823 1047 925 813 350 105 -263 -571 -767 -804 -617 -239 98 350 396 448 378 439 551 612 719 441 39 -522 -1048 -1076 -800 -91 292 592 492 190 20 -199 -171 -22 281 538 722 759 718 583 228 -165 -764 -1186 -1398 -1236 -704 -184 286 498 519 542 542 679 701 725 590 334 151 -174 -439 -600 -735 -642 -513 -254 -130 -291 -428 -739 -904 -953 -888 -765 -646 -492 -340 -310 -281 -291 -279 -200 29 366 601 678 580 457 74 -383 -840 -1230 -1261 -963 -468 -51 129 22 -238 -474 -455 -305 29 400 683 958 1139 1239 1141 1007 738 444 253 175 216 345 397 440 396 410 355 244 20 -240 -572 -750 -596 -379 -118 125 363 335 354 390 395 397 264 -72 -499 -917 -1149 -1001 -1629 -1142 -423 89 463 765 987 878 767 326 97 -240 -514 -685 -716 -555 -219 92 336 382 431 365 426 538 601 706 431 38 -490 -957 -979 -738 -87 289 593 491 188 20 -194 -166 -22 280 543 738 781 739 594 230 -163 -739 -1120 -1303 -1157 -676 -182 291 512 534 556 557 702 730 757 613 343 153 -174 -437 -590 -716 -626 -507 -255 -132 -293 -424 -721 -876 -923 -865 -751 -638 -489 -341 -314 -286 -295 -282 -203 30 381 637 721 612 474 75 -382 -822 -1185 -1211 -941 -473 -53 138 23 -249 -491 -470 -319 31 434 697 987 1180 1292 1199 1060 771 457 257 176 218 352 407 452 408 421 366 252 20 -241 -568 -741 -592 -385 -122 132 386 360 382 423 431 435 289 -78 -529 -952 -1183 -1039 SQ 20221219T141717.5Z 8388.595342 347 1 50.02 2 87.72 233.48 117.40 277.82 0.00 1 0.00 1 0 22 1.420e-02 1.384e-02 2.387e-03 2.345e-03 2.699416e+03 1.270794e+02 1.414816e+02 -6.000000e+00 7.363433e+02 273.57 -2.00 20.00 -2.00 505.51 0.00e+00 5.652312e-03 6.340997e-05 2.910832e-04 -6.000000e+00 3.854935e-04 273.57 -2.00 20.00 -2.00 501.69 0.00e+00 9.298103e-01 4.335914e-02 6.017376e-02 -6.000000e+00 2.762161e-01 222.26 -2.00 3.58 -2.00 503.07 0.00e+00 15.399 0.072 489.00 513.00 1 6.1673e+00 3.2894e-03 2.5526e-02 -6.0000e+00 2.1089e-02 -4.5904e-02 8.4679e-04 1.5920e-03 -6.0000e+00 5.5180e-03 4.0734e-04 9.0188e-04 2.0660e-04 9.1053e-05 -6.0000e+00 1.2547e-03 -5.2272e-03 1.9910e-05 12 8 8 10 0 0 49.95 14.88 10.39 -0.00523 0.00103 20 4 1 2 0 2500.650 14.99 14.88 5.62 22.40 21.98 42.1 1031.71 -4659 -3563 -1794 -616 629 869 1525 1913 1181 636 85 -625 -1633 -2218 -2067 -1701 -477 157 1039 907 1153 920 1206 1500 1600 1435 880 161 -1130 -2254 -2255 -1636 148 961 1721 1302 665 163 -121 9 121 988 1839 2138 1836 1922 1473 911 -41 -1403 -2485 -2796 -2324 -1003 -212 1131 1222 1249 1170 1345 1613 2097 1734 1362 684 145 -229 -907 -1429 -1441 -1530 -938 -465 -96 -395 -1041 -1752 -1905 -1957 -2417 -1498 -1349 -1220 -877 -506 -444 -755 -574 -424 130 832 1625 1519 1359 926 -271 -682 -1671 -2691 -2594 -1887 -634 196 963 -357 -491 -1363 -1030 -672 -20 1050 1535 2028 2038 2405 2332 2318 1407 753 507 114 490 875 925 668 858 533 691 367 -92 -719 -1544 -1734 -1597 -925 -434 -199 656 463 590 658 658 972 588 -87 -1096 -1826 -2562 -2018 -3928 -3067 -1584 -511 528 745 1321 1659 1024 545 72 -527 -1371 -1855 -1724 -1429 -405 136 915 803 1019 817 1075 1349 1449 1304 801 144 -994 -1947 -1943 -1430 132 887 1603 1213 616 150 -111 9 113 931 1751 2064 1795 1876 1432 878 -39 -1319 -2301 -2565 -2140 -942 -206 1114 1219 1247 1169 1345 1626 2120 1764 1381 690 145 -226 -894 -1399 -1398 -1491 -921 -463 -96 -394 -1026 -1709 -1849 -1900 -2370 -1469 -1329 -1212 -878 -508 -446 -759 -575 -426 132 854 1692 1594 1415 952 -276 -681 -1644 -2621 -2517 -1854 -639 204 1010 -378 -512 -1416 -1068 -708 -22 1140 1691 2253 2284 2711 2653 2639 1600 848 564 126 544 980 1044 759 975 608 790 421 -104 -810 -1727 -1933 -1795 -1056 -505 -217 719 512 655 735 740 1095 664 -97 -1207 -1977 -2772 -2200