File name: Pandora162s1_Brussels-Uccle_20220204_L2Fit_fwvt1c15d20210909p1-8.txt File generation date: 20250630T111357.4Z Data description: Level 2 spectral fitting results file Data file version: fwvt1c15p1-8 Local principal investigator: Michel Van Roozendael Network principal investigator: Alexander Cede Instrument type: Pandora Instrument number: 162 Spectrometer number: 1 Processing software version used: BlickP v1.8.75 Instrument operation file used: Pandora162_OF_v8d20191031.txt Instrument calibration file used: Pandora162s1_CF_v15d20210909.txt Level 1 file used: Pandora162s1_Brussels-Uccle_20220204_L1_smca1c15d20210909p1-8.txt Full location name: Uccle (BIRA-IASB, KMI-IRM) Short location name: Brussels-Uccle Country of location: Belgium Location latitude [deg]: 50.7980 Location longitude [deg]: 4.3580 Location altitude [m]: 107 Local noon date: 20220204 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: 1648 1798 Nominal wavelengths inside fitting window [nm]: 492.0476 492.1683 492.2890 492.4097 492.5303 492.6510 492.7716 492.8922 493.0128 493.1334 493.2539 493.3745 493.4950 493.6155 493.7360 493.8564 493.9769 494.0973 494.2178 494.3382 494.4585 494.5789 494.6993 494.8196 494.9399 495.0602 495.1805 495.3008 495.4211 495.5413 495.6615 495.7817 495.9019 496.0221 496.1422 496.2624 496.3825 496.5026 496.6227 496.7427 496.8628 496.9828 497.1029 497.2229 497.3428 497.4628 497.5828 497.7027 497.8226 497.9425 498.0624 498.1823 498.3021 498.4220 498.5418 498.6616 498.7814 498.9012 499.0209 499.1407 499.2604 499.3801 499.4998 499.6194 499.7391 499.8587 499.9783 500.0980 500.2175 500.3371 500.4567 500.5762 500.6957 500.8152 500.9347 501.0542 501.1736 501.2931 501.4125 501.5319 501.6513 501.7707 501.8900 502.0094 502.1287 502.2480 502.3673 502.4865 502.6058 502.7250 502.8443 502.9635 503.0827 503.2018 503.3210 503.4401 503.5593 503.6784 503.7974 503.9165 504.0356 504.1546 504.2736 504.3926 504.5116 504.6306 504.7496 504.8685 504.9874 505.1063 505.2252 505.3441 505.4630 505.5818 505.7006 505.8194 505.9382 506.0570 506.1758 506.2945 506.4132 506.5320 506.6506 506.7693 506.8880 507.0066 507.1253 507.2439 507.3625 507.4810 507.5996 507.7181 507.8367 507.9552 508.0737 508.1921 508.3106 508.4291 508.5475 508.6659 508.7843 508.9027 509.0210 509.1394 509.2577 509.3760 509.4943 509.6126 509.7308 509.8491 509.9673 --------------------------------------------------------------------------------------- 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-252: Unweighted slant column residuals for each pixel inside the fitting window multiplied by 1e5, 9e5=pixel was not used for fitting Columns 253-403: 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 20220204T121009.9Z 8070.507059 158 1 1.79 2 67.01 183.57 65.48 136.26 0.00 1 0.00 1 0 15 3.266e-03 3.287e-03 2.699e-03 2.695e-03 1.510676e+03 1.516456e+02 1.232851e+02 -6.000000e+00 1.751364e+02 267.18 -2.00 20.00 -2.00 505.26 0.00e+00 4.386170e-04 7.205238e-05 2.461910e-05 -6.000000e+00 8.906712e-05 267.18 -2.00 20.00 -2.00 500.46 0.00e+00 5.678868e-02 5.380677e-02 1.144963e-02 -6.000000e+00 6.562684e-02 223.14 -2.00 3.79 -2.00 502.05 0.00e+00 2.540 0.000 489.00 513.00 1 4.4123e+00 4.1129e-03 3.7466e-03 -6.0000e+00 5.0172e-03 9.0851e-03 1.0541e-03 3.4810e-04 -6.0000e+00 1.2897e-03 1.6477e-02 8.1175e-04 2.8871e-04 2.0564e-05 -6.0000e+00 3.5963e-04 -1.2581e-02 -1.4166e-05 12 1 1 12 1 1 999.00 18.49 1.90 -0.01258 -0.00695 1 15 2 2 0 673.890 15.05 18.49 14.78 20.96 21.17 45.7 998.80 -268 -125 -36 126 -512 277 48 -54 312 300 487 -99 -92 -57 -364 92 141 270 -177 -164 -606 676 -249 391 316 221 133 177 76 -861 -159 178 -308 -32 498 -84 -90 580 -254 -331 -156 302 -111 171 342 697 -404 330 544 -415 -758 -918 193 -70 280 -134 -433 293 -679 -131 510 551 216 171 747 149 -214 -90 -399 -530 -610 -80 -285 40 140 -228 255 -176 -56 59 -258 -162 -118 -403 -87 -32 -152 -154 269 469 31 10 111 -8 43 -354 795 92 -97 -168 -182 -330 -450 -619 -109 -654 30 36 -66 243 -107 187 -150 111 30 657 449 -290 597 162 449 -48 205 171 18 -78 -151 74 95 150 270 145 4 47 -117 -159 150 63 -16 28 166 -395 173 -96 413 -23 -208 256 -332 -450 -198 -260 -122 -36 127 -525 288 50 -57 327 311 498 -102 -93 -57 -365 92 142 279 -183 -173 -632 713 -259 407 329 231 138 185 78 -871 -160 178 -313 -33 519 -88 -93 594 -259 -342 -161 310 -113 176 354 722 -421 338 558 -420 -765 -906 189 -69 274 -133 -441 294 -692 -133 514 556 220 174 751 149 -216 -89 -392 -526 -601 -79 -281 39 137 -225 248 -172 -54 56 -249 -157 -115 -391 -84 -32 -151 -153 267 465 30 10 112 -9 43 -354 786 90 -93 -161 -173 -316 -432 -602 -108 -647 30 36 -64 237 -106 187 -150 112 30 665 457 -292 605 163 444 -47 202 169 17 -76 -148 72 92 148 265 141 4 45 -112 -151 144 61 -16 27 166 -390 171 -95 408 -23 -208 250 -319 -436 -191 SQ 20220204T124638.4Z 8070.532388 173 1 49.93 2 67.83 193.01 61.78 145.28 0.00 1 0.00 1 0 15 1.888e-03 1.889e-03 1.805e-04 1.801e-04 1.890277e+03 1.039042e+01 1.491797e+02 -6.000000e+00 1.012738e+02 267.18 -2.00 20.00 -2.00 505.27 0.00e+00 4.800494e-04 4.806558e-06 2.903625e-05 -6.000000e+00 5.150524e-05 267.18 -2.00 20.00 -2.00 500.40 0.00e+00 9.558488e-02 3.605838e-03 1.497038e-02 -6.000000e+00 3.795442e-02 223.14 -2.00 3.79 -2.00 502.01 0.00e+00 2.628 0.000 489.00 513.00 1 2.9627e+00 2.7505e-04 3.9644e-03 -6.0000e+00 2.9016e-03 5.5242e-03 7.0626e-05 4.1929e-04 -6.0000e+00 7.4602e-04 6.9063e-02 1.0823e-03 1.9214e-05 2.4111e-05 -6.0000e+00 2.0797e-04 -9.4748e-03 -4.4006e-05 10 0 0 10 0 0 99.98 18.50 1.89 -0.00947 -0.00197 831 167 1 2 0 59.570 14.94 18.50 15.47 20.96 21.17 46.1 1000.30 -61 -86 45 72 -155 203 409 333 279 230 183 181 -36 -70 -40 -164 -135 46 151 -107 -101 134 98 185 306 157 67 194 -15 -289 -43 -74 -227 143 141 -307 -245 31 0 -103 14 30 -22 45 143 150 94 52 80 -108 -356 -261 -378 -410 -177 -127 -274 -266 -155 -15 144 273 170 129 170 109 -10 -51 -15 -124 -191 -66 -211 -172 50 -99 -344 -185 26 -85 -190 -119 -310 -269 -31 -82 -182 -86 -36 35 2 71 304 174 10 115 53 -65 -74 -153 -298 -451 -409 -64 -124 -484 -382 -140 -122 81 152 -13 39 114 153 196 309 444 271 158 275 163 94 321 249 -17 131 95 -20 225 294 -48 -13 49 -160 -16 280 6 -206 74 120 15 131 184 103 100 26 -156 -219 -206 -142 -60 -86 44 72 -160 213 437 359 301 246 193 188 -37 -71 -40 -166 -139 48 161 -111 -106 140 103 193 320 165 71 200 -16 -300 -43 -75 -233 150 146 -322 -257 32 0 -105 14 31 -23 47 149 158 98 54 81 -108 -362 -259 -370 -402 -177 -124 -274 -271 -158 -15 147 278 174 132 173 110 -10 -50 -15 -119 -190 -65 -202 -166 48 -96 -331 -175 26 -83 -186 -112 -293 -257 -29 -80 -178 -85 -36 34 2 71 307 176 11 115 52 -62 -70 -141 -273 -416 -386 -63 -124 -490 -385 -140 -120 79 151 -13 40 112 152 196 310 445 271 156 268 164 94 319 248 -17 126 91 -19 224 291 -47 -12 47 -153 -15 271 5 -198 72 117 15 130 183 103 99 25 -150 -208 -192 -133 SQ 20220204T161936.0Z 8070.680278 236 1 49.83 2 88.03 241.25 58.27 205.90 0.00 1 0.00 1 0 15 1.410e-02 1.414e-02 1.146e-02 1.142e-02 4.912424e+03 7.121098e+02 4.081799e+02 -6.000000e+00 7.561057e+02 267.18 -2.00 20.00 -2.00 505.39 0.00e+00 2.461428e-03 3.069743e-04 1.186309e-04 -6.000000e+00 3.846695e-04 267.18 -2.00 20.00 -2.00 500.73 0.00e+00 1.546166e+00 2.265156e-01 3.526773e-02 -6.000000e+00 2.836840e-01 223.15 -2.00 3.79 -2.00 502.29 0.00e+00 16.298 0.094 489.00 513.00 1 6.8742e+00 1.7280e-02 2.6714e-02 -6.0000e+00 2.1683e-02 2.1119e-02 4.3384e-03 1.6924e-03 -6.0000e+00 5.5830e-03 1.7792e-04 5.2610e-03 1.2158e-03 7.3017e-05 -6.0000e+00 1.5538e-03 1.0027e-02 -2.4950e-04 12 8 8 10 0 0 99.64 18.51 2.65 0.01003 0.04652 23 5 2 2 0 2165.920 14.94 18.51 10.07 20.96 21.06 45.3 993.38 963 -454 1290 271 -433 433 -848 1699 1509 -2191 -2794 -718 -1032 802 1301 -329 -707 837 394 854 -559 1078 1053 514 1973 2297 2251 -1877 737 -553 446 -1837 -909 123 -1217 1787 -197 -1860 690 580 2670 -1348 -1262 -840 1050 223 -1720 110 -1049 -1933 -742 -844 -1151 553 -432 -576 -42 -457 400 -1638 570 -303 -912 2191 739 1239 2366 50 418 1713 1611 -2183 647 -708 -1627 -941 -742 -1469 1833 -733 1685 -600 -1689 -719 703 -1647 477 -93 2857 -2506 1415 1058 2434 614 1531 214 -168 -1470 -2283 1560 -1930 924 -2525 -1059 1078 -1224 -186 -895 -1185 280 -1885 -2045 219 971 911 930 3735 807 1835 -1587 -399 -1826 -201 599 1180 923 -870 678 -220 -172 1398 -1188 -41 1553 -210 -2104 1487 -1369 -938 324 -555 3474 567 -457 391 3654 1761 -688 427 -2250 -3419 838 -406 1170 255 -422 429 -865 1705 1517 -2250 -2848 -707 -1006 760 1220 -315 -691 823 395 864 -578 1098 1077 528 2006 2332 2279 -1958 734 -544 427 -1802 -900 124 -1271 1832 -206 -1958 702 587 2651 -1395 -1320 -886 1099 236 -1841 115 -1086 -1974 -730 -813 -1101 523 -422 -580 -43 -473 412 -1726 588 -317 -963 2244 763 1259 2343 49 408 1635 1528 -2150 624 -701 -1626 -932 -724 -1423 1710 -694 1558 -570 -1637 -697 677 -1627 463 -91 2709 -2512 1373 1043 2394 614 1508 210 -161 -1398 -2140 1383 -1758 827 -2387 -1022 1047 -1232 -185 -889 -1166 272 -1897 -2095 222 990 938 959 3770 843 1894 -1680 -412 -1886 -203 601 1183 927 -891 684 -224 -175 1395 -1204 -41 1496 -206 -2110 1464 -1414 -982 339 -589 3590 603 -493 420 3776 1833 -721 435 -2334 -3589